殷海濤

研究背景

（一）開(kāi)展防震減災(zāi)工作的現(xiàn)實(shí)需要

沂沭斷裂帶是中國(guó)東部著名的郯（城）廬（江）深大斷裂在山東境內(nèi)的區(qū)段，也是整個(gè)郯廬斷裂帶上出露最好、規(guī)模最大、第四紀(jì)活動(dòng)最強(qiáng)烈的地區(qū)。該斷裂帶北起渤海萊州灣，向南經(jīng)昌邑、濰坊、安丘、莒縣、臨沂、郯城，延入江蘇新沂一帶，縱貫山東省中部，長(zhǎng)達(dá)330多公里。沂沭斷裂帶主要由第四紀(jì)活動(dòng)的鄌郚—葛溝斷裂（F4）、沂水—湯頭斷裂（F3）、白芬子—浮來(lái)山斷裂（F2）、安丘—莒縣斷裂（F5）和昌邑—大店斷裂（F1）5條主干斷裂組成（圖1），這五條主斷裂在魯中沂、沭河谷地構(gòu)成了20—40公里的“兩塹夾一壘”的構(gòu)造。第四紀(jì)以來(lái)，該斷裂帶新構(gòu)造運(yùn)動(dòng)強(qiáng)烈，據(jù)統(tǒng)計(jì)沂沭帶地區(qū)歷史上曾發(fā)生過(guò)5級(jí)以上地震25次，其中8級(jí)以上地震1次，為1668年郯城8.5級(jí)地震；7級(jí)以上地震4次，如公元前70年安丘7.0級(jí)地震和1969年渤海7.4級(jí)地震（圖1）。此外，與之相關(guān)聯(lián)的益都—無(wú)棣斷裂、雙山—李家莊斷裂、蒼尼斷裂、沂水—沂源等北西向分支斷裂，第四紀(jì)時(shí)期也有明顯活動(dòng)，發(fā)生過(guò)多次中強(qiáng)地震，如1995年蒼山5.2級(jí)地震。

最近幾年，沂沭斷裂帶及附近地區(qū)也發(fā)生了多次強(qiáng)有感地震和震群，如2013年11月23日萊州發(fā)生的ML5.0級(jí)地震，此次地震是繼蒼山5.2級(jí)地震后，山東內(nèi)陸發(fā)生的最大地震。所以，沂沭斷裂帶作為山東乃至華北地區(qū)的主要地震活動(dòng)區(qū)之一，不僅在歷史上發(fā)生過(guò)多次強(qiáng)烈地震，近期地震活動(dòng)也相對(duì)活躍，對(duì)其運(yùn)動(dòng)狀態(tài)進(jìn)行系統(tǒng)的分析研究具有重要的現(xiàn)實(shí)意義。

圖1 研究區(qū)地震活動(dòng)分布圖

（二）重要的理論研究?jī)r(jià)值

對(duì)長(zhǎng)期演化、構(gòu)造復(fù)雜、新構(gòu)造運(yùn)動(dòng)強(qiáng)的沂沭斷裂帶來(lái)說(shuō)，其現(xiàn)階段的構(gòu)造活動(dòng)如何，局部地區(qū)的形變演化如何，反映了怎樣的地球動(dòng)力學(xué)過(guò)程，這些問(wèn)題多年來(lái)一直備受國(guó)內(nèi)外地質(zhì)、地震學(xué)家的關(guān)注（張國(guó)民，1988；魏光興，1993；鄧起東，2002；張培震，2003）。孕育地震的能量主要來(lái)源于地殼差異運(yùn)動(dòng)產(chǎn)生的應(yīng)變能積累，因而從地殼形變觀測(cè)資料所獲取現(xiàn)今地殼運(yùn)動(dòng)微動(dòng)態(tài)定量信息，對(duì)探討地震的孕育和發(fā)生是十分必要的。近些年，GPS技術(shù)在地殼形變監(jiān)測(cè)和地震中期預(yù)測(cè)方面發(fā)揮了重要的作用。GPS空間大地測(cè)量技術(shù)在大、中空間尺度上把水平觀測(cè)精度提高到了10-8-10-9水平，使得大區(qū)域地殼運(yùn)動(dòng)監(jiān)測(cè)能力有了根本性提高，為地球動(dòng)力學(xué)和地震模型的建立提供了很好的約束條件，如用于板塊邊界帶的運(yùn)動(dòng)、同震錯(cuò)動(dòng)、斷層系統(tǒng)動(dòng)力學(xué)以及孕震過(guò)程研究等方面的工作，不僅促使動(dòng)力學(xué)研究更逼近實(shí)際、同時(shí)也為地震的動(dòng)力學(xué)預(yù)報(bào)逐步創(chuàng)造條件（周碩愚，2006）。

