川東南東溪地區龍馬溪組裂縫分布預測及頁巖氣保存條件評價

謝佳彤，付小平，秦啟榮，李 虎

川東南東溪地區龍馬溪組裂縫分布預測及頁巖氣保存條件評價

謝佳彤1，付小平1，秦啟榮2，李 虎3

(1. 中國石化勘探分公司，四川 成都 610041；2. 西南石油大學 地球科學與技術學院，四川 成都 610500；3. 四川建筑職業技術學院，四川 成都 610399)

頁巖氣的商業性開發證實成熟頁巖具有較大的勘探潛力，川東南地區是我國南方頁巖氣勘探開發的戰略先導區。裂縫發育情況是頁巖氣保存和開發部署的關鍵因素，對于構造環境復雜特殊的盆緣地區，在加大頁巖裂縫研究難度的同時，也給頁巖氣勘探開發提供了新的方向。頁巖氣勘探作為我國油氣勘探的新領域，尚未形成成熟的裂縫評價體系。以川東南東溪地區為例，基于三軸巖石力學實驗結果，首先采用有限元數值模擬手段進行應力場模擬，而后利用巖石力學參數進行模型構建和模擬結果的反復試算，以獲取東溪地區最大、最小主應力和應力差分布圖，最后運用莫爾–庫倫準則計算巖體破裂系數預測裂縫分布規律。結果表明：研究區裂縫發育情況主要分為4個級別，其中，Ⅰ級裂縫呈條帶狀分布在東西兩側斷裂附近，Ⅱ級裂縫在Ⅰ級裂縫周圍發育，而Ⅲ級裂縫分布在北部背斜核部和西部背斜翼部；優選頁巖氣埋深、距剝蝕區距離、距齊岳山斷裂距離、斷裂作用、應力差及壓力系數等10個保存條件參數，采用組合權重法確定保存條件參數權重，其中一級參數權重分別為埋深(0.2)、距剝蝕區距離(0.1)、距齊岳山斷裂距離(0.1)、斷裂作用(0.25)、應力差(0.15)、壓力系數(0.2)，以此建立相對完善的頁巖氣評價體系與評價標準。優選出2類頁巖氣勘探目標有利區，其中A類、B類有利區分別分布在研究區南部背斜核部及翼部和中部寬緩褶皺兩翼。研究成果為川東南地區的頁巖氣勘探提供重要的參考價值。

川東南；東溪地區；頁巖儲層；構造應力場模擬；裂縫分布預測；頁巖氣保存

對于具有自生自儲性質的頁巖而言，裂縫發育在烴類氣體儲集及橫縱向運移中均起重要作用。當今頁巖氣勘探的側重點之一主要在突破裂縫分布精準預測方法及其與保存條件之間的評價。目前常用方法有曲率法、模擬法、分形法等[1-6]，但是這些裂縫預測方法都具有一定的局限性，在裂縫預測中比較單一，裂縫預測準確性不高。構造應力場分布與裂縫發育密切相關，是控制裂縫發育的最主要因素之一，構造變形和構造活動性是頁巖氣保存的重要影響因素[7]。裂縫預測的地質研究中主要根據巖性組合以及巖相劃分等參數直接對裂縫密度進行統計進而直接進行定性推導，更多者傾向于對單一的應力場模擬進行計算[8-9]。四川盆地龍馬溪組頁巖氣已取得商業性開發。由于川東南含有盆內、盆緣和盆外三區，其受到的構造運動不同，保存條件也不同，為了準確評價，蓋層條件、壓力系數、構造條件、物質基礎、地層水條件和氣體組分等6個主要參數被挑選用于評價[10]，其中構造運動規模、頂底板密封性和巖層壓力系數3個參數已被用于四川盆地頁巖氣保存條件評價，取得了較好的效果[11]；構造部位與地層傾角已作為頁巖氣評價體系的重要指標[12-13]。但受控于頁巖儲層的各向異性及不同地區區域構造、地質條件差異，針對特定儲層評價時，其頁巖氣保存評價參數以及權重系數的確定不能直接復制他區現有參數[10-16]。

四川盆地東溪地區第1口高產井DYS1井是國內首口埋深大于4 200 m頁巖氣井，具有較強的代表性。為了深入剖析其高產經驗，筆者詳細剖析其裂縫分布預測方法和頁巖氣保存條件，利用構造應力場模擬方法預測應力場，并在巖體破裂準則的基礎上預測裂縫展布規律；進而劃分頁巖氣保存條件的指標，運用組合權重法確定評價指標的權重系數，以預測和優選有利區帶。研究方法為下一步區域頁巖氣勘探開發提供參考和借鑒。

1 區域地質特征

