裂縫性油藏天然裂縫動(dòng)靜態(tài)綜合預(yù)測(cè)方法
——以鄂爾多斯盆地華慶油田三疊系長(zhǎng)63儲(chǔ)集層為例

蘇皓，雷征東，張荻萩，李俊超，張澤人，鞠斌山，李治平

（1. 中國(guó)地質(zhì)大學(xué)（北京）能源學(xué)院，北京 100083；2. 中國(guó)石油勘探開(kāi)發(fā)研究院，北京 100083；3. 中國(guó)石油經(jīng)濟(jì)技術(shù)研究院，北京 100724；4. 中國(guó)石油東方地球物理勘探有限公司研究院，河北涿州 072750）

裂縫性油藏天然裂縫動(dòng)靜態(tài)綜合預(yù)測(cè)方法
——以鄂爾多斯盆地華慶油田三疊系長(zhǎng)63儲(chǔ)集層為例

蘇皓1,2，雷征東2，張荻萩3，李俊超2，張澤人4，鞠斌山1，李治平1

（1. 中國(guó)地質(zhì)大學(xué)（北京）能源學(xué)院，北京 100083；2. 中國(guó)石油勘探開(kāi)發(fā)研究院，北京 100083；3. 中國(guó)石油經(jīng)濟(jì)技術(shù)研究院，北京 100724；4. 中國(guó)石油東方地球物理勘探有限公司研究院，河北涿州 072750）

針對(duì)已有裂縫預(yù)測(cè)方法適用范圍小、精度低、對(duì)高成本資料需求程度高的問(wèn)題，以鄂爾多斯盆地華慶油田三疊系長(zhǎng)63儲(chǔ)集層為例，綜合考慮地質(zhì)靜態(tài)資料和生產(chǎn)動(dòng)態(tài)資料，提出了一種預(yù)測(cè)新方法?；谟吞锒喾N地質(zhì)靜態(tài)資料，得到巖性、沉積相、地層厚度、巖石破裂指數(shù)、裂縫發(fā)育強(qiáng)度 5個(gè)約束裂縫發(fā)育的控制條件，應(yīng)用多元線性回歸分析方法，建立了裂縫發(fā)育的 5種約束條件與裂縫密度關(guān)系的定量表達(dá)式，計(jì)算得到全區(qū)的裂縫密度體；根據(jù)井史、示蹤劑、干擾試井及吸水剖面等生產(chǎn)動(dòng)態(tài)資料，通過(guò)油藏工程分析方法綜合判斷出裂縫在平面上和縱向上的方向及分布范圍，并結(jié)合數(shù)值模擬技術(shù)對(duì)裂縫密度體進(jìn)行檢驗(yàn)和定量校正，建立了既符合地質(zhì)靜態(tài)認(rèn)識(shí)又符合生產(chǎn)動(dòng)態(tài)的三維離散裂縫地質(zhì)模型。數(shù)值模擬擬合驗(yàn)證表明，該裂縫模型擬合度高，具較高的可靠性和適用性。圖15表1參32

裂縫性油藏；天然裂縫；裂縫地質(zhì)模型；裂縫預(yù)測(cè)；華慶油田；鄂爾多斯盆地

0 引言

對(duì)于裂縫性油田水驅(qū)開(kāi)發(fā)，大量注入水沿裂縫水竄嚴(yán)重影響了油田的水驅(qū)開(kāi)發(fā)效果，成為制約油田開(kāi)發(fā)的一個(gè)重要因素。裂縫作為油氣滲流的重要通道，對(duì)特低、超低滲透儲(chǔ)集層的開(kāi)發(fā)至關(guān)重要，研究裂縫分布規(guī)律、定量預(yù)測(cè)裂縫發(fā)育程度是油氣勘探開(kāi)發(fā)中的難題[1-2]。天然裂縫的預(yù)測(cè)方法目前主要有3類：①基于構(gòu)造本身的結(jié)構(gòu)特征來(lái)探討構(gòu)造體與裂縫發(fā)育的關(guān)系[3-4]，如構(gòu)造主曲率法預(yù)測(cè)裂縫、分析斷層與裂縫位置關(guān)系等，該種方法主要適用于預(yù)測(cè)斷層、褶皺構(gòu)造伴生的裂縫；②基于裂縫的構(gòu)造成因，應(yīng)用構(gòu)造應(yīng)力有限元模擬技術(shù)和破裂法則預(yù)測(cè)裂縫分布[5-6]，由于數(shù)學(xué)模型的局限性，此方法預(yù)測(cè)結(jié)果與實(shí)際仍有一定出入；③基于高精度的地震資料[7-11]，提取裂縫發(fā)育相關(guān)屬性參數(shù)進(jìn)行分析預(yù)測(cè)，如地震相干體方法、多波地震、螞蟻?zhàn)粉櫟确椒?，該方法建立的裂縫模型精度較高，是目前的主流方法，但由于地震資料的數(shù)量和分辨率問(wèn)題，以及地震解釋的人為因素，仍難以得到反映實(shí)際儲(chǔ)集層的裂縫模型。

綜上所述，裂縫發(fā)育控制因素的多樣性使得裂縫發(fā)育具有很強(qiáng)的隨機(jī)性和非均質(zhì)性。目前的裂縫預(yù)測(cè)方法大多建立在靜態(tài)地質(zhì)資料基礎(chǔ)上，但由于地質(zhì)資料的數(shù)量少、精度低、獲取成本高等問(wèn)題，得到的裂縫模型具有很大的不確定性和不適用性，不能很好地用于指導(dǎo)油田開(kāi)發(fā)，對(duì)構(gòu)造不發(fā)育、沒(méi)有地震數(shù)據(jù)資料的地區(qū)更是如此。本文針對(duì)上述問(wèn)題，以鄂爾多斯盆地華慶油田長(zhǎng) 63儲(chǔ)集層為例，基于該區(qū)多種地質(zhì)上和生產(chǎn)上的靜、動(dòng)態(tài)資料，在常規(guī)裂縫建模方法的基礎(chǔ)上，考慮多種裂縫發(fā)育控制因素，實(shí)現(xiàn)多條件相互約束，并結(jié)合油田開(kāi)發(fā)中的認(rèn)識(shí)和經(jīng)驗(yàn)，綜合得到該區(qū)裂縫模型，使之既滿足地質(zhì)認(rèn)識(shí)又符合開(kāi)發(fā)動(dòng)態(tài)，為后期油田水驅(qū)調(diào)整政策提供有效依據(jù)。

1 區(qū)域地質(zhì)背景

1.1 區(qū)域構(gòu)造特征