華北作為中國(guó)大陸開(kāi)展GPS測(cè)量最早的地區(qū)之一，積累了較為豐富的觀測(cè)資料，多位專家學(xué)者已經(jīng)做出了大量有價(jià)值的研究成果（李延興，2003；張培震，2003；江在森，2003；牛之俊，2005；郭良遷，2007；楊國(guó)華，2007；王敏，2011），為研究華北地區(qū)的地殼運(yùn)動(dòng)奠定了基礎(chǔ)。但這些成果所采用的數(shù)據(jù)大多以流動(dòng)站為主，其時(shí)-空分辨率和可靠性等都具有一定的局限性，而且針對(duì)地殼運(yùn)動(dòng)垂向變化的研究相對(duì)較少。針對(duì)這個(gè)問(wèn)題，在前人研究的基礎(chǔ)上，本研究利用沂沭斷裂帶及其附近地區(qū)130個(gè)山東CORS站點(diǎn)（Continuously Operating Reference Stations）連續(xù)5年（2011—2016）的觀測(cè)數(shù)據(jù)，結(jié)合16個(gè)跨斷層水準(zhǔn)測(cè)量和154個(gè)流動(dòng)重力測(cè)點(diǎn)資料，構(gòu)建沂沭斷裂帶構(gòu)造活動(dòng)高時(shí)-空分辨率的三維形變場(chǎng)，并通過(guò)分析不同區(qū)段的斷層運(yùn)動(dòng)速率、應(yīng)變積累水平等信息研究斷裂帶近期的地殼運(yùn)動(dòng)演化過(guò)程和方式，為全面深入認(rèn)識(shí)沂沭斷裂帶現(xiàn)今的三維形變特征以及地震危險(xiǎn)區(qū)的研判提供定量支持。

研究目標(biāo)、內(nèi)容和技術(shù)路線

（一） 研究目標(biāo)

利用高精度、高密度的GPS觀測(cè)資料獲取沂沭斷裂帶及其北西向分支斷裂的三維地殼活動(dòng)圖像，并結(jié)合水準(zhǔn)、重力等資料分析其運(yùn)動(dòng)特征；通過(guò)分析不同斷層區(qū)段的構(gòu)造形變態(tài)勢(shì)、應(yīng)變積累水平以及斷層蠕滑速率等參數(shù)，為未來(lái)地震危險(xiǎn)性判定提供定量依據(jù)，進(jìn)一步結(jié)合實(shí)際地震活動(dòng)特征，為地震的動(dòng)力學(xué)預(yù)報(bào)提供支持。

（二）研究?jī)?nèi)容

本項(xiàng)目在整理沂沭斷裂帶及其周邊地區(qū)的高質(zhì)量、高密度GPS觀測(cè)數(shù)據(jù)的基礎(chǔ)上，獲取高精度、高時(shí)-空分辨率的地殼形變場(chǎng)圖像，并結(jié)合垂向觀測(cè)手段來(lái)分析斷裂帶的時(shí)-空變化特征，并為區(qū)域地震預(yù)測(cè)提供定量依據(jù)。為達(dá)到該研究目標(biāo)，需要開(kāi)展以下幾方面的研究工作：

1.批量GPS觀測(cè)數(shù)據(jù)資料的分析整理與高精度綜合處理。主要涉及批量CORS站數(shù)據(jù)的篩選、高精度數(shù)據(jù)的解算、GPS三維坐標(biāo)時(shí)間序列分析、多源非構(gòu)造形變干擾分析及剔除等內(nèi)容，最終得到高精度的GPS地殼形變解算結(jié)果。

2.沂沭斷裂帶地區(qū)現(xiàn)今的三維形變特征及時(shí)-空演化過(guò)程分析。主要涉及參考框架的選取、區(qū)域速度場(chǎng)的計(jì)算（最優(yōu)空間插值算法）、跨斷層基線分析以及跨斷層水準(zhǔn)和流動(dòng)重力觀測(cè)資料分析等內(nèi)容，最終確定沂沭斷裂帶及其主要分支斷裂的運(yùn)動(dòng)特征。