東溪地區位于四川盆地東南部，隸屬于重慶市綦江縣石壕鎮，研究區主要為東溪構造，分布在川東南盆內及盆緣拗褶區，面積約240 km2。四川盆地主體受到大斷裂擠壓形成一個大構造單元，四面均受到逆沖斷裂和走滑斷裂的擠壓，北部以米倉山–大巴山斷褶帶為界，西部為龍門山大斷裂，東部為齊岳山斷裂帶，南部主要為黔中隆起的北部邊緣帶(圖1)。志留系龍馬溪組為研究區目的層位，從上到下可分為2種巖性段，整體厚度為100～300 m，龍馬溪組上段厚度為20～80 m，巖性為深灰色泥巖夾粉砂質、灰質泥頁巖，其TOC含量相對較低；龍馬溪組下段厚度為100～180 m，巖性為富含筆石化石的黑色頁巖，局部夾微粒黃鐵礦條帶，下段由上至下TOC含量逐漸增大(3%～6%)，是構成頁巖氣的主力產氣層段。截至2020年，東溪地區頁巖氣測試產量已突破100萬m3/d，潛力非常大，即將成為第二個焦石壩頁巖氣田[17-18]。

2 構造應力場模擬

常見的構造應力場模擬主要運用有限元法，大型的有限元軟件將基礎的計算模式鑲嵌在算法中，可以快速得到模擬結果且更加準確，有限元模擬具體步驟如圖2。

2.1 地質模型

地質模型是地質體的一種反射，地質模型的構建應充分考慮研究區實際地質特征，特別是由構造作用產生的穿層性斷層，同時對部分斷層做相應的合并處理。考慮到模型的復雜性、準確性及資料掌握完整性，本次建模采用二維方法，以降低多期構造運動帶來的復雜構造格局的影響。

圖1 四川盆地區域構造及巖性柱狀圖

圖2 構造應力場數值模擬流程

2.2 力學模型

根據斷層、斜坡、背斜、向斜及正常沉積區域等對力學模型進行構造單元劃分，對不同構造單元進行參數賦值，主要參數包括巖石彈性模量()、泊松比()、黏聚力()、內摩擦角()和密度()等。本次研究采用彈性靜力學模型，主要分析研究區在邊界區域應力作用下，斷裂發生、斷層應力迅速釋放、構造形成前后的巖體應力–應變關系。

2.2.1 邊界條件確定

東溪地區力學模型邊界設置、構造力作用方位、約束條件如圖3所示。根據三軸巖石力學實驗，在目的層選取DYS1井的不同樣品進行實驗，確定模型邊界條件：右上施加最小主應力3=80 MPa，右下施加最大主應力1=120 MPa。位移邊界條件：NW、SW向兩邊固定，設置位移和轉角為零，其他單元均可自由移動[19]。

圖3 東溪地區構造應力場模擬邊界條件

2.2.2 巖石力學參數

在建立的研究區模型中，對所劃分的斷層、斜坡、背斜、向斜及正常沉積區域等不同的構造區進行材料參數的差異賦值。對斷層而言，根據斷距大小及規模進行三級劃分處理，其中，齊岳山斷裂作為基底斷層劃分為一級斷裂，二級及三級斷裂劃分的主要依據是延伸規模。根據巖石類型差異對五峰–龍馬溪組進行單元劃分，主要巖石類型包括黑色含粉砂質泥巖、黑色炭質泥巖和灰白色灰巖。

假設斷層與斷層所圍限的塊體為連續介質，斷層處材料屬性設置不同。DYS1井樣品巖石力學參數見表1。根據前人研究經驗[20]并依據模型的反復試算結果，將斷層的材料參數處理如下：構造形成前將斷層參數強化，其彈性模量比正常區域大，泊松比比正常區域小，密度不變。

2.2.3 斷層處理方案

在有限元數值模擬中，國內對斷層處理主要方法分為弱化法和接觸面法[21]，通過對比分析，本次采用斷層弱化法。對斷層、構造高點、斜坡等不同部位巖石彈性參數進行差異化賦值處理，斷層區域內巖石彈性模量小于正常地層，泊松比大于正常地層，構造高區巖石彈性模量大于正常地層，泊松比則相對較小。黏聚力及內摩擦角在實驗數據的基礎上，結合模擬情況，并在地質認識的基礎上進行相應的調整，構造上易破碎區黏聚力可調小，相反則調大，直至模擬結果與實際測試結果一致(表2)。

表1 DYS1井巖石力學參數

表2 不同單元類型巖石力學參數

2.2.4 反演標準

