劉琦　李臘月　胡樂銀

摘要：為掌握北京及周邊地區(qū)跨斷層形變觀測場地目前的效能情況，從而實(shí)現(xiàn)準(zhǔn)確研判該區(qū)域地震危險(xiǎn)性，降低未來地震可能造成的災(zāi)害風(fēng)險(xiǎn)，綜合考慮跨斷層場地監(jiān)測條件和地震預(yù)測效能兩方面設(shè)計(jì)了觀測站網(wǎng)評(píng)估方案，分別從場地監(jiān)測環(huán)境、監(jiān)測有效性、地下水位變化及地表沉降影響、場地預(yù)測配套性、最優(yōu)預(yù)測指標(biāo)效能和區(qū)域變形配套性等方面開展研究，給出了北京及周邊地區(qū)跨斷層形變站網(wǎng)的綜合評(píng)估結(jié)果及后續(xù)站網(wǎng)優(yōu)化建議。結(jié)果表明：①目前北京及周邊地區(qū)的跨斷層場地對區(qū)域內(nèi)大多數(shù)主要斷層具有較好的監(jiān)測能力和地震預(yù)測效能，但八寶山—黃莊—高麗營斷裂上多數(shù)觀測場地監(jiān)測條件相對較差，預(yù)測效能普遍偏低，建議選擇該斷裂的其它位置補(bǔ)充新建若干觀測場地。②除了斷裂構(gòu)造，區(qū)域地質(zhì)構(gòu)造、水文地質(zhì)條件等也對觀測場地效能有重要影響，后續(xù)新場地選址時(shí)需要重點(diǎn)考慮這些因素。

關(guān)鍵詞：跨斷層觀測；地震預(yù)測效能；監(jiān)測條件；地下水位；構(gòu)造環(huán)境；北京

中圖分類號(hào)：P315.727文獻(xiàn)標(biāo)識(shí)碼：A文章編號(hào)：1000-0666（2023）04-0491-12

doi：10.20015/j.cnki.ISSN1000-0666.2023.0055

0引言

北京是中國的政治、文化、國際交往和科技創(chuàng)新中心，經(jīng)濟(jì)發(fā)達(dá)、人口稠密，同時(shí)北京及周邊地區(qū)也是強(qiáng)震活動(dòng)較為頻繁的區(qū)域之一，歷史上曾發(fā)生過1679年三河—平谷8級(jí)地震及多次6級(jí)以上地震（蔣銘，馬宗晉，1985；張國民等，1988），特別是新中國成立以來該區(qū)域發(fā)生了1978年唐山7.8級(jí)和灤縣7.1級(jí)地震，造成了24萬多人死亡和巨大經(jīng)濟(jì)損失。因此，針對北京及周邊地區(qū)具備強(qiáng)震孕震能力的主要斷層，開展有效地?cái)鄬舆\(yùn)動(dòng)監(jiān)測，捕捉地震孕育過程中可能的異?，F(xiàn)象，對準(zhǔn)確研判區(qū)域的地震危險(xiǎn)性具有重要意義，針對現(xiàn)有的斷層形變觀測開展監(jiān)測、預(yù)測效能評(píng)估是尤為重要的環(huán)節(jié)。

早在20世紀(jì)60年代前后，北京及周邊地區(qū)已開始陸續(xù)布設(shè)跨斷層的形變觀測場地，經(jīng)過數(shù)十年不斷的發(fā)展，目前已經(jīng)有20多個(gè)場地，積累了數(shù)十年的觀測數(shù)據(jù)，針對這些場地的監(jiān)測環(huán)境以及實(shí)際的地震預(yù)測效能，前人已開展了大量的研究（陸明勇等，2011，2013，2018），中國地震局也曾組織相關(guān)單位進(jìn)行評(píng)估分析和站網(wǎng)改造方案設(shè)計(jì)。由于觀測環(huán)境處于動(dòng)態(tài)變化之中，特別是隨著城市化發(fā)展的穩(wěn)步推進(jìn)，以往研究中沒有充分考慮的一些影響因素也變得突出，例如華北地區(qū)相對匱乏的水量補(bǔ)給和不斷提升的用水需求之間的矛盾劇增，造成了超量開采地下水的問題。受此影響，華北平原成為中國沉降速率和影響面積最大的地區(qū)，局部地區(qū)的累積沉降量甚至超過了3.4 m，北京地區(qū)的最大累積沉積量超過1.4 m（何慶成等，2006；楊艷，2015；郭海朋等，2017），因此在分析北京及周邊地區(qū)跨斷層觀測條件時(shí)，需要考慮這種地下水位變化及其引發(fā)的區(qū)域性沉降所帶來的影響。地下水儲(chǔ)量、水位變化等的動(dòng)態(tài)監(jiān)測以及區(qū)域水準(zhǔn)、GNSS、InSAR等觀測資料的積累為深入開展相關(guān)影響的分析提供了重要基礎(chǔ)（梁芳等，2013；孫啟凱等，2017；Zhang et al，2019；史珉等，2021；Bai et al，2022；Li et al，2022）。

以往的地震異常判識(shí)和預(yù)測效能評(píng)估大多基于定性認(rèn)識(shí)，專家經(jīng)驗(yàn)對結(jié)果有重要影響。近幾年在中國地震局監(jiān)測預(yù)報(bào)司和地殼形變分析預(yù)報(bào)技術(shù)管理組的統(tǒng)一組織下，相關(guān)信息提取、異常確定、效能評(píng)估等的定量方法得到了長遠(yuǎn)發(fā)展（中國地震局監(jiān)測預(yù)報(bào)司，2020），并且已經(jīng)開始在業(yè)務(wù)中進(jìn)行推廣應(yīng)用，為開展跨斷層場地的預(yù)測效能評(píng)估提供了重要可靠的方法基礎(chǔ)。

