基于有限差分法的三峽庫(kù)首區(qū)蓄水前后應(yīng)力場(chǎng)變化與地震活動(dòng)研究

摘要： 水庫(kù)蓄水前后構(gòu)造應(yīng)力場(chǎng)的變化在水庫(kù)誘發(fā)地震研究中具有重要意義。為了探究三峽水庫(kù)蓄水后應(yīng)力場(chǎng)變化及蓄水與水庫(kù)地震的相關(guān)性，以三峽水庫(kù)庫(kù)首區(qū)九畹溪斷裂、仙女山斷裂和高橋斷裂為研究對(duì)象，選取自2003年開(kāi)始蓄水之后72個(gè)月的水位變化數(shù)據(jù)，模擬研究不同蓄水期庫(kù)水荷載對(duì)三峽庫(kù)首區(qū)有效應(yīng)力及剪應(yīng)變?cè)隽康挠绊??；赗HINO軟件實(shí)現(xiàn)真三維地層建模，并應(yīng)用有限差分軟件FLAC^3D模擬不同蓄水期庫(kù)首區(qū)有效應(yīng)力及地下500 m處剪應(yīng)變?cè)隽康淖兓?。從模擬結(jié)果可知，庫(kù)首區(qū)的最大主應(yīng)力和最小主應(yīng)力都隨著水位的上升而增大;與蓄水前相比，蓄水深度為135、145和177 m時(shí)，最大主應(yīng)力增幅分別為5.8%、13.1%和16.8%，最小主應(yīng)力增幅分別為20.5%、20.5%和32.6%;主應(yīng)力數(shù)值與蓄水位存在正相關(guān)性，蓄水位的增大對(duì)地震的發(fā)生起到了催化和推動(dòng)作用，提高了地震風(fēng)險(xiǎn)性。研究結(jié)果可為三峽水庫(kù)誘發(fā)型地震預(yù)測(cè)提供參考。

關(guān)鍵詞： 三峽水庫(kù);FLAC^3D;有效應(yīng)力;剪應(yīng)變?cè)隽?/p>

中圖分類號(hào)： P315.72 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)： 1000-0844（2024）06-1475-09

DOI：10.20000/j.1000-0844.20230207002

Stress field changes and seismic activities before and after the

impoundment in the head area of the Three Gorges

Reservoir based on the finite-difference method

WANG Renlong^1，2， ZHANG Lili^1，3， ZHANG Meidong¹， QIN Weibing²，

YAO Yunsheng¹， WANG Yun¹， DAI Yiming¹， HU Caixiong¹

（1. Institute of Disaster Prevention， Sanhe 065201， Hebei， China;

2. China Three Gorges Corporation， Beijing 430010， China;

3. Hebei Provincial Key Laboratory of Earthquake Dynamics， Sanhe 065201， Hebei， China）

Abstract： The changes in the structural stress field before and after reservoir impoundment are of considerable importance in the study of reservoir-induced seismicity. This study focuses on the Jiuwanxi， Xiannvshan， and Gaoqiao faults in the head area of the Three Gorges Reservoir to investigate the stress field changes and the correlation between reservoir impoundment and earthquakes in the Three Gorges Reservoir. The water level change over 72 months since the impoundment in 2003 is selected to simulate the effect of reservoir water loads during different impoundment periods on the effective stress and shear strain increment in the study area. The true 3D geological modeling is implemented using the software RHINO， and the changes in the effective stress of the study area and the shear strain increment at a depth of 500 m underground are simulated using the finite-difference software FLAC^3D. Simulation results show that the maximum and minimum principal stresses in the head area increase with the rise in water levels. Compared with those before impoundment， the maximum principal stress is increased by 5.8%， 13.1%， and 16.8% when the impoundment depth is 135， 145， and 177 m， respectively， and the minimum principal stress is increased by 20.5%， 20.5%， and 32.6%， respectively. A positive correlation exists between the principal stress value and the water level， and the change in water level plays a catalytic role in the occurrence of earthquakes. This study can provide a reference for earthquake predictions in the Three Gorges Reservoir.

Keywords： Three Gorges Reservoir; FLAC^3D; effective stress; shear strain increment

