滄州河間構(gòu)造地裂縫應(yīng)力場(chǎng)二維平面模型研究

宋 偉，馬學(xué)軍，呂鳳蘭，王兵虎，邵長(zhǎng)慶

（中國(guó)地質(zhì)調(diào)查局水文地質(zhì)環(huán)境地質(zhì)調(diào)查中心，河北·保定 071051）

滄州河間構(gòu)造地裂縫應(yīng)力場(chǎng)二維平面模型研究

宋 偉，馬學(xué)軍，呂鳳蘭，王兵虎，邵長(zhǎng)慶

（中國(guó)地質(zhì)調(diào)查局水文地質(zhì)環(huán)境地質(zhì)調(diào)查中心，河北·保定 071051）

通過(guò)求解平面應(yīng)力的問(wèn)題的有限單元方法反演了滄州河間區(qū)域構(gòu)造應(yīng)力場(chǎng)，并確定了最大住應(yīng)力場(chǎng)和最大剪應(yīng)力集中區(qū)，從而探索分析構(gòu)造地裂縫的潛在危險(xiǎn)區(qū)；建立了地裂縫二維平面模型，分析河間地區(qū)區(qū)域構(gòu)造應(yīng)力場(chǎng)以及區(qū)域應(yīng)力集中區(qū)和地裂縫的對(duì)應(yīng)關(guān)系。

地裂縫；構(gòu)造應(yīng)力場(chǎng)；有限元分析；二維平面模型

滄州是華北地裂縫分布較為密集的地區(qū)，地下水開(kāi)采造成地下水位下降是地裂縫形成主要誘因之一[1～3]，其與區(qū)域內(nèi)大陸動(dòng)力學(xué)基本背景特征尤其是與盆地內(nèi)深部隱伏活動(dòng)斷裂的持續(xù)活動(dòng)也密切相關(guān)[4,5]。揭示該地區(qū)地裂縫分布與構(gòu)造應(yīng)力場(chǎng)關(guān)系是地裂縫形成機(jī)理研究、實(shí)現(xiàn)地裂縫預(yù)測(cè)預(yù)報(bào)的關(guān)鍵因素之一[6]。

通過(guò)數(shù)值模擬方法，綜合利用地質(zhì)、地球物理、地球化學(xué)等研究成果，從數(shù)學(xué)方法上證實(shí)或證偽已有的地質(zhì)模型，在滄州地裂縫分布區(qū)建立符合區(qū)域地質(zhì)特征的、可說(shuō)明地裂縫分布與構(gòu)造應(yīng)力場(chǎng)關(guān)系的數(shù)學(xué)力學(xué)模型。對(duì)所建立的模型的力學(xué)機(jī)制和動(dòng)力學(xué)過(guò)程進(jìn)行分析探討與驗(yàn)證，以合理再現(xiàn)復(fù)雜地質(zhì)構(gòu)造和事件的演化過(guò)程，并與地質(zhì)—地球物理—地球化學(xué)方面的資料進(jìn)行聯(lián)系對(duì)比，從而獲得地裂縫分布動(dòng)力學(xué)的定量化理解[7]。

1 滄州河間地裂縫概況

經(jīng)實(shí)地調(diào)查表明，河間市地裂縫共有87處，其中地面塌陷9處。河間地區(qū)的地裂縫長(zhǎng)度多在2km以內(nèi)，個(gè)別超過(guò)3km。其中最長(zhǎng)的為留古寺鎮(zhèn)艾家莊村西南發(fā)育的地裂縫，推測(cè)長(zhǎng)度大于3km；其次為留古寺鎮(zhèn)北王莊東約1.5km處發(fā)育的另一條地裂縫，推測(cè)長(zhǎng)度約3km。但在調(diào)查期間，大部分地裂縫已不可見(jiàn)，只在地表有低洼痕跡。多數(shù)地裂縫長(zhǎng)度在100～200m之間，且地表多可見(jiàn)，長(zhǎng)度500m以上地裂縫多數(shù)已不可見(jiàn)。

2 模型區(qū)選擇與建設(shè)方法

2.1 模型區(qū)選擇