3.區(qū)域地震危險(xiǎn)性研判。主要是通過(guò)分析斷裂構(gòu)造變形態(tài)勢(shì)、區(qū)域應(yīng)變積累水平以及斷層閉鎖程度等參數(shù)，為區(qū)域地震危險(xiǎn)性判定提供定量依據(jù)，并結(jié)合沂沭帶地區(qū)實(shí)際的地震活動(dòng)特點(diǎn)，為地震的動(dòng)力學(xué)預(yù)報(bào)提供支持。涉及內(nèi)容包括主應(yīng)變率場(chǎng)、最大剪切應(yīng)變率場(chǎng)、面應(yīng)變率場(chǎng)的計(jì)算，斷層蠕滑深度和蠕滑速率的反演等工作。

（三）技術(shù)路線

綜合考慮本項(xiàng)目實(shí)施的目標(biāo)、內(nèi)容以及現(xiàn)有的技術(shù)手段和資金預(yù)算情況，按照?qǐng)D2所示的技術(shù)路線開(kāi)展研究。

1.高時(shí)-空分辨率GPS觀測(cè)網(wǎng)點(diǎn)的篩選

圖2 技術(shù)路線流程圖

沂沭斷裂帶現(xiàn)今的運(yùn)動(dòng)整體表現(xiàn)為右旋走滑，但由于其斷層的復(fù)雜性，其不同分段也表現(xiàn)出不同的運(yùn)動(dòng)特性。為了最有效地捕獲和分辨斷裂運(yùn)動(dòng)特征，我們對(duì)沂沭斷裂帶及其周邊現(xiàn)存的連續(xù)觀測(cè)基準(zhǔn)站資料進(jìn)行分析研究。由于山東CORS站資料存在站點(diǎn)不穩(wěn)定、連續(xù)性不高等問(wèn)題，采用TEQC程序和自主開(kāi)發(fā)的Fortran程序?qū)⒉荒軡M足地殼運(yùn)動(dòng)觀測(cè)的站點(diǎn)數(shù)據(jù)進(jìn)行剔除。剔除涉及兩方面內(nèi)容：一是利用TEQC軟件對(duì)每天的原始觀測(cè)數(shù)據(jù)進(jìn)行觀測(cè)時(shí)長(zhǎng)和觀測(cè)質(zhì)量（觀測(cè)噪聲、多路徑效應(yīng)等誤差）的分析，將質(zhì)量差的站點(diǎn)進(jìn)行剔除；二是對(duì)解算結(jié)果中存在的奇異點(diǎn)進(jìn)行篩選，剔除精度較低、序列不穩(wěn)定或明顯異于周邊站點(diǎn)的資料?？紤]到站點(diǎn)的分布均勻性，對(duì)站點(diǎn)較為密集的區(qū)域進(jìn)行站點(diǎn)抽稀處理確保整個(gè)結(jié)果的正確性和可靠性。

2. GPS觀測(cè)數(shù)據(jù)的高精度處理

目前山東GPS連續(xù)觀測(cè)網(wǎng)的密度較大，站點(diǎn)的平均距離小于50km，在GPS的數(shù)據(jù)處理過(guò)程中，可克服多種誤差干擾，獲得高精度的地殼形變與構(gòu)造運(yùn)動(dòng)結(jié)果。我們使用行業(yè)通用的高精度GPS數(shù)據(jù)處理軟件GAMIT（V10.6）和GIPSY（V6.2）分別進(jìn)行解算，兩套軟件的處理策略和優(yōu)勢(shì)弱點(diǎn)有所不同，可進(jìn)行相互檢核和查錯(cuò)，確保處理結(jié)果的精度和可靠性。GPS數(shù)據(jù)后期處理，包括站點(diǎn)時(shí)間序列分析、單日解聯(lián)合平差獲取速度場(chǎng)等工作采用QOCA（V1.35）來(lái)完成。

3.參考框架的選擇

為了更好地繪制現(xiàn)階段沂沭斷裂帶3D速度場(chǎng)，挖掘沂沭斷裂帶最顯著的運(yùn)動(dòng)特征，采用自主編寫的matlab程序，分別計(jì)算該區(qū)域相對(duì)于歐亞塊體、華南塊體、鄂爾多斯塊體（張培震院士建議）的水平速度場(chǎng)，發(fā)現(xiàn)其整體運(yùn)動(dòng)特征；然后再建立沂沭帶地區(qū)的“局域自洽”參考框架（山東16個(gè)站，選取最穩(wěn)定的站點(diǎn)，CASH、RIZH、TAIN、YATI、WUDI等），發(fā)現(xiàn)其局域內(nèi)部細(xì)微的運(yùn)動(dòng)特征。

