張建勇 崔振東 周 健 韓偉歌 于 慧 王燚釗

1.中國科學(xué)院頁巖氣與地質(zhì)工程重點實驗室·中國科學(xué)院地質(zhì)與地球物理研究所 2.中國科學(xué)院地球科學(xué)研究院3.中國科學(xué)院大學(xué)地球與行星科學(xué)學(xué)院 4.北京石油化工學(xué)院機械工程學(xué)院5.中國礦業(yè)大學(xué)（北京）地球科學(xué)與測繪工程學(xué)院 6.北京科技大學(xué)土木與資源工程學(xué)院

1 流體注入誘發(fā)地震簡介

增強型地?zé)崮艿拈_發(fā)、二氧化碳地質(zhì)封存和頁巖氣開采等深部地質(zhì)工程都需要向深部地層中注入高壓流體。此過程可能引起斷層活化，誘發(fā)低級別有感地震[1-7]。這些斷層活化風(fēng)險除了可能影響深部地質(zhì)工程的正常施工和安全運營外，還可能對人民群眾生命財產(chǎn)安全造成不利影響，降低社會公眾對于此類地下注入流體工程的信任度。因此，研究此類地震的機理、影響因素、風(fēng)險評估方法以及應(yīng)對措施，平復(fù)學(xué)界各方爭議，給民眾以科學(xué)、清晰的認(rèn)識迫在眉睫又至關(guān)重要。

20世紀(jì)60年代，Healy等[8]在研究洛基山Arsenal地震序列時首次提出：向地下深部注入流體可能誘發(fā)地震（圖1[9]）。隨后，Raleigh等[10]在Rangely油田開展現(xiàn)場注水試驗。首次證明了流體注入壓力增加會導(dǎo)致地震的產(chǎn)生，并用Hubbert—Rubey有效應(yīng)力原理進行了解釋。后來的工程實踐中有諸多注入流體誘發(fā)地震案例，美國及加拿大南部地區(qū)發(fā)生過幾次震級較大的注入流體誘發(fā)的地震。例如：2010年Arkansas州的地震矩震級（MW）4.0級地震[11-12]，2011年Oklahoma州的MW5.7級地震[13-14]、2011年Colorado州的MW5.3級地震[15-16]、2011年Ohio州的MW3.9級地震[16]，2012年Texas州的MW4.8級地震[17]。

我國也有流體注入誘發(fā)地震的案例發(fā)生[18]。四川盆地作為我國主要的油氣產(chǎn)區(qū)之一，有多處注水采氣井及廢水回注井。比如，四川盆地某國家級頁巖氣開采示范區(qū)歷史上地震活動頻率較低，隨著頁巖氣的開采，近年來在該地區(qū)監(jiān)測到的地震活動頻率顯著增加；四川盆地西南部某地區(qū)1988—2006年向地層深部注入廢水超過100×104m3，該項注水工作已引發(fā)了超過32 000次有地表記錄的地震。

在流體注入活動誘發(fā)地震獲得越來越多關(guān)注的同時，斷層活化的機理、影響因素、判別標(biāo)準(zhǔn)以及斷層活化的致災(zāi)模型[19]等方面的研究也取得了諸多進展。此類研究有助于探明適于進行地下注入活動的場地工程地質(zhì)條件，研究儲層各項物理力學(xué)性質(zhì)對斷層活化的影響，并產(chǎn)生用于指導(dǎo)實際地下流體注入工程的指南和規(guī)范，從而科學(xué)平衡風(fēng)險與效益，讓地下流體注入工程為人類造福。

2 斷層活化影響因素及關(guān)鍵評價指標(biāo)

注入活動誘發(fā)地震受諸多因素共同控制，比如：流體注入速率及注入量、儲層巖石的物理力學(xué)性質(zhì)、斷層的展布與幾何形態(tài)、斷層所在區(qū)域地應(yīng)力狀態(tài)等。針對不同影響因素，國內(nèi)外學(xué)者們運用了諸多手段加以研究，不同因素對斷層活化風(fēng)險的影響存在較大差異。

2.1 流體注入速率及注入量

