大華北地區(qū)地下流體資料趨勢(shì)性轉(zhuǎn)折時(shí)空演化特征分析

劉冬英，孫小龍，邵志剛

(1.中國地震局地殼應(yīng)力研究所，北京　100085；2.中國地震局地震預(yù)測(cè)研究所，北京　100026)

大華北地區(qū)地下流體資料趨勢(shì)性轉(zhuǎn)折時(shí)空演化特征分析

劉冬英1，孫小龍1，邵志剛2

(1.中國地震局地殼應(yīng)力研究所，北京100085；2.中國地震局地震預(yù)測(cè)研究所，北京100026)

摘要：對(duì)大華北地區(qū)地下流體資料進(jìn)行梳理，篩選連續(xù)性好、相對(duì)穩(wěn)定且存在趨勢(shì)轉(zhuǎn)折的資料進(jìn)行趨勢(shì)轉(zhuǎn)折變化分析，根據(jù)轉(zhuǎn)折時(shí)間點(diǎn)和形態(tài)進(jìn)行時(shí)空特征分析，并提出了趨勢(shì)轉(zhuǎn)折的定量分析方法。結(jié)果顯示，大華北地區(qū)流體趨勢(shì)轉(zhuǎn)折的測(cè)項(xiàng)數(shù)在2006年和2011年出現(xiàn)明顯增多的現(xiàn)象，此現(xiàn)象與2006年文安5.1級(jí)地震和2011年日本9.0級(jí)地震的影響有關(guān)。

關(guān)鍵詞：大華北；地下流體；趨勢(shì)性轉(zhuǎn)折；時(shí)空特征

0引言

地下流體是地殼中十分活躍而又重要的組成部分。由于地下流體的普遍性、流動(dòng)性、靈敏性以及物質(zhì)和能量的可交換性,使其成為重要的前兆觀測(cè)手段之一。大量觀測(cè)資料和研究表明,在地震前地下流體出現(xiàn)十分豐富的趨勢(shì)性異常。陸明勇等[1]分析了大華北地區(qū)自1976年以來MS≥6.0地震前地下流體長趨勢(shì)變化的形態(tài)、演化特征，認(rèn)為地下流體的長趨勢(shì)變化出現(xiàn)時(shí)間早、范圍廣，具有協(xié)調(diào)性、重復(fù)性、遷移性等特征，主要沿構(gòu)造帶遷移，遷移方向?yàn)橛赏鈬騼?nèi)遷移或內(nèi)部向外遷移；形態(tài)主要為趨勢(shì)轉(zhuǎn)折、破年變、加速等，其中絕大多數(shù)為趨勢(shì)轉(zhuǎn)折變化。在相關(guān)的孕震理論模型中，也有地震理論[2]提到：當(dāng)應(yīng)力應(yīng)變積累到一定程度，這種趨勢(shì)性異常表現(xiàn)為趨勢(shì)性轉(zhuǎn)折，這些趨勢(shì)性轉(zhuǎn)折對(duì)于地震預(yù)測(cè)具有中期尺度的指示意義。在一些強(qiáng)震之前確實(shí)也觀測(cè)到這種趨勢(shì)性轉(zhuǎn)折，例如1964年日本新潟7.5級(jí)地震前4年，觀測(cè)到地形變(長水準(zhǔn))觀測(cè)曲線群體性趨勢(shì)性轉(zhuǎn)折；焦青等[3]總結(jié)出，文安地震前,在華北北部、特別是首都圈地區(qū)前兆觀測(cè)資料出現(xiàn)了大范圍、同步轉(zhuǎn)折變化。大量震例表明，地下流體轉(zhuǎn)折是大華北地區(qū)最重要的地下流體前兆現(xiàn)象之一。張國民等[4]提出的構(gòu)造塊體成組孕震模型表明，成組強(qiáng)震之間的相互關(guān)聯(lián)性會(huì)導(dǎo)致構(gòu)造應(yīng)力場(chǎng)整體增強(qiáng)。同一構(gòu)造塊體會(huì)出現(xiàn)源兆和場(chǎng)兆等2種前兆變化，源兆直接反映震源發(fā)展過程，場(chǎng)兆是由構(gòu)造應(yīng)力場(chǎng)整體增強(qiáng)所造成的，所以大范圍的構(gòu)造應(yīng)力調(diào)整時(shí)，會(huì)出現(xiàn)大范圍的前兆變化。根據(jù)震例分析，遠(yuǎn)場(chǎng)前兆異常的數(shù)量與強(qiáng)震震級(jí)不存在正比關(guān)系，但同一構(gòu)造區(qū)前兆異常數(shù)量相對(duì)較多的時(shí)段與強(qiáng)震連發(fā)的時(shí)段相對(duì)應(yīng)。梅世蓉等[5]、車用太等[6]、杜方等[7]對(duì)大量的觀測(cè)資料研究后發(fā)現(xiàn)，地震前特別是強(qiáng)震前地下流體觀測(cè)資料經(jīng)常出現(xiàn)長趨勢(shì)變化，并且分布范圍廣、出現(xiàn)時(shí)間早。如1976年7月28日唐山7.8級(jí)大地震前，該地區(qū)地下流體觀測(cè)資料在以往多年基本平穩(wěn)的背景下，從1973年初開始出現(xiàn)大范圍的趨勢(shì)性轉(zhuǎn)折變化[8]。也有研究顯示，地震后研究區(qū)異常增多，可能與區(qū)域應(yīng)力的充分釋放有關(guān)[9]。