華慶油田位于鄂爾多斯盆地伊陜斜坡中南部（見(jiàn)圖1），其構(gòu)造演化與華北板塊的關(guān)系十分密切[12]。盆地在區(qū)域上表現(xiàn)為向西南傾斜，由于應(yīng)力強(qiáng)度的差異，盆地外部邊緣構(gòu)造作用活躍，受力強(qiáng)度大，斷裂、破裂較為發(fā)育，而盆地內(nèi)部則相反。研究區(qū)塊處于盆地內(nèi)部，地勢(shì)平緩，傾角不足1°，構(gòu)造相對(duì)簡(jiǎn)單，斷裂、褶皺發(fā)育較少，局部地區(qū)發(fā)育有鼻狀隆起構(gòu)造帶[13]。

圖1 華慶油田位置圖

1.2 地應(yīng)力場(chǎng)特征

根據(jù)構(gòu)造運(yùn)動(dòng)的強(qiáng)度及時(shí)間，將鄂爾多斯盆地主要構(gòu)造應(yīng)力場(chǎng)分為 4期：印支期、燕山期、喜馬拉雅期、新構(gòu)造期[14]。其中燕山期和喜馬拉雅期是裂縫形成的主要時(shí)期，燕山期構(gòu)造地應(yīng)力水平擠壓作用造成斷裂發(fā)育，喜馬拉雅期構(gòu)造運(yùn)動(dòng)主要造成地層的抬升，對(duì)延長(zhǎng)組裂縫進(jìn)行改造[15]。劉格云應(yīng)用赤平投影方法[16]，得出燕山期構(gòu)造應(yīng)力場(chǎng)最大主應(yīng)力方向?yàn)楸蔽魑鳌蠔|東向，優(yōu)勢(shì)方位為 116°～296°，最小主應(yīng)力優(yōu)勢(shì)方位則為26°～206°；喜馬拉雅期最大主應(yīng)力方向北東—南西向，方位為45°～225°，最小主應(yīng)力優(yōu)勢(shì)方位為 135°～315°。

1.3 儲(chǔ)集層沉積特征和物性特征

研究區(qū)長(zhǎng)6段自上而下分為3個(gè)亞段：長(zhǎng)63、長(zhǎng)62、長(zhǎng) 61，各亞段巖性相同，主要為砂巖、粉砂巖、泥巖等沉積巖，其中長(zhǎng)63亞段為主力層位。長(zhǎng)63亞段以砂質(zhì)碎屑流沉積為主，伴有滑塌和濁流沉積的重力流沉積體系，主要巖性為粉細(xì)—細(xì)砂巖。研究區(qū)主要發(fā)育砂質(zhì)碎屑流、濁流和半深湖—深湖3種沉積亞相，鑒于砂質(zhì)碎屑流分布面積較廣，為了進(jìn)一步研究裂縫發(fā)育規(guī)律，將砂質(zhì)碎屑流亞相按厚度和物性的差異進(jìn)一步劃分為主體和邊緣兩個(gè)相帶。

長(zhǎng)63儲(chǔ)集層整體上表現(xiàn)為低孔、低滲的特征，孔滲分布范圍較大，分選差，具有很強(qiáng)的非均質(zhì)性。全區(qū)孔隙度為6%～16%，平均10.8%；全區(qū)覆壓滲透率分布范圍跨度較大，為（0.04～0.60）×10-3μm2，平均0.34×10-3μm2，為超低滲儲(chǔ)集層。

1.4 天然裂縫發(fā)育特征

在前人認(rèn)識(shí)的基礎(chǔ)上[17-18]，根據(jù)研究區(qū)延河剖面露頭及成像測(cè)井等資料分析發(fā)現(xiàn)，長(zhǎng) 63儲(chǔ)集層裂縫發(fā)育優(yōu)勢(shì)方位主要為北東東向、北東向，少量裂縫為近南北向，與裂縫形成期的古應(yīng)力場(chǎng)方向基本一致。

對(duì)研究區(qū)66口取心井進(jìn)行巖心觀察統(tǒng)計(jì)，發(fā)現(xiàn)其中有60口取心井的巖心中都發(fā)育天然裂縫，存在天然裂縫的巖心井占90.9%，裂縫開(kāi)度為30～100 μm，平均值為60 μm，主要發(fā)育小開(kāi)度的微裂縫系統(tǒng)。

從成因來(lái)看，研究區(qū)天然裂縫可分為構(gòu)造裂縫和成巖裂縫。構(gòu)造裂縫占絕大多數(shù)，主要分布于砂巖中，伴有少量礦物充填，裂縫面較為平整，可以觀察到擦痕、階步等天然裂縫特征，構(gòu)造裂縫以剪切縫為主，張應(yīng)力和剪應(yīng)力聯(lián)合作用下形成的張扭縫少量發(fā)育；研究區(qū)還發(fā)育平行層面分布的成巖縫，這種裂縫的規(guī)模較小，橫向連續(xù)性差，且在上覆巖層壓力的作用下通常為閉合狀態(tài)，因此本文在預(yù)測(cè)井間裂縫時(shí)不考慮成巖裂縫。

從裂縫傾角來(lái)看，天然裂縫以近垂直的高角度縫為主，斜交裂縫與低角度裂縫發(fā)育較少，高角度縫所占比例高達(dá)90%。

進(jìn)一步對(duì)取心井的裂縫密度進(jìn)行分析計(jì)算，可知全區(qū)裂縫平均密度為0.76條/m，其中61%的井裂縫發(fā)育密度達(dá)到0.5條/m以上，因此認(rèn)為該區(qū)天然裂縫較為發(fā)育，但由于儲(chǔ)集層的強(qiáng)非均質(zhì)性，不同取心井的裂縫密度差異較大，在區(qū)域上北區(qū)總體較南區(qū)裂縫密度大。

2 基于靜態(tài)地質(zhì)資料預(yù)測(cè)裂縫

華慶油田地勢(shì)平緩，構(gòu)造簡(jiǎn)單，斷層、復(fù)雜褶皺等構(gòu)造地質(zhì)體發(fā)育較少，且油田缺乏高精度地震資料，因此只能從裂縫的構(gòu)造成因著手，對(duì)裂縫分布進(jìn)行預(yù)測(cè)。構(gòu)造縫的發(fā)育主要受構(gòu)造應(yīng)力、巖性、地層厚度、沉積相等地質(zhì)因素控制，可通過(guò)分析每個(gè)因素與裂縫發(fā)育的關(guān)系，找出控制裂縫發(fā)育的約束條件，并綜合多種因素實(shí)現(xiàn)對(duì)裂縫在三維空間分布的約束和預(yù)測(cè)。

2.1 裂縫發(fā)育的約束條件

2.1.1 裂縫與巖性、地層厚度、沉積相的關(guān)系