對模型進行網格節點劃分，模擬應力場大小及分布，不斷調試模型，使得模擬計算得到的應力場大小、應力方向與施加的約束條件相近或一致。模型計算過程中，調整最大主應力、最小主應力的大小及方向，通過不斷試算，將計算結果與實驗結果進行比較，直到試算結果與實驗結果一致。

2.3 數學模型

利用有限元原理建立相應的數學模型，根據劃分模型單元節點與未知點之間的平衡關系，搭建二者之間的數學函數方程式。在限定邊界條件之后，求解方程式，得到未知點結果，計算結果的可靠性與網格單元劃分的精細程度呈正相關。

由前述分析可知，研究區內斷裂極為發育，在模型構造分析中，僅就規模最大的15條斷裂進行模擬。該模型經離散后，東溪模型共有7 500個單元和15 352個節點，離散后的模型如圖4所示。

圖4 三維有限元模型

2.4 構造應力場分布

以研究區的主要地質構造特征為基礎，建立構造應力場有限元分析力學模型[22]。經過大量模擬試算，東溪地區構造形成后初步模擬結果如圖5所示。在持續受到NE向構造擠壓力的作用下，在DYS1井斷層附近及構造脊部附近應力釋放，其應力值較小，最大主應力值分布比較平均，介于–62.80～ 157.3 MPa。斷層附近出現應力集中現象，最大、最小主應力均分布于大斷層附近，DYS1井附近最大主應力值為–95.6 MPa，向DYS1井附近應力值呈逐漸減小的趨勢，向NE方向斷層減少埋深增加，應力值呈逐漸增加趨勢。

通過最大、最小主應力分布求取二者的應力差值，應力差分布由南向北呈逐漸降低趨勢，整體為60～102 MPa。大斷裂附近應力差值較大，破裂程度高，DYS1井應力差在40 MPa左右(圖5c)。

3 裂縫分布預測方法

3.1 巖體破裂系數

巖石破裂系數表示在應力應變中巖石受到應力大于巖石所承受的強度時產生破壞的邊界值，一般可用該值表示巖石變形程度，在構造地質領域中用破裂系數來表示巖體的裂縫發育情況，由于地層埋深設為彈性變形，為了使模擬預測高效便捷，主要采用莫爾–庫倫準則進行判別，部分地區采用格里菲斯準則。根據莫爾理論再結合格里菲斯理論，作巖石綜合破裂破壞接近程度圖解(圖6)，巖體破裂系數表示如下：

式中：f、k分別為剪應力和最大抗剪能力，MPa；σ1為最大主應力，MPa；σ3為最小主應力，MPa；c為黏聚力，MPa；φ為內摩擦角，(°)。

圖6 巖石綜合破裂準則破壞接近程度圖解

在巖石力學理論中巖石破裂系數一般具有以下判斷標準：當<1時巖體比較穩定，沒有裂縫產生；當≥1時巖體受到的應力使巖體失去穩定性發生破壞而產生裂縫，在實際情況中巖體的破裂往往比較復雜，理論上<1時巖體不破裂，但實際中即使<1的巖體，其內部也會因為受到應力產生細小的微裂縫。因此，在運用破裂理論進行裂縫破裂程度模擬時，不同研究區其標準不同，破裂系數一般可以定性表示裂縫的相對發育程度，值越大，裂縫就越發育。

3.2 巖體破裂特征

通過式(1)，對東溪地區進行分區破裂特征模擬預測，從平面上得到巖體破裂程度分布特征。研究區主要形成2個方向的裂縫，破裂程度分布范圍廣，巖體破裂系數最大值為1.921，最小值0.658。除邊界條件影響綜合分析，研究區一般巖體的主體破壞程度系數為0.895～1.763，大部分接近或大于破裂臨界值(圖7)。

圖7 東溪地區巖體破裂程度

從圖7可知，巖體破壞程度最高的地區主要分布在東溪地區東南部，具體位于近NW向斷層及NE向斷層周圍，呈條帶狀展布，斷層與斷層周圍破裂帶、斷層外圍以及構造高部位，破裂系數達1.447～1.605；距離斷層較遠且由西翼變形區、東翼變形區過渡到東部向斜地區破壞程度逐漸降低，最終破壞度普遍小于1.053；DYS1井由于埋深較大且遠離斷層，破裂程度較低，破裂系數為1左右。

3.3 裂縫分布定量預測

裂縫發育的定量預測，實際上是利用裂縫發育的力學特性，將地質學、物理學及數理統計學相結合，對裂縫發育分布的理論性探索[23]，目前裂縫分布的預測方法主要有：①傳統地質預測方法，通過尋求具體地區地表露頭與地下巖心裂縫發育之間的對應關系，達到裂縫預測的目的；②利用數學統計學原理進行分布預測；③利用現代地球物理勘探技術，通過地震波屬性處理實現裂縫預測；④數值模擬法。利用巖石破裂準則，通過古構造應力場的反演可實現裂縫預測的目的[24]。本次主要根據現有資料采用巖體破裂系數法進行裂縫初步預測。