流體的總注入量是影響斷層活化極為重要的因素。在注入過程中，隨著總注入量的增加，注入壓力增大，巖層的孔隙水壓力增大，使得斷層面上有效應(yīng)力降低，導(dǎo)致斷層活化風(fēng)險增大。Nicol等[20]對北美地區(qū)30處注水誘發(fā)地震場址的注水資料進行統(tǒng)計分析，得出誘發(fā)地震震級大小與流體總注入量呈正相關(guān)，其相關(guān)方程為：

式中y表示地震震級，級；x表示總注入量，1 000 m3。

在現(xiàn)場施工中，流體注入速率與注入量是工程技術(shù)人員關(guān)注的重點。張致偉等[21]在四川某地區(qū)注水誘發(fā)地震的研究中，2009年1—7月，地震活動頻次隨注水壓力升高而增大，譜振幅相關(guān)系數(shù)由0.84增加至0.92，誘發(fā)地震頻次與注水壓力呈現(xiàn)良好的相關(guān)性。

由上述分析可知，當(dāng)其他條件相同時，注入速率越大，斷層活化的可能性越高；另一方面，長時間以低速率注入流體也會導(dǎo)致斷層活化，流體沿儲層向斷層區(qū)域擴散，在斷層處積累了較強的孔隙壓力，此時斷層活化所觸發(fā)的地震震級往往高于高速率注入所誘發(fā)的地震。

2.2 斷層的地應(yīng)力條件

斷層區(qū)域的地應(yīng)力條件是斷層受力狀態(tài)的主控因素，也是斷層活化最為敏感的因素。對于注入流體誘發(fā)斷層活化，斷層發(fā)生剪切破壞的潛力要遠大于拉張破壞。因此，斷層處的水平應(yīng)力越大，越能起到限制斷層剪切滑移的側(cè)限作用[22]。

基于上述結(jié)論，學(xué)者們研究了斷層的橫縱應(yīng)力比對斷層活化的影響。在其他條件保持不變的情況下，改變斷層區(qū)域的橫縱應(yīng)力比，對比斷層滑移距離的變化。斷層所在位置的橫縱應(yīng)力比越大，水平應(yīng)力越高，斷層活化的可能性越小，即水平地應(yīng)力起到限制斷層剪切滑移作用[23]。

2.3 斷層巖石滲透率

不同滲透率的斷層在注水過程中表現(xiàn)的差別很大。流體沿著巖層內(nèi)部的節(jié)理裂隙等通道從注水井向周圍擴散，斷層滲透率越低，流體越容易在斷層處富集，導(dǎo)致斷層區(qū)域孔隙水壓力增大，進而誘發(fā)斷層活化。

Figueiredo等[23]通過TOUGH-FLAC水力耦合軟件分析了滲透率對斷層活化的影響，在兩條斷層中間注水，改變斷層巖石的滲透率，統(tǒng)計不同滲透率下斷層活化情況（表1）。分析表明，在其他條件相同時，滲透率越低（小于10－3mD）的斷層，流體更容易在斷層處富集，誘發(fā)斷層活化的風(fēng)險更高。

表1 不同初始滲透率對斷層活化影響分析數(shù)據(jù)表[23]

2.4 斷層的類型及產(chǎn)狀

斷層的類型及產(chǎn)狀影響地應(yīng)力與斷層的相互作用方式，進而影響斷層活化風(fēng)險。Swierczek等[24-25]以澳大利亞Gippsland盆地不同性質(zhì)的斷層為模型建模，在相同條件下注入CO2，監(jiān)測3種不同類型斷層的滑移情況。由于不同類型的斷層與地應(yīng)力相互作用方式不同，因此，注入CO2誘發(fā)斷層活化的臨界孔隙壓力也不同。監(jiān)測結(jié)果顯示，走滑逆斷層臨界孔隙壓力為2 MPa，走滑斷層為4 MPa，逆斷層為6 MPa。同時，由受力分析可知，高傾角斷層發(fā)生活化的可能性也要比低傾角斷層大得多。

3 斷層活化風(fēng)險評估方法

3.1 基于摩爾—庫倫準(zhǔn)則的基本理論解釋方法