巖性、地層厚度、沉積相等地質(zhì)因素是影響裂縫發(fā)育的內(nèi)因。巖性不同，巖石的成分、孔隙結(jié)構(gòu)也不同，在相同的構(gòu)造作用下，巖性不同是裂縫發(fā)育程度存在差異的主要因素[19]；室內(nèi)實(shí)驗(yàn)、巖心觀察及野外調(diào)查統(tǒng)計(jì)表明，裂縫發(fā)育受地層厚度的控制[20]，相對(duì)于厚巖層，薄巖層更容易產(chǎn)生密集的裂縫；沉積相則通過(guò)控制儲(chǔ)集層的巖性、層厚及巖石組合來(lái)控制裂縫發(fā)育程度。

本文基于露頭及巖心觀察資料，統(tǒng)計(jì)了裂縫在不同巖性、不同厚度地層和不同沉積相中的發(fā)育頻率，并繪制了裂縫發(fā)育頻率與各因素的關(guān)系圖。由圖 2可見(jiàn)，裂縫主要發(fā)育在石英、方解石、白云石等脆性礦物含量較高的砂巖中，脆性礦物含量低的泥巖中裂縫發(fā)育較差；研究區(qū)裂縫平均間距與地層厚度呈較好的線性關(guān)系（見(jiàn)圖 3），隨著地層厚度增大，裂縫平均間距增加，即裂縫密度變??；裂縫在砂質(zhì)碎屑流邊緣相帶最發(fā)育，砂質(zhì)碎屑流主體相帶次之（見(jiàn)圖 4），濁流和湖相沉積中泥質(zhì)巖類比例較高，裂縫發(fā)育程度明顯低于砂質(zhì)碎屑流。

圖2 研究區(qū)長(zhǎng)63儲(chǔ)集層裂縫發(fā)育頻率與巖性關(guān)系圖

圖3 研究區(qū)長(zhǎng)63儲(chǔ)集層地層厚度與裂縫間距關(guān)系圖

2.1.2 裂縫與構(gòu)造應(yīng)力的關(guān)系

構(gòu)造應(yīng)力是構(gòu)造縫形成的外因，是巖石破裂的觸發(fā)機(jī)制、構(gòu)造縫發(fā)育的主控因素。為定量表征構(gòu)造應(yīng)力與裂縫發(fā)育的關(guān)系，運(yùn)用有限元數(shù)值模擬技術(shù)[21]，對(duì)裂縫形成期的構(gòu)造應(yīng)力場(chǎng)分布進(jìn)行模擬。

圖4 研究區(qū)長(zhǎng)63儲(chǔ)集層裂縫發(fā)育頻率與沉積相關(guān)系圖

2.1.2.1 模型的建立

構(gòu)造應(yīng)力場(chǎng)模型主要包括地質(zhì)模型、力學(xué)模型和數(shù)學(xué)模型[22]。

地質(zhì)模型的建立主要包括地質(zhì)體的選取和邊界條件的確定。由于研究區(qū)無(wú)明顯斷層發(fā)育，裂縫發(fā)育主要受巖性、地層厚度及巖性組合的影響，而這些因素主要受沉積作用的控制，因此地質(zhì)模型中地質(zhì)體主要依據(jù)沉積相分布選取。初始受力邊界條件則根據(jù)裂縫形成期（燕山期和喜馬拉雅期）的最大主應(yīng)力和最小主應(yīng)力的方向、大小分別設(shè)定，最終的邊界條件需要經(jīng)過(guò)多次的正反演來(lái)確定。

力學(xué)模型主要根據(jù)地質(zhì)體力學(xué)性質(zhì)和巖石力學(xué)參數(shù)確定。根據(jù)巖石力學(xué)實(shí)驗(yàn)，目的層巖石力學(xué)性質(zhì)表現(xiàn)為脆性，因此將其作為彈脆性材料處理；巖石力學(xué)參數(shù)主要考慮了不同地質(zhì)單元體的彈性模量、泊松比、黏聚力、內(nèi)摩擦角等，這些參數(shù)由巖石力學(xué)實(shí)驗(yàn)得到。

數(shù)學(xué)模型建立的依據(jù)為有限元原理，其基本思路是將連續(xù)的地質(zhì)體離散為有限個(gè)通過(guò)節(jié)點(diǎn)相連的連續(xù)單元，給每個(gè)單元賦予實(shí)際的巖石力學(xué)參數(shù)，再根據(jù)邊界受力條件和節(jié)點(diǎn)平衡條件，即可計(jì)算出每個(gè)單元內(nèi)應(yīng)力與應(yīng)變的近似值。

2.1.2.2 模擬結(jié)果

模擬結(jié)果主要為最大主應(yīng)力、最小主應(yīng)力和差應(yīng)力。研究區(qū)長(zhǎng)63儲(chǔ)集層在燕山期構(gòu)造應(yīng)力場(chǎng)的最大主應(yīng)力方向?yàn)楸蔽魑鳌蠔|東向，最小主應(yīng)力為北北東—南南西向，最大主應(yīng)力為70～85 MPa，最小主應(yīng)力為20～30 MPa，差應(yīng)力為50～65 MPa；喜馬拉雅期構(gòu)造應(yīng)力場(chǎng)的最大主應(yīng)力方向?yàn)楸睎|—南西向，最小主應(yīng)力為北西—南東向，其中最大主應(yīng)力為60～80 MPa，最小主應(yīng)力為13～20 MPa，差應(yīng)力為40～58 MPa。該區(qū)最大和最小主應(yīng)力方位分布比較穩(wěn)定，主要受遠(yuǎn)場(chǎng)應(yīng)力和邊界條件控制，應(yīng)力變化較大。

根據(jù)研究區(qū)最大和最小主應(yīng)力，結(jié)合格里菲斯、庫(kù)倫-莫爾巖石破裂準(zhǔn)則[23]，能夠定性地判斷巖石中能否發(fā)生張剪破裂，并引入巖石破裂指數(shù)I，其定義為：

式中I——巖石破裂指數(shù)，無(wú)因次；σ——張應(yīng)力，MPa；σ′——巖石抗張強(qiáng)度，MPa；τ——剪應(yīng)力，MPa；τ′——巖石抗剪強(qiáng)度，MPa。

由（1）式、（2）式計(jì)算出研究區(qū)各部位的巖石破裂指數(shù)（見(jiàn)圖5），該指數(shù)可定量表征巖石受構(gòu)造應(yīng)力作用后達(dá)到破裂狀態(tài)的程度，即全區(qū)各部位形成裂縫的趨勢(shì)。

圖5 研究區(qū)長(zhǎng)63儲(chǔ)集層巖石破裂指數(shù)分布圖

2.1.3 裂縫發(fā)育強(qiáng)度