在巖體破裂程度分布預測的基礎上，結合東溪地區的構造地質背景、構造地質特征以及現有的實驗測試結果及生產數據等資料，依據東溪地區的巖體破裂程度對裂縫發育的影響制定判別標準，除斷裂帶以外將東溪地區裂縫發育劃分4個級別(表3)。

表3 東溪地區龍馬溪組巖石裂縫預測的η值標準

根據表3，并以巖石力學及構造地質學理論為指導，綜合現場鉆錄井及實際生產資料分析，預測東溪地區龍馬溪組裂縫發育情況(圖8)，除斷裂帶以外，分為Ⅰ—Ⅳ級。

圖8 東溪地區裂縫發育預測結果

1) 斷裂帶

斷裂帶巖體破裂程度最高，發育在巖體破裂最明顯的斷裂及其附近周圍區域。研究區NE向斷裂為早期形成，近SN向以及NW向斷層晚期形成，晚期形成的破裂對早期有疊加改造作用，使得形成的復合斷裂破裂值更大，在復合斷裂附近，值普遍在1.684以上，研究區巖體破裂最強區域值高達2.000。

2) Ⅰ級裂縫發育區

主要分布在斷裂帶周圍破裂帶、斷層外圍以及局部構造高點附近，整體沿斷層兩側呈NW和NE向條帶狀分布，該區域內值均在1.368～1.684。Ⅰ級裂縫發育區分布受區域應力場直接控制。

3) Ⅱ級裂縫較發育區

褶皺脊部頂部、Ⅰ級發育區的外圍及相鄰次級斷層的周圍為Ⅱ級裂縫發育區，分布面積比較廣，占東溪地區的30%，圍繞Ⅰ級裂縫區呈NE向條帶狀分布，該區巖體較高，為1.053～1.368。

4) Ⅲ級裂縫發育臨界區

主要位于除斷裂帶、Ⅰ級和Ⅱ級裂縫發育區以外的構造低部位區，分布在東溪地區的北部及西北方向，分布面積小于Ⅱ級裂縫發育區，整體呈近南北及北西向延伸的帶狀，巖石破壞接近程度較小，為0.737～1.053。

5) Ⅳ級裂縫不發育區

分布面積較小，主要分布在東溪地區的北部，巖石破壞程度非常小，基本不破壞，<0.737。

4 頁巖氣保存條件評價

4.1 影響因素

1) 構造與埋深

東溪構造主體呈“兩凹夾一隆”式構造格局，可進一步細化為5個次級構造(圖9)。東溪斷背斜區(埋深4 000～4 500 m)，東西兩側分別與向斜接觸，南部與桃子蕩斷洼接觸；西部向斜區(埋深4 500～5 500 m)，與東溪斷背斜呈斷層接觸；桃子蕩斷洼區(埋深3 500～4 500 m)，斷層夾持下降盤；東部抬升斷塊區(埋深1 500～3 500 m)，齊岳山及其派生斷裂切割；東部向斜區(埋深4 500～5 500 m)，西翼與斷背斜接觸，東翼部與東部抬升斷塊區呈斷洼接觸。東部向斜區可分為東翼變形區、向斜核部區、西翼變形區，埋深范圍分別為4 500～5 000 m、5 000～6 000 m、4 000～5 000 m。

2) 內部斷裂和裂縫

東溪地區發育多級多組斷層：NNW向延伸較長，斷距較大，高陡褶皺帶控制斷裂，與現今最大水平主應力近垂直，封堵性較好；NE向斷層，主要為齊岳山斷裂帶斷層，切穿二疊系，與現今最大水平主應力近垂直，封堵性較好(圖9、圖10)。DX-F1斷層為三級斷裂，對保存條件影響較小。DX-F1斷層具有相似特征，為向斜內部斷裂，延伸長度短(32 km)，斷距有向兩端明顯減小的趨勢，與現今最大主應力近垂直，無走滑性質，對保存條件影響不大(表4)。

圖9 東溪地區構造分區

圖10 東溪地區斷裂剖面解釋結果

3) 壓力系數

對東溪龍馬溪組優質泥頁巖二維地層進行壓力系數預測，結果發現，東溪地區壓力系數變化范圍不大，整體較高，均在1以上，遠離剝蝕區的壓力系數較大，靠近大斷裂附近壓力系數變小(圖11)，DYS1井區壓力系數約1.2～1.4，東溪地區整體表現為異常高壓，利于頁巖氣保存；靠近剝蝕區，受到齊岳山大斷裂的影響，壓力系數逐漸降低(<1)，頁巖氣易逸散，不利于保存。