0 引言三峽水庫(kù)自2003年5月19日開(kāi)始蓄水，至2003年9月停止，期間共記錄大小地震2 216次。水位達(dá)到125 m時(shí)觸發(fā)了2.1級(jí)地震;水位達(dá)到135 m時(shí)出現(xiàn)了小震群活動(dòng);當(dāng)水位保持在135 m時(shí)，出現(xiàn)了不同間隔時(shí)間的小震。6月7—15日，地震日頻次呈上升趨勢(shì)，15—24日地震頻次呈衰減趨勢(shì)。顯然，地震活動(dòng)與水庫(kù)蓄水進(jìn)程^［1-4］有關(guān)，水庫(kù)地震頻度、強(qiáng)度和庫(kù)水位綜合統(tǒng)計(jì)如圖1所示。目前世界公認(rèn)最早的水庫(kù)地震是發(fā)生于1931年的希臘馬拉松水庫(kù)地震^［5］，兩年后阿爾及利亞又發(fā)生富達(dá)水庫(kù)地震。自此，水庫(kù)地震逐漸被重視起來(lái)，在之后的幾十年里，國(guó)內(nèi)外學(xué)者對(duì)此進(jìn)行了廣泛研究，提出了多種研究方法，如地震監(jiān)測(cè)、地質(zhì)勘探、實(shí)驗(yàn)室模擬和遙感技術(shù)等。近年來(lái)，數(shù)值模擬在水庫(kù)誘發(fā)地震的研究中得到了長(zhǎng)足發(fā)展。20世紀(jì)80年代，數(shù)值模擬^［6］主要采用有限元法和有限差分法模擬水庫(kù)蓄水對(duì)地下應(yīng)力場(chǎng)和地下水流場(chǎng)的影響，進(jìn)而推測(cè)對(duì)地震活動(dòng)產(chǎn)生影響的可能性;20世紀(jì)90年代，數(shù)值模擬開(kāi)始由靜態(tài)模型向動(dòng)態(tài)模型轉(zhuǎn)變，在考慮水庫(kù)蓄水的時(shí)間序列并分析水位升降對(duì)地震活動(dòng)影響的同時(shí)，開(kāi)始考慮地震對(duì)地下流體壓力的反饋?zhàn)饔?21世紀(jì)初，學(xué)者們開(kāi)始將地質(zhì)學(xué)、地球物理學(xué)和水文學(xué)等多學(xué)科知識(shí)融入數(shù)值模擬中，并建立更為復(fù)雜和真實(shí)的地下模型，以提高模擬結(jié)果的準(zhǔn)確性。隨著計(jì)算技術(shù)的不斷發(fā)展和高性能計(jì)算設(shè)備的普及，數(shù)值模擬的計(jì)算效率和模擬精度得到了大幅提升。近年來(lái)，學(xué)者們逐漸將不同類型的模型進(jìn)行耦合^［7］，并建立更為復(fù)雜的多物理場(chǎng)模型，進(jìn)一步提高數(shù)值模擬的真實(shí)性和可靠性。

傳統(tǒng)的水庫(kù)地震模擬主要關(guān)注土體的動(dòng)力響應(yīng)，忽視了水體和結(jié)構(gòu)體的相互作用。本文基于有限差分法^［8］，優(yōu)化模型網(wǎng)格，建立更加復(fù)雜的地質(zhì)模型;同時(shí)，結(jié)合地質(zhì)學(xué)和水文學(xué)等學(xué)科知識(shí)，將水土結(jié)構(gòu)相互作用考慮在內(nèi)，采用多物理場(chǎng)耦合，從多個(gè)角度提高模擬的精確性，并從數(shù)值解的角度分析蓄水后研究區(qū)應(yīng)力場(chǎng)及剪應(yīng)變?cè)隽康淖兓?，以反映地震與蓄水的相關(guān)性。

1 地質(zhì)背景及相關(guān)理論

1.1 地質(zhì)背景