因此，本文在充分考慮前人已有研究基礎(chǔ)上，針對北京及周邊地區(qū)跨斷層形變觀測站網(wǎng)設(shè)計(jì)了評(píng)估方案，從跨斷層場地的監(jiān)測條件以及地震預(yù)測效能兩方面開展綜合評(píng)估，給出綜合評(píng)估結(jié)果，針對關(guān)鍵性的影響因素進(jìn)行了分析討論，提出了未來北京及周邊地區(qū)跨斷層形變站網(wǎng)的優(yōu)化建議。

北京及周邊地區(qū)的跨斷層形變測量始于20世紀(jì)60年代末，目前，共有跨斷層觀測場地27處，其中流動(dòng)場地25處、定點(diǎn)臺(tái)站2個(gè)，觀測場地遍及NW走向的南口—孫河斷裂、施莊斷裂、淶水—高碑店斷裂以及NE走向延礬盆地北緣斷裂、南口山前斷裂、紫荊關(guān)斷裂、黃莊—高麗營斷裂、夏墊斷裂、唐山斷裂等多條主要活動(dòng)斷裂（圖1）。為給出當(dāng)前北京及周邊地區(qū)跨斷層場地監(jiān)測條件較為客觀的評(píng)估結(jié)果，本文將從場地的監(jiān)測環(huán)境、監(jiān)測有效性、地下水位變化及地表沉降影響等幾個(gè)方面開展評(píng)估分析。

1.1監(jiān)測環(huán)境評(píng)估

參考中國地震局相關(guān)單位前期工作基礎(chǔ)以及前人的研究成果，針對北京及周邊地區(qū)的跨斷層形變場地的監(jiān)測環(huán)境開展定性評(píng)估，主要考慮測樁類型、是否跨過目標(biāo)斷層、周圍是否存在顯著干擾等因素。場地監(jiān)測環(huán)境評(píng)估的具體準(zhǔn)則為：當(dāng)同時(shí)滿足測樁類型為基巖、觀測跨過目標(biāo)斷層、周邊無顯著干擾3個(gè)條件時(shí)，定為A類；當(dāng)上述條件同時(shí)滿足2個(gè)時(shí)，定為B類；當(dāng)上述條件只滿足一個(gè)或都不滿足時(shí)，定為C類，評(píng)估結(jié)果見表1。由于不同研究所考慮的因素存在一定差異，并且近幾年部分場地經(jīng)過了升級(jí)改造，周邊環(huán)境也有一定變化，因此本文與前人的結(jié)果存在部分差異，但總體較為一致。北京及周邊地區(qū)跨斷層場地的監(jiān)測環(huán)境整體較好，其中A類場地14處，B類場地8處，合計(jì)占總場地?cái)?shù)的81.5%。

1.2監(jiān)測有效性評(píng)估

跨斷層形變觀測的目的是監(jiān)測斷層運(yùn)動(dòng)的動(dòng)態(tài)演化，從而為地震預(yù)測提供可靠參考。因此，能否有效反映真實(shí)的斷層運(yùn)動(dòng)，是評(píng)估跨斷層場地監(jiān)測條件的一項(xiàng)重要指標(biāo)。針對跨斷層場地監(jiān)測有效性的評(píng)估，一方面可以直接基于原始觀測曲線定性判斷斷層活動(dòng)特征，另一方面可以基于幾何關(guān)系公式定量計(jì)算斷層活動(dòng)參數(shù)，從而分析斷層活動(dòng)特征（薄萬舉等，1998；李臘月等，2016）。為此，筆者設(shè)定了跨斷層場地監(jiān)測有效性的評(píng)估準(zhǔn)則：從長趨勢背景運(yùn)動(dòng)角度，若原始觀測曲線及計(jì)算的活動(dòng)參數(shù)結(jié)果反映的斷層活動(dòng)方式與地質(zhì)學(xué)給出的斷層背景活動(dòng)一致，則定為A類；若原始觀測曲線及計(jì)算的活動(dòng)參數(shù)結(jié)果反映的斷層活動(dòng)方式有小幅波動(dòng)，但整體上與地質(zhì)學(xué)給出的斷層背景活動(dòng)一致，則定為B類；若原始觀測曲線及計(jì)算的活動(dòng)參數(shù)結(jié)果反映的斷層活動(dòng)方式不明確，或與地質(zhì)學(xué)、其它手段給出的斷層背景活動(dòng)不一致，則定為C類。

針對水準(zhǔn)數(shù)據(jù)，選用起始年的觀測數(shù)據(jù)作為參考基準(zhǔn)，利用每年某個(gè)月份的測值減去基準(zhǔn)年對應(yīng)月份的測值，計(jì)算每年在基準(zhǔn)年基礎(chǔ)上的累積水準(zhǔn)觀測變化量Δh，Δh>0表示正斷層活動(dòng)，上盤下降；Δh<0表示逆斷層活動(dòng)，上盤上升。經(jīng)過該處理，基本可以消除季節(jié)變化等因素引起的年變影響，得到的斷層垂直升降累積量能相對客觀地描述斷層在觀測期內(nèi)的活動(dòng)狀態(tài)。