表4 東溪地區斷裂系統及斷裂特征統計

注釋：TS表示時間域志留系；TP表示時間域二疊系；T?表示時間域寒武系。

圖11 東溪地區龍馬溪組壓力系數預測平面圖

4.2 評價方法

運用主客觀組合法確定評價參數的權重，其中，專家經驗法為主觀權重法，灰色關聯度為客觀權重法，二者評價過程中各有優缺點，其中，主觀權重法能充分體現樣本屬性，但是客觀性較差，而客觀權重法不能體現專家對不同屬性權重賦值差異，甚至有時會出現賦值權重與實際相悖的情況[25]，因此，為了兼顧2種評價方法的優點同時又盡量減少權重的主觀隨意性，使評價參數權重值的主、客觀評價結果近統一，使評價結果真實可靠。因此，本次評價綜合考慮評價指標數據的真實客觀性與主觀經驗的權重賦值。

2種常用的主客觀組合權重計算方法：“乘法”集成法和“加法”集成法。不同方法適用情況存在差異，其中“乘法”集成法適用于評價參數多、多參數權重分配均勻的情況，“加法”集成法應用過程中依賴使用者的主觀分析。其公式如下：

式中：w為第個參數因子的組合權重，共計有個；a、b分別為第個參數因子的客觀權重和主觀權重；為相關系數。

為了降低主觀賦權與客觀賦權的誤差，將適用“加法”集成法對專家經驗法和灰色關聯度法進行組合賦權[26]，運用組合賦權法中的“加法”集成法對研究區龍馬溪組進行分析統計，為了使權重相關性較大，使其誤差降低到最小，取相關系數為0.6，最終計算得出油氣地質條件子項的各參數值的權重(表5)。

4.3 有利區帶預測

在有利區預測過程中，參考前人關于參數及權重分析方法，重點考慮構造要素，結合東溪地區實際地質特征及勘探現狀[27-29]，優選埋深、距離剝蝕區距離、距離大斷裂距離、斷裂作用、應力差及壓力系數6個主要參數及次級斷裂規模、級別與期次、傾角大小、距大規模次級斷裂距離及發育程度等子參數作為關鍵參數進行頁巖氣保存條件評價。將地質方法和數學理論相結合，基于勘探實踐經驗，總結出適合東溪地區的海相頁巖氣綜合評價指標體系，并對各參數賦予權重系數(表5)。

表5 東溪地區龍馬溪組頁巖氣評價參數

依據權重系數及評價參數對東溪地區進行參數疊加，并進行區帶排序，同時采用均值聚類分析法對區帶進行分類，劃分為2類有利區(圖12)。評價結果表明，頁巖氣綜合評價指標為0.652 5～1。其中，當綜合指標值介于0.9～1.0時為A類有利區，位于東溪背斜區，具有埋深大(4 500～5 000 m)、遠離剝蝕區、斷裂規模小、發育時期早的特點，整體處于Ⅱ級裂縫發育區，次級斷裂傾角小于45°，應力差較低，預測壓力系數大于1.4；當綜合評價值介于0.8～0.9時為B類有利區，是較有利勘探目標區，B類有利區位于東翼變形區和西翼變形區，頁巖埋藏深度介于4 000～4 500 m，整體處于Ⅲ級裂縫發育區，地層應力差小，斷裂規模小，預測地層壓力系數在1.2以上。

5 結論

a. 四川盆地東溪地區應力場模擬結果表明，不同地區的破裂程度差異較大，由于研究區埋深較大，且構造比較穩定，最大最小主應力差值相對較小，差值較大區域一般分布在大斷裂附近或埋深較小地區。

b. 研究區裂縫發育程度劃分為Ⅰ—Ⅳ級，其中Ⅱ級裂縫發育區和Ⅲ級裂縫發育區對頁巖氣的運移及保存效果最好，屬于優勢發育區。其中，Ⅱ級分布面積比較廣，巖體較高，為1.053～1.368；Ⅲ級裂縫發育區整體呈近SN和NW向延伸的帶狀，巖石破壞接近程度較小，為0.737～1.053。

圖12 東溪地區龍馬溪組頁巖氣有利區評價

c. 運用組合權重法對東溪地區頁巖氣保存評價參數體系進行權重的劃分和計算，將研究區保存評價體系進行半定量分析和各參數權重賦值，并優選埋深、距離剝蝕區距離、距離齊岳大斷裂距離、應力差及壓力系數等10項作為關鍵參數進行頁巖氣保存條件評價。

