劉 建 惠 晨 樊建明 呂文雅 王繼偉 尹 陳 王浩南
1. 中國(guó)石油長(zhǎng)慶油田分公司勘探開(kāi)發(fā)研究院, 陜西 西安 710018;2. 中國(guó)石油大學(xué) (北京) 地球科學(xué)學(xué)院, 北京 102249;3. 中國(guó)石油大學(xué) (北京) 油氣資源與探測(cè)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 北京 102249
空氣滲透率小于1 mD 的儲(chǔ)層被稱為致密儲(chǔ)層(鄒才能等, 2015), 致密油氣是繼頁(yè)巖氣之后全球非常規(guī)油氣勘探開(kāi)發(fā)的新熱點(diǎn) (賈承造等,2012)。 致密儲(chǔ)層物性差, 利用常規(guī)開(kāi)采技術(shù)開(kāi)發(fā)難度大, 針對(duì)該類儲(chǔ)層, 需采用大型壓裂改造進(jìn)行投產(chǎn), 而地應(yīng)力對(duì)壓裂增產(chǎn)效果具有重要影響(張斌等, 2012), 同時(shí), 地應(yīng)力還影響致密油藏井網(wǎng)部署、 注水管理和壓裂改造方案等方面 (曾聯(lián)波, 2008; 萬(wàn)曉龍等, 2009; 翁劍橋, 2020)。因此地應(yīng)力的研究在致密油氣開(kāi)發(fā)過(guò)程中扮演著十分重要的角色。
現(xiàn)今地應(yīng)力主要由重力、 構(gòu)造應(yīng)力、 孔隙流體壓力和熱應(yīng)力耦合而成 (廖新武等, 2015), 是多期構(gòu)造作用和新構(gòu)造活動(dòng)共同影響而成 (張浩等, 2020)。 文中的地應(yīng)力指現(xiàn)今地應(yīng)力狀態(tài), 通常泛指第四紀(jì)中更新世以來(lái)的地應(yīng)力 (曾聯(lián)波和田崇魯, 1998)。 地應(yīng)力的研究已經(jīng)有上百年的歷史, 國(guó)內(nèi)外學(xué)者在地應(yīng)力測(cè)量和計(jì)算方面開(kāi)展了大量的研究工作 (Fairhurst, 1964; Hast, 1969;Bell and Gough, 1979; 許忠淮, 1990), 提出了許多測(cè)量和計(jì)算地應(yīng)力的方法。 每種方法都有其優(yōu)缺點(diǎn), 應(yīng)根據(jù)實(shí)際情況選擇較合理的方法。 除此之外, 數(shù)值模擬也在地應(yīng)力研究中被廣泛應(yīng)用,應(yīng)力場(chǎng)數(shù)值模擬開(kāi)始于20 世紀(jì)中期, 隨著計(jì)算機(jī)技術(shù)的發(fā)展, 中國(guó)在20 世紀(jì)70 年代將該方法應(yīng)用在地學(xué)領(lǐng)域 (張勝利, 2011)。 數(shù)值模擬常用的方法為離散元與有限元方法, 常用的軟件有Petrel,Ansys, PFC, FLAC, ABAQUS 等(毛哲等, 2018)。其中, 有限元數(shù)值模擬法具有很高的適應(yīng)性, 可以模擬地質(zhì)演化過(guò)程、 階段等, 對(duì)受力復(fù)雜、 加載特殊的模擬問(wèn)題也可以實(shí)現(xiàn) (權(quán)凱, 2014)。 由于三維地質(zhì)建模復(fù)雜繁瑣 (劉愛(ài)華等, 2013), 目前大多學(xué)者進(jìn)行二維應(yīng)力場(chǎng)數(shù)值模擬, 但其只能展示平面的應(yīng)力分布, 三維應(yīng)力場(chǎng)數(shù)值模擬相比之下對(duì)空間的應(yīng)力變化就有更好地展示, 但需要考慮的因素和參數(shù)也更多, 其精度主要取決于地質(zhì)模型的建立和巖石力學(xué)參數(shù)的確定, 目前三維應(yīng)力場(chǎng)數(shù)值模擬多存在將同一地層賦予相同的巖石物理參數(shù), 對(duì)于層內(nèi)存在透鏡體、 砂泥巖互層等現(xiàn)象, 無(wú)法體現(xiàn)單層內(nèi)垂向上非均質(zhì)性及平面非均質(zhì)性對(duì)地應(yīng)力的影響等問(wèn)題 (劉愛(ài)華等,2013; 權(quán)凱, 2014; 劉洪濤等, 2016; 毛哲等,2018; 徐珂等, 2018)。
眾多學(xué)者從單井應(yīng)力測(cè)量、 二維應(yīng)力場(chǎng)數(shù)值模擬并結(jié)合生產(chǎn)動(dòng)態(tài)等方面對(duì)鄂爾多斯盆地的地應(yīng)力特征、 地應(yīng)力的應(yīng)用進(jìn)行了研究 (周新桂等,2009; 杜瑋暄等, 2010; 徐磊, 2016)。 研究表明,由于鄂爾多斯盆地的構(gòu)造變化較弱, 地應(yīng)力方向的變化較弱; 同時(shí), 水平最大主應(yīng)力的方向與差應(yīng)力的大小對(duì)壓裂縫的展布和復(fù)雜程度有很大的影響, 從而對(duì)壓裂生產(chǎn)有重要的影響。 總體來(lái)說(shuō),目前對(duì)鄂爾多斯盆地地應(yīng)力分布規(guī)律研究以二維地應(yīng)力數(shù)值模擬為主, 對(duì)地應(yīng)力分布的影響因素研究較少, 同時(shí)對(duì)合水地區(qū)長(zhǎng)6 致密儲(chǔ)層現(xiàn)今地應(yīng)力研究工作開(kāi)展較少, 缺少對(duì)該地層現(xiàn)今地應(yīng)力分布的認(rèn)識(shí), 一定程度上制約了該地區(qū)井網(wǎng)部署及壓裂施工等工作的進(jìn)程。
文章以鄂爾多斯盆地合水油田長(zhǎng)6 致密砂巖儲(chǔ)層為研究對(duì)象, 在對(duì)單井地應(yīng)力分析和區(qū)域地應(yīng)力特點(diǎn)綜合分析的基礎(chǔ)上, 結(jié)合研究區(qū)實(shí)際構(gòu)造特點(diǎn)、 沉積特點(diǎn)、 巖相特點(diǎn)建立三維非均質(zhì)地質(zhì)模型; 通過(guò)高溫高壓三軸巖石力學(xué)實(shí)驗(yàn)和研究區(qū)壓裂施工的實(shí)際數(shù)據(jù), 確定不同巖相的巖石物理參數(shù), 建立符合研究區(qū)實(shí)際的三維力學(xué)模型; 利用Ansys 有限元數(shù)值模擬軟件開(kāi)展三維現(xiàn)今應(yīng)力場(chǎng)數(shù)值模擬。 以研究區(qū)實(shí)測(cè)的地應(yīng)力為驗(yàn)證依據(jù),通過(guò)不斷調(diào)整模型和邊界條件, 進(jìn)行反復(fù)模擬,得到研究區(qū)的地應(yīng)力分布規(guī)律。 在此基礎(chǔ)上, 探討了研究區(qū)地應(yīng)力分布的影響因素及地應(yīng)力對(duì)研究區(qū)井網(wǎng)部署和壓裂施工的影響, 為致密油藏的高效開(kāi)發(fā)提供地質(zhì)依據(jù)。
鄂爾多斯盆地地處中國(guó)中部, 是一個(gè)疊合的克拉通盆地, 即中生代發(fā)育的盆地疊加在古生代盆地之上 (呂文雅等, 2020)。 自三疊系地層沉積以后, 該區(qū)中新生代構(gòu)造應(yīng)力場(chǎng)可以分為4 期, 分別為印支期、 燕山期、 喜馬拉雅期、 新構(gòu)造期,其中, 新構(gòu)造期為中更新世以來(lái)的構(gòu)造應(yīng)力場(chǎng),即現(xiàn)今應(yīng)力場(chǎng)。 該區(qū)以現(xiàn)今最大主應(yīng)力為北東東—南西西方向擠壓和最小主應(yīng)力方向?yàn)楸北蔽鳌夏蠔|方向拉張為主, 現(xiàn)今最大主應(yīng)力方位約為NE70°—NE80°, 且隨著深度增加, 現(xiàn)今最大、 最小主應(yīng)力呈線性增加 (曾聯(lián)波, 2008)。 研究區(qū)合水油田位于鄂爾多斯盆地的西南部 (圖1), 是三疊系延長(zhǎng)組和侏羅系延安組兩個(gè)含油層組的疊合發(fā)育區(qū)。 構(gòu)造上位于伊陜斜坡 (圖1), 總體為向西傾的單斜, 傾角不足1° (吳松濤等, 2015)。 三疊系延長(zhǎng)組自下而上分為10 個(gè)油層組: 長(zhǎng)10—長(zhǎng)1。 研究區(qū)主要發(fā)育長(zhǎng)10—長(zhǎng)2 儲(chǔ)層 (鄧勝徽等,2018), 此次的研究主力層為三疊系延長(zhǎng)組長(zhǎng)6 油層, 屬于湖泊環(huán)境下的辮狀河三角洲前緣亞相沉積, 主要發(fā)育水下分流河道、 河口壩、 分流間灣和席狀砂微相 (白薷等, 2012)。 研究區(qū)長(zhǎng)6 儲(chǔ)層巖性主要為灰色細(xì)—中粒巖屑長(zhǎng)石砂巖夾泥巖、細(xì)—中粒長(zhǎng)石巖屑砂巖和細(xì)—中粒長(zhǎng)石砂巖等,平均孔隙度為8.47%, 平均滲透率為0.11 mD (鞠瑋等, 2014), 屬于典型的致密砂巖儲(chǔ)層。
圖1 研究區(qū)所處構(gòu)造位置示意圖 (據(jù)趙向原等, 2015修改)Fig.1 Structural location of the study area (modified after Zhao et al. , 2015)
地應(yīng)力方向?qū)W(wǎng)的部署具有重要影響, 判斷地應(yīng)力方向的方法包括波速各向異性法、 差應(yīng)變法、 井壁崩落法和微地震法等 (戴俊生等,2014)。 其中, 微地震監(jiān)測(cè)人工裂縫技術(shù)常用于測(cè)試深部的地應(yīng)力信息, 測(cè)試結(jié)果可以反映地下較大范圍的地應(yīng)力信息 (圖2), 被認(rèn)為是目前現(xiàn)今地應(yīng)力測(cè)試的有效方法 (戴俊生等, 2016)。 根據(jù)研究區(qū)兩口井的微地震監(jiān)測(cè)結(jié)果, 得到人工裂縫的延伸方向大致為NE70°—NE80°, 即可推斷水平最大主應(yīng)力的方位為NE70°—NE80°。
單井地應(yīng)力大小的測(cè)量方法有很多種, 包括水壓致裂法、 應(yīng)力解除法、 聲發(fā)射法、 差應(yīng)變法等 (景鋒等, 2008; 解東亮, 2013)。 其中, 水壓致裂法的應(yīng)用十分普遍, 其顯著的優(yōu)點(diǎn)為操作簡(jiǎn)單、 原位測(cè)量、 測(cè)量深度可達(dá)數(shù)千米等 (鄔立等,2009; 朱琳琳等, 2018)。 利用水壓致裂法來(lái)判斷地應(yīng)力的原理為: 在注水壓裂過(guò)程中, 壓裂段孔壁會(huì)出現(xiàn)裂隙, 停止注水后, 裂隙會(huì)立即停止延伸, 在地應(yīng)力場(chǎng)作用下被高壓液體漲破的裂隙會(huì)趨于閉合, 把保持裂隙張開(kāi)時(shí)的平衡壓力稱為瞬時(shí)關(guān)閉壓力, 它與垂直裂隙面的最小水平主應(yīng)力大小相同 (李振剛等, 2018; 曹金鳳等, 2012;印興耀等, 2018), 由此可得到最小水平主應(yīng)力的大小。
根據(jù)研究區(qū)3 口井的水力壓裂資料分析, 通過(guò)水力壓裂得到該井點(diǎn)長(zhǎng)6 儲(chǔ)層的單井水平最大主應(yīng)力值在34 ~37 MPa 之間, 平均值為35.29 MPa,水平最小主應(yīng)力在27 ~32 MPa 之間, 平均值為29.33 MPa, 水平差應(yīng)力在3 ~7 MPa 之間, 平均值為5.96 MPa (表1)。
圖2 研究區(qū)G1 井、 X1 井微地震監(jiān)測(cè)示意圖Fig.2 DiagramofthemicroseismicmonitoringonthewellsG1andX1inthestudyarea
表1 水力壓裂法測(cè)得單井應(yīng)力值Table 1 Values of the present-day in-situ stress in the single wells measured by the hydrofracturing method
應(yīng)力場(chǎng)數(shù)值模擬是研究構(gòu)造應(yīng)力場(chǎng)的一種重要手段, 包括了有限元法、 有限差分法、 邊界元法、 離散元法等 (李玉江等, 2009)。 其基本原理為: 將地質(zhì)體離散成有限個(gè)連續(xù)的單元, 單元之間以結(jié)點(diǎn)相連, 用離散的結(jié)構(gòu)模型去逼近實(shí)際的油氣儲(chǔ)層結(jié)構(gòu), 建立以結(jié)點(diǎn)位移或單元內(nèi)力為未知量, 以整體剛度矩陣為系數(shù)的聯(lián)立方程組, 可求取各節(jié)點(diǎn)主應(yīng)力值的大小和方向, 即儲(chǔ)層應(yīng)力場(chǎng)(劉廣峰等, 2009)。 有限元數(shù)值模擬在其發(fā)展過(guò)程中由二維模擬發(fā)展到三維模擬, 相比較于二維模擬, 三維模擬能更好地體現(xiàn)區(qū)域地應(yīng)力場(chǎng)在垂向上的特征。
為了明確研究區(qū)低滲透儲(chǔ)層地應(yīng)力的空間特征, 文章利用Ansys 軟件對(duì)合水地區(qū)延長(zhǎng)組長(zhǎng)63-1小層進(jìn)行了三維有限元數(shù)值模擬。 首先根據(jù)研究區(qū)實(shí)際構(gòu)造特點(diǎn)、 沉積特點(diǎn)、 巖相特點(diǎn)建立三維非均質(zhì)地質(zhì)模型, 研究區(qū)發(fā)育湖泊環(huán)境下的辮狀河三角洲前緣亞相沉積 (白薷等, 2012), 主要發(fā)育水下分流河道和少部分分流間灣微相, 構(gòu)造上研究區(qū)整體表現(xiàn)為一個(gè)西傾的單斜 (吳松濤等,2015), 構(gòu)造變形較弱, 且無(wú)斷層, 因此研究區(qū)三維非均質(zhì)地質(zhì)模型主要體現(xiàn)在巖性和砂體厚度的差異上 (圖3), 在建立三維非均質(zhì)地質(zhì)模型的基礎(chǔ)上, 通過(guò)高溫高壓三軸巖石力學(xué)實(shí)驗(yàn)和研究區(qū)壓裂施工的實(shí)際數(shù)據(jù), 確定不同巖性及不同砂體厚度的巖石力學(xué)參數(shù) (表2), 研究區(qū)的主要巖性為厚層細(xì)砂巖、 中層粉砂巖、 薄層粉砂質(zhì)泥巖和泥巖, 其中, 對(duì)厚層細(xì)砂巖賦值泊松比0.28, 楊氏模量35.70 GPa; 中層細(xì)砂巖泊松比為0.30, 楊氏模量賦值為32.80 GPa; 薄層泥質(zhì)粉砂巖/粉砂質(zhì)泥巖賦值泊松比0.32, 楊氏模量30.00 GPa; 對(duì)泥巖賦值泊松比0.35, 楊氏模量28.00 GPa; 屬性賦值結(jié)束后建立適合研究區(qū)實(shí)際的三維巖石力學(xué)模型。
圖3 研究區(qū)巖性分布圖 (據(jù)長(zhǎng)慶勘探開(kāi)發(fā)研究院, 2020修改)Fig.3 Diagram of the lithologic distribution in the study area(modified after Petroleum Exploration and Production Research Institute of Changqing Oilfield Company, 2020)
表2 研究區(qū)巖石力學(xué)參數(shù)施加值Table 2 Assignment of values to rock mechanics parameters in the study area
在三維非均質(zhì)地質(zhì)模型和三維力學(xué)模型建立的基礎(chǔ)上, 對(duì)模型進(jìn)行網(wǎng)格劃分, 得到三維有限元網(wǎng)格模型 (圖4)。 由震源機(jī)制解研究可得, 研究區(qū)水平最大主應(yīng)力方向?yàn)楸睎|東—南西西向(盛書(shū)中等, 2015), 同時(shí)結(jié)合單井地應(yīng)力大小和方向的研究, 對(duì)三維力學(xué)模型施加邊界條件 (表3), 其中, 將水平最大主應(yīng)力方向設(shè)為NE70°, 大小為36 MPa; 將最小水平主應(yīng)力方向設(shè)為NW20°,大小為28 MPa; 垂向應(yīng)力在鉛錘方向隨著深度的增加應(yīng)力大小也會(huì)增加, 由于研究區(qū)為一向西傾的單斜, 構(gòu)造變形及起伏較弱, 且?guī)r性主要以細(xì)—粉砂巖為主, 因此對(duì)于垂向應(yīng)力采用依照密度施加, 根據(jù)密度資料得到研究區(qū)巖石的平均密度為2.3 g/cm3。 在邊界條件施加完成之后, 利用Ansys 有限元數(shù)值模擬軟件開(kāi)展三維現(xiàn)今應(yīng)力場(chǎng)數(shù)值模擬。 以研究區(qū)實(shí)測(cè)的地應(yīng)力為驗(yàn)證依據(jù), 通過(guò)不斷調(diào)整模型和邊界條件, 進(jìn)行多次反復(fù)模擬,得到研究區(qū)的地應(yīng)力分布規(guī)律。
圖4 研究區(qū)三維有限元網(wǎng)格模型Fig.4 Diagram of the 3D finite element mesh model of the study area
表3 邊界條件施加值Table 3 Assignment of values to boundary conditions
通過(guò)三維應(yīng)力場(chǎng)數(shù)值模擬, 得到研究區(qū)長(zhǎng)63-1層的地應(yīng)力分布圖 (圖5—圖7), 從模擬結(jié)果來(lái)看, 水平最大主應(yīng)力范圍為34 ~42 MPa; 水平最小主應(yīng)力范圍為25 ~36 MPa; 水平差應(yīng)力范圍為3 ~10 MPa。 平面上水平最大、 最小主應(yīng)力的高值主要分布在細(xì)砂巖區(qū)域, 而粉砂巖區(qū)域和粉砂質(zhì)泥巖區(qū)域應(yīng)力值較低, 研究區(qū)邊緣泥巖區(qū)域應(yīng)力值最低。 垂向上呈現(xiàn)出與平面上較一致的規(guī)律,在厚度較大的細(xì)砂巖區(qū)域, 水平最大主應(yīng)力值較高, 而在厚度小的細(xì)砂巖和薄層的粉砂質(zhì)泥巖及泥巖區(qū)域, 應(yīng)力值降低。 將模擬結(jié)果與壓裂所得的N1 井、 N2 井和Z1 井的地應(yīng)力結(jié)果進(jìn)行對(duì)比(表4), 結(jié)果表明誤差在10.00%以內(nèi), 說(shuō)明模擬的可信度相對(duì)較可好。
表4 長(zhǎng)63-1 層單井應(yīng)力測(cè)量值與應(yīng)力場(chǎng)數(shù)值模擬值對(duì)比表Table 4 Comparison of the simulation results and the measured values of the stress field in the Chang 63-1 Formation
圖5 長(zhǎng)63-1 層最大水平主應(yīng)力分布圖Fig.5 Distribution diagram of the maximum horizontal principle stress in the Chang 63-1 Formation
圖6 長(zhǎng)63-1 層最小水平主應(yīng)力分布圖Fig.6 Distribution diagram of the minimum horizontal principle stress in the Chang 63-1 Formation
圖7 長(zhǎng)63-1 層水平差應(yīng)力分布圖Fig.7 Distribution diagram of the horizontal differential stress in the Chang 63-1 Formation
一般來(lái)說(shuō), 地應(yīng)力分布主要受三個(gè)因素的影響: 一是構(gòu)造格架, 包括斷層產(chǎn)狀和地層厚度變化; 二是區(qū)域應(yīng)力場(chǎng); 三是巖石物理力學(xué)性質(zhì)(譚成軒等, 2006); 由于研究區(qū)構(gòu)造變形較弱,斷層不發(fā)育, 構(gòu)造格架對(duì)現(xiàn)今地應(yīng)力的影響并不大。 主要影響研究區(qū)現(xiàn)今地應(yīng)力分布的因素為巖性與砂體厚度。 巖性對(duì)地應(yīng)力的影響主要體現(xiàn)在巖石的礦物組分、 結(jié)構(gòu)及膠結(jié)程度等方面。 其中,細(xì)砂巖主要為厚層、 中層發(fā)育, 泥質(zhì)粉砂巖/粉砂質(zhì)泥巖為薄層發(fā)育。 由于巖性的不同, 其巖石力學(xué)性質(zhì)也不同, 主要表現(xiàn)為細(xì)砂巖的楊氏模量值最大, 泊松比最小, 而泥巖的楊氏模量最小, 泊松比最大, 進(jìn)而造成了研究區(qū)現(xiàn)今地應(yīng)力分布的差異, 結(jié)果顯示為細(xì)砂巖區(qū)域應(yīng)力值都高于粉砂巖、 粉砂質(zhì)泥巖及泥巖區(qū)。 這表明影響研究區(qū)地應(yīng)力分布的主要因素為巖性, 高強(qiáng)度的巖石較低強(qiáng)度的巖石應(yīng)力值高。
在低滲透儲(chǔ)層壓裂施工過(guò)程中, 地應(yīng)力對(duì)體積壓裂的影響主要體現(xiàn)在: ①最大主應(yīng)力方向影響了單一壓裂縫的延伸方向和長(zhǎng)度 (張志強(qiáng)等,2016); ②水平差應(yīng)力大小是壓裂的關(guān)鍵因素, 其值較小時(shí), 才能形成復(fù)雜的縫網(wǎng), 否則易形成延伸方向與最大主應(yīng)力方向一致的單一壓裂縫 (徐珂等, 2018)。 因此, 結(jié)合研究區(qū)布井情況分析,Z1 井、 Z2 井、 N5 井所處位置的水平差應(yīng)力大致在2 ~4 MPa, 對(duì)這幾口井進(jìn)行體積壓裂更容易形成復(fù)雜的縫網(wǎng), 而N1 井、 N6 井所處位置水平差應(yīng)力較高, 為8 ~10 MPa, 更容易形成單一的壓裂縫, 這些單一的壓裂縫延伸的方向受到水平最大主應(yīng)力方向的控制, 推斷其大致朝著NE70°的方向延伸。 但同時(shí), 壓裂縫的延伸也受到天然裂縫的影響, 因此在布井時(shí), 也要考慮到研究區(qū)天然裂縫的發(fā)育情況。
另一方面, 在水平應(yīng)力差相等的情況下, 所處地區(qū)的水平最大應(yīng)力值越低, 壓裂縫形成所需的起裂壓力就越低 (張來(lái)功, 2016); 因此為盡可能地降低開(kāi)發(fā)成本, 在差應(yīng)力相等的區(qū)域, 油氣工業(yè)井一般部署在應(yīng)力值低的地方。 如N7 井與N1井所在地區(qū), 差應(yīng)力大致相等, 而N1 井所在地區(qū)的最大主應(yīng)力值較低, 因此開(kāi)發(fā)成本相對(duì)N7 井較低一些。
(1) 根據(jù)井下微地震監(jiān)測(cè)分析表明, 研究區(qū)水平最大主應(yīng)力的優(yōu)勢(shì)方向?yàn)镹E70°—NE80°。
(2) 通過(guò)應(yīng)力場(chǎng)三維數(shù)值模擬得到, 水平最大、 最小主應(yīng)力的高值主要分布在研究區(qū)的中部,即砂巖區(qū)域, 而研究區(qū)邊緣區(qū)為泥巖區(qū)域, 應(yīng)力值較低, 也證明了研究區(qū)的應(yīng)力分布主要受巖石力學(xué)性質(zhì)的影響。
(3) 在差應(yīng)力相等的區(qū)域, 盡量選擇水平最大主應(yīng)力值低的地區(qū), 其所需起裂壓力較低, 可以一定程度上降低開(kāi)發(fā)成本。
(4) 在實(shí)際進(jìn)行布井時(shí), 不僅需要考慮現(xiàn)今地應(yīng)力的方向和大小, 同時(shí)要結(jié)合研究區(qū)天然裂縫的發(fā)育情況綜合考慮。