三峽水庫(kù)庫(kù)首區(qū)^［9］位于揚(yáng)子準(zhǔn)地臺(tái)上揚(yáng)子臺(tái)褶帶上，范圍從長(zhǎng)江兩岸南北向外延20～30 km，東至宜昌市夷陵區(qū)三斗坪下游5 km，西至重慶市巫山縣境內(nèi)的巴東上游20 km。研究區(qū)分布有多條斷裂，分別是九畹溪斷裂、仙女山斷裂和高橋斷裂等。九畹溪斷裂帶位于黃陵背斜的南側(cè)，總長(zhǎng)31 km，由東、西兩條平行斷裂組成;斷層帶寬2～5 m，穿切寒武系至下二疊系及下白堊系地層;斷裂破碎帶主要由糜棱巖、構(gòu)造巖、構(gòu)造角礫巖、碎裂巖和斷層泥等構(gòu)成，斷層泥結(jié)構(gòu)上相對(duì)不均勻，粗砂含量約40％，黏粒含量14％，變形強(qiáng)度較大。仙女山斷裂帶位于黃陵背斜的西南側(cè)，全長(zhǎng)超過(guò)80 km;其地形呈右行左階展布，形成將軍山和花橋場(chǎng)兩個(gè)拉分性質(zhì)的階區(qū);斷裂帶切割長(zhǎng)陽(yáng)復(fù)背斜和清江河谷，整體走向?yàn)镹W340°～350°，傾角為60°～80°;斷裂帶有明顯分段性，從衛(wèi)星影像上可以看出明顯的線性構(gòu)造，與九畹溪斷裂帶在老林河水電站附近交匯。高橋斷裂帶位于秭歸盆地的西北緣，全長(zhǎng)近40 km，切割上古生界至三疊系地層;斷裂主體走向?yàn)镹E50°，傾向SE，傾角50°～80°;該斷裂帶具有明顯的分段現(xiàn)象，可大致分為北東段和南西段兩個(gè)部分;高橋斷裂帶從地貌角度分析，反差明顯，形成嶺埡、凹槽或溝谷，且沿?cái)嗔褞в斜婪e物發(fā)育或水系斷錯(cuò)。三峽庫(kù)首區(qū)區(qū)域地質(zhì)構(gòu)造如圖2所示。

1.2 相關(guān)理論

目前，國(guó)內(nèi)外學(xué)者大多采用太沙基有效應(yīng)力原理來(lái)研究庫(kù)水誘發(fā)地震^［10］。根據(jù)固體力學(xué)原理，當(dāng)莫爾應(yīng)力圓位于莫爾強(qiáng)度曲線的右側(cè)時(shí)，固體不會(huì)產(chǎn)生破裂，即不會(huì)產(chǎn)生地震;當(dāng)莫爾應(yīng)力圓觸碰到莫爾強(qiáng)度曲線時(shí)，固體就會(huì)產(chǎn)生破裂引發(fā)地震。如圖3所示的應(yīng)力圓中，其右側(cè)是安全區(qū)，左側(cè)是危險(xiǎn)區(qū)。假設(shè)地層巖石為多孔介質(zhì)，孔隙中充滿了流體，地層巖石受外應(yīng)力（地應(yīng)力）和內(nèi)應(yīng)力（孔隙壓力）的共同作用。地層巖石受到的有效應(yīng)力就是地應(yīng)力減去一部分孔隙壓力。其中，雙重有效應(yīng)力^［11］中的結(jié)構(gòu)有效應(yīng)力為：

σ^s_eff=σ-?_cp （1）

式中：σ^s_eff為結(jié)構(gòu)有效應(yīng)力（單位：MPa）;σ為應(yīng)力（單位：MPa）;?_c為觸點(diǎn)孔隙度，其數(shù)值大于本體孔隙度，但小于1;p為介質(zhì)內(nèi)壓（單位：MPa）。

用有效應(yīng)力繪制的應(yīng)力圓與用地應(yīng)力繪制的應(yīng)力圓，大小相同，位置不同。由于有效應(yīng)力通常小于地應(yīng)力，因此有效應(yīng)力的應(yīng)力圓會(huì)位于地應(yīng)力的應(yīng)力圓左側(cè)。

如果向地層注水，即提高孔隙壓力，有效應(yīng)力會(huì)減小，如圖3所示，應(yīng)力圓會(huì)左移使得地層巖石傾向于破裂，地震發(fā)生的危險(xiǎn)性增強(qiáng)，這就是注水誘發(fā)地震^［12］的基本原理。當(dāng)然，注水能否誘發(fā)地震還要看注水增壓的幅度和應(yīng)力圓的初始位置，如果應(yīng)力圓的初始位置遠(yuǎn)離莫爾強(qiáng)度曲線且十分靠右，而注水增壓的幅度較小，因而不會(huì)誘發(fā)地震;如果地層應(yīng)力圓的初始位置在莫爾強(qiáng)度曲線附近，而注水增壓的幅度較大，則會(huì)誘發(fā)地震。

2 模型建立

根據(jù)三峽庫(kù)首區(qū)1∶50 000的數(shù)字地圖，運(yùn)用RHINO軟件建立涵蓋地表單元的三維模型，如圖4所示。模型厚度取自水庫(kù)蓄水前的水位至地殼深度12 km處，共分為6層。模型采用六面體^［13］的計(jì)算單元，共計(jì)247 160個(gè)單元。與傳統(tǒng)的四面體計(jì)算單元相比，六面體計(jì)算單元更為規(guī)則，擁有更好外觀的同時(shí)不易發(fā)生畸變，且精度更高。此外，使用六面體計(jì)算的單元數(shù)遠(yuǎn)小于四面體計(jì)算的單元數(shù)，它在斷裂處使用非結(jié)構(gòu)化網(wǎng)格，從計(jì)算的角度來(lái)看具有更大的優(yōu)勢(shì)。