d. 將研究區頁巖保存條件劃分為2類有利區。A類：主要位于深埋平緩區或寬緩褶皺核部，即東溪南部背斜核部及翼部，遠離主干斷裂，頂底板致密且連續，Ⅱ級裂縫發育區，早期以及三級或四級裂縫發育，壓力系數大，物質基礎好；B類：主要分布于寬緩褶皺兩翼，發育三級或四級裂縫，Ⅲ級裂縫發育區，物質基礎較好，壓力系數較大。

[1] 鄒才能，楊智，張國生，等. 常規–非常規油氣“有序聚集”理論認識及實踐意義[J]. 石油勘探與開發，2014，41(1)：14–27.

ZOU Caineng，YANG Zhi，ZHANG Guosheng，et al. Conventional and unconventional petroleum“orderly accumulation”：Concept and practical significance[J]. Petroleum Exploration and Development，2014，41(1)：14–27.

[2] CURTIS J B. Fractured shale–gas systems[J]. AAPG Bulletin，2002，86(11)：1921–1938.

[3] FENG Weiping，WANG Feiyu，GUAN Jing，et al. Geologic structure controls on initial productions of Lower Silurian Longmaxi shale in south China[J]. Marine and Petroleum Geology，2018，91：163–178.

[4] 周文. 裂縫性油氣儲集層評價方法[M]. 成都：四川科學技術出版社，1998.

ZHOU Wen. Evaluation methods of fracture reservoir in oil and gas pool[M]. Chengdu：Sichuan Science and Technology Press，1998.

[5] 徐旭輝，申寶劍，李志明，等. 頁巖氣實驗地質評價技術研究現狀及展望[J]. 油氣藏評價與開發，2020，10(1)：1–8.

XU Xuhui，SHEN Baojian，LI Zhiming，et al. Status and prospect of experimental technologies of geological evaluation for shale gas[J]. Reservoir Evaluation and Development，2020，10(1)：1–8.

[6] 魏志紅. 四川盆地及其周緣五峰組–龍馬溪組頁巖氣的晚期逸散[J]. 石油與天然氣地質，2015，36(4)：659–665.

WEI Zhihong. Late fugitive emission of shale gas from Wufeng–Longmaxi Formation in Sichuan Basin and its periphery[J]. Oil and Gas Geology，2015，36(4)：659–665.

[7] 蔡周榮，夏斌，萬志峰. 下揚子蕪湖地區后期構造活動特征及其對古生界頁巖氣保存的影響[J]. 煤炭學報，2013，38(5): 890–895.

CAI Zhourong，XIA Bin，WAN Zhifeng. Characteristics of late tectonic activity in Wuhu and its influence on shale gas conservation in Paleozoic boundary[J]. Journal of China Coal Society，2013，38(5)：890–895.

[8] 任浩林，劉成林，劉文平，等. 四川盆地富順–永川地區五峰組—龍馬溪組應力場模擬及裂縫發育區預測[J]. 地質力學學報，2020，26(1): 74–83.

REN Haolin，LIU Chenglin，LIU Wenping，et al. Stress field simulation and fracture development prediction of the Wufeng Formation–Longmaxi Formation in the Fushun-Yongchuan Block，Sichuan Basin[J]. Journal of Geomechanics，2020，26(1)：74–83.

[9] 張守仁，萬天豐，陳建平. 川西坳陷孝泉–新場地區須家河組二—四段構造應力場模擬及裂縫發育區帶預測[J]. 石油與天然氣地質，2004(1)：70–74.

ZHANG Shouren，WAN Tianfeng，CHEN Jianping. Tectonic stress field modeling and fracture prediction in T3x2-4 strata in Xiaoquan-Xinchang area, western Sichuan depression[J]. Oil & Gas Geology，2004(1)：70–74.

[10] 聶海寬，何治亮，劉光祥，等. 中國頁巖氣勘探開發現狀與優選方向[J]. 中國礦業大學學報，2020，49(1)：13–35.

NIE Haikuan，HE Zhiliang，LIU Guangxiang，et al. Status and direction of shale gas exploration and development in China[J]. Journal of China University of Mining & Technology，2020，49(1)：13–35.

[11] 胡東風，張漢榮，倪楷，等. 四川盆地東南緣海相頁巖氣保存條件及其主控因素[J]. 天然氣工業，2014，34(6)：17–23.

HU Dongfeng，ZHANG Hanrong，NI Kai，et al. Main controlling factors for gas preservation conditions of marine shales in southeastern margins of the Sichuan Basin[J]. Natural Gas Industry，2014，34(6)：17–23.

[12] 王濡岳，胡宗全，劉敬壽，等. 中國南方海相與陸相頁巖裂縫發育特征及主控因素對比：以黔北岑鞏地區下寒武統為例[J]. 石油與天然氣地質，2018，39(4)：631–640.