本文通過收集、整理、分析大華北地區(qū)的地下流體觀測(cè)資料，確定各個(gè)觀測(cè)曲線的轉(zhuǎn)折點(diǎn)，統(tǒng)計(jì)出現(xiàn)轉(zhuǎn)折的觀測(cè)點(diǎn)的數(shù)量、趨勢(shì)變化持續(xù)時(shí)間、變化形態(tài)及幅度，進(jìn)而討論趨勢(shì)性轉(zhuǎn)折測(cè)點(diǎn)的時(shí)空演化特征。

1資料篩選

以國家地震前兆臺(tái)網(wǎng)中心數(shù)據(jù)庫中的流體資料為基礎(chǔ)，本文收集了大華北地區(qū)連續(xù)性好、相對(duì)穩(wěn)定且存在趨勢(shì)轉(zhuǎn)折的地下流體資料72條觀測(cè)數(shù)據(jù),占目前大華北地區(qū)流體觀測(cè)項(xiàng)數(shù)的比例為5%。比例較小的原因主要是剔除了存在明顯干擾因素的數(shù)據(jù)，如受華北地區(qū)地下水開采影響較大的部分水位數(shù)據(jù)。72條觀測(cè)數(shù)據(jù)中，水位數(shù)據(jù)40項(xiàng)、水溫?cái)?shù)據(jù)29項(xiàng)、其它化學(xué)及離子數(shù)據(jù)3項(xiàng)。化學(xué)測(cè)項(xiàng)相對(duì)較少，主要因?yàn)榛瘜W(xué)量異常多屬于閾值類，存在趨勢(shì)轉(zhuǎn)折現(xiàn)象的異常較少。圖1為收集整理的地下流體觀測(cè)資料的觀測(cè)時(shí)間進(jìn)程圖。從圖中可以看出，篩選出存在趨勢(shì)轉(zhuǎn)折的數(shù)據(jù)其觀測(cè)時(shí)間基本都在5年以上，其中數(shù)據(jù)觀測(cè)時(shí)間超過10年的數(shù)據(jù)約有30項(xiàng)；另外，觀測(cè)資料數(shù)量的增加存在2個(gè)明顯的時(shí)間點(diǎn)，分別為2002年和2007年。眾所周知，這2個(gè)時(shí)間點(diǎn)分別是中國地震局前兆觀測(cè)臺(tái)網(wǎng)進(jìn)行“九五”和“十五”觀測(cè)技術(shù)數(shù)字化改造的時(shí)段。在這2個(gè)時(shí)間點(diǎn)，新增了部分觀測(cè)點(diǎn)或測(cè)項(xiàng)。

圖1　地下流體資料觀測(cè)時(shí)間進(jìn)程圖

2流體觀測(cè)資料趨勢(shì)轉(zhuǎn)折分析方法