4.斷裂帶垂直運(yùn)動(dòng)特征的量化分析

山東GPS連續(xù)觀測(cè)網(wǎng)的密度較大，可以利用其垂直時(shí)間序列建立垂直速度場(chǎng)。我們選取研究區(qū)2011至2016共6年、38個(gè)連續(xù)站的時(shí)間序列數(shù)據(jù)，獲取了山東地區(qū)高精度、高時(shí)空分辨率的地殼垂直形變速率（杜存鵬等，2018，圖3），結(jié)果表明：1）山東大部分地區(qū)呈現(xiàn)差異性下降狀態(tài)，只有中部泰山山脈地區(qū)出現(xiàn)上升現(xiàn)象，與研究區(qū)繼承性構(gòu)造活動(dòng)相一致。2）山東半島區(qū)域的基準(zhǔn)站呈現(xiàn)正常的沉降趨勢(shì)，而山東省的西北部地區(qū)和西部地區(qū)由于過(guò)度開(kāi)采地下水，沉降較大。

利用跨斷層GPS基線可量化分析沂沭帶及其北西向分支斷裂現(xiàn)今的運(yùn)動(dòng)狀態(tài)。一方面，分別利用年周期法，小波分析法，發(fā)現(xiàn)其年變特性和低頻特性，并對(duì)多源非構(gòu)造形變干擾進(jìn)行剔除，得到更加純凈的構(gòu)造形變信息，進(jìn)而發(fā)現(xiàn)斷裂活動(dòng)的演化規(guī)律。另一方面，利用垂直形變速率合成法（李杰等，2010）計(jì)算跨斷層水準(zhǔn)資料不同時(shí)段的變化速率，對(duì)斷層活動(dòng)方式進(jìn)行分析。結(jié)果表明（李杰等，2018）：沂沭斷裂帶及其北西向分支斷裂近期垂直形變速率均較小，沒(méi)有明顯的速率異常和受力狀態(tài)改變；1995年蒼山5.2級(jí)、2011年日本9.0級(jí)地震（圖4）前后斷裂垂直形變速率和受力方式發(fā)生了變化。

圖3 山東區(qū)域2011—2016年垂直形變場(chǎng)

圖4 日本9.0級(jí)地震前后垂直形變特征

5.區(qū)域重力場(chǎng)分析

重力資料處理采用《LGADJ》程序進(jìn)行自由網(wǎng)平差，以獲得測(cè)網(wǎng)中心基準(zhǔn)下的重力變化。為了體現(xiàn)起算點(diǎn)的穩(wěn)定性，以泰安、日照兩個(gè)絕對(duì)重力觀測(cè)點(diǎn)為約束，采用魯南測(cè)網(wǎng)、膠東測(cè)網(wǎng)及魯中測(cè)區(qū)合并計(jì)算整體平差方法，再利用最小二乘配置方法對(duì)重力平差結(jié)果進(jìn)行擬合推估，以便突出顯示構(gòu)造因素的重力效應(yīng)。對(duì)區(qū)域重力場(chǎng)不同時(shí)段（半年尺度、年尺度等）以及多年累積變化量進(jìn)行分析，發(fā)現(xiàn)區(qū)域地殼物質(zhì)的轉(zhuǎn)化狀態(tài)。以2016年8月至2017年8月一年尺度重力場(chǎng)特征為例，顯示（圖5）重力零值線在濟(jì)南、淄博與濱州之間劇烈扭曲，形成明顯的重力梯度帶，沂沭斷裂帶北段重力場(chǎng)差異變化較大。

圖5 山東地區(qū)重力場(chǎng)變化等值線

圖6 研究區(qū)總釋放能量折算震級(jí)M與兩側(cè)塊體相對(duì)運(yùn)動(dòng)的關(guān)系

6.斷層蠕滑深度及蠕滑速率的反演