3 數(shù)值模擬

基于RHINO的建模功能，將涵蓋高程信息的三維地層數(shù)據(jù)轉(zhuǎn)化為三維立體模型。使用Griddle 插件生成體網(wǎng)格，并將網(wǎng)格輸入到FLAC^3D中，利用FLAC^3D內(nèi)置的分組功能賦予模型不同的地層、河流及斷層參數(shù)。以蓄水后庫(kù)區(qū)重力場(chǎng)、形變場(chǎng)、介質(zhì)變化及地震活動(dòng)為約束，模擬分析不同蓄水期（蓄水前、135 m、145 m和177 m）庫(kù)水荷載對(duì)九畹溪斷裂、仙女山斷裂和高橋斷裂區(qū)域地震活動(dòng)的影響。通過(guò)FLAC^3D軟件模擬不同水位時(shí)地下500 m處剪應(yīng)變?cè)隽?sup>［14^］，深入探討研究區(qū)典型斷裂誘震機(jī)制^［15］及地震成因。對(duì)模型表面進(jìn)行分組并設(shè)置邊界條件，將四周及底部設(shè)置為不透水邊界，在不同蓄水深度下對(duì)庫(kù)區(qū)施加相對(duì)應(yīng)的荷載。

當(dāng)模型達(dá)到平衡時(shí)，分析不同水位荷載下主應(yīng)力與剪應(yīng)變?cè)隽考捌浞植继攸c(diǎn)，探討蓄水誘發(fā)地震^［16］的可能性。

地質(zhì)分層及巖層力學(xué)參數(shù)如表1所列。值得注意的是FLAC^3D中的滲透系數(shù)k與一般土力學(xué)中滲透系數(shù)K的概念不同，F(xiàn)LAC^3D中k的國(guó)際單位是m²/（Pa·s），與土力學(xué)中滲透系數(shù)K（cm/s）之間存在如下?lián)Q算關(guān)系：

k［m²/（Pa·s）］=K（cm /s）×1.02×10^-6（2）

因此，在FLAC^3D中需要將實(shí)驗(yàn)獲得的土體滲透系數(shù)參數(shù)乘以1.02×10^-6，才能用于計(jì)算。

斷裂內(nèi)部的力學(xué)參數(shù)比較復(fù)雜，所以根據(jù)《巖石力學(xué)參數(shù)手冊(cè)》^［17］，結(jié)合庫(kù)首區(qū)的地質(zhì)和水文地質(zhì)條件，確定了表2所列的斷裂力學(xué)參數(shù)。

3.1 計(jì)算模式選擇

在進(jìn)行流固分析^［18］時(shí)，有兩種計(jì)算模式：一種是無(wú)滲流模式，在這種模式下孔隙水壓力（孔壓）保持不變，土體單元的屈服判斷由有效應(yīng)力決定;另外一種是滲流模式，這種模式求解方法較多，孔隙水壓力會(huì)隨著浸潤(rùn)線的改變而改變。滲流模式和無(wú)滲流模式不是互相獨(dú)立的，而是可以同時(shí)考慮的。本文模擬水位變化對(duì)巖土體的影響時(shí)，先使用滲流模式計(jì)算水流的變化，再使用非滲流模式計(jì)算巖土體的力學(xué)響應(yīng)。

3.2 模型加載方式

選取2003年開(kāi)始蓄水之后72個(gè)月的水位變化數(shù)據(jù)做研究（圖5）。根據(jù)不同蓄水期，分三個(gè)階段分析：第一階段為2003年6月—2006年9月，水庫(kù)開(kāi)始蓄水至135 m并保持該水位運(yùn)行，此階段為穩(wěn)定水位荷載;第二階段為2006年9月—2008年9月，水位從145 m蓄水至156 m，每年在高水位運(yùn)行半年之后，逐漸下降至145 m的低水位，在156 m蓄水末期曾試驗(yàn)性蓄水至172 m，為175 m試驗(yàn)性蓄水做準(zhǔn)備;第三階段為2008年9月21日—2020年12月，可以視作175 m試驗(yàn)性蓄水。因此，每年的水位變化過(guò)程是從145 m的低水位試驗(yàn)性蓄水至接近175 m的高水位，在高水位運(yùn)行半年之后，再將水位下降至145 m的循環(huán)過(guò)程。為了探究三峽水庫(kù)蓄水位的最大范圍，在175 m的基礎(chǔ)上再增加2 m的蓄水高度。177 m水位具有一定的研究意義，因此對(duì)177 m的蓄水情況進(jìn)行了模擬。