針對各跨斷層形變場地開展相關(guān)斷層活動(dòng)參數(shù)計(jì)算，結(jié)合已有地質(zhì)研究結(jié)果（高文學(xué)，馬瑾，1993；郝書儉，尤惠川，2001；徐錫偉等，2002；尤惠川等，2002；張世民等，2007；查明等，2011；楊海盟等，2014），依據(jù)上述準(zhǔn)則進(jìn)行跨斷層場地監(jiān)測有效性的評(píng)估，具體評(píng)估結(jié)果見表2，其中A、B類場地共計(jì)14處，占總場地?cái)?shù)的51.9%。密云、墻子路、八寶山、狼山、淶水、南孟、寧河、小寺沖等場地監(jiān)測有效性相對較差，主要和場地監(jiān)測環(huán)境相對較差、所處斷層活動(dòng)性相對較弱等有關(guān)，導(dǎo)致觀測數(shù)據(jù)容易出現(xiàn)與地質(zhì)運(yùn)動(dòng)背景不一致的情況。另外，古北口、鎮(zhèn)羅營、南口、百善、京西實(shí)驗(yàn)場等場地所跨斷裂以走滑運(yùn)動(dòng)為主，但目前僅有水準(zhǔn)觀測，無法有效監(jiān)測斷層的主要運(yùn)動(dòng)方式。

1.3地下水位變化及地表沉降的影響評(píng)估

地下水位下降一方面可以造成區(qū)域局部發(fā)生顯著的地表沉降，進(jìn)而可能對周邊的觀測產(chǎn)生影響；另一方面，水位變化本身也可能對相關(guān)的觀測直接造成影響，因此在分析過程中筆者同時(shí)考慮了上述兩種影響過程。前人研究已經(jīng)給出了1959—2016年研究區(qū)的淺層地下水位累積變化分布結(jié)果（Pang et al，2020；竇甜甜等，2022），在此基礎(chǔ)上，筆者基于該區(qū)域覆蓋相對完整的InSAR監(jiān)測結(jié)果（史珉等，2021），同時(shí)提取了2016—2018年地表沉降速率相對較高（大于50 mm/a）的區(qū)域（圖2）。由于淺層地下水位變化結(jié)果是近60年的累積變化量，地表沉降速率僅為近幾年較短時(shí)間內(nèi)的結(jié)果，因此二者的高值區(qū)在細(xì)節(jié)上存在一定差異，但整體分布具有一致性，顯示了超量開采地下水對區(qū)域地表快速沉降現(xiàn)象的控制作用?？傮w上來看，淺層地下水位下降和地表沉降最為突出的區(qū)域主要集中在研究區(qū)的中部及南部。

根據(jù)空間分布關(guān)系初步分析，密云、德勝口、南口、百善、八寶山、淶水、南孟、唐山等場地位于淺層地下水位累積下降最為顯著的區(qū)域附近，可能會(huì)受其影響，其中百善、南孟等場地還同時(shí)位于近幾年地表沉降較快的區(qū)域附近，受影響程度可能更加嚴(yán)重。在此基礎(chǔ)上，筆者進(jìn)一步分析相關(guān)影響的時(shí)間過程，考慮到地表沉降產(chǎn)生和發(fā)展過程與地下水開采過程基本保持同步且略滯后（何慶成等，2006；楊艷等，2010，2013），因此可以利用地下水儲(chǔ)量時(shí)序變化來表示地下水位、地表沉降的變化過程，并通過與跨斷層觀測數(shù)據(jù)的對比來分析地下水位、地表沉降等因素對跨斷層場地的實(shí)際影響。在初步判定為受影響的場地中，密云、德勝口、南口、百善、八寶山、淶水、唐山等場地的觀測數(shù)據(jù)并未表現(xiàn)出與地下水儲(chǔ)量較為同步的變化過程，即未受其顯著控制（圖3a），僅南孟場地受影響較為嚴(yán)重（圖3b），這可能主要與場地局部地質(zhì)環(huán)境特征以及測樁類型是否為基巖有關(guān)。進(jìn)一步研究發(fā)現(xiàn)，除了南孟場地之外，京西實(shí)驗(yàn)場、大灰廠、上萬、?？谟葓龅氐挠^測背景變化與水儲(chǔ)量的變化一致性也較好（圖3c）。上述場地周邊雖然沒有淺層地下水位顯著下降區(qū)和地表快速沉降區(qū)，但是明顯受到水儲(chǔ)量變化的直接影響，造成這種現(xiàn)象的可能原因在后文做進(jìn)一步討論。綜合上述地下水位變化及地表沉降對跨斷層觀測的實(shí)際影響程度，將南孟、京西實(shí)驗(yàn)場、大灰廠、上萬、?？谟?處場地定為C類，密云、德勝口、南口、百善、八寶山、淶水、唐山7處場地定為B類，其它15處場地定為A類，A類場地?cái)?shù)占總場地?cái)?shù)的55.6%。

2地震預(yù)測效能評(píng)估

2.1基于單測項(xiàng)單方法的預(yù)測效能定量評(píng)估

跨斷層觀測能否監(jiān)測到可靠的地震前兆異常變化是場地評(píng)估中需要重點(diǎn)考慮的因素。歷史震例資料的統(tǒng)計(jì)結(jié)果顯示，一些跨斷層觀測在多次地震前都曾記錄到顯著的異常變化，并且異常主要以趨勢轉(zhuǎn)折、破年變、速率變化、超閾值變化等為主（曹建玲等，2016；陳長云等，2016；李瑞莎等，2016）。因此，針對單個(gè)場地的單個(gè)測項(xiàng)，分別利用趨勢轉(zhuǎn)折、破年變、速率、閾值等方法的不同算法（洪敏等，2019）進(jìn)行相應(yīng)特征信息的提取，并利用R值檢驗(yàn)方法進(jìn)行最優(yōu)預(yù)測指標(biāo)的自動(dòng)確定以及預(yù)測效能的定量評(píng)估（中國地震局監(jiān)測預(yù)報(bào)司，2020）。