將CO2等流體注入地層內(nèi)部的脆性巖石中，會極大地提高誘發(fā)地震的風(fēng)險。脆性介質(zhì)內(nèi)部的斷層和裂隙為液體在其內(nèi)部的流動提供了高滲透性的通道。因此，在已有斷層存在的地質(zhì)條件下，孔隙流體壓力的變化成為引起地震活動的關(guān)鍵[26]。

基于此，斷層活化的力學(xué)機制可以用有效應(yīng)力理論解釋，公式為：

式中τn表示斷層面上抗剪強度，MPa；μ表示斷層面的內(nèi)摩擦系數(shù)，無量綱；σn表示斷層面上的正應(yīng)力，MPa；p表示孔隙水壓力，MPa；c表示黏聚力，MPa。

引入活化趨勢因子（Reactivation Tendency Factor）公式為：式中fR表示活化趨勢因子，無量綱；τ0表示斷層固有抗剪強度，MPa。

斷層活化的力學(xué)機制可用圖2[27]三維應(yīng)力莫爾圓解釋，注入流體使得斷層面上的孔隙水壓力p增大，有效正應(yīng)力(σn－p)降低，導(dǎo)致斷層面的抗剪強度[μ(σn－p)＋c]降低。當(dāng)作用于斷層面上的剪應(yīng)力高于其抗剪強度，即fR＞1時，斷層發(fā)生剪切滑動，導(dǎo)致斷層活化進而誘發(fā)地震。

圖2 用三維應(yīng)力莫爾圓判斷斷層活化過程力學(xué)機理圖[27]

美國俄亥俄州Perry地區(qū)地震事件是應(yīng)用摩爾—庫倫原理解釋注水誘發(fā)地震的典型案例。圖3是結(jié)合該地區(qū)注水資料以及地應(yīng)力資料繪制的應(yīng)力莫爾圓[28]，隨著流體不斷注入，斷層面的孔隙水壓力增大，有效應(yīng)力減小，莫爾圓沿著橫軸向左移動。當(dāng)莫爾圓觸及摩爾—庫倫破壞包線，斷層發(fā)生滑移。

圖3 美國俄亥俄州Perry地區(qū)地震事件的應(yīng)力莫爾圓解釋圖[28]

3.2 臨界孔隙壓力預(yù)測法

由有效應(yīng)力原理可知，觸發(fā)斷層活化的關(guān)鍵在于最大流體壓力增量。若能確定某一斷層發(fā)生活化的臨界孔隙壓力值，則可對實際流體注入工程起到指導(dǎo)作用。巖層所處的地應(yīng)力體系是影響斷層受力狀態(tài)的外因，斷層的黏聚力以及內(nèi)摩擦角等性質(zhì)是影響斷層最大孔隙壓力增量的內(nèi)因。將兩因素結(jié)合運用到斷層活化臨界孔隙壓力的分析中，是評判斷層活化風(fēng)險的新思路。Wiprut和Zoback等[29-31]提出臨界孔隙壓力法（以下簡稱CPP法）評估斷層活化風(fēng)險（圖4-a）。在已知斷層區(qū)域應(yīng)力場、孔隙壓力、斷層傾向以及斷層巖石的摩爾包線的情況下，綜合考慮斷層的持水度（σ3－p）、滑移趨勢（τ/σn'）、膨脹趨勢（σ1'－σn'）/（σ1'－σn'）以及摩爾—庫倫破壞準(zhǔn)則等一系列方法[32]，計算出臨界孔隙壓力。此方法沒有考慮斷層的黏聚力，故誘發(fā)斷層活化的最大孔隙壓力增量（?p）過于保守。因此，Mildren等[33-34]提出考慮斷層黏聚力的最大流體壓力增量方法（以下簡稱FAST法）。由于考慮黏聚力，因此摩爾庫倫破壞曲線上移，臨界孔隙壓力（?p）增大，如圖4-b所示。