蓄水期時(shí)，三峽水庫(kù)的水位變化波動(dòng)較大，每年水位落差可達(dá)幾十米，并且具有周期性和循環(huán)性的特點(diǎn)?？紤]到水庫(kù)水位短時(shí)間內(nèi)不會(huì)發(fā)生突變，所以設(shè)置加載方式為靜荷載。在靜荷載的加載方式下，默認(rèn)水庫(kù)水位不變，即荷載穩(wěn)定。無(wú)滲流模式下，在水庫(kù)上方施加庫(kù)水的重力荷載，這時(shí)需要考慮構(gòu)造應(yīng)力場(chǎng)，將水位荷載直接作用于網(wǎng)格上，并通過(guò)設(shè)定材料參數(shù)來(lái)考慮水與土體之間的相互作用。在上述基礎(chǔ)上設(shè)置孔隙度、滲透系數(shù)等參數(shù)，模擬庫(kù)水對(duì)土體的荷載。

3.3 模擬結(jié)果分析

在靜荷載的加載方式下，不同蓄水期庫(kù)首區(qū)的最大主應(yīng)力^［19］如圖6所示。由圖6可得，蓄水前庫(kù)首區(qū)的最大主應(yīng)力在-1.208 6～-95.038 MPa之間;當(dāng)蓄水深度為135 m時(shí)，最大主應(yīng)力有所增大，為-2.232 4～-100.58 MPa;當(dāng)蓄水深度為145 m時(shí)，最大主應(yīng)力較135 m深度增幅不大，為-2.407 6～-107.52 MPa;當(dāng)蓄水深度達(dá)到177 m時(shí)，最大主應(yīng)力變?yōu)?3.056 5～-111.04 MPa之間。壓應(yīng)力的最大值較蓄水前增大了16.002 MPa，增幅為16.8%。同時(shí)，不同蓄水期庫(kù)首區(qū)的最小主應(yīng)力也會(huì)發(fā)生變化，如圖7所示。蓄水前庫(kù)區(qū)的最小主應(yīng)力范圍在-4.005 2～-263.54 MPa;當(dāng)蓄水深度為135 m時(shí)，最小主應(yīng)力范圍為-0.430～-317.62 MPa，增大了54.08 MPa；在145 m水位，最小主應(yīng)力范圍為-0.432～-328.64 MPa，增大了11.02 MPa;當(dāng)蓄水深度達(dá)到177 m時(shí)，最小主應(yīng)力范圍為-0.478 69～-349.57 MPa。最小主應(yīng)力的增幅達(dá)到了86.01 MPa，增長(zhǎng)率為32.6%，并且最大主應(yīng)力及最小主應(yīng)力的分布遵循從上到下依次增大的規(guī)律。這說(shuō)明蓄水之后，最小主應(yīng)力絕對(duì)值的增幅大于最大主應(yīng)力，即在蓄水之后區(qū)域的剪應(yīng)力會(huì)明顯增大。也就是說(shuō)，蓄水會(huì)導(dǎo)致庫(kù)區(qū)原有地層特性發(fā)生變化^［20］，這種變化會(huì)打破原有的平衡狀態(tài)，并達(dá)到新的應(yīng)力平衡狀態(tài)^［21］。