R值檢驗(yàn)是目前在地震預(yù)報(bào)領(lǐng)域應(yīng)用最為廣泛的一種效能評(píng)估方法，綜合考慮了有震報(bào)準(zhǔn)率、無震報(bào)準(zhǔn)率、漏報(bào)率、虛報(bào)率等評(píng)估因素，通過不同因素的組合和折中給出相應(yīng)的效能評(píng)估結(jié)果（國家地震局科技監(jiān)測司，1990）。由于上述因素并非完全獨(dú)立，所以R值檢驗(yàn)有多種等價(jià)的表達(dá)形式，常用的一種形式如下：R=有震報(bào)準(zhǔn)率-預(yù)報(bào)占時(shí)率=

（報(bào)對的地震次數(shù)/應(yīng)預(yù)報(bào)的地震總次數(shù)）-（預(yù)報(bào)占用時(shí)間/研究的總時(shí)間）（3）當(dāng)指定了目標(biāo)震級(jí)范圍、空間范圍及預(yù)報(bào)時(shí)窗后，通過調(diào)整不同的異常閾值就可以對應(yīng)不同的有震報(bào)準(zhǔn)率和預(yù)報(bào)占時(shí)率組合，進(jìn)而獲得不同的R值，從中選擇R值最高所對應(yīng)的異常閾值作為該指定預(yù)報(bào)規(guī)則下（特定的震級(jí)、空間、時(shí)間范圍）的相對最優(yōu)閾值。綜合所有預(yù)報(bào)規(guī)則下的評(píng)估結(jié)果，就可以確定最優(yōu)的異常閾值、相應(yīng)的預(yù)報(bào)規(guī)則及預(yù)測效能。

選取了1970—2020年北京及周邊地區(qū)（38.5°～42°N，112°～119°E）5級(jí)以上地震用于效能檢驗(yàn)，共計(jì)16次地震事件。首次為1976年4月6日內(nèi)蒙古和林格爾6.2級(jí)地震，末次為2020年7月12日河北唐山5.1級(jí)地震，1976年7月28日河北唐山7.8級(jí)、灤縣7.1級(jí)以及震后短期內(nèi)5級(jí)以上余震作為1次地震事件，1989年10月18—19日山西大同5級(jí)震群同樣作為1次地震事件。

R值同時(shí)考慮了有震報(bào)準(zhǔn)率和預(yù)報(bào)占時(shí)率，表示扣除了隨機(jī)概率后的預(yù)報(bào)成功率，當(dāng)R=0時(shí)表示預(yù)報(bào)沒有起作用。一般除了考慮R值大小外，還需要考慮其信度。同樣的R值，當(dāng)其基于的地震次數(shù)不一樣時(shí)，信度也不一樣，因此需要考慮R0值，即保證97.5%的置信度的最低R值。當(dāng)R≥R0時(shí)，可以認(rèn)為這個(gè)R值至少有97.5%的置信度（國家地震局科技監(jiān)測司，1990）。筆者依據(jù)通過效能及置信度檢驗(yàn)（R>0且R≥R0）的測項(xiàng)-方法數(shù)來評(píng)估場地預(yù)測配套性，通過數(shù)越多表明預(yù)測配套性越好，越有利于在實(shí)際業(yè)務(wù)應(yīng)用中綜合評(píng)估地震危險(xiǎn)性：通過數(shù)≤4的場地定義其預(yù)測配套性為C類；5≤通過數(shù)≤10的場地定義為B類，通過數(shù)＞10的場地定義為A類。

眾多通過檢驗(yàn)的測項(xiàng)-方法中最優(yōu)預(yù)測指標(biāo)的效能同樣具有重要的評(píng)估意義，這代表了“王牌”手段的預(yù)測能力：所有測項(xiàng)-方法均未通過檢驗(yàn)，或最優(yōu)測項(xiàng)-方法的R值通過檢驗(yàn)，但R≤0.3的場地定義其最優(yōu)預(yù)測指標(biāo)效能為C類；最優(yōu)測項(xiàng)-方法的R值通過檢驗(yàn)，并且0.3＜R≤0.5的場地定義為B類；最優(yōu)測項(xiàng)-方法的R值通過檢驗(yàn)，并且R＞0.5的場地定義為A類。

表3給出了基于單測項(xiàng)單方法的場地預(yù)測效能評(píng)估結(jié)果，場地預(yù)測配套性A類有6處場地，B類5處，合計(jì)占總場地?cái)?shù)40.7%；最優(yōu)預(yù)測指標(biāo)效能A類有8處場地，B類9處，合計(jì)占總場地?cái)?shù)的63.0%，比預(yù)測配套性的比例更高，表明部分場地雖然異常配套性較差，但是具備一些“王牌”的預(yù)測手段，仍然可以在實(shí)際的預(yù)報(bào)業(yè)務(wù)中發(fā)揮重要作用。

2.2基于區(qū)域變形配套性的預(yù)測效能評(píng)估

在區(qū)域統(tǒng)一構(gòu)造應(yīng)力場作用下，關(guān)聯(lián)的斷層系統(tǒng)往往會(huì)出現(xiàn)協(xié)同性的運(yùn)動(dòng)，顯示出對區(qū)域應(yīng)力場變化的敏感響應(yīng)。通過監(jiān)測關(guān)聯(lián)斷層系統(tǒng)整體的斷層運(yùn)動(dòng)變化，有助于了解區(qū)域大范圍、長時(shí)間的應(yīng)力場調(diào)整過程，對于準(zhǔn)確研判未來一段時(shí)期區(qū)域整體的強(qiáng)震趨勢具有重要意義。