WANG Ruyue，HU Zongquan，LIU Jingshou，et al. Comparative analysis of characteristics and controlling factors of fractures in marine and continental shales：A case study of the Lower Cambrian in Cengong area，northern Guizhou Province[J]. Oil & Gas Geology，2018，39(4)：631–640.

[13] 魏祥峰，李宇平，魏志紅，等. 保存條件對四川盆地及周緣海相頁巖氣富集高產的影響機制[J]. 石油實驗地質，2017，39(2)：147–153.

WEI Xiangfeng，LI Yuping，WEI Zhihong，et al. Effects of preservation conditions on enrichment and high yield of shale gas in Sichuan Basin and its periphery[J]. Petroleum Geology and Experiment，2017，39(2)：147–153.

[14] SURPLESS B，WIGGINTON S S. The impact of inter–bed cohesion on fold–related fracture development，still well anticline，west Texas(USA)[J]. Journal of Structural Geology，2020，134：103974.

[15] GAO Jian，ZHANG Jiankun，HE Sheng，et al. Overpressure generation and evolution in Lower Paleozoic gas shales of the Jiaoshiba region，China：Implications for shale gas accumulation[J]. Marine and Petroleum Geology，2019，102：844–859.

[16] GUO Peng，REN Desheng，XUE Yanhui. Simulation of multi–period tectonic stress fields and distribution prediction of tectonic fractures in tight gas reservoirs：A case study of the Tianhuan depression in western Ordos Basin，China[J]. Marine and Petroleum Geology，2019，109：530–546.

[17] JU Wei，WANG Jilin，FANG Huihuang，et al. Paleotectonic stress field modeling and prediction of natural fractures in the Lower Silurian Longmaxi shale reservoirs，Nanchuan region，South China[J]. Marine and Petroleum Geology，2019，100：20–30.

[18] 耿一凱. 川東南地區龍馬溪組頁巖儲層質量控制因素研究[D]. 北京：中國石油大學(北京)，2017.

GENG Yikai. Controlling factors on quality of shale reservoir of Longmaxi Formation in the southeastern Sichuan Basin[D]. Beijing：China University of Petroleum(Beijing)，2017.

[19] 尤廣芬. 儲層裂縫預測的三維有限元數值模擬研究[D]. 成都：成都理工大學，2010.

YOU Guangfen. Reservoir fracture prediction of three–dimensional finite element numerical simulation[D]. Chengdu：Chengdu University of Technology，2010.

[20] 丁文龍，尹帥，王興華，等. 致密砂巖氣儲層裂縫評價方法與表征[J]. 地學前緣，2015，22(4)：173–187.

DING Wenlong，YIN Shuai，WANG Xinghua，et al. Assessment method and characterization of tight sandstone gas reservoir fractures[J]. Earth Science Frontiers，2015，22(4)：173–187.

[21] 曾海容，宋惠珍. 斷裂端部特殊單元模型[J]. 地震地質，1999，21(3)：215–220.

ZENG Hairong，SONG Huizhen. A special model for the tip of faults[J]. Seismology and Geology，1999，21(3)：215–220.

[22] 曾維特. 渝東南地區下古生界頁巖儲層裂縫發育特征與分布預測研究[D]. 北京：中國地質大學(北京)，2014.

ZENG Weite. The characteristics and distribution prediction of fractures in Lower Palaeozoic shale，southeast of Chongqing area[D]. Beijing：China University of Geosciences(Beijing)，2014.

[23] LIU Jingshou，DING Wenlong，WANG Ruyue，et al. Simulation of paleo tectonic stress fields and quantitative prediction of multi–period fractures in shale reservoirs：A case study of the Niutitang Formation in the Lower Cambrian in the Cen’gong block，South China[J]. Marine and Petroleum Geology，2017，84：289–310.

[24] GUO Xusheng. Major factors controlling the shale gas accumulations in Wufeng–Longmaxi Formation of the Pingqiao shale gas field in Fuling area，Sichuan Basin，China[J]. Journal of Natural Gas Geoscience，2019，4(3)：129–138.

[25] 謝佳彤，秦啟榮，李斌，等. 塔里木盆地中央隆起帶良里塔格組預探階段有利區帶評價[J]. 特種油氣藏，2018，25(6)：25–31.

XIE Jiatong，QIN Qirong，LI Bin，et al. Favorable zone evaluation of Lianglitage Formation within pre–exploration stage in the central uplift of Tarim Basin[J]. Special Oil & Gas Reservoirs，2018，25(6)：25–31.