隨著庫(kù)首區(qū)主應(yīng)力的變化，剪應(yīng)變?cè)隽糠植家舶l(fā)生了變化。不同蓄水期庫(kù)首區(qū)地表以下500 m處的剪應(yīng)變?cè)隽糠植既鐖D8所示。由圖8可知，蓄水前剪應(yīng)變?cè)隽繛?.012 7；當(dāng)蓄水深度在135 m時(shí)，剪應(yīng)變?cè)隽繛?.015 5;在蓄水深度為145 m時(shí)，剪應(yīng)變?cè)隽繛?.016 1;在蓄水深度達(dá)177 m時(shí)，應(yīng)變?cè)隽繛?.021 2，即剪應(yīng)變?cè)隽侩S著水位的上升而逐漸增大。從切片云圖看，地表以下500 m深度處剪應(yīng)變?cè)隽繋ы樓衅缑嫜由?sup>［22^］，大體集中分布在水庫(kù)、仙女山斷裂、九畹溪斷裂及高橋斷裂處，表現(xiàn)為多個(gè)分支。隨著蓄水深度不斷增加，庫(kù)水滲透加劇，導(dǎo)致土壤以及巖層的含水量增多，水的滲透和擴(kuò)散會(huì)使周?chē)鷰r層發(fā)生一定程度的軟化，最終導(dǎo)致原本不穩(wěn)定的斷層體系更加不穩(wěn)定。隨著斷層注水的增多，剪應(yīng)變?cè)隽恐饾u增大。如果斷層原本處于一種平衡狀態(tài)，注水后增加的水壓力可能會(huì)打破這種平衡，使得斷層重新運(yùn)動(dòng)，這樣就會(huì)導(dǎo)致斷層內(nèi)部結(jié)構(gòu)平衡被打破，從而誘發(fā)地震。借鑒剪應(yīng)變?cè)隽吭诙逊e體邊坡穩(wěn)定性研究中的應(yīng)用^［23］，剪應(yīng)變?cè)隽繋г跀鄬蛹八畮?kù)下方連續(xù)性較好，而在遠(yuǎn)離壩區(qū)區(qū)域的連續(xù)性較差，剪應(yīng)變?cè)隽繋У姆植挤秶∏∈亲畈环€(wěn)定的區(qū)域，也是最容易發(fā)生失穩(wěn)破壞的。水庫(kù)庫(kù)區(qū)失穩(wěn)，一般都是沿著剪應(yīng)變最大的部位發(fā)生，大量工程實(shí)例也驗(yàn)證了這一點(diǎn)。

4 結(jié)論

通過(guò)分析三峽庫(kù)首區(qū)蓄水前后應(yīng)力場(chǎng)變化與地震活動(dòng)的關(guān)系，得到以下結(jié)論：

（1） 從模擬結(jié)果分析可知，庫(kù)首區(qū)的最大主應(yīng)力和最小主應(yīng)力都會(huì)隨著水位的上升而增大。與蓄水前相比，蓄水深度分別為135、145和177 m時(shí)，最大主應(yīng)力增幅分別為5.8%、13.1%和16.8%，最小主應(yīng)力增幅分別為20.5%、20.5%和32.6%。主應(yīng)力數(shù)值與蓄水位存在正相關(guān)性，蓄水位的變化對(duì)地震的發(fā)生起到了催化作用。

（2） 不同蓄水期的剪應(yīng)變?cè)隽吭诰蓬迪獢嗔?、仙女山斷裂和高橋斷裂較為集中，庫(kù)水荷載對(duì)斷層產(chǎn)生的剪應(yīng)力的直接作用極大地增加了地震發(fā)生的可能性。

（3） 庫(kù)首區(qū)受到滲流的影響增大，考慮到水庫(kù)的荷載增多等因素，使得水庫(kù)周?chē)牡孛姘l(fā)生彈性變形。蓄水導(dǎo)致庫(kù)首區(qū)應(yīng)力場(chǎng)發(fā)生變化，體現(xiàn)在地面垂直應(yīng)力增加的同時(shí)水平應(yīng)力減小，對(duì)地下水和地質(zhì)構(gòu)造產(chǎn)生影響;庫(kù)水荷載產(chǎn)生的附加應(yīng)力疊加在構(gòu)造應(yīng)力場(chǎng)的間接作用，對(duì)誘發(fā)三峽庫(kù)首區(qū)水庫(kù)地震也起到了推動(dòng)作用。

（4） 本研究沒(méi)有考慮溫度的影響。實(shí)際上，溫度也會(huì)影響地下流體的滲流，從而影響庫(kù)區(qū)的應(yīng)力分布。因此，在以后的研究中考慮溫度的影響會(huì)使研究結(jié)果更精確。

參考文獻(xiàn)（References）

［1］王孔偉，路永強(qiáng)，聶進(jìn)，等.三峽庫(kù)區(qū)仙女山和九畹溪斷裂帶水庫(kù)地震變化規(guī)律［J］.吉林大學(xué)學(xué)報(bào)（地球科學(xué)版），2021，51（2）：624-637.

WANG Kongwei，LU Yongqiang，NIE Jin，et al.Earthquake variation law of Xiannüshan and Jiuwanxi fault zones in Three Gorges Reservoir area［J］.Journal of Jilin University （Earth Science Edition），2021，51（2）：624-637.

［2］何超楓.三峽庫(kù)首區(qū)仙女山—九畹溪斷裂帶活化特征研究［D］.武漢：中國(guó)地震局地震研究所，2016.

HE Chaofeng.Study on fault reactivation feature of Xiannvshan—Jiuwanxi fault zone in the head area of the Three Gorges［D］.Wuhan：Institute of seismology，China Earthquake Administration，2016.