跨斷層形變觀測在監(jiān)測斷層運(yùn)動(dòng)方面一直發(fā)揮著重要作用，目前主要是基線和水準(zhǔn)觀測兩類?；€觀測由于所跨斷層走向、基線方位、斷層運(yùn)動(dòng)方式等的不同，直接進(jìn)行有物理含義的區(qū)域綜合難度較大，而水準(zhǔn)觀測因?yàn)橹环从硵鄬哟瓜蛏系倪\(yùn)動(dòng)變化，所以是進(jìn)行區(qū)域綜合比較好的選擇。

北京及周邊地區(qū)跨斷層場地都具有水準(zhǔn)觀測，為開展基于區(qū)域變形配套性的預(yù)測效能評(píng)估提供了基礎(chǔ)。筆者使用垂直形變速率合成方法進(jìn)行計(jì)算（焦青，周俊萍，1998；蘇琴，2005；張晶等，2012），按照跨斷層場地的區(qū)域集中程度、地質(zhì)背景以及構(gòu)造相關(guān)性等，主要分成首都圈西北部、首都圈東北部、八寶山斷裂相關(guān)區(qū)域3個(gè)區(qū)域（圖2）開展計(jì)算分析。首都圈西北部包含狼山、施莊、土木、小水峪、沿河城、燕家臺(tái)、張山營等場地（小寺沖場地觀測時(shí)間較短，未參與合成計(jì)算）；首都圈東北部包含北石城、德勝口、古北口、密云、南口、墻子路、張家臺(tái)、鎮(zhèn)羅營等場地（百善場地?cái)?shù)據(jù)質(zhì)量較差，未參與合成計(jì)算）；八寶山斷裂相關(guān)區(qū)域包含八寶山、大灰廠、?？谟?、京西實(shí)驗(yàn)場、上萬等場地，區(qū)域南部相對零散的場地未進(jìn)行考慮。

3個(gè)區(qū)域的垂直形變速率合成結(jié)果表現(xiàn)出明顯差異，首都圈西北部的計(jì)算結(jié)果呈現(xiàn)2個(gè)顯著的活動(dòng)時(shí)段（圖4a）：一是1970—1979年，在該時(shí)段末期發(fā)生了1976年唐山7.8級(jí)、灤縣7.1級(jí)地震及多次5級(jí)以上余震，該時(shí)段區(qū)域斷層系統(tǒng)整體相對活動(dòng)反映了唐山、灤縣大地震前區(qū)域孕震應(yīng)力環(huán)境變化以及后續(xù)的震后調(diào)整；二是1985—1999年，持續(xù)時(shí)間相對較長，期間區(qū)域內(nèi)雖未發(fā)生大地震，但是中強(qiáng)地震持續(xù)活躍，發(fā)生了9次5級(jí)以上地震，主要集中在山西帶北段和張渤帶西段。之后區(qū)域斷層運(yùn)動(dòng)相對平穩(wěn)，期間地震活動(dòng)也相對平靜，僅發(fā)生1次2006年文安5.1級(jí)地震；隨后2011年日本9.0級(jí)地震對該區(qū)域產(chǎn)生了持續(xù)深遠(yuǎn)的影響，期間區(qū)域內(nèi)5級(jí)地震長時(shí)間平靜可能與日本地震震后引起的區(qū)域應(yīng)力、斷層運(yùn)動(dòng)狀態(tài)的調(diào)整過程有關(guān)（張貝等，2015）；2020年唐山5.1級(jí)地震前后，區(qū)域斷層運(yùn)動(dòng)系統(tǒng)再次表現(xiàn)出一定的活躍跡象，可能指示了日本地震的影響正趨于結(jié)束。整體來看，該區(qū)域的斷層垂直形變速率合成結(jié)果對于區(qū)域應(yīng)力狀態(tài)調(diào)整較為敏感，能夠相對準(zhǔn)確地反映區(qū)域內(nèi)地震活動(dòng)趨勢，各觀測場地的變形配套性相對較好。

首都圈東北部的計(jì)算結(jié)果同樣呈現(xiàn)兩個(gè)顯著的活動(dòng)時(shí)段（圖4b）：一是1976—1979年，主要反映了唐山、灤縣大地震后的調(diào)整過程，并且受震后影響的程度明顯高于首都圈西北部的結(jié)果，而對于震前應(yīng)力環(huán)境的演化過程反而不是特別敏感；二是1985—1999年，顯示了與首都圈西北部類似的演化過程。之后該區(qū)域結(jié)果同樣顯示了相對平穩(wěn)的運(yùn)動(dòng)狀態(tài)以及日本地震對于區(qū)域應(yīng)力、斷層運(yùn)動(dòng)狀態(tài)的影響；該區(qū)域垂直形變速率合成結(jié)果在2021年下半年表現(xiàn)為高值變化，可能與2021年北京市出現(xiàn)的近幾十年內(nèi)最強(qiáng)降雨這類罕見氣象條件有關(guān)，屬于干擾因素的影響。因此，對比首都圈西北部的結(jié)果發(fā)現(xiàn)，本區(qū)域結(jié)果對個(gè)別區(qū)域應(yīng)力調(diào)整過程不敏感，并且易受干擾因素影響，表明該區(qū)域各場地的區(qū)域變形配套性略差。

與前兩個(gè)區(qū)域相比，八寶山斷裂相關(guān)區(qū)域的計(jì)算結(jié)果從1990年前后開始整體波動(dòng)加劇，并未表現(xiàn)出與區(qū)域應(yīng)力場調(diào)整以及周邊中強(qiáng)地震明顯的對應(yīng)關(guān)系（圖4c）。從針對各場地受地下水位變化及地表沉降影響的分析結(jié)果可知，該區(qū)參與計(jì)算的多個(gè)場地受影響嚴(yán)重（圖2），計(jì)算結(jié)果的不穩(wěn)定可能與地下水儲(chǔ)量快速變化，導(dǎo)致不同觀測場地出現(xiàn)相應(yīng)的調(diào)整運(yùn)動(dòng)有關(guān)（圖3c）。另外，該區(qū)域結(jié)果同樣受到2021年極端天氣的影響，并且受影響程度比首都圈東北部的結(jié)果更嚴(yán)重。整體來看，八寶山斷裂相關(guān)區(qū)域的觀測場地其區(qū)域變形配套性所反映的效能最差。