為發(fā)現(xiàn)沂沭斷裂帶淺層蠕滑特征，我們以其周邊GPS地殼形變資料為約束，基于OKADA斷裂位錯(cuò)模型，采用網(wǎng)格搜索和模擬退火方法反演確定蠕滑深度和蠕滑速率。斷層位錯(cuò)模型需要的參數(shù)，包括斷裂的長(zhǎng)度、走向、傾角參數(shù)等信息，采用山東省地震工程研究院的最新研究成果“山東省1∶50萬(wàn)地震構(gòu)造圖”中的結(jié)論。脆性層深度通過(guò)區(qū)域小震精定位深度的統(tǒng)計(jì)大致確定。這樣，我們可以利用GPS速度場(chǎng)資料作為約束，更好地反演獲得斷層的蠕滑深度及其速率，為地震中長(zhǎng)期預(yù)測(cè)提供支持。

7.區(qū)域地震預(yù)測(cè)的動(dòng)力學(xué)分析

在獲取高精度、高密度三維GPS地殼形變觀測(cè)資料的基礎(chǔ)上，將全球參考框架的站速度進(jìn)行內(nèi)插并且轉(zhuǎn)換為區(qū)域無(wú)旋轉(zhuǎn)基準(zhǔn)，計(jì)算應(yīng)變場(chǎng)（主應(yīng)變、剪應(yīng)變、面應(yīng)變等）；結(jié)合實(shí)際的震情形勢(shì)，建立了區(qū)域構(gòu)造運(yùn)動(dòng)演化過(guò)程與地震活動(dòng)特征之間的相關(guān)關(guān)系，并通過(guò)定量研究孕震構(gòu)造的形變態(tài)勢(shì)、斷裂運(yùn)動(dòng)轉(zhuǎn)換方式、斷層分段閉鎖程度、深部滑動(dòng)速率和應(yīng)力-應(yīng)變積累水平等多項(xiàng)研究結(jié)果的綜合分析，力爭(zhēng)發(fā)現(xiàn)其動(dòng)力學(xué)機(jī)制，為區(qū)域動(dòng)力學(xué)預(yù)報(bào)提供幫助。圖6是以研究區(qū)ML0以上地震擬合建立的沂沭斷裂帶兩側(cè)地區(qū)總的能量釋放（折算為震級(jí)，時(shí)間窗長(zhǎng)為6個(gè)月）與沂沭斷裂帶兩側(cè)塊體相對(duì)運(yùn)動(dòng)速率V的線性關(guān)系，可以為區(qū)域動(dòng)力學(xué)預(yù)報(bào)提供技術(shù)支持。

成果總結(jié)與應(yīng)用

本項(xiàng)目收集了研究區(qū)內(nèi)112個(gè)連續(xù)GPS觀測(cè)站、88個(gè)非連續(xù)GPS觀測(cè)站、20條跨斷層水準(zhǔn)測(cè)線、154個(gè)流動(dòng)重力觀測(cè)點(diǎn)數(shù)據(jù)及地震目錄、地質(zhì)構(gòu)造、斷層模型數(shù)據(jù)。在此基礎(chǔ)上，解決了項(xiàng)目技術(shù)路線中的幾個(gè)關(guān)鍵問(wèn)題，表現(xiàn)在：1）通過(guò)GPS數(shù)據(jù)篩選、高精度數(shù)據(jù)解算、序列分析及非構(gòu)造形變干擾剔除等過(guò)程，獲得了研究區(qū)高精度、高時(shí)-空分辨率的地殼形變場(chǎng)。2）利用GPS垂直向及水準(zhǔn)資料研究了目標(biāo)區(qū)垂直形變特征及其與蒼山地震、日本地震及地形構(gòu)造的關(guān)系，在此基礎(chǔ)上建立了三維速度場(chǎng)。3）以日本地震為時(shí)間節(jié)點(diǎn)，研究了地震前后區(qū)域應(yīng)變積累水平、沂沭斷裂帶分段形變特征等參數(shù)，為區(qū)域地震危險(xiǎn)性判定提供了定量依據(jù)。4）通過(guò)塊體模型建立了區(qū)域構(gòu)造運(yùn)動(dòng)演化過(guò)程與地震活動(dòng)特征之間的相關(guān)關(guān)系，探討了區(qū)域地震活動(dòng)形勢(shì)的動(dòng)力學(xué)機(jī)制，結(jié)合重力觀測(cè)結(jié)果為區(qū)域動(dòng)力學(xué)預(yù)測(cè)提供幫助。研究成果及時(shí)應(yīng)用于山東省地震趨勢(shì)會(huì)商中，為區(qū)域地震趨勢(shì)研判提供了重要技術(shù)支撐。