本文采用的流體觀測(cè)資料趨勢(shì)轉(zhuǎn)折分析方法類似于司學(xué)蕓等[10]所提出的線性擬合法，但在其基礎(chǔ)上做了適當(dāng)?shù)母倪M(jìn)，即：在趨勢(shì)轉(zhuǎn)折時(shí)間點(diǎn)提取時(shí)采用線性擬合法，但在趨勢(shì)幅度提取時(shí)，并不是采用斜率的對(duì)比，而是將斜率轉(zhuǎn)化為轉(zhuǎn)折角度，進(jìn)而來探討趨勢(shì)轉(zhuǎn)折的幅度大小、形態(tài)等。

2.1轉(zhuǎn)折點(diǎn)的提取

趨勢(shì)變化異常分析時(shí)，為了更精確地確定趨勢(shì)轉(zhuǎn)折時(shí)間點(diǎn)，采用滑動(dòng)計(jì)算殘差的方法：首先，以數(shù)據(jù)起始時(shí)間和第2個(gè)預(yù)估轉(zhuǎn)折時(shí)間點(diǎn)為數(shù)據(jù)控制始末點(diǎn)，擬定其中間任意時(shí)刻為第1個(gè)趨勢(shì)轉(zhuǎn)折時(shí)間點(diǎn)，利用線性擬合法可得到2條趨勢(shì)變化線，分別計(jì)算由每一個(gè)擬定轉(zhuǎn)折點(diǎn)所得到的趨勢(shì)線與原始曲線的擬合殘差均值，擬合殘差均值最小者即為最佳趨勢(shì)轉(zhuǎn)折點(diǎn),得到第1個(gè)趨勢(shì)轉(zhuǎn)折時(shí)間點(diǎn)的精確值。

同樣，利用第1個(gè)趨勢(shì)轉(zhuǎn)折時(shí)間點(diǎn)的精確值和第3個(gè)預(yù)估趨勢(shì)轉(zhuǎn)折時(shí)間點(diǎn)，滑動(dòng)求得第2個(gè)趨勢(shì)轉(zhuǎn)折時(shí)間點(diǎn)的精確值。以此類推，可搜索得到每一個(gè)趨勢(shì)轉(zhuǎn)折時(shí)間點(diǎn)，并確定其精確值。

利用以上方法，可得到每一條觀測(cè)曲線不同時(shí)段內(nèi)的趨勢(shì)轉(zhuǎn)折線。

2.2趨勢(shì)性轉(zhuǎn)折的歸一化

為了分析對(duì)比同一區(qū)域內(nèi)不同觀測(cè)點(diǎn)或不同觀測(cè)項(xiàng)目趨勢(shì)變化之間的差異性，我們對(duì)各觀測(cè)曲線進(jìn)行時(shí)間和幅度上的歸一化，使各觀測(cè)曲線在相同的時(shí)間尺度和變化幅度上具有可比性，進(jìn)而分析同一地區(qū)各觀測(cè)點(diǎn)趨勢(shì)轉(zhuǎn)折的群體性特征。

為了更好地從時(shí)間和空間上對(duì)比分析趨勢(shì)轉(zhuǎn)折異常的持續(xù)時(shí)間和異常幅度，我們首先對(duì)觀測(cè)資料進(jìn)行了坐標(biāo)歸一化處理。處理后的觀測(cè)數(shù)據(jù)，其時(shí)間數(shù)據(jù)和觀測(cè)數(shù)據(jù)均分布在0～1之間(圖2)。

圖2　大灰廠水位趨勢(shì)轉(zhuǎn)折異常識(shí)別圖

2.3趨勢(shì)轉(zhuǎn)折幅度的提取