[26] 漆艷茹. 確定指標權重的方法及應用研究[D]. 沈陽：東北大學，2010.

QI Yanru. The method of determining index weight and its application research[D]. Shenyang：Northeastern University，2010.

[27] 翟常博. 川東南綦江–仁懷地區頁巖氣成蔵條件及有利目標區研究[D]. 北京：中國地質大學(北京)，2013.

ZHAI Changbo. Reservoir forming condition and exploration prospect of shale–gas in Qijiang-Renhuai area of southeast Sichuan[D]. Beijing：China University of Geosciences(Beijing)，2013.

[28] 殷啟春，方朝剛，鄭紅軍，等. 下揚子地區奧陶紀頁巖氣地質條件及遠景區優選[J]. 華東地質，2020，41(1)：70–78.

YIN Qichun，FANG Chaogang，ZHENG Hongjun，et al. Geological conditions of Ordovician shale gas and optimization for prospective areas in the Lower Yangtze region[J]. East China Geology，2020，41(1)：70–78.

[29] XIE Jiatong，QIN Qirong，FAN Cunhui. Quantitative prediction of fracture distribution of the Longmaxi Formation in the Dingshan area，China using FEM numerical simulation[J]. Acta Geologica Sinica(English Edition)，2019，93(6)：1662–1672.

Prediction of fracture distribution and evaluation of shale gas preservation conditions in Longmaxi Formation in Dongxi area

XIE Jiatong1, FU Xiaoping1, QIN Qirong2, LI Hu3

(1. Sinopec Exploration Company, Chengdu 610041, China; 2. School of Geoscience and Technology, Southwestern Petroleum University, Chengdu 610500, China; 3. Sichuan College of Architectural Technology, Chengdu 610399, China)

Commercial development of shale gas confirmed that mature shale has great potential for exploration, southeast Sichuan is the strategic pilot area of shale gas exploration and development in southern China, cracks are the key factor of shale gas preservation and development deployment, for the structural environment complex and special basin edge area, increasing the difficulty of shale crack research, but also provides a new direction for shale gas exploration and development. As a new field of oil and gas exploration in China, shale gas exploration has not yet developed a mature crack evaluation system. Taking Dongxi area of the southeast of Sichuan as an example, the stress field was simulated in the region by finite element numerical simulation. The model was constructed and calculated by rock mechanical parameters, the diagrams of maximum principal stress, minimum principal stress and differential stress distribution in Dongxi area were obtained through repeated debugging of the simulation results. Mohr-Coulomb criterion was used to calculate rock fracture coefficient and to predict the law of crack distribution in the study area, The prediction shows that the crack distribution of the study area is divided into four levels, GradeⅠcracks are distributed in bands near the east and west sides, grade cracks develop around gradeⅠcracks, the grade cracks are distributed in the northern dorsal nucleus and the western dorsal wing, respectively. At the same time, 10 preservation condition parameters, such as shale gas burial depth, distance from the ablative area, fracture distance from Qiyun Mountain, fracture action, differential stress and pressure coefficient, were clarified, and were used to determine the primary and secondary parameter weight by the combined weight method. Among them, the first stage parameters are buried depth(0.2), distance from the exfoliation area(0.1), fracture distance(0.1), fracture action(0.25), differential stress(0.15), and pressure coefficient(0.2). A relatively sound shale gas evaluation system and evaluation standards have been established. The two types of shale gas exploration targets are preferred, the favorable areas of the typeⅠand the typeⅡare located in the core and wings of the anticline core in the south of the study area and wings of the gentle folds in the middle of the study area, This study provides an important reference for shale gas exploration in southeast Sichuan.

southeast Sichuan; Dongxi area; shale reservoir; tectonic stress field simulation; fracture distribution prediction; shale gas conservation

語音講解

TE132.2

A

1001-1986(2021)06-0035-11

2021-06-03；

2021-10-13

國家科技重大專項任務(2017ZX05036-003-003)

謝佳彤，1992年生，女，黑龍江大慶人，博士，工程師，從事頁巖氣地質綜合研究工作. E-mail：963955769@qq.com

謝佳彤，付小平，秦啟榮，等. 川東南東溪地區龍馬溪組裂縫分布預測及頁巖氣保存條件評價[J]. 煤田地質與勘探，2021，49(6)：35–45. doi: 10.3969/j.issn.1001-1986.2021.06.004

XIE Jiatong，FU Xiaoping，QIN Qirong，et al.Prediction of fracture distribution and evaluation of shale gas preservation conditions in Longmaxi Formation in Dongxi area[J]. Coal Geology & Exploration，2021，49(6)：35–45. doi: 10.3969/ j.issn.1001-1986.2021.06.004

移動閱讀

(責任編輯 范章群)