在上述分析基礎(chǔ)上，筆者采用類似2.1節(jié)的方式進(jìn)行區(qū)域垂直形變速率合成結(jié)果的預(yù)測效能定量評(píng)估（表4），定量評(píng)估結(jié)果與定性分析的認(rèn)識(shí)一致。綜合上述結(jié)果，將首都圈西北部參與計(jì)算場地的區(qū)域變形配套性定為A類，首都圈東北部參與計(jì)算場地定為B類，八寶山斷裂相關(guān)區(qū)域參與計(jì)算場地以及未參與上述3個(gè)區(qū)域計(jì)算的其它場地均定為C類。A類場地7個(gè)、B類場地8個(gè)，合計(jì)占總場地?cái)?shù)的55.6%。

3討論

圖5給出了綜合場地監(jiān)測條件以及預(yù)測效能的北京及周邊地區(qū)跨斷層場地最終評(píng)估結(jié)果。該結(jié)果表現(xiàn)出幾個(gè)顯著特點(diǎn)：

（1）跨斷層場地的監(jiān)測條件對于其地震預(yù)測效能具有較大影響，監(jiān)測條件相對較差的場地其預(yù)測效能各項(xiàng)評(píng)估指標(biāo)普遍偏低。

（2）區(qū)域地質(zhì)構(gòu)造環(huán)境的差異對場地監(jiān)測條件有較大影響，進(jìn)而影響其預(yù)測效能。整體來看，位于山地基巖區(qū)的場地觀測條件相對較好，預(yù)測效能普遍較高，位于平原區(qū)沉積層之上的場地則容易受各種干擾因素影響，導(dǎo)致監(jiān)測條件較差，預(yù)測效能普遍偏低。

（3）區(qū)域局部的水文地質(zhì)條件差異對其附近場地受地下水超采、強(qiáng)降雨等影響的顯著程度可能具有重要控制作用。實(shí)際觀測數(shù)據(jù)受地下水儲(chǔ)量變化顯著影響的京西實(shí)驗(yàn)場、大灰廠、上萬、牛口峪等場地，在空間上并未處于淺層地下水位下降和地表沉降最為突出的區(qū)域（圖2）?；趨^(qū)域水文地質(zhì)研究結(jié)果分析（雷坤超等，2016；Lei et al，2021），造成這種現(xiàn)象的原因可能與上述場地位于單一砂礫含水層上方有關(guān)，這同時(shí)造成了上述場地對于強(qiáng)降雨影響也較為敏感。

前人研究已經(jīng)揭示了不同觀測臺(tái)站或場地的監(jiān)測能力及其對地震前兆異常的敏感程度存在顯著差異，臺(tái)站或場地所處的斷裂構(gòu)造位置可能起到主要作用，特別是在斷裂端點(diǎn)、拐點(diǎn)與交匯點(diǎn)等構(gòu)造復(fù)雜、應(yīng)力較易集中的區(qū)域，觀測臺(tái)站較為敏感（車用太，1990；焦明若，張國民，1998；劉琦等，2016），而且不同運(yùn)動(dòng)性質(zhì)的斷層優(yōu)勢的敏感觀測區(qū)域可能也不同（馬棟等，2019），這些認(rèn)識(shí)也得到了巖石破裂實(shí)驗(yàn)等的支持（馬瑾等，2007）。本文的研究進(jìn)一步顯示了區(qū)域地質(zhì)構(gòu)造、水文地質(zhì)條件等對觀測臺(tái)站或場地的監(jiān)測能力及預(yù)測效能等的影響，因此在后續(xù)針對北京及周邊地區(qū)觀測站網(wǎng)布局的優(yōu)化改造中，除了考慮斷裂構(gòu)造等因素，區(qū)域地質(zhì)構(gòu)造、水文地質(zhì)條件等同樣需要重點(diǎn)考慮。

4結(jié)論

本文針對北京及周邊地區(qū)的跨斷層觀測場地嘗試綜合考慮監(jiān)測條件、預(yù)測效能進(jìn)行評(píng)估，分別從監(jiān)測環(huán)境、監(jiān)測有效性、地下水位變化/地表沉降影響以及場地預(yù)測配套性、最優(yōu)預(yù)測指標(biāo)效能、區(qū)域變形配套性等方面開展了詳細(xì)的研究和分析討論?；诰C合評(píng)估結(jié)果，主要獲得以下結(jié)論：

（1）北京及周邊地區(qū)的跨斷層場地布設(shè)對區(qū)域內(nèi)大多數(shù)主要斷層具有較好的監(jiān)測能力，并有較好的地震預(yù)測效能。

（2）八寶山—黃莊—高麗營斷裂多數(shù)場地（例如八寶山、京西實(shí)驗(yàn)場、大灰廠、上萬、?？谟龋┠壳氨O(jiān)測條件相對較差，預(yù)測效能普遍偏低，建議針對該斷裂選擇其它合適的斷層段位置補(bǔ)充新建若干跨斷層觀測場地。

（3）基于評(píng)估結(jié)果，建議停測百善、南孟、寧河、淶水、京西實(shí)驗(yàn)場、大灰廠等預(yù)測效能較低且監(jiān)測條件一般的場地。