建立模型時(shí)充分考慮了地殼淺部結(jié)構(gòu)、介質(zhì)條件的均勻性等多種因素，特別是地質(zhì)構(gòu)造運(yùn)動(dòng)中對(duì)斷裂活動(dòng)的決定作用（表1）[8,9]。數(shù)值模擬研究區(qū)內(nèi)的活動(dòng)斷裂反映了該地區(qū)的構(gòu)造運(yùn)動(dòng)的發(fā)生和發(fā)展歷史及其強(qiáng)度，從而成為不同的地貌和構(gòu)造單元邊界，控制了該地區(qū)區(qū)域構(gòu)造應(yīng)力場(chǎng)的分布和變化特征[10]。

表1 數(shù)值模擬研究區(qū)主要活動(dòng)斷裂簡(jiǎn)表Table1 The main active faults in the study area of numerical simulation

模型區(qū)內(nèi)有兩個(gè)新構(gòu)造活動(dòng)區(qū)，即滄縣斷壟和冀中坳陷，且位于唐山—河間—磁縣地震帶中段，新構(gòu)造活動(dòng)比較活躍。唐山—河間—磁縣地區(qū)為北東向地震帶，具有統(tǒng)一的現(xiàn)代構(gòu)造應(yīng)力場(chǎng)，其主壓應(yīng)力（P）軸和主張應(yīng)力（T）軸近于水平，方向分別是NEE—SWW和NNW—SSE向，中間應(yīng)力（T）軸近直立[11]。該地震帶的地震機(jī)制解譯等資料確定的震源斷層多為N15～35E，傾角都在80以上，以右旋走滑為主，并且有共軛破裂的特點(diǎn)。該地震帶是一條現(xiàn)代右旋剪切錯(cuò)動(dòng)的地殼破裂帶，但其不是一條貫通連續(xù)的破裂帶，而是有一系列具有相同運(yùn)動(dòng)性質(zhì)的斷裂組成。數(shù)值模擬研究區(qū)主要受滄縣斷隆西緣斷裂帶的控制，該地區(qū)斷層總體走向呈NNW走向，該斷裂帶發(fā)育有兩條NNW向和NWW向的次級(jí)斷裂。因此，NE—NNE走向的活動(dòng)斷裂可能是宏觀上影響或控制該地區(qū)地裂縫發(fā)育的主要因素。

2.2 模型建立方法

本次模擬計(jì)算首先以研究區(qū)內(nèi)基本構(gòu)造格局為基礎(chǔ)，建立數(shù)值分析平面模型，分析研究區(qū)內(nèi)的區(qū)域構(gòu)造應(yīng)力場(chǎng)的平面分布特征和規(guī)律（圖1）。研究二維狀態(tài)條件下，區(qū)域構(gòu)造應(yīng)力場(chǎng)的總體特征和規(guī)律，以及活動(dòng)斷裂對(duì)地裂縫成生、擴(kuò)展和發(fā)育分布的控制關(guān)系[12]。

圖1 平面數(shù)值模擬研究區(qū)位置示意圖Fig.1 Location map of the study area of plane numerical simulation

3 模型的建設(shè)

3.1 二維平面模型的建立

（1）基本假設(shè)

在建立模型時(shí)的四個(gè)基本假設(shè)：平面應(yīng)變狀態(tài)；同一地層巖石為均勻的彈性各向同性材料；斷層面假定為軟弱帶，并按照有限元模型中的接觸面設(shè)定；巖土體服從摩爾—庫(kù)侖破壞準(zhǔn)則。

（2）模型幾何參數(shù)

通過(guò)對(duì)地質(zhì)體結(jié)構(gòu)特征和斷裂帶的幾何形狀進(jìn)行簡(jiǎn)化，建立平面二維有限元計(jì)算模型（圖2）。該模型是參照河間及其周邊地區(qū)的新構(gòu)造運(yùn)動(dòng)方面的資料建立的[13,14]。模型NW—SE向長(zhǎng)115km，設(shè)置為X向，NE—SW向長(zhǎng)95km，設(shè)置為Y向。本次計(jì)算采用摩爾庫(kù)侖本構(gòu)方程，共剖分14593個(gè)計(jì)算單元，14377個(gè)計(jì)算節(jié)點(diǎn)。采用四節(jié)點(diǎn)四邊形平面應(yīng)變單元?jiǎng)澐謳缀文Ｐ停▓D3）。