在逐一搜索出最佳轉(zhuǎn)折時(shí)間點(diǎn)和轉(zhuǎn)折速率后，為了更好地對(duì)比不同觀測(cè)點(diǎn)之間的異常幅度，采用了轉(zhuǎn)折角度來定量分析其異常幅度。如果觀測(cè)曲線在出現(xiàn)轉(zhuǎn)折前的擬合斜率為Ka,轉(zhuǎn)折后的擬合斜率為Kb,那么其轉(zhuǎn)折角度θ=arctan(Kb)-arctan(Ka)。以河北無極靜水位觀測(cè)數(shù)據(jù)為例 (圖3)，如果觀測(cè)曲線出現(xiàn)順時(shí)針轉(zhuǎn)折，那么其轉(zhuǎn)折形態(tài)應(yīng)為轉(zhuǎn)折下降型，其轉(zhuǎn)折角度為負(fù)值；如果觀測(cè)曲線出現(xiàn)逆時(shí)針轉(zhuǎn)折，那么其轉(zhuǎn)折形態(tài)應(yīng)為轉(zhuǎn)折上升型，其轉(zhuǎn)折角度為正值。

圖3　趨勢(shì)轉(zhuǎn)折幅度識(shí)別方法示意圖

3大華北地區(qū)流體趨勢(shì)性轉(zhuǎn)折的時(shí)空特征

圖4為大華北地區(qū)流體新增趨勢(shì)轉(zhuǎn)折異常數(shù)量統(tǒng)計(jì)圖，從圖4可以看出：

1)2009年后的趨勢(shì)轉(zhuǎn)折異常數(shù)量明顯增多，這是由于流體臺(tái)網(wǎng)數(shù)字化改造 ，新增觀測(cè)點(diǎn)或測(cè)項(xiàng)，流體觀測(cè)資料數(shù)量增多所致；

2)新增趨勢(shì)轉(zhuǎn)折異常數(shù)量在2006年和2011年明顯多于臨近的其它年份，可能與2006年7月4日文安5.1級(jí)地震和2011年3月11日日本9.0級(jí)地震的影響有關(guān)。

圖4　大華北地區(qū)流體觀測(cè)資料新增趨勢(shì)轉(zhuǎn)折異常數(shù)量統(tǒng)計(jì)

圖5為大華北地區(qū)2001—2012年流體新增趨勢(shì)轉(zhuǎn)折異常時(shí)空分布圖，從圖5可以看出，2001—2005年期間，研究區(qū)域內(nèi)新增趨勢(shì)轉(zhuǎn)折點(diǎn)數(shù)量均維持在4項(xiàng)左右，2006年突增至9項(xiàng)，2007—2008年基本恢復(fù)至正常水平(5項(xiàng))，這應(yīng)與2006年文安5.1級(jí)地震的發(fā)生有一定關(guān)系；2009年后數(shù)量明顯增多，正如前文所述，與流體臺(tái)網(wǎng)數(shù)字化改造有關(guān)；2011年數(shù)量顯明高于其前后2年，且增加的趨勢(shì)轉(zhuǎn)折異常點(diǎn)多集中分布于東北地區(qū)，這應(yīng)與2011年日本9.0級(jí)地震的發(fā)生有一定關(guān)系。

2006年大華北地區(qū)流體趨勢(shì)轉(zhuǎn)折異常點(diǎn)數(shù)量為9項(xiàng)，明顯高于其當(dāng)時(shí)的平均數(shù)量。從趨勢(shì)轉(zhuǎn)折異常出現(xiàn)的時(shí)間點(diǎn)來看，這9項(xiàng)異常中有4項(xiàng)異常是出現(xiàn)在2006年7月4日文安5.1級(jí)地震發(fā)生前，有5項(xiàng)是出現(xiàn)在地震之后；2011年大華北地區(qū)流體趨勢(shì)轉(zhuǎn)折異常點(diǎn)數(shù)量為25項(xiàng)，其中有5項(xiàng)異常是出現(xiàn)在2011年3月11日日本9.0級(jí)地震發(fā)生前，有20項(xiàng)是出現(xiàn)在地震之后。因此，從趨勢(shì)轉(zhuǎn)折異常點(diǎn)數(shù)量來看，2006年文安5.1級(jí)地震和2011年日本9.0級(jí)地震發(fā)生的年份，大華北地區(qū)趨勢(shì)轉(zhuǎn)折異常數(shù)量顯明高于其它年份，但從異常出現(xiàn)的時(shí)間點(diǎn)來看，趨勢(shì)轉(zhuǎn)折異常多出現(xiàn)于地震發(fā)生之后。分析認(rèn)為，大華北地區(qū)的趨勢(shì)轉(zhuǎn)折數(shù)量增多，是由于地震發(fā)生后區(qū)域應(yīng)力場(chǎng)調(diào)整所致[11]。