（4）考慮所跨斷層的主要運(yùn)動(dòng)方式，建議古北口、鎮(zhèn)羅營、南口等場地增加基線觀測，以有效監(jiān)測所跨斷層的斷層運(yùn)動(dòng)。

（5）后續(xù)若要補(bǔ)充新建跨斷層形變觀測場地，在場地選址時(shí)除了考慮所跨斷裂基本活動(dòng)特性之外，還需考慮區(qū)域地質(zhì)構(gòu)造、水文地質(zhì)條件等因素。

感謝北京市地震局劉桂萍研究員、中國地震局地震預(yù)測研究所孫漢榮研究員、薛兵研究員、云南省地震局洪敏正高級(jí)工程師、山西省地震局劉瑞春高級(jí)工程師、河北省地震局馬棟高級(jí)工程師、中國地震局第一監(jiān)測中心陳長云高級(jí)工程師、中國地震局第二監(jiān)測中心李瑞莎高級(jí)工程師在本項(xiàng)工作開展及文章撰寫過程中給予的意見建議。

參考文獻(xiàn)：

薄萬舉，謝覺民，郭良遷.1998.八寶山斷裂帶形變分析與探討［J］.地震，18（1）：63-68.

曹建玲，賈曉東，李臘月，等.2016.華北地區(qū)跨斷層前兆異常特征綜述［J］.地震，36（4）：61-75.

車用太.1990.地下水動(dòng)態(tài)觀測井映震靈敏條件研究［J］.中國地震，6（4）：69-75.

陳長云，鄭智江，李臘月，等.2016.南北地震帶中南段典型強(qiáng)震震前跨斷層形變特征分析［J］.地震，36（4）：47-60.

竇甜甜，程惠紅，周元澤，等.2022.華北平原地下水開采對區(qū)域地震活動(dòng)性的影響［J］.地球物理學(xué)報(bào)，65（8）：2931-2944.

高文學(xué)，馬瑾.1993.首都圈地震地質(zhì)環(huán)境與地震災(zāi)害［M］.北京：地震出版社.

郭海朋，白晉斌，張有全，等.2017.華北平原典型地段地面沉降演化特征與機(jī)理研究［J］.中國地質(zhì)，44（6）：1115-1127.

國家地震局科技監(jiān)測司.1990.地震分析預(yù)報(bào)方法程式指南［M］.北京：地震出版社.

郝書儉，尤惠川.2001.唐山活動(dòng)斷層的淺層地震精細(xì)探測［J］.地震地質(zhì)，23（1）：93-97.

何慶成，劉文波，李志明.2006.華北平原地面沉降調(diào)查與監(jiān)測［J］.高校地質(zhì)學(xué)報(bào)，12（2）：195-209.

洪敏，張晶，邵德盛，等.2019.跨斷層形變智能化資料處理與定量化預(yù)測指標(biāo)［C］//2019年中國地球科學(xué)聯(lián)合學(xué)術(shù)年會(huì)論文集（八）——專題21：活動(dòng)地塊邊界帶強(qiáng)震機(jī)理與預(yù)測、專題22、中國地震科學(xué)實(shí)驗(yàn)場.北京，26.

蔣銘，馬宗晉.1985.華北第三、四地震活躍期的對比［J］.地震，5（6）：5-11.

焦明若，張國民.1998.地震前兆復(fù)雜性成因機(jī)理研究的討論（二）——地震前兆復(fù)雜性成因機(jī)理［J］.地震，18（2）：112-118.

焦青，周俊萍.1998.跨斷層垂直位移速率的動(dòng)態(tài)演化特征與強(qiáng)震的關(guān)系［J］.地震，18（3）：265-273.

雷坤超，羅勇，陳蓓蓓，等.2016.北京平原區(qū)地面沉降分布特征及影響因素［J］.中國地質(zhì)，43（6）：2216-2225.

李臘月，季靈運(yùn)，李玉江，等.2016.基于跨斷層測量的鮮水河斷裂帶現(xiàn)今活動(dòng)特征及其與強(qiáng)震關(guān)系研究［J］.地震學(xué)報(bào)，38（5）：739-750.

李瑞莎，張希，唐紅濤，等.2016.青藏塊體東北緣跨斷層形變典型震例總結(jié)與預(yù)測判據(jù)指標(biāo)分析［J］.地震，36（4）：35-46.

梁芳，孫建寶，沈正康，等.2013.L波段InSAR數(shù)據(jù)觀測的北京及其周邊2007—2010年間累計(jì)地殼形變基本特征［J］.地震，33（4）：43-54.

劉琦，閆偉，李智蓉，等.2016.南北地震帶定點(diǎn)形變前兆異常指標(biāo)初建［J］.地震，36（4）：76-88.

陸明勇，黃寶森，劉天海，等.2011.首都圈跨斷層流動(dòng)形變觀測資料映震能力及場地監(jiān)測能力分析［J］.地震，31（3）：121-134.

陸明勇，牛安福，黃寶森，等.2013.首都圈地區(qū)跨斷層流動(dòng)形變監(jiān)測場地優(yōu)化布設(shè)探討［J］.大地測量與地球動(dòng)力學(xué)，33（S1）：1-4，15.

陸明勇，于海生，李志雄，等.2018.華北地區(qū)跨斷層流動(dòng)形變監(jiān)測場地優(yōu)化初探［J］.中國地震，34（4）：837-844.

馬棟，馬利軍，侯曉真，等.2019.華北地區(qū)地殼形變連續(xù)觀測靈敏臺(tái)站研究［J］.地震研究，42（3）：405-410.

馬瑾，劉力強(qiáng)，劉培洵，等.2007.斷層失穩(wěn)錯(cuò)動(dòng)熱場前兆模式：雁列斷層的實(shí)驗(yàn)研究［J］.地球物理學(xué)報(bào)，50（4）：1141-1149.