［3］ZHANG Z X，ZHANG H J，HAN L Q，et al.Multi-slip surfaces searching method for earth slope with weak interlayer based on local maximum shear strain increment［J］.Computers and Geotechnics，2022，147：104760.

［4］朱家正，孫玉軍.四川大崗山水庫(kù)蓄水對(duì)地震活動(dòng)影響的數(shù)值模擬研究［J］.地球物理學(xué)報(bào)，2022，65（10）：3930-3943.

ZHU Jiazheng，SUN Yujun.Numerical simulation of the effect of water storage on seismic activity in Dagangshan Reservoir，Sichuan，China［J］.Chinese Journal of Geophysics，2022，65（10）：3930-3943.

［5］楊清源，胡毓良，陳獻(xiàn)程，等.國(guó)內(nèi)外水庫(kù)誘發(fā)地震目錄［J］.地震地質(zhì)，1996，18（4）：453-461.

YANG Qingyuan，HU Yuliang，CHEN Xiancheng，et al.Catalogue of reservoir-induced earthquakes at home and abroad［J］.Seismology and Geology，1996，18（4）：453-461.

［6］龔鋼延.孔隙壓力擴(kuò)散與水庫(kù)誘發(fā)地震活動(dòng)性的初步研究［J］.西北地震學(xué)報(bào)，1988，10（1）：52-59.

GONG Gangyan.Primary research on pore pressure diffusion with reservoir-induced seismicity［J］.Northwestern Seismological Journal，1988，10（1）：52-59.

［7］虞金鑫，王秋良，趙艷南.基于流固耦合理論的三峽庫(kù)首區(qū)蓄水后附加應(yīng)力場(chǎng)研究［J］.大地測(cè)量與地球動(dòng)力學(xué)，2020，40（9）：928-930，941.

YU Jinxin，WANG Qiuliang，ZHAO Yannan.Stress field of the Three Gorges based on fluid-solid coupling simulation［J］.Journal of Geodesy and Geodynamics，2020，40（9）：928-930，941.

［8］李保民.基于有限差分法的水庫(kù)透水性能研究［J］.海河水利，2023（5）：77-80.

LI Baomin.Research on permeability performance of reservoir based on finite difference method［J］.Haihe Water Resources，2023（5）：77-80.

［9］王秋良，張麗芬，廖武林，等.三峽庫(kù)首區(qū)斷裂構(gòu)造與地震活動(dòng)特征［J］.大地測(cè)量與地球動(dòng)力學(xué)，2013，33（5）：29-33.

WANG Qiuliang，ZHANG Lifen，LIAO Wulin，et al.Fault tectonics and seismic activity characteristics of Three Gorges Reservoir［J］.Journal of Geodesy and Geodynamics，2013，33（5）：29-33.

［10］許忠淮，閻明，趙仲和.由多個(gè)小地震推斷的華北地區(qū)構(gòu)造應(yīng)力場(chǎng)的方向［J］.地震學(xué)報(bào)，1983，5（3）：268-279.

XU Zhonghuai，YAN Ming，ZHAO Zhonghe.Evaluation of the direction of tectonic stress in North China from recorded data of a large number of small earthquakes［J］.Acta Seismologica Sinica，1983，5（3）：268-279.

［11］李傳亮，朱蘇陽(yáng).再談雙重有效應(yīng)力：對(duì)《雙重有效應(yīng)力再認(rèn)識(shí)及其綜合作用》一文的討論與分析［J］.石油科學(xué)通報(bào)，2019，4（4）：414-429.

LI Chuanliang，ZHU Suyang.Discussion of dual effective stresses in porous media：analysis on the paper of re-recognition of dual effective stresses and the comprehensive effect［J］.Petroleum Science Bulletin，2019，4（4）：414-429.

［12］周永利.基于AutoCAD與ANSYS建立FLAC^3D模型的邊坡穩(wěn)定性分析［J］.中國(guó)科技信息，2009（19）：90-92.

ZHOU Yongli.Slope stability analysis of FLAC^3D model based on AutoCAD and ANSYS［J］.China Science and Technology Information，2009（19）：90-92.

［13］張玉新，鐘毅，王華昆，等.塑性成形中的六面體網(wǎng)格劃分與重劃分技術(shù)［J］.鍛壓技術(shù)，2004，29（6）：33-35，61.

ZHANG Yuxin，ZHONG Yi，WANG Huakun，et al.Meshing and remeshing techniques of hexahedral elements in plastic forming［J］.Forging amp; Stamping Technology，2004，29（6）：33-35，61.