（3）模型材料力學(xué)參數(shù)

數(shù)值模擬所采用的材料是根據(jù)已有的地質(zhì)資料和數(shù)值模擬資料，用經(jīng)驗(yàn)類比法確定的，主要模擬華北平原第四系地層在斷裂活動(dòng)作用下的水平應(yīng)力場(chǎng)分布特征。模型中將地質(zhì)結(jié)構(gòu)介質(zhì)材料劃分為非斷裂帶（plain）和斷裂帶（fault）兩大類，且將斷裂帶按照軟弱帶處理[15]。通過(guò)上述假設(shè)和處理，就可以通過(guò)控制斷裂物理力學(xué)參數(shù)實(shí)現(xiàn)對(duì)模擬區(qū)變形和應(yīng)力分布的有效控制[16]。除此之外，模型中忽略了各條斷層的斷裂面力學(xué)參數(shù)的差異，模型各介質(zhì)單元具體物理力學(xué)參數(shù)見(jiàn)表2與圖4所示。

圖2 數(shù)值模擬研究區(qū)平面幾何模型圖Fig.2 The plane geometry model of numerical simulation study area

圖3 數(shù)值模擬研究區(qū)平面模型網(wǎng)格剖分圖Fig.3 The mesh generation of plane model in numerical simulation study area

表2 數(shù)值模擬研究區(qū)平面模型材料參數(shù)表Table2 Table of material parameters of plane model in numerical simulation study area

圖4 數(shù)值模擬研究區(qū)平面材料模型圖Fig.4 The plane material model of numerical simulation study area

（4）模型邊界條件

該平面模型中的邊界條件為：左邊界X向固定邊界，約束模型的X向位移；右邊界為-X向拉張應(yīng)力，以模擬數(shù)值模擬區(qū)內(nèi)的NWW—NEE向的拉張區(qū)域構(gòu)造應(yīng)力場(chǎng)的作用；上部邊界為Y向固定邊界，約束模型的Y向位移；下部邊界在Y向自由，不約束位移；同時(shí)上下邊界還有一對(duì)逆時(shí)針的剪力作用，以模擬研究區(qū)內(nèi)東西向的左旋剪切作用（圖5）。

圖5 數(shù)值模擬研究區(qū)平面模型邊界條件圖Fig.5 The boundary conditions of plane model in numerical simulation study area

3.2 計(jì)算結(jié)果及其分析

為了研究方便，我們?cè)谄矫婢W(wǎng)格上布置了兩條剖面線Ⅰ—Ⅰ 和Ⅱ—Ⅱ（圖3），通過(guò)對(duì)剖面Ⅰ—Ⅰ 和Ⅱ—Ⅱ上的應(yīng)力、應(yīng)變和位移曲線變化規(guī)律的分析，揭示模型計(jì)算分析區(qū)相應(yīng)的應(yīng)力、應(yīng)變和位移場(chǎng)特征。

（1）移場(chǎng)特征

通過(guò)對(duì)平面模型X、Y向位移等值線云圖和等值線圖（圖6～圖9）以及總位移位移矢量圖（圖10）的分析，可以得出如下認(rèn)識(shí)。

X向位移等值線圖顯示了模型南部區(qū)域具有向X向正方向的位移的特征，而其北部區(qū)域則具有相對(duì)向X負(fù)向移動(dòng)的趨勢(shì)。該結(jié)果顯示了對(duì)應(yīng)的研究區(qū)北西西向斷裂具有左旋特點(diǎn)（河間—完縣斷裂f8），北東向斷裂具有右旋特點(diǎn)（里坦斷凹東緣斷裂f（5）滄縣斷隆西緣斷裂F17）。Y向位移等值線圖同樣顯示了北西西向斷層具有左旋的特點(diǎn)，北東向斷裂具有右旋特點(diǎn)。