史珉，宮輝力，陳蓓蓓，等.2021.Sentinel-1A京津冀平原區(qū)2016—2018年地面沉降InSAR監(jiān)測［J］.自然資源遙感，33（4）：55-63.

蘇琴.2005.川滇菱形塊體各邊界垂直形變速率與斷層活動(dòng)特性［J］.四川地震，（1）：43-47.

孫啟凱，池國民，徐東卓.2017.首都圈地區(qū)地殼垂直形變特征及剖面分析［J］.大地測量與地球動(dòng)力學(xué)，37（5）：497-501.

徐錫偉，吳衛(wèi)民，張先康，等.2002.首都圈地區(qū)地殼最新構(gòu)造變動(dòng)與地震［M］.北京：科學(xué)出版社.

楊海盟，戴俊生，汪必峰，等.2014.牛駝鎮(zhèn)凸起發(fā)育過程［J］.東北石油大學(xué)學(xué)報(bào)，38（6）：22-29，91.

楊艷，賈三滿，王海剛，等.2013.北京規(guī)劃新城地面沉降影響分析［J］.城市規(guī)劃，37（11）：67-11.

楊艷，賈三滿，王海剛.2010.北京平原區(qū)地面沉降現(xiàn)狀及發(fā)展趨勢分析［J］.上海地質(zhì)，31（4）：23-28.

楊艷.2015.京津冀區(qū)域地面沉降災(zāi)害防治思考［J］.城市地質(zhì)，10（1）：1-7.

尤惠川，徐錫偉，吳建平，等.2002.唐山地震深淺構(gòu)造關(guān)系研究［J］.地震地質(zhì)，24（4）：571-582.

查明，尉亞民，高長海，等.2011.牛駝鎮(zhèn)凸起南段潛山勘探潛力分析［J］.巖性油氣藏，23（2）：10-14.

張貝，張懷，石耀霖.2015.有限元模擬彈性位錯(cuò)的等效體力方法［J］.地球物理學(xué)報(bào)，58（5）：1666-1674.

張國民，馬宗晉，蔣銘.1988.華北強(qiáng)震規(guī)律的研究［J］.中國地震，4（3）：66-70.

張晶，劉峽，方穎，等.2012.日本9.0級(jí)地震前后華北形變場變化特征與地震趨勢分析［J］.地震，32（2）：64-72.

張世民，王丹丹，劉旭東，等.2007.北京南口—孫河斷裂帶北段晚第四紀(jì)活動(dòng)的層序地層學(xué)研究［J］.地震地質(zhì)，29（4）：730-743.

中國地震局監(jiān)測預(yù)報(bào)司.2020.地震危險(xiǎn)性判定技術(shù)方法系列叢書——形變分析預(yù)測技術(shù)方法工作手冊［M］.北京：地震出版社.

Bai Z，Wang Y，Balz T.2022.Beijing land subsidence revealed using PS-InSAR with long time series TerraSAR-X SAR data［J］.Remote Sensing，14（11）：2529.

Lei K，Ma F，Chen B，et al.2021.Three－dimensional surface deformation characteristics based on time series InSAR and GPS technologies in Beijing，China［J］.Remote Sensing，13（19）：3964.

Li Z，Yu C，Xiao R，et al.2022.Entering a new ear of InSAR：advanced techniques and emerging applications［J］.Journal of Geodesy and Geoinformation Science，5（1）：1-4.

Pang Y，Zhang H，Cheng H，et al.2020.The modulation of groundwater exploitation on crustal stress in the North China Plain，and its implications on seismicity［J］.Journal of Asian Earth Sciences，189：104141.

Zhang T，Shen W，Wu W，et al.2019.Recent surface deformation in the Tianjin area revealed by Sentinel-1A data［J］.Remote Sensing，11（2）：130.

Evaluation of the Cross－fault Observation Network

in Beijing and Its Adjacent Areas

LIU Qi LI Layue HU Leyin

（1.Institute of Earthquake Forecasting，China Earthquake Administration，Beijing 100036，China）

（2.The First Monitoring and Application Center，China Earthquake Administration，Tianjin 300180，China）

（3.Beijing Earthquake Agency，Beijing 100080，China）

Abstract

Beijing and its adjacent areas are characterized by developed economy，dense population and frequent historical strong earthquakes.In order to accurately study and judge the seismic risk in Beijing and its adjacent areas and reduce the disaster risk that may be caused by earthquakes，it is necessary to first understand the effectiveness of the current observational sites in this region.For this reason，this paper comprehensively considers the monitoring conditions and the effectiveness of earthquake forecasting of the cross－fault sites to design the evaluation scheme of the observational network.This paper studies the monitoring environment，the monitoring effectiveness，the effect of the groundwater level change and the surface subsidence，the matching of forecasting，the effectiveness of the optimal forecasting index and the matching of regional deformation.This paper proposes the comprehensive evaluation results and suggestions for the optimization of the observational network.The distribution of the current cross－fault sites in Beijing and its adjacent areas have good monitoring capacity for most of the major faults in the region and have good earthquake forecasting efficiency.However，the monitoring conditions of most observational sites on the Babaoshan－Huangzhuang－Gaoliying fault are relatively poor，and their forecasting efficiency is generally low.It is advisable to select other locations on the fault to supplement the observational sites.In addition to the fault structure，regional geological settings and hydrogeological conditions have an important impact on the effectiveness of the observational site，and need to be taken into account when selecting new sites.

Keywords：cross－fault observation；earthquake forecasting effectiveness；monitoring conditions；groundwater；tectonic conditions；Beijing