圖6為大華北地區(qū)當(dāng)前各觀測(cè)點(diǎn)趨勢(shì)轉(zhuǎn)折異常的時(shí)空演化特征圖。從圖6a圖可以看出，當(dāng)前大華北地區(qū)以趨勢(shì)轉(zhuǎn)折下降異常為主，且多集中分布于張渤地震帶和郯廬地震帶北段；從圖6b圖可以看出，當(dāng)前大華北地區(qū)異常幅度較大、持續(xù)時(shí)間較長的趨勢(shì)轉(zhuǎn)折異常主要集中分布在首都圈及環(huán)渤海地區(qū)和山西地震帶北端，而東北地區(qū)的趨勢(shì)轉(zhuǎn)折多以持續(xù)時(shí)間短、異常幅度大為主要特征。這些趨勢(shì)異常集中分布的區(qū)域，可能為當(dāng)前或今后中強(qiáng)地震發(fā)生的主體活動(dòng)區(qū)域。

4結(jié)論與討論

本文基于國家地震前兆臺(tái)網(wǎng)中心數(shù)據(jù)庫中的流體資料，篩選了大華北地區(qū)連續(xù)性好、相對(duì)穩(wěn)定且存在趨勢(shì)轉(zhuǎn)折的地下流體資料72條觀測(cè)數(shù)據(jù)，利用線性擬合法定量提取了各觀測(cè)資料中的趨勢(shì)轉(zhuǎn)折時(shí)間點(diǎn)和轉(zhuǎn)折幅度，在此基礎(chǔ)上統(tǒng)計(jì)分析了研究區(qū)域內(nèi)的流體趨勢(shì)轉(zhuǎn)折異常點(diǎn)時(shí)空分布特征，得出以下認(rèn)識(shí)：

1)從趨勢(shì)轉(zhuǎn)折異常點(diǎn)數(shù)量來看，大華北地區(qū)2006年和2011年新增趨勢(shì)轉(zhuǎn)折異常數(shù)量明顯多于臨近的其它年份，可能與2006年7月4日文安5.1級(jí)地震和2011年3月11日日本9.0級(jí)地震的影響有關(guān)；

圖5　大華北地區(qū)2001—2012年流體新增趨勢(shì)轉(zhuǎn)折異常時(shí)空分布圖

注：a　圖中藍(lán)色圓圈表示趨勢(shì)轉(zhuǎn)折下降點(diǎn)，紅色圓圈表示趨勢(shì)轉(zhuǎn)折上升點(diǎn)，圓圈大小表示其趨勢(shì)轉(zhuǎn)折幅度的大??；b　圖中圓圈大小表示趨勢(shì)轉(zhuǎn)折幅度的大小，顏色表示異常持續(xù)時(shí)間的長短。圖6　大華北地區(qū)流體趨勢(shì)轉(zhuǎn)折異常時(shí)空演化特征圖

2)從異常出現(xiàn)的時(shí)間點(diǎn)來看，2006年和2011年大華北地區(qū)新增趨勢(shì)轉(zhuǎn)折異常點(diǎn)并非出現(xiàn)于地震發(fā)生之前，而是多出現(xiàn)于地震發(fā)生之后，可能與地震發(fā)生后區(qū)域應(yīng)力場(chǎng)調(diào)整有關(guān)，也即地震的發(fā)生引起了區(qū)域應(yīng)力場(chǎng)的調(diào)整，進(jìn)而引起了地下流體觀測(cè)資料的趨勢(shì)變化；

3)大華北地區(qū)異常幅度較大、持續(xù)時(shí)間較長的趨勢(shì)轉(zhuǎn)折異常主要集中分布在首都圈及環(huán)渤海地區(qū)和山西地震帶北端，這些趨勢(shì)異常集中分布的區(qū)域，可能為當(dāng)前或今后中強(qiáng)地震發(fā)生的主體活動(dòng)區(qū)域；