研究區(qū)成像測(cè)井資料較少，但大多數(shù)井都有常規(guī)測(cè)井資料。常規(guī)測(cè)井資料包含地層巖性、物性、導(dǎo)電性等多種信息數(shù)據(jù)，大多數(shù)據(jù)都能反映地層裂縫的發(fā)育特征，但不同測(cè)井曲線對(duì)裂縫的響應(yīng)特征不同，因此具很強(qiáng)的多解性。綜合概率指數(shù)法[24-25]能夠提取出有效反映裂縫特征的曲線信息來(lái)綜合識(shí)別裂縫。該方法在分析多種測(cè)井參數(shù)曲線與裂縫的相關(guān)性后，根據(jù)各個(gè)參數(shù)的影響程度確定權(quán)系數(shù)，并將各參數(shù)綜合計(jì)算得到反映裂縫在縱向上相對(duì)發(fā)育程度的綜合指標(biāo)，即裂縫發(fā)育強(qiáng)度，其值越大則說(shuō)明裂縫越發(fā)育，反之，裂縫相對(duì)不發(fā)育?？梢杂昧芽p發(fā)育強(qiáng)度來(lái)定量分析裂縫分布特征。

利用該方法對(duì)研究區(qū) 166口井的常規(guī)測(cè)井曲線（雙側(cè)向電阻率、中子孔隙度、自然電位、聲波時(shí)差、孔隙結(jié)構(gòu)指數(shù)等）進(jìn)行分析計(jì)算，得到每口井的裂縫發(fā)育強(qiáng)度曲線。與實(shí)際巖心資料對(duì)比，其符合率達(dá)到80%以上（見(jiàn)圖 6）。采用協(xié)同克里金方法，對(duì)井間區(qū)域進(jìn)行插值計(jì)算，得到全區(qū)裂縫發(fā)育強(qiáng)度屬性體（見(jiàn)圖7），也將其作為控制該區(qū)裂縫發(fā)育程度的約束條件。

圖6 部分井裂縫發(fā)育強(qiáng)度曲線與巖心裂縫對(duì)比圖

圖7 研究區(qū)長(zhǎng)63儲(chǔ)集層裂縫發(fā)育強(qiáng)度分布圖

2.2 裂縫密度體的計(jì)算

綜上，研究區(qū)長(zhǎng)63儲(chǔ)集層裂縫的發(fā)育與巖性、地層厚度、沉積相、巖石破裂指數(shù)以及裂縫發(fā)育強(qiáng)度 5種屬性有關(guān)。采用數(shù)學(xué)分析方法，分析出裂縫密度與這 5種相關(guān)屬性的定量關(guān)系，即可通過(guò)相關(guān)屬性計(jì)算出裂縫密度體，從而實(shí)現(xiàn)裂縫在三維空間的預(yù)測(cè)。定量相關(guān)性分析方法有多種，本文采用較為快捷、分析結(jié)果較為準(zhǔn)確的多元線性回歸分析預(yù)測(cè)方法[26-27]。

結(jié)合地質(zhì)建模，將巖性（x1）、地層厚度（x2）、沉積相（x3）、巖石破裂指數(shù)（x4）、裂縫發(fā)育強(qiáng)度（x5）5種約束條件屬性體網(wǎng)格化到同一模型中，并將由巖心等資料得出的單井裂縫密度曲線也重采樣到同一套網(wǎng)格，因此裂縫井所經(jīng)過(guò)的每個(gè)網(wǎng)格中就包含以上 5種約束條件以及裂縫密度（y）這6套屬性體，即可對(duì)單井上的網(wǎng)格進(jìn)行裂縫密度與 5種屬性的關(guān)系分析，以實(shí)現(xiàn)對(duì)井間網(wǎng)格裂縫密度的計(jì)算。從60口裂縫井中篩選出上述6套屬性值均大于0的網(wǎng)格作為分析樣品數(shù)據(jù)，這樣既保證了樣品的數(shù)量又保證了樣品的可靠性。將樣品的巖性、地層厚度、沉積相、破裂指數(shù)、裂縫發(fā)育強(qiáng)度作為自變量，裂縫密度作為因變量，對(duì)數(shù)據(jù)進(jìn)行多元線性回歸，得到裂縫發(fā)育密度與相關(guān)屬性的線性回歸模型為：

回歸模型的復(fù)相關(guān)系數(shù)為 0.918，決定系數(shù)為0.843，修正決定系數(shù)0.721，表明模型擬合程度較好；檢驗(yàn)回歸方程的統(tǒng)計(jì)量值為20.273，顯著性值為0，說(shuō)明該模型高度顯著，可信度高。（3）式中巖石破裂指數(shù)的回歸系數(shù)最高，也反映了構(gòu)造應(yīng)力是影響該區(qū)裂縫發(fā)育的主控因素。利用以上模型計(jì)算可得到全區(qū)裂縫密度屬性體（見(jiàn)圖8）。

圖8 研究區(qū)長(zhǎng)63儲(chǔ)集層裂縫密度分布圖

3 基于生產(chǎn)動(dòng)態(tài)資料校正裂縫密度體

由于裂縫密度體的計(jì)算完全從地質(zhì)規(guī)律的角度出發(fā)，所以往往與生產(chǎn)資料不能較好地吻合，如有些區(qū)域計(jì)算的裂縫密度很小，而實(shí)際此區(qū)域的油井在生產(chǎn)上卻表現(xiàn)出裂縫性水淹的特征，因此需要利用生產(chǎn)動(dòng)態(tài)資料對(duì)該密度體進(jìn)行檢驗(yàn)和校正。

3.1 利用動(dòng)態(tài)資料識(shí)別裂縫的原理

裂縫性油藏在進(jìn)行注水開(kāi)發(fā)時(shí)，由于裂縫滲透率遠(yuǎn)大于基質(zhì)，注入水將首先沿著高滲透裂縫通道流動(dòng)。因此可根據(jù)這一特性，結(jié)合生產(chǎn)動(dòng)態(tài)資料，運(yùn)用油藏工程方法判斷裂縫在平面上和縱向上的發(fā)育情況，從而對(duì)裂縫密度屬性體進(jìn)行檢驗(yàn)校正。

平面上，裂縫性水竄造成對(duì)應(yīng)油井層內(nèi)方向性見(jiàn)水。結(jié)合井史、干擾試井、示蹤劑等資料，根據(jù)油水井注水動(dòng)態(tài)響應(yīng)特征及示蹤劑采出曲線綜合判斷油井產(chǎn)出水的來(lái)源，即是源于注入水、原生水亦或是裂縫性水竄等，識(shí)別出裂縫性水竄在平面上的方向和分布范圍。