Lee等[35]實地取樣測得斷層巖石的物理力學(xué)參數(shù)及所處地應(yīng)力狀態(tài)數(shù)據(jù)，利用CPP法對韓國Pohang盆地注入CO2誘發(fā)斷層活化風(fēng)險進行了數(shù)值模擬預(yù)測。在不考慮斷層的黏聚力時，數(shù)值模擬的實驗結(jié)果表明：在μ=0.25時，該斷層的臨界孔隙壓力為3 MPa；μ=0.6時，臨界孔隙壓力為10 MPa。Ruth等[36]使用FAST方法對澳大利亞Gippsland盆地注入CO2誘發(fā)斷層活化進行評估，結(jié)果顯示，此處無黏聚力斷層的臨界孔隙壓力為3.8 MPa，愈合斷層的臨界孔隙壓力為15.6 MPa。上述兩種方法需要實地測量斷層的各項物理力學(xué)參數(shù)，從而得到誘發(fā)斷層活化的大致臨界孔隙壓力預(yù)測值。此預(yù)測值可以用來判斷不同斷層活化風(fēng)險的相對大小，由于黏聚力和內(nèi)摩擦角難以準(zhǔn)確測量等問題，若將預(yù)測值用于指導(dǎo)實際流體注入工程，其準(zhǔn)確性尚存疑。

3.3 三維評估法

以往對斷層活化風(fēng)險評估方法多是針對某一具體斷層，對包含諸多斷層的場地區(qū)塊安全性分析不具有普適性。對場地內(nèi)諸多斷層的活化風(fēng)險性分級，可從中找出危險區(qū)域。Sibson[37]最早提出根據(jù)斷層橫縱有效應(yīng)力比（SR=σ1'/σ3'）的大小，將場地內(nèi)的斷層活化的危險性等級分為高中低3類，然而此方法評價指標(biāo)過于單一，且只考慮等圍壓情況。Leclère和Fabbri[27]在此基礎(chǔ)上提出了一種三向受力狀態(tài)下，基于橫縱有效應(yīng)力比的斷層活化風(fēng)險三維定量評價方法。聯(lián)合應(yīng)力形狀比Φ=（σ2'－σ3'）/（σ1'－σ3'）[38]與SR，從而將x、y兩個方向的有效正應(yīng)力σ2'與σ3'用z方向的有效正應(yīng)力σ1'表示，結(jié)合斷層的內(nèi)摩擦系數(shù)、黏聚力等因素，經(jīng)公式（4）得出Q值。

其中：

圖4 臨界孔隙壓力法判斷斷層活化過程力學(xué)機理圖[29-34]

式中l(wèi)、m、n分別表示斷層面的法線在z、x、y三軸正應(yīng)力方向上的余弦；σ1、σ2、σ3分別表示z、x、y方向的正應(yīng)力，MPa；σ1'、σ2'、σ3'分別表示z、x、y方向的有效正應(yīng)力，MPa；Φ表示應(yīng)力形狀比，無量綱。

將斷層Q值與界限值Qo比較，其中Qo由斷層的內(nèi)摩擦角和黏聚力得出。評價標(biāo)準(zhǔn)如表2所示。

表2 基于Q值的斷層活化風(fēng)險評估標(biāo)準(zhǔn)表

三維評估法對孔隙水壓力升高的場地內(nèi)所有可能發(fā)生活化的斷層進行評估，從而標(biāo)定區(qū)域內(nèi)各個斷層發(fā)生斷層活化的風(fēng)險性高低。Leclère和Fabbri[27]將此方法應(yīng)用于日本Tohoku-Oki地區(qū)主震后，地震區(qū)域沿板塊界面產(chǎn)生超孔隙水壓力，導(dǎo)致周圍145處斷層出現(xiàn)斷層活化。根據(jù)測得區(qū)域內(nèi)所有可能發(fā)生活化的斷層主應(yīng)力、應(yīng)力形狀比、斷層產(chǎn)狀以及各項物理力學(xué)參數(shù)，代入公式計算出Q值，將全部斷層的Q值繪制成圖5。結(jié)果顯示，實際地球物理方法所探測到的斷層活化位置與三維評估法所得藍色區(qū)域大致重疊。

圖5 日本Tohoku-Oki地區(qū)斷層活化風(fēng)險評估圖[27]

3.4 數(shù)值模擬方法