4)日本9.0級(jí)地震之后，東北地區(qū)的趨勢(shì)轉(zhuǎn)折多以持續(xù)時(shí)間短、異常幅度大為主要特征，在這些趨勢(shì)異常集中區(qū)于2013年10月31日、11月23日相繼發(fā)生了前郭5.5級(jí)和5.8級(jí)地震，表明這些趨勢(shì)異常對(duì)區(qū)域中強(qiáng)地震的發(fā)生具有一定的指示意義。

基于司學(xué)蕓等[10]的分析方法，本文在利用線性擬合法提取趨勢(shì)異常變化幅度過程中，采用轉(zhuǎn)折角度的方法來表征趨勢(shì)轉(zhuǎn)折幅度，比原有采用斜率差的方法更具一定的科學(xué)性和準(zhǔn)確性。由數(shù)學(xué)三角函數(shù)知識(shí)可知，當(dāng)觀測(cè)曲線與時(shí)間軸的夾角小于45°時(shí)，其斜率值小于1，當(dāng)夾角大于45°時(shí)，其斜率值大于1。如果2組觀測(cè)曲線轉(zhuǎn)折相同的角度，但若轉(zhuǎn)折前后的角度在45°～90°之間，其斜率差值會(huì)遠(yuǎn)遠(yuǎn)大于轉(zhuǎn)折前后角度在0°～45°之間的斜率差值，這就會(huì)造成轉(zhuǎn)折幅度計(jì)算的不準(zhǔn)確性。

分析結(jié)果顯示，大華北地區(qū)趨勢(shì)轉(zhuǎn)折異常數(shù)量于2011年日本9.0級(jí)地震后明顯增多，2012年異常數(shù)量雖然有所減少，但減小的異常點(diǎn)均分布于首都圈及其附近區(qū)域，東北地區(qū)減少不明顯，并且東北地區(qū)的異常點(diǎn)有向黑吉蒙交界地區(qū)遷移的特征(圖5)。分析認(rèn)為，日本地震的發(fā)生對(duì)中國大陸東北地區(qū)的影響最為顯著，這種影響一方面會(huì)引起區(qū)域應(yīng)力場(chǎng)的調(diào)整，反映在流體觀測(cè)資料上，就表現(xiàn)為流體趨勢(shì)性、群體性的異常變化；另一方面，這種區(qū)域應(yīng)力場(chǎng)的調(diào)整有可能會(huì)誘發(fā)受影響區(qū)域內(nèi)地震的發(fā)生，如前郭2013年10月31日5.5級(jí)和2013年11月23日5.8級(jí)地震的發(fā)生。

因此，東北地區(qū)趨勢(shì)轉(zhuǎn)折異常數(shù)量在2011年日本9.0級(jí)地震后明顯增多，反映了東北地區(qū)區(qū)域構(gòu)造應(yīng)力場(chǎng)在日本9.0級(jí)地震后出現(xiàn)了明顯的調(diào)整作用。與本文分析結(jié)果類似，GPS觀測(cè)資料的分析結(jié)果也顯示出，2011年日本9.0級(jí)地震的發(fā)生導(dǎo)致了東北和華北地區(qū)一系列NNE走向的斷裂產(chǎn)生不同程度的張性應(yīng)變(圖7)，且震后受影響較大的區(qū)域主要為東北地區(qū)[11]。

a　1999—2007年　　b　2011年8月—2012年8月圖7　日本9.0級(jí)地震前后華北與東北地區(qū)GPS主應(yīng)變率場(chǎng)

在地下流體資料的使用中，由于趨勢(shì)變化持續(xù)時(shí)間長，在短時(shí)間內(nèi)對(duì)觀測(cè)點(diǎn)附近的中強(qiáng)以上地震的反映能力并不突出，因此在震情年度會(huì)商及短臨跟蹤中發(fā)揮的作用有限。而大范圍、群體性的趨勢(shì)變化可能反映了一定區(qū)域的應(yīng)力變化特征，從而能夠?yàn)榇笮蝿?shì)、中長期預(yù)測(cè)提供一定的參考。

參考文獻(xiàn)：

[1]陸明勇, 范雪芳, 周偉, 等. 華北強(qiáng)震前地下流體長趨勢(shì)變化特征及其產(chǎn)生機(jī)理的研究[J]. 西北地震學(xué)報(bào), 2010, 32(2): 129-138.