縱向上，注入水會(huì)沿裂縫發(fā)育的高滲透層段突進(jìn)，而低滲透層段則較少吸水或不吸水，因此裂縫較發(fā)育層位的相對(duì)吸水量就大，在吸水剖面上形成尖峰狀，這種具有尖峰狀且大規(guī)模吸水的層段很有可能就是裂縫發(fā)育段或水淹層位。因此可利用吸水剖面測(cè)試資料判斷裂縫在縱向上的發(fā)育情況。

3.2 裂縫密度體的校正

綜合上述裂縫識(shí)別方法，對(duì)研究區(qū)所有水淹井組進(jìn)行動(dòng)態(tài)分析，識(shí)別出平面上裂縫水竄方向和分布范圍，并將裂縫性水竄的水線標(biāo)定在平面井位圖上（見(jiàn)圖 9）；并對(duì)每條水線對(duì)應(yīng)注水井的吸水剖面資料進(jìn)行分析，判斷出縱向上裂縫性吸水的層段，得到注水井的裂縫性吸水剖面統(tǒng)計(jì)表（見(jiàn)表 1）。利用平面和縱向上分析結(jié)果即可判斷出裂縫在儲(chǔ)集層三維空間發(fā)育的范圍。

圖9 研究區(qū)長(zhǎng)63儲(chǔ)集層裂縫性水竄水線平面分布圖

結(jié)合生產(chǎn)動(dòng)態(tài)分析結(jié)果，利用數(shù)值模擬法可對(duì)地質(zhì)方法得到的裂縫密度體進(jìn)行校正[28-29]?；舅悸窞椋豪么Ｕ牧芽p密度體進(jìn)行建模得到裂縫地質(zhì)模型后，根據(jù)裂縫開(kāi)度、密度、延伸長(zhǎng)度等參數(shù)，計(jì)算出裂縫的滲透率、孔隙度、形狀因子等屬性，將這些屬性參數(shù)與基質(zhì)模型結(jié)合，即可建立油藏?cái)?shù)值模擬模型進(jìn)行數(shù)模研究，通過(guò)調(diào)整裂縫密度值，擬合生產(chǎn)歷史數(shù)據(jù)，確定符合實(shí)際生產(chǎn)動(dòng)態(tài)的裂縫密度值和分布，實(shí)現(xiàn)對(duì)密度體的定量校正。在超低滲透油藏中，裂縫密度的大小對(duì)單井產(chǎn)量、含水率等生產(chǎn)指標(biāo)的影響通常十分顯著（見(jiàn)圖10）。據(jù)此，重點(diǎn)對(duì)動(dòng)態(tài)分析標(biāo)定的裂縫性水竄區(qū)域進(jìn)行小范圍裂縫密度擬合校正，最終得到一個(gè)擬合度較高的裂縫密度體（見(jiàn)圖11）。

表1 研究區(qū)長(zhǎng)63儲(chǔ)集層部分裂縫井吸水剖面統(tǒng)計(jì)表

圖10 裂縫密度對(duì)生產(chǎn)參數(shù)的影響

將校正前后的裂縫密度體進(jìn)行對(duì)比，發(fā)現(xiàn)在研究區(qū)北區(qū)中部、西南部及南區(qū)個(gè)別區(qū)域存在顯著的差別，這些區(qū)域內(nèi)的裂縫性水竄特征在校正前的密度體中均未能體現(xiàn)，而校正后的密度體不僅在其他區(qū)域仍保持基于地質(zhì)認(rèn)識(shí)的計(jì)算結(jié)果，在這些區(qū)域內(nèi)也體現(xiàn)了由生產(chǎn)動(dòng)態(tài)認(rèn)識(shí)得到的水竄范圍，提高了裂縫模型的可靠性和適用性。

3.3 三維離散裂縫地質(zhì)模型

3.3.1 離散裂縫模型

以巖心觀察等資料統(tǒng)計(jì)分析的單井裂縫分布作為單井上的約束條件，以校正后的三維裂縫密度屬性體作為預(yù)測(cè)井間裂縫的約束條件，并結(jié)合不同部位裂縫的產(chǎn)狀，進(jìn)行離散裂縫建模，即可得到研究區(qū)天然裂縫的三維模型（見(jiàn)圖12）。模型中裂縫主要走向?yàn)榻鼥|西向，其次為北東方向，近南北向裂縫較少，3組天然裂縫均具有很強(qiáng)的各向異性，受儲(chǔ)集層非均質(zhì)性影響，裂縫分布規(guī)律性較差，在研究區(qū)的北部地區(qū)裂縫發(fā)育情況較好。

圖11 校正后的研究區(qū)長(zhǎng)63儲(chǔ)集層裂縫密度分布圖

3.3.2 裂縫滲透率模型

由Oda公式可知[30-31]，裂縫滲透率的大小由裂縫尺寸、裂縫密度等參數(shù)計(jì)算得到。利用建模軟件對(duì)研究區(qū)裂縫滲透率進(jìn)行計(jì)算，可得到全區(qū)的裂縫滲透率分布模型（見(jiàn)圖13）。研究區(qū)裂縫滲透率為（500～2 000）×10-3μm2，局部地區(qū)滲透率大于 2 000×10-3μm2。

3.3.3 裂縫模型的驗(yàn)證

含水率擬合是裂縫性水驅(qū)油藏中重要的歷史擬合指標(biāo)之一。由于常規(guī)的儲(chǔ)集層建模方法通常沒(méi)有考慮裂縫等因素，利用建立的模型進(jìn)行數(shù)值模擬時(shí)，常通過(guò)修改基質(zhì)滲透率、傳導(dǎo)率、相滲曲線或者表皮系數(shù)等方法來(lái)擬合含水率等歷史數(shù)據(jù)，但這種方法往往擬合效果不佳[32]，難以模擬裂縫性水驅(qū)油藏的復(fù)雜見(jiàn)水情況。

圖12 研究區(qū)長(zhǎng)63儲(chǔ)集層天然裂縫離散網(wǎng)絡(luò)模型

圖13 研究區(qū)長(zhǎng)63儲(chǔ)集層裂縫滲透率模型

利用數(shù)值模擬方法對(duì)比常規(guī)模型、校正前及校正后裂縫模型的全區(qū)含水率擬合曲線（見(jiàn)圖14）。通過(guò)設(shè)置相同模擬參數(shù)對(duì)比發(fā)現(xiàn)，常規(guī)方法含水率擬合的效果最差，這是由于該模型不能很好地考慮和模擬裂縫，使注入水失去了竄流通道，整體含水率偏低；校正前的裂縫模型僅從地質(zhì)角度考慮了裂縫分布，與生產(chǎn)實(shí)際仍存在矛盾，擬合效果與實(shí)際仍有一定偏差；本文建立的裂縫模型根據(jù)實(shí)際資料綜合考慮了地質(zhì)上的認(rèn)識(shí)規(guī)律和生產(chǎn)動(dòng)態(tài)分析的結(jié)果，全區(qū)含水率擬合度最高，與實(shí)際值偏差最小。