目前針對斷層活化風(fēng)險性分析，使用頻率較高的數(shù)值模擬方法有離散元法、有限差分法、有限元法以及邊界元法等。Yoon等[39]利用水力耦合離散元方法對斷層帶進行建模，研究注水過程中流體注入與斷層帶破裂的關(guān)系；Rutqvist[40]基于TOUGH-FLAC水力耦合軟件，開發(fā)了流體注入誘發(fā)斷層活化的數(shù)值計算模型；Keranen等[13-14]利用有限元方法建立了地震定位與孔隙壓力關(guān)系模型；Pereira等[41]利用非線性有限元方法開發(fā)了二維地質(zhì)力學(xué)模型；龔發(fā)雄等[42]介紹了Thomas所開發(fā)的基于邊界元方法的Poly 3D軟件。這些方法在斷層活化的模擬上得到很好的應(yīng)用，其斷層活化評估方法的關(guān)系、原理以及評判標(biāo)準(zhǔn)如圖6所示。

圖6 斷層活化評價方法的使用流程圖

4 問題討論與對策

4.1 問題討論

注入流體誘發(fā)斷層活化研究已有近60年的歷史，但仍有許多值得研究的領(lǐng)域，也存在諸多問題。

1）易被誘發(fā)活化斷層的界定較為困難。斷層深埋于地下，其物理力學(xué)性質(zhì)、周圍地應(yīng)力狀態(tài)等信息難以準(zhǔn)確測得，對斷層孔隙壓力實施監(jiān)測存在難度和現(xiàn)實可行性。目前的監(jiān)測多是現(xiàn)場微震監(jiān)測，雖然具有一定的研究意義，但是很難做到提前預(yù)知。

2）地下注入超臨界CO2存在水巖化學(xué)反應(yīng)、流固耦合以及相態(tài)轉(zhuǎn)化等諸多問題，增加斷層活化的不確定性。在地下注入超臨界CO2過程中，CO2與原巖的水化學(xué)作用會對巖層滲透率、斷層所能承受的最大孔隙壓力等性質(zhì)進行改造，會增加斷層活化的不確定性。除此之外，超臨界CO2在誘發(fā)斷層活化過程中，由于壓力以及溫度的變化存在相態(tài)轉(zhuǎn)化問題，其誘震過程及作用機制尚不明晰，影響斷層活化判斷。目前的CO2地質(zhì)封存相關(guān)工程案例中，并未監(jiān)測到震級超過3級的地震[43]。但可以預(yù)知的是，在誘發(fā)斷層活化過程中，超臨界CO2由于其壓縮性強，且自身相態(tài)轉(zhuǎn)化要釋放能量，相比于水等介質(zhì)，對斷層活化的破壞性更強，危險性更大。目前此方面的研究欠缺，研究手段也面臨固—液—氣三相耦合等諸多問題。

3）數(shù)值模擬方法缺乏物理實驗驗證。目前斷層活化的研究多是基于數(shù)值模擬，研究人員從現(xiàn)場獲取巖樣，在室內(nèi)測得其物理力學(xué)參數(shù)，通過計算機對注入—滲流—誘震過程進行數(shù)值模擬。然而該方面的物理實驗經(jīng)驗匱乏，若能通過物理模型實驗以及現(xiàn)場注入流體誘發(fā)地震實驗對數(shù)值模擬結(jié)果進行檢驗，可以為數(shù)值模擬所獲得的規(guī)律性認(rèn)識提供強有力的證據(jù)。目前來看，物理模型試驗的設(shè)計思路較為復(fù)雜，如何準(zhǔn)確有效的設(shè)計模型也是一大難題；在現(xiàn)場注入誘震實驗中，如何確保安全性也是值得研究的問題。

4）斷層活化存在滯后響應(yīng)效應(yīng)。流體注入初期的地震事件大多震級較小，產(chǎn)生原因多與流體注入導(dǎo)致巖層及裂隙破壞有關(guān)。震級較大的地震往往在流體注入一段時間后發(fā)生[44]。流體從注入井向周圍經(jīng)歷一定時間的擴散，在斷層位置富集。隨著注入量的增加，斷層處孔隙水壓力不斷增加，導(dǎo)致斷層活化。目前的研究鮮有分析流體需要多長時間的運移與積累，故無法給出誘發(fā)斷層活化的大致時間。

4.2 對策及建議

面臨眾多問題，未來注入流體誘發(fā)地震的研究方向以及應(yīng)對措施可從以下幾個方面取得突破。