[2]米亞奇金B(yǎng) N,布霍斯W F.地震理論與實(shí)驗(yàn)譯文集[M]. 國家地震局地球物理研究所第三研究室譯. 北京:地震出版社, 1979.

[3]焦青, 劉耀煒, 楊選輝, 等. 文安5.1級(jí)地震前兆變化特征分析[J]. 大地測(cè)量與地球動(dòng)力學(xué), 2007, 27(6): 96-99.

[4]張國民, 耿魯明, 張永仙, 等. 構(gòu)造塊體的成組孕震模型和前兆場(chǎng)某些特征的分析[J]. 地震學(xué)報(bào), 1995, 17(1): 1-10.

[5]梅世蓉, 馮德益, 張國民, 等. 中國地震預(yù)報(bào)概論[M]: 北京: 地震出版社, 1993.

[6]車用太, 劉五洲, 魚金子, 等. 板內(nèi)強(qiáng)震的中地殼硬夾層孕震與流體促震假設(shè)[J]. 地震學(xué)報(bào), 2000, 22(1): 93-101.

[7]杜方, 吳江.昆侖山口西8.1級(jí)地震引起的地下流體遠(yuǎn)場(chǎng)響應(yīng)研究[J]. 地震研究, 2008, 31(3): 197-202.

[8]陸明勇, 劉耀煒, 范雪芳, 等. 地下流體長趨勢(shì)異常變化與強(qiáng)震預(yù)測(cè)的初步研究[J].地震研究, 2009, 32(4): 323-332.

[9]余丹, 魏星, 紀(jì)壽文. 文安5.1級(jí)地震前后地殼形變特征研究[J]. 華北地震科學(xué), 2014, 32(2): 62-67, 72.

[10]司學(xué)蕓, 孫小龍, 邵志剛, 等. 南北地震帶流體資料趨勢(shì)性轉(zhuǎn)折與強(qiáng)震的關(guān)系[J].中國地震, 2013, 29(1): 148-156.

[11]王敏, 李強(qiáng), 王凡, 等. 全球定位系統(tǒng)測(cè)定的2011年日本宮城MW9.0級(jí)地震遠(yuǎn)場(chǎng)同震位移[J].科學(xué)通報(bào), 2011, 56(20): 1593-1596.

[12]孫麗娟, 冀林旺, 金鎮(zhèn)洪, 等. 山龍峪井水位異常與地震關(guān)系初探[J].東北地震研究, 2004, 20(1): 42-47.

Temporal-spatial Evolution Characteristics of Trend Turning Data of Underground Fluid in the Great North China Area

LIU Dong-ying1, SUN Xiao-long1, SHAO Zhi-gang2

(1.The Institute of Crustal Dynamic, Beijing 100085, China;2. Institute of Earthquake Science, CEA, Beijing 100026, China)

Abstract:We select out continue and relatively stable underground fluid data that has trend turning in great North China Area to analyze its changing characteristics .According to the turning point and form, we find out the temporal-spatial characteristics and then put forward a quantitative method for tread turning analysis. The results show that the fluid trend turning changes increased rapidly in the great North China Area in 2006 and 2011, this phenomenon may be related to with the 2006 Wen’an MS5.1 earthquake and the 2011 Japan MW9.0 earthquake respectively.

Key words:great North China Area; underground fluid; trend turning; the temporal-spatial characteristics

收稿日期:2016-01-15

基金項(xiàng)目：中國地震局震情跟蹤合同制定向任務(wù)(2016020306)

作者簡(jiǎn)介：劉冬英(1971—)，女，北京人，高級(jí)工程師，主要從事地震活動(dòng)性及地下流體研究．E-mail:ldy7117@163.com

中圖分類號(hào):P315. 723

文獻(xiàn)標(biāo)志碼:A

文章編號(hào):1003-1375(2016)02-0001-06

doi:10.3969/j.issn.1003-1375.2016.02.001

劉冬英 ,孫小龍 ,邵志剛．大華北地區(qū)地下流體資料趨勢(shì)性轉(zhuǎn)折時(shí)空演化特征分析[J]．華北地震科學(xué)，2016，34(2):1-6.