圖14 不同方法擬合區(qū)塊含水率效果對(duì)比

此外，將歷史擬合后得到的全區(qū)含水飽和度分布圖、裂縫分布圖以及單井日產(chǎn)量擬合圖（以區(qū)塊中部典型水淹井組為例）進(jìn)行綜合對(duì)比，發(fā)現(xiàn)非裂縫性水淹井（低含水井）和裂縫性水淹井（高含水井）的單井日產(chǎn)量數(shù)據(jù)擬合程度均較高，且高含水井的來(lái)水方向均與動(dòng)態(tài)分析結(jié)論一致（見(jiàn)圖15）。綜上所述，認(rèn)為本文建立的裂縫模型與地質(zhì)上及開(kāi)發(fā)上的認(rèn)識(shí)吻合，能更準(zhǔn)確地表征裂縫分布，很好地應(yīng)用于實(shí)際生產(chǎn)。

圖15 研究區(qū)中部含水飽和度、裂縫分布和單井日產(chǎn)量擬合曲線綜合對(duì)比圖

4 結(jié)論

華慶油田長(zhǎng)63儲(chǔ)集層的天然裂縫較為發(fā)育，大部分裂縫為近垂直的高角度縫，裂縫開(kāi)度小，發(fā)育方向主要為北東東向，部分近北東方向以及少量近南北向，平均裂縫密度為0.76條/m，北區(qū)總體較南區(qū)裂縫相對(duì)更為發(fā)育。

采用動(dòng)靜態(tài)結(jié)合的方法對(duì)天然裂縫進(jìn)行預(yù)測(cè)。根據(jù)油田多種地質(zhì)靜態(tài)資料，分析得到約束裂縫發(fā)育的5個(gè)控制條件：巖性、沉積相、地層厚度、巖石破裂指數(shù)、裂縫發(fā)育強(qiáng)度；通過(guò)多元線性回歸分析，得到 5種約束屬性與裂縫密度的定量關(guān)系式，計(jì)算出全區(qū)的裂縫密度體；根據(jù)井史、示蹤劑、干擾試井及吸水剖面等生產(chǎn)動(dòng)態(tài)資料，利用油藏工程動(dòng)態(tài)分析方法綜合判斷出裂縫在平面上和縱向上的方向及分布范圍，在動(dòng)態(tài)分析的裂縫分布范圍內(nèi)結(jié)合數(shù)值模擬方法對(duì)裂縫密度體進(jìn)行擬合校正，確定符合實(shí)際生產(chǎn)動(dòng)態(tài)的裂縫密度值，并以校正后的裂縫密度體構(gòu)建天然裂縫模型。動(dòng)靜態(tài)結(jié)合預(yù)測(cè)方法得到的裂縫模型既滿足地質(zhì)靜態(tài)認(rèn)識(shí)又符合生產(chǎn)動(dòng)態(tài)，提高了裂縫預(yù)測(cè)的可靠性和適用性。

[1]袁士義, 宋新民, 冉啟全. 裂縫性油藏開(kāi)發(fā)技術(shù)[M]. 北京: 石油工業(yè)出版社, 2004.YUAN Shiyi, SONG Xinmin, RAN Qiquan. Exploitation technology of fractured reservoirs[M]. Beijing: Petroleum Industry Press, 2004.

[2]曾聯(lián)波, 柯式鎮(zhèn), 劉洋. 低滲透油氣儲(chǔ)層裂縫研究方法[M]. 北京:石油工業(yè)出版社, 2010.ZENG Lianbo, KE Shizhen, LIU Yang. Research method of fracture in low permeability reservoir[M]. Beijing: Petroleum Industry Press,2010.

[3]彭紅利, 熊鈺, 孫良田, 等. 主曲率法在碳酸鹽巖氣藏儲(chǔ)層構(gòu)造裂縫預(yù)測(cè)中的應(yīng)用研究[J]. 天然氣地球科學(xué), 2005, 16(3):343-346.PENG Hongli, XIONG Yu, SUN Liangtian, et al. The research on main curvature method application in carbonate gas reservoir’s structural fracture prediction[J]. Natural Gas Geoscience, 2005,16(3): 343-346.

[4]曾錦光, 羅元華, 陳太源. 應(yīng)用構(gòu)造面主曲率研究油氣藏裂縫問(wèn)題[J]. 力學(xué)學(xué)報(bào), 1982, 18(2): 94-98.ZENG Jinguang, LUO Yuanhua, CHEN Taiyuan. Oil and gas reservoir fracture study by principal curvature of structural surface[J].Chinese Journal of Theoretical and Applied Mechanics, 1982, 18(2):94-98.

[5]張守仁, 萬(wàn)天豐, 陳建平. 川西坳陷孝泉—新場(chǎng)地區(qū)須家河組二—四段構(gòu)造應(yīng)力場(chǎng)模擬及裂縫發(fā)育區(qū)帶預(yù)測(cè)[J]. 石油與天然氣地質(zhì), 2004, 25(1): 70-74, 80.ZHANG Shouren, WAN Tianfeng, CHEN Jianping. Tectonic stress field modeling and fracture prediction in T3x2-4strata in Xiaoquan-Xinchang area, western Sichuan Depression[J]. Oil and Gas Geology, 2004, 25(1): 70-74, 80.

[6]孫曉慶. 古構(gòu)造應(yīng)力場(chǎng)有限元數(shù)值模擬的應(yīng)用及展望[J]. 斷塊油氣田, 2008, 15(3): 31-33.SUN Xiaoqing. Present situation and prospect of application for finite element numerical simulation of paleotectonic stress fields[J].Fault-Block Oil & Gas Field, 2008, 15(3): 31-33.

[7]李陽(yáng), 侯加根, 李永強(qiáng). 碳酸鹽巖縫洞型儲(chǔ)集體特征及分類分級(jí)地質(zhì)建模[J]. 石油勘探與開(kāi)發(fā), 2016, 43(4): 600-606.LI Yang, HOU Jiagen, LI Yongqiang. Features and hierarchical modeling of carbonate fracture-cavity reservoirs[J]. Petroleum Exploration and Development, 2016, 43(4): 600-606.

[8]劉振峰, 曲壽利, 孫建國(guó), 等. 地震裂縫預(yù)測(cè)技術(shù)研究進(jìn)展[J]. 石油物探, 2012, 51(2): 191-198.LIU Zhenfeng, QU Shouli, SUN Jianguo, et al. Progress of seismic fracture characterization technology[J]. Geophysical Prospecting for Petroleum, 2012, 51(2): 191-198.