由區(qū)域總位移矢量圖可以看出，模擬區(qū)域整體運(yùn)動(dòng)方向以近北西西—東西向?yàn)橹鳌Ｔ跀嗔褞Ц浇?，位移方向發(fā)生偏轉(zhuǎn)，且斷層兩盤之間的位移量值有明顯的區(qū)別。

圖6 平面模型X向位移云圖Fig.6 The X direction displacement nephogram of plane model

圖7 平面模型X向位移等值線圖Fig.7 The X direction displacement contour of plane model

圖8 平面模型Y向位移云圖Fig.8 The Y direction displacement nephogram of plane model

圖9 平面模型Y向位移等值線圖Fig.9 The Y direction displacement contour of plane model

圖10 平面模型總位移矢量圖Fig.1 0 The total displacement vector of plane model

（2）應(yīng)力分布特征

圖11 平面模型X向應(yīng)力等值線云圖Fig.1 1 The X direction stress curve contour nephogram of plane model

圖12 平面模型X向應(yīng)力等值線圖Fig.1 2 The X direction stress curve contour of plane model

圖13 平面模型Y向應(yīng)力等值線云圖Fig.1 3 The Y direction stress curve contour nephogram of plane model

圖14 平面模型Y向應(yīng)力等值線圖Fig.1 4 The Y direction stress curve contour of plane model

圖15 平面模型最大主應(yīng)力等值線云圖Fig.1 5 The maximum principal stress contour nephogram of plane model

圖16 平面模型最大主應(yīng)力等值線圖Fig.1 6 The maximum principal stress contour of plane model

圖17 最大主應(yīng)力矢量圖Fig.1 7 The maximum principal stress vector of plane model

圖18 最大主應(yīng)力矢量圖(局部放大)Fig.1 8 The maximum principal stress vector of plane model (Enlargement of partial positions)

通過(guò)分析圖11～圖18，可以得到以下認(rèn)識(shí)：

X和Y向應(yīng)力等值線云圖或等值線圖顯示了在斷裂帶的拐點(diǎn)、端點(diǎn)和交匯點(diǎn)易發(fā)生應(yīng)變集中現(xiàn)像，且應(yīng)力變化梯度。例如，徐水-安新斷裂f1與雄縣斷裂交匯區(qū)f2，以及里坦斷凹北緣斷裂（f3）與里坦斷凹西緣斷裂（f4）、里坦斷凹東緣斷裂（f5）交匯區(qū)、斷裂f3的端部、滄縣斷隆西緣斷裂（F17）在文安一帶的拐彎處、等地段，均存在應(yīng)力集中現(xiàn)象。而這些應(yīng)力集中區(qū)域正是實(shí)際當(dāng)中地裂縫的高發(fā)區(qū)域；

模擬區(qū)域最大主應(yīng)力方向總體上呈NW—SE方向，與已有的地震機(jī)制解所得的模擬區(qū)的主拉張應(yīng)力方向基本一致[17]。在斷裂帶交匯處，應(yīng)力方向變化比較復(fù)雜，出現(xiàn)了應(yīng)力的高值異常。此外，在任丘斷裂兩側(cè)是應(yīng)力方向的轉(zhuǎn)變點(diǎn)，任丘斷層（f7）的上盤主拉應(yīng)力是東西向，而在斷層的下盤，主拉張應(yīng)力逐漸變?yōu)镾E—SSE向。最大主拉應(yīng)壓應(yīng)力方向與板塊的絕對(duì)運(yùn)動(dòng)方向存在著相關(guān)性，基本上就是板塊內(nèi)部的縮短方向、地塊的運(yùn)動(dòng)方向或者裂縫帶的走向。

（3）應(yīng)變特征

通過(guò)圖19～圖26，可以得到以下認(rèn)識(shí)：

X向應(yīng)變和等值線圖均顯示了在斷裂帶上應(yīng)變梯度變化大，斷裂帶上的應(yīng)變集中所形成的應(yīng)變量幾乎集中了擬區(qū)內(nèi)的總形變量，而且主要以拉張應(yīng)變?yōu)橹鱗18,19]。這一現(xiàn)象與實(shí)際情況基本相符，活動(dòng)斷裂作為構(gòu)造單元的軟弱帶，是新構(gòu)造運(yùn)動(dòng)最活躍的地帶，幾乎控制了區(qū)域構(gòu)造構(gòu)造應(yīng)力場(chǎng)的應(yīng)變[20,21]。通過(guò)剖面I應(yīng)變曲線可以更清楚地看到這一特征。Y向應(yīng)變特征與此類似，不再贅述。剖面II的應(yīng)力應(yīng)變特征與剖面I類似，其應(yīng)力應(yīng)變集中現(xiàn)象同樣出現(xiàn)在斷裂帶上。