［14］李井岡，張麗芬，廖武林，等.褶皺構(gòu)造中的地震：2017年三峽庫(kù)區(qū)巴東M4.3地震序列成因討論［J］.地球物理學(xué)報(bào)，2018，61（9）：3701-3712.

LI Jinggang，ZHANG Lifen，LIAO Wulin，et al.Earthquakes in the fold structure：the genesis of the M4.3 earthquake sequence in the Three Gorges Reservoir area in 2017［J］.Chinese Journal of Geophysics，2018，61（9）：3701-3712.

［15］劉港，趙海軍，馬鳳山，等.裂隙巖體滲流應(yīng)力耦合研究現(xiàn)狀與展望＼[C＼]//2017年全國(guó)工程地質(zhì)學(xué)術(shù)年會(huì)論文集.桂林：《工程地質(zhì)學(xué)報(bào)》編輯部，2017：235-247.

LIU Gang，ZHAO Haijun，MA Fengshan，et al.Present situation and prospect of seepage stress coupling research in fractured rock mass［C］//Proceedings of the 2017 National Academic Conference on Engineering Geology.Guilin：Editorial Department of Journal of Engineering Geology，2017：235-247.

［16］孟慶筱，姚運(yùn)生，廖武林，等.三峽蓄水進(jìn)程中庫(kù)首區(qū)地震活動(dòng)與庫(kù)水位的關(guān)聯(lián)性研究［J］.大地測(cè)量與地球動(dòng)力學(xué)，2021，41（7）：714-720.

MENG Qingxiao，YAO Yunsheng，LIAO Wulin，et al.Analysis of the cross-correlation between seismicity and water level at the head area of the Three Gorges Reservoir in the impoundment process［J］.Journal of Geodesy and Geodynamics，2021，41（7）：714-720.

［17］水利水電科學(xué)研究院，水利水電規(guī)劃設(shè)計(jì)院，水利水電情報(bào)研究所，等.巖石力學(xué)參數(shù)手冊(cè)＼[M＼].北京：水利電力出版社，1991.

Institute of Water Resources and Hydropower Research，Renewable Energy Engineering Institute，Water Resources and Hydropower Information Research Institute，et al.Handbook of rock mechanics parameters＼[M＼].Beijing：Water Resources and Electric Power Press，1991.

［18］謝向東.基于流固耦合理論的丹江口水庫(kù)誘發(fā)地震機(jī)理研究［D］.武漢：武漢大學(xué)，2010.

XIE Xiangdong.Study on mechanism of reservoir-induced seismicity based on theory of fluid-solid coupling in the Danjiangkou Reservoir area［D］.Wuhan：Wuhan University，2010.

［19］YAO Y S，WANG Q L，LIAO W L，et al.Influences of the Three Gorges Project on seismic activities in the reservoir area［J］.Science Bulletin，2017，62（15）：1089-1098.

［20］張紅艷，謝富仁，荊振杰.京西北盆嶺構(gòu)造區(qū)現(xiàn)代構(gòu)造應(yīng)力場(chǎng)的非均勻特征［J］.地球物理學(xué)報(bào)，2009，52（12）：3061-3071.

ZHANG Hongyan，XIE Furen，JING Zhenjie.Research on heterogeneity of the present tectonic stress field in the basin-and-range province northwest of Beijing［J］.Chinese Journal of Geophysics，2009，52（12）：3061-3071.

［21］常婷.三峽秭歸M4.5、吉林松原M5.5地震序列庫(kù)侖應(yīng)力變化的對(duì)比研究［D］.武漢：中國(guó)地震局地震研究所，2021.

CHANG Ting.Comparative study on Coulomb stress change of Three Gorges Zigui M4.5 and Jilin Songyuan M5.5 seismic sequence［D］.Wuhan：Institute of seismology，China Earthquake Administration，2021.

［22］CUI G J，ZHANG C Q，YE J P，et al.Influences of dynamic normal disturbance and initial shear stress on fault activation characteristics［J］.Geomechanics and Geophysics for Geo-Energy and Geo-Resources，2022，8（5）：159.

［23］丁秀美，劉光士，黃潤(rùn)秋，等.剪應(yīng)變?cè)隽吭诙逊e體邊坡穩(wěn)定性研究中的應(yīng)用［J］.地球科學(xué)進(jìn)展，2004，19（增刊1）：318-323.

DING Xiumei，LIU Guangshi，HUANG Runqiu，et al.Application of shear strain increment in stability study of accumulation slope［J］.Advances in Earth Science，2004，19（Suppl01）：318-323.

（本文編輯：賈源源）