[9]侯加根, 馬曉強(qiáng), 劉鈺銘, 等. 縫洞型碳酸鹽巖儲(chǔ)層多類多尺度建模方法研究: 以塔河油田四區(qū)奧陶系油藏為例[J]. 地學(xué)前緣,2012, 19(2): 59-66.HOU Jiagen, MA Xiaoqiang, LIU Yuming, et al. Modeling of carbonate fracture-vuggy reservoir: A case study of Ordovician reservoir of 4th Block in Tahe Oilfield[J]. Earth Science Frontiers,2012, 19(2): 59-66.

[10]郎曉玲, 郭召杰. 基于DFN離散裂縫網(wǎng)絡(luò)模型的裂縫性儲(chǔ)層建模方法[J]. 北京大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 49(6): 964-972.LANG Xiaoling, GUO Zhaojie. Fractured reservoir modeling method based on discrete fracture network model[J]. Acta Scientiarum Naturalium Universitatis Pekinensis, 2013, 49(6): 964-972.

[11]龍旭, 武林芳. 螞蟻?zhàn)粉檶傩泽w提取參數(shù)對(duì)比試驗(yàn)及其在塔河四區(qū)裂縫建模中的應(yīng)用[J]. 石油天然氣學(xué)報(bào), 2011, 33(5): 76-81.LONG Xu, WU Linfang. Parameter comparative experiments on ant tracking attribute extraction and its application in fracture modeling in Region 4 of Tahe Oilfield[J]. Journal of Oil and Gas Technology,2011, 33(5): 76-81.

[12]楊華, 陳洪德, 付金華. 鄂爾多斯盆地晚三疊世沉積地質(zhì)與油藏分布規(guī)律[M]. 北京: 科學(xué)出版社, 2012.YANG Hua, CHEN Hongde, FU Jinhua. Law of sedimentary geology and petroleum distribution in late Triassic Ordos Basin[M]. Beijing:Science Press, 2012.

[13]曾聯(lián)波. 低滲透砂巖儲(chǔ)層裂縫的形成與分布[M]. 北京: 科學(xué)出版社, 2008.ZENG Lianbo. Formation and distribution of fractures in low permeability sandstone reservoirs[M]. Beijing: Science Press,2008.

[14]張泓, 孟召平, 何宗蓮. 鄂爾多斯煤盆地構(gòu)造應(yīng)力場(chǎng)研究[J]. 煤炭學(xué)報(bào), 2000, 25(增刊): 1-5.ZHANG Hong, MENG Zhaoping, HE Zonglian. Study on the tectonic stress fields in the Ordos Coal Basin[J]. Journal of China Coal Society, 2000, 25(Supp.): 1-5.