圖19 X向應(yīng)力等值線云圖Fig.1 9 The X direction strain contour nephogram of plane model

圖20 X向應(yīng)變等值線圖Fig.2 0 The X direction strain contour of plane model

圖21 X向應(yīng)變?cè)茍D（極?。〧ig.2 1 The X direction strain nephogram of plane model(minimum)

圖22 X向應(yīng)變等值線圖（極?。〧ig.2 2 The X direction strain contour of plane model(minimum)

圖23 I剖面X向應(yīng)變曲線圖Fig.2 3 The X direction strain curve of I profle

圖24 I剖面Y向應(yīng)變曲線圖Fig.2 4 The Y direction strain curve of I profle

圖25 II剖面應(yīng)變曲線圖Fig.2 5 The strain curve of II profle

圖26 II剖面應(yīng)力曲線圖Fig.2 6 The stress curve of II profle

4 初步結(jié)論

從上述計(jì)算分析可以得到如下結(jié)論：

（1）研究區(qū)發(fā)育的構(gòu)造型地裂縫有兩種類型，一種與斷裂活動(dòng)直接相關(guān)，一種是地殼拉張應(yīng)力作用下的表層變形開(kāi)裂。華北地區(qū)現(xiàn)代地殼處于一致性較好的應(yīng)力場(chǎng)中，其主壓應(yīng)力軸為NEE向，主張應(yīng)力軸為NNW向，主應(yīng)力軸的傾向都不大。在這個(gè)統(tǒng)一的應(yīng)力場(chǎng)作用下，NNE和NWW是剪切應(yīng)力最大的方向，也是本區(qū)主要構(gòu)造線展布方向。而地裂縫的走向也多為NNE、NE向；

（2）活動(dòng)斷裂作為構(gòu)造單元的軟弱帶，是新構(gòu)造運(yùn)動(dòng)最活躍的地帶，幾乎控制了區(qū)域構(gòu)造構(gòu)造應(yīng)力場(chǎng)的應(yīng)變。其中，在斷裂帶的拐點(diǎn)、端點(diǎn)和交匯點(diǎn)易發(fā)生應(yīng)變集中，而這些應(yīng)力集中區(qū)域正是實(shí)際當(dāng)中地裂縫的高發(fā)區(qū)域；

（3）在區(qū)域拉張應(yīng)力場(chǎng)與重力場(chǎng)的聯(lián)合作用下，活動(dòng)斷裂帶沿線區(qū)段出現(xiàn)拉張與擠壓緊密伴生的應(yīng)力集中，并在斷裂帶之外的其他區(qū)域，還會(huì)在多個(gè)局部地段發(fā)生次級(jí)的拉應(yīng)力集中現(xiàn)象。因此，在區(qū)域構(gòu)造應(yīng)力場(chǎng)作用下，不僅在斷裂帶分布地段容易發(fā)育地裂縫，在斷裂帶分布地段之外的次級(jí)拉應(yīng)力集中區(qū)段內(nèi)，也會(huì)導(dǎo)致地裂縫的出現(xiàn)；

（4）平面模型可較好地反映出斷裂帶上的應(yīng)力應(yīng)變集中現(xiàn)象，針對(duì)斷裂帶附近區(qū)域，地表地裂縫的發(fā)育特征可有效預(yù)測(cè)評(píng)價(jià)。

References)

[1] 王景明,王江帥,劉金峰,等.地裂縫及其災(zāi)害的理論與應(yīng)用[M].西安:陜西科學(xué)技術(shù)出版社,2000.

Wang J M, Wang J S, Liu J F, et al.Theory of ground fissures hazards and its disasters[M].Shaanxi Science and Technology Press, 2000.