1）加強監(jiān)測。加強地震流動測震臺站和數(shù)字測震臺網(wǎng)的監(jiān)測，利用地震學(xué)方法定位震源震中位置，獲得震源機制解。同時與合成孔徑雷達干涉所監(jiān)測到的地表變形相結(jié)合，在注入流體區(qū)域建立地震預(yù)警機制?？蒲腥藛T應(yīng)加強對監(jiān)測數(shù)據(jù)的研究，結(jié)合震源機制解、地震波形分析等對其機理深入探討。

2）科學(xué)選址。在實際工程中，流體注入場地的選址應(yīng)密切關(guān)注活動斷裂帶、深大斷裂帶的封閉性和活動性[45-46]，避免在地殼中的脆性巖石中注入流體。資料表明，地殼脆性巖體中的老斷層對孔隙壓力變化極為敏感，注入流體活動極易誘發(fā)斷層活化。因此流體注入選址時應(yīng)對現(xiàn)場的地質(zhì)條件、地質(zhì)構(gòu)造以及地震歷史等因素進行綜合考慮。

3）重視歷史注入流體資料的研究。流體注入量以及注入速率是影響斷層活化極為重要的因素，也是工程技術(shù)人員在實際工程中的可控制因素。對于注入流體誘發(fā)地震區(qū)域的歷史注入資料應(yīng)妥善保存并加以詳細研究。將有資料記載的斷層活化區(qū)的斷層物理力學(xué)參數(shù)、儲層性質(zhì)以及所處地應(yīng)力狀態(tài)與注入資料對比，分析其在流體壓力逐漸增大時的活化特征，探索斷層活化的臨界流體壓力值。從而在實際注入過程中根據(jù)儲層實際情況合理調(diào)節(jié)注入量與注入速率。

4）盡快制訂完善地下流體注入技術(shù)規(guī)范及指南。2015年3月，美國聯(lián)邦政府頒布了水力壓裂法細則，該細則于同年6月開始實施，這是美國在水力壓裂方面首個國家層面上的技術(shù)細則。我國也發(fā)布了各項地下流體注入指南及規(guī)范《水力壓裂安全技術(shù)要求》《二氧化碳捕集、利用與封存環(huán)境風(fēng)險評估技術(shù)指南（試行）》等。然而，此類規(guī)范及指南沒有對注入流體誘發(fā)斷層活化的風(fēng)險評估做具體定量化的規(guī)定。我國應(yīng)盡快建立一套流體注入誘發(fā)斷層活化的評價體系，并制定相關(guān)的技術(shù)標(biāo)準(zhǔn)及規(guī)范，平衡經(jīng)濟效益與地震風(fēng)險，從而對現(xiàn)場注入流體提供指導(dǎo)。

5 結(jié)論

1）影響斷層活化的因素眾多，其中地應(yīng)力是主控因素。不同傾向、傾角的斷層是否活化受所處的地應(yīng)力狀態(tài)控制，當(dāng)?shù)貞?yīng)力對斷層滑移起約束作用時，斷層活化的風(fēng)險相對較低。斷層的滲透率、內(nèi)摩擦角、黏聚力等物理力學(xué)參數(shù)是斷層活化的內(nèi)因。流體注入速率及注入量是引發(fā)斷層活化的誘發(fā)因素。

2）目前的斷層活化風(fēng)險評估方法多是摩爾—庫倫有效應(yīng)力原理的拓展和延伸。在實際運用中，可根據(jù)實際工程地質(zhì)情況靈活運用。各種評估方法的絕對準(zhǔn)確性尚存疑，在實際使用以及研究過程中需不斷改進和創(chuàng)新。

3）增強型地?zé)崮荛_發(fā)、二氧化碳地質(zhì)封存以及頁巖氣開采等向深部注入高壓流體的新型地質(zhì)工程領(lǐng)域，人類的工程實踐經(jīng)驗尚不充足，對斷層活化機理的認(rèn)識遠未成熟，理論研究和技術(shù)水平也有待提高，目前仍存在大量亟待解決的問題。怎樣降低深部地質(zhì)工程的誘震風(fēng)險，讓新技術(shù)造福人類，仍需廣大科研人員和技術(shù)人員不斷探索。