[15]周新桂, 張林炎, 黃臣軍, 等. 華慶地區(qū)長(zhǎng) 63儲(chǔ)層裂縫分布模型與裂縫有效性[J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 2012, 43(3):689-697.ZHOU Xingui, ZHANG Linyan, HUANG Chenjun, et al. Distraction network conceptual model and validity of fractures in Chang 63low permeable reservoir in Huaqing area[J]. Journal of Jilin University(Earth Science Edition), 2012, 43(3): 689-697.

[16]劉格云. 長(zhǎng)慶油田華慶地區(qū)延長(zhǎng)組長(zhǎng) 63儲(chǔ)層裂縫發(fā)育特征研究[D]. 蘭州: 蘭州大學(xué), 2010: 7-14.LIU Geyun. Research on fracture developing characteristics of Chang 63reservoir, Yanchang group in Huaqing area, Changqing petroleum field[D]. Lanzhou: Lanzhou University, 2010: 7-14.

[17]周新桂, 張林炎, 黃臣軍. 華慶探區(qū)長(zhǎng) 63儲(chǔ)層破裂壓力及裂縫開(kāi)啟壓力估測(cè)與開(kāi)發(fā)建議[J]. 中南大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013,44(7): 2812-2818.ZHOU Xingui, ZHANG Linyan, HUANG Chenjun. Estimation of formation breakdown pressure and fracture open pressure of Chang 63low permeable reservoir in Huaqing area and development suggestions[J]. Journal of Central South University (Science and Technology), 2013, 44(7): 2812-2818.

[18]劉格云, 黃臣軍, 周新桂, 等. 鄂爾多斯盆地三疊系延長(zhǎng)組裂縫發(fā)育程度定量評(píng)價(jià)[J]. 石油勘探與開(kāi)發(fā), 2015, 42(4): 444-453.LIU Geyun, HUANG Chenjun, ZHOU Xingui, et al. Quantitative evaluation of fracture development in Triassic Yanchang Formation,Ordos Basin, NW China[J]. Petroleum Exploration and Development,2015, 42(4): 444-453.

[19]周永勝, 張流. 裂縫性儲(chǔ)集層巖心裂縫統(tǒng)計(jì)分析[J]. 世界地質(zhì),2000, 19(2): 117-124.ZHOU Yongsheng, ZHANG Liu. Statistics and analysis of drilling core fractures of fractured reservoir[J]. World Geology, 2000, 19(2):117-124.

[20]WU H. POLLARD D D. An experiment study of the relationship between joint spacing and layer thickness[J]. Journal of Structural Geology, 1995, 17(6): 887-905.

[21]夏建芳, 葉南海. 有限元法原理與 ANSYS應(yīng)用[M]. 北京: 國(guó)防工業(yè)出版社, 2011.XIA Jianfang, YE Nanhai. Principle of finite element method and application of ANSYS[M]. Beijing: National Defense Industry Press,2011.

[22]曾聯(lián)波, 文世鵬, 肖淑蓉. 低滲透油氣儲(chǔ)集層裂縫空間分布的定量預(yù)測(cè)[J]. 勘探家, 1998(2): 24-26.ZENG Lianbo, WEN Shipeng, XIAO Shurong. Quantitative prediction of fracture space distribution in low permeability reservoir[J]. Petroleum Explorationist, 1998(2): 24-26.

[23]鄧攀, 魏國(guó)齊, 楊泳. 儲(chǔ)層構(gòu)造裂縫定量預(yù)測(cè)中地質(zhì)數(shù)學(xué)模型的建立與應(yīng)用研究[J]. 天然氣地球科學(xué), 2006, 17(4): 480-484.DENG Pan, WEI Guoqi, YANG Yong. To establish and research of three-dimensional geological and mathematical model for quantitative prediction of structural fracture[J]. Natural Gas Geoscience, 2006, 17(4): 480-484.

[24]王兆生, 曾聯(lián)波, 張振國(guó), 等. 巴喀油田裂縫發(fā)育特征及常規(guī)測(cè)井識(shí)別方法[J]. 新疆地質(zhì), 2012, 30(3): 359-362.WANG Zhaosheng, ZENG Lianbo, ZHANG Zhenguo, et al. Feature of fracture development and evaluate fractures by conventional logging data in Baka Oilfield[J]. Xinjiang Geology, 2012, 30(3):359-362.

[25]LYU Wenya, ZENG Lianbo, LIU Zhongqun, et al. Fracture responses of conventional logs in tight-oil sandstones: A case study of the upper Triassic Yanchang formation in southwest Ordos Basin,China[J]. AAPG Bulletin, 2016, 100(9): 1399-1417.

[26]張景陽(yáng), 潘光友. 多元線性回歸與BP神經(jīng)網(wǎng)絡(luò)預(yù)測(cè)模型對(duì)比與運(yùn)用研究[J]. 昆明理工大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 38(6): 61-67.ZHANG Jingyang, PAN Guangyou. Comparison and application of multiple regression and BP neural network prediction model[J].Journal of Kunming University of Science and Technology (Natural Science Edition), 2013, 38(6): 61-67.

[27]王惠文, 孟潔. 多元線性回歸的預(yù)測(cè)建模方法[J]. 北京航空航天大學(xué)學(xué)報(bào), 2007, 33(4): 500-504.WANG Huiwen, MENG Jie. Predictive modeling on multivariate linear regression[J]. Journal of Beijing University of Aeronautics and Astronautics, 2007, 33(4): 500-504.

[28]樊建明, 魏忠元, 何永宏, 等. 超低滲透油藏考慮天然裂縫的數(shù)值模擬擬合新思路[J]. 低滲透油氣田, 2009, 14 (3): 83-85.FAN Jianming, WEI Zhongyuan, HE Yonghong, et al. Consider the ultra-low permeability reservoir numerical simulation of the natural cracks in fitting new ideas[J]. Low Permeability Oil & Gas Fields,2009, 14(3): 83-85.

[29]樊建明, 屈雪峰, 王沖, 等. 鄂爾多斯盆地致密儲(chǔ)集層天然裂縫分布特征及有效裂縫預(yù)測(cè)新方法[J]. 石油勘探與開(kāi)發(fā), 2016, 43(5):740-748.FAN Jianming, QU Xuefeng, WANG Chong, et al. Natural fracture distribution and a new method predicting effective fractures in tight oil reservoirs of Ordos Basin, NW China[J]. Petroleum Exploration and Development, 2016, 43(5): 740-748.

[30]ODA M. Permeability tensor for discontinuous rock mass[J].Geotechnique, 1985, 35(4): 483-495.

[31]DERSHOWITZ B, LAPOINTE P, EIBEN T, et al. Integration of discrete feature network methods with conventional simulator approaches[J]. SPE Reservoir Evaluation & Engineering, 1998, 3(2):165-170.

[32]池建萍, 鄭強(qiáng), 祁軍, 等. 復(fù)雜裂縫性油藏歷史擬合中的特殊做法[J]. 新疆石油地質(zhì), 2004, 25(5): 517-519.CHI Jianping, ZHENG Qiang, QI Jun, et al. The special approach for history match of a complex fractured reservoir[J]. Xinjiang Petroleum Geology, 2004, 25(5): 517-519.

Dynamic and static comprehensive prediction method of natural fractures in fractured oil reservoirs: A case study of Triassic Chang 63reservoirs in Huaqing Oilfield, Ordos Basin, NW China

SU Hao1,2, LEI Zhengdong2, ZHANG Diqiu3, LI Junchao2, ZHANG Zeren4, JU Binshan1, LI Zhiping1
(1.School of Energy Resource,China University of Geosciences,Beijing100083,China; 2.PetroChina Research Institute of Petroleum Exploration & Development,Beijing100083,China; 3.PetroChina Research Institute of Economics & Technology,Beijing100724,China; 4.Geophysical Research Institute,Bureau of Geophysical Prospecting,PetroChina,Zhuozhou072750,China)

In consideration of the limited adaptability scope, low accuracy and high demand of great cost data of existent fracture prediction methods, a new fracture predicting method was advanced by implementing geological static data and production dynamic data from the Triassic Chang 63reservoirs in the Huaqing Oilfield. Five constraints, lithology, sedimentary facies, thickness, rock rupture index and fracture intensity controlling the development of fractures were sorted out based on the static geological data. The multiple linear regression method was adopted to work out the quantitative relationships between the five constraints and fracture density, and the fracture density property of the whole area was calculated. Based on production dynamic data of well history, tracer, well interference test and intake profile test, the direction and distribution of fracture horizontally and vertically were figured out by reservoir engineering analysis method. The fracture density property was verified and quantitatively corrected with numerical simulation, and a 3D discrete fracture geological model in agreement with both geological cognition and dynamic production performance was built. The numerical simulation shows that the fracture model has higher fitting consistency, high reliability and adaptability.

fractured oil reservoir; natural fracture; fracture geological model; fracture prediction; Huaqing Oilfield; Ordos Basin

國(guó)家科技重大專項(xiàng)（2017ZX05013-002）

TE344

A

1000-0747(2017)06-0919-11

10.11698/PED.2017.06.09

蘇皓, 雷征東, 萩張荻 , 等. 裂縫性油藏天然裂縫動(dòng)靜態(tài)綜合預(yù)測(cè)方法: 以鄂爾多斯盆地華慶油田三疊系長(zhǎng)63儲(chǔ)集層為例[J]. 石油勘探與開(kāi)發(fā), 2017, 44(6): 919-929.

SU Hao, LEI Zhengdong, ZHANG Diqiu, et al. Dynamic and static comprehensive prediction method of natural fractures in fractured oil reservoirs: A case study of Triassic Chang 63reservoirs in Huaqing Oilfield, Ordos Basin, NW China[J]. Petroleum Exploration and Development, 2017, 44(6): 919-929.

蘇皓（1990-），男，廣西南寧人，碩士，中國(guó)地質(zhì)大學(xué)（北京）博士研究生，主要從事低滲透油藏開(kāi)發(fā)方面的研究工作。地址：北京市海淀區(qū)學(xué)院路29號(hào)，中國(guó)地質(zhì)大學(xué)（北京）能源學(xué)院，郵政編碼：100083。E-mail：suhao0912@qq.com

聯(lián)系作者簡(jiǎn)介：雷征東（1979-），男，重慶市人，博士，中國(guó)勘探開(kāi)發(fā)研究院高級(jí)工程師，主要從事油藏?cái)?shù)值模擬、低滲透油氣田開(kāi)發(fā)等領(lǐng)域的研究工作。地址：北京市海淀區(qū)學(xué)院路20號(hào)，中國(guó)石油勘探開(kāi)發(fā)研究院油田開(kāi)發(fā)研究所，郵政編碼：100083。E-mail：leizhengdong@petrochina.com.cn

2016-12-26

2017-09-18

（編輯 郭海莉）