[2] 白永輝,張麗.河北省滄州市地質(zhì)災(zāi)害與地下水關(guān)系研究[J].中國(guó)地質(zhì)災(zāi)害與防治學(xué)報(bào),2005,16(3).

Bai Y H, Zhang L.Study on relationship between geological hazard and groundwater in Cangzhou, Hebei Province[J].The Chinese Journal of Geological Hazard and Control, 2005,16(3).

[3] 韓占成,韓彥霞.滄州市地下水環(huán)境地質(zhì)問(wèn)題與防治措施[J].地下水,2006,28(3).

Han Z C, Han Y X.Groundwater environmental geological problems and prevention measures in Cangzhou[J].Underground Water, 2006,28(3).

[4] 馬學(xué)軍,呂鳳蘭,宋偉,等.華北平原地裂縫影響因素分析[J].上海國(guó)土資源,2016,37(1):74-76.

Ma X J, Lü F L, Song W, et al.Major infuence factors of ground fssure in the North China Plain[J].Shanghai Land & Resources, 2016,37(1):74-76.

[5] 呂鳳蘭,馬學(xué)軍,邵長(zhǎng)慶,等.河北平原地裂縫分布特征及成因分析[J].上海國(guó)土資源,2014,35(4):49-52,57.

Lü F L, Ma X J, Shao C Q, et al.The distribution and causes of ground fissures on the Hebei plain[J].Shanghai Land & Resources, 2014,35(4):49-52,57.

[6] 李昌存.河北平原地裂縫研究[D].中國(guó)地質(zhì)大學(xué)博士論文,2003.

Li C C.Study on ground fissures in Hebei plain[D].Doctoral dissertation, China University of Geosciences.2003.

[7] 林舸, 趙重斌,張晏華,等.地質(zhì)構(gòu)造變形數(shù)值模擬研究的原理、方法及相關(guān)進(jìn)展[J].地球科學(xué)進(jìn)展,2005,20(5).

Lin G, Zhao C B, Zhang Y H, et al.The principle, method and related research progress on the numerical modeling of geological structural deformation[J].Advances in Earth Science, 2005,20(5).

[8] 鄭貴州,申永利.地質(zhì)特征三維分析及三維模擬現(xiàn)狀研究[J].地球科學(xué)進(jìn)展,2004,19(2).

Zheng G Z, Shen Y L.3D analysis of geological characteristics and status research of 3D geology modelling[J].Advances in Earth Sciences, 2004,19(2).

[9] 陳志新,伍素蘭.大同市地裂縫成因模型及數(shù)值模擬//2l世紀(jì)環(huán)境地學(xué)論文[C].石油出版社,1996.

Chen Z X, Wu S L.Genetic model and numerical simulation of ground fissures in Datong//2l century environmental geology papers[C].Petroleum Press, 1996.

[10] 彭建兵,范文,李喜安,等.汾渭盆地地裂縫成因研究中的若干關(guān)鍵問(wèn)題[J].工程地質(zhì)學(xué)報(bào),2007,15(4).

Peng J B, Fan W, Li X A, et al.Some key questions in the formation of ground fssures in the Fen-Wei Basin[J].Journal of Engineering Geology, 2007,15(4).

[11] 國(guó)家地震局地質(zhì)研究所.河北平原及其周圍地區(qū)新構(gòu)造圖及說(shuō)明書(shū)[M].地震出版社,1985.

Geological Institute of State Seismological Bureau.New structures and specifcations of the Hebei plain and its surrounding areas[M].Seismological Press, 1985.

[12] 馬潤(rùn)勇.華北平原典型地裂縫發(fā)育地段活動(dòng)斷裂調(diào)查與構(gòu)造應(yīng)力場(chǎng)數(shù)值模擬研究[R].長(zhǎng)安大學(xué),2009.

Ma R Y.Investigation of active faults and numerical simulation of tectonic stress feld in typical fracture development areas of NorthChina Plain[R].Chang'an University, 2009.

[13] 王凱.鮮水河斷裂帶關(guān)鍵地段構(gòu)造應(yīng)力場(chǎng)數(shù)值模擬研究[D].長(zhǎng)安大學(xué)碩士畢業(yè)論文,2011.

Wang K.The numerical simulation research on tectonic stress feld in key areas of Xianshuihe fault belt[D].Chang'an University, 2011.

[14] 宋衛(wèi)華,張宏偉,徐秀茹.區(qū)域構(gòu)造應(yīng)力場(chǎng)的數(shù)值模擬與應(yīng)用[J].遼寧工程技術(shù)大學(xué)學(xué)報(bào),2006,25(1).

Song W H, Zhang H W, Xu X R.Numerical simulation and application of regional tectonic stress feld[J].Journal of Liaoning Technical University, 2006,25(1).

[15] 陳望和,倪明云.河北第四紀(jì)地質(zhì)[M].地質(zhì)出版社,1987.

Chen W H, Ni M Y.Quaternary Geology of Hebei[M].Geology Press, 1987.

[16] 王瑜.中生代以來(lái)華北地區(qū)造山帶與盆地的演化及動(dòng)力學(xué)[M].地震出版社,1998.

Wang Y.Evolutional processes and dynamics of the orogenic belts and basins in North China since the Mesozoic[M].Seismological Press, 1988.

[17] 鄧起東,程紹平,閔偉,等.鄂爾多斯塊體新生代構(gòu)造活動(dòng)和動(dòng)力學(xué)討論[J].地質(zhì)力學(xué)學(xué)報(bào),1999,5(3).

Deng Q D, Cheng S P, Min W, et al.Discussion on cenozoic tectonics and dynamics of Ordos block[J].Journal of Geomechanics, 1999,5(3).

[18] 馬廷著,黃佩玉.唐山7.0級(jí)地震前后的斷裂運(yùn)動(dòng)[J].地震地質(zhì), 1982,4(4).

Ma T Z, Huang P Y.Fracture movements before and after the Tangshan 7 earthquake[J].Earthquake Geology, 1982,4(4).

[19] 鄧起東,張培震,冉勇康,等.中國(guó)活動(dòng)構(gòu)造基本特征[J].中國(guó)科學(xué)(D輯),2002,32(14).

Deng Q D, Zhang P Z, Ran Y K, et al.Basic characteristics of active tectonics in China[J].Chinese Science (D), 2002,32(14).

[20] 丁國(guó)瑜.中國(guó)巖石圈動(dòng)力學(xué)概論[M].地質(zhì)出版社,1991.

Ding G Y.An introduction to lithospheric dynamics in China[M].Geology Press, 1991.

[21] 徐繼山.華北陸緣盆地地裂縫成因機(jī)理研究[D].長(zhǎng)安大學(xué)博士論文,2012.

Xu J S.Study on genetic mechanism of ground fssures in North China continental margin basin[D].Chang'an University.2012.

The two-dimensional plane model study of tectonic ground fssure stress feld in Hejian County of Cangzhou city, Hebei Province

SONG Wei, MA Xue-Jun, Lü Feng-Lan, WANG Bing-Hu, SHAO Chang-Qing
(Center for Hydrogeology and Environmental Geology, China Geological Survey, Hebei Baoding 071051, China)

By solving the plane should force the fnite element inversion method the Cangzhou Hejian regional tectonic stress feld, and determine the maximum live should force feld and the maximum shear stress concentration zone, so as to explore the analysis of potential danger area of tectonic ground fssures; establish the two-dimensional plane model of ground fssure and calculation analysis in the area of Hejian regional tectonic stress feld and regional stress concentration zone and fssure in the relationship.

ground fssure; tectonic stress feld; fnite element analysis; two dimensional plane model

P642.27

A

2095-1329(2017)03-0083-07

10.3969/j.issn.2095-1329.2017.03.019

2017-05-31

修回日期: 2017-07-09

宋偉(1985-),男,學(xué)士,工程師,主要從事地質(zhì)災(zāi)害研究.

電子郵箱: sw262914595@126.com

聯(lián)系電話: 0312-8908251

中國(guó)地質(zhì)調(diào)查局地質(zhì)調(diào)查項(xiàng)目“京津冀地面沉降地裂縫調(diào)查及地質(zhì)環(huán)境監(jiān)測(cè)”(DD20160236)