王 宏，王蓓蕾，王成俊

（陜西延長(zhǎng)石油（集團(tuán)）有限責(zé)任公司研究院）

“兩步法”相控隨機(jī)建模在延長(zhǎng)油田儲(chǔ)層描述中的應(yīng)用

王 宏，王蓓蕾，王成俊

（陜西延長(zhǎng)石油（集團(tuán)）有限責(zé)任公司研究院）

延長(zhǎng)油田低滲透砂巖儲(chǔ)層的地質(zhì)建模具有儲(chǔ)層描述資料少、精度低、不確定性因素多的特點(diǎn)。文中采用“兩步法”相控隨機(jī)建模技術(shù)，以密井網(wǎng)測(cè)井解釋資料為基礎(chǔ)，結(jié)合巖心分析、生產(chǎn)動(dòng)態(tài)等已知信息，先建立相控巖相模型，再以此為約束條件進(jìn)行孔、滲屬性建模，來(lái)減少儲(chǔ)層描述的不確定性。以延長(zhǎng)油田BLC油區(qū)為例，應(yīng)用“兩步法”相控隨機(jī)建模技術(shù)，建立了符合該油區(qū)實(shí)際生產(chǎn)狀態(tài)的低滲透砂巖地質(zhì)模型，為區(qū)塊合理開(kāi)發(fā)決策和改善油田開(kāi)發(fā)效果提供了地質(zhì)依據(jù)。

相控建模；隨機(jī)建模；儲(chǔ)層描述；延長(zhǎng)油田

延長(zhǎng)油田BLC油區(qū)位于鄂爾多斯盆地的東北部，屬于安塞三角洲沉積體系，物源方向?yàn)楸睎|—南西向。該區(qū)油藏主要受構(gòu)造控制，局部也受到巖性的影響，多是以構(gòu)造控制為主的巖性-構(gòu)造復(fù)合油藏。長(zhǎng)22儲(chǔ)層碎屑巖體系表現(xiàn)為碎屑成分以長(zhǎng)石為主、石英次之，顆粒大小不均勻，層內(nèi)、層間滲透性差異較大，顆粒粒度偏細(xì)。由于孔、滲參數(shù)受沉積相帶、成巖壓實(shí)和后生成巖變化的影響，導(dǎo)致砂巖粗細(xì)分布不均，儲(chǔ)層非均質(zhì)性嚴(yán)重[1］。三維地質(zhì)建模技術(shù)可以對(duì)儲(chǔ)層的非均質(zhì)性進(jìn)行高密度的定量表征，提供井間儲(chǔ)層參數(shù)的米級(jí)或10米級(jí)的變化[2］?！皟刹椒ā毕嗫仉S機(jī)建模技術(shù)能夠降低隨機(jī)模擬算法的不確定性，所建立的相模型與孔、滲模型在空間上相互匹配，與所研究的地質(zhì)現(xiàn)象相吻合，可以輸入到油藏?cái)?shù)值模擬中，以提高油藏動(dòng)態(tài)預(yù)測(cè)的可靠性。

1 “兩步法”相控隨機(jī)建模

1.1 隨機(jī)模擬

隨機(jī)模擬是指利用已知的信息，以隨機(jī)函數(shù)為算法，應(yīng)用隨機(jī)模擬方法，人工合成可選、等概率和高精度的反映變量空間分布的模型，其能夠反映出參數(shù)的細(xì)微變化[3］。隨機(jī)模擬算法對(duì)優(yōu)化方案設(shè)計(jì)、降低風(fēng)險(xiǎn)提供了重要的參考資料，但不能給出最優(yōu)的地質(zhì)模型結(jié)果[4］。

1.2 相控隨機(jī)建模

為了降低隨機(jī)模擬的不確定性，應(yīng)給隨機(jī)模擬算法加上約束條件，如沉積微相平面分布范圍等，這樣就會(huì)大大降低儲(chǔ)層預(yù)測(cè)的不確定性[5］。沉積微相圖代表沉積單元在不同沉積環(huán)境時(shí)的平面分布狀況，平面上沉積環(huán)境不同會(huì)導(dǎo)致儲(chǔ)層性質(zhì)差異較大。相控建模就是基于此發(fā)展起來(lái)的一項(xiàng)盡可能降低儲(chǔ)層預(yù)測(cè)的不確定性的建模技術(shù)[6］。

相控隨機(jī)建模是以密井網(wǎng)測(cè)井解釋資料為基礎(chǔ)，以沉積微相平面分布圖作為邊界條件，采用隨機(jī)建模算法建立儲(chǔ)層巖相及其內(nèi)部孔、滲屬性的模型[7］。但由于同一微相內(nèi)部，垂向上的巖相差異較大且孔、滲屬性值域較寬，導(dǎo)致相控條件下的巖相模型與孔、滲屬性模型匹配性較差，即好的巖相內(nèi)部往往出現(xiàn)孔、滲屬性低值，而差的巖相內(nèi)部則出現(xiàn)孔、滲屬性高值。

1.3 “兩步法”相控隨機(jī)建模

為了解決相控隨機(jī)建模存在的巖相模型與孔、滲屬性模型匹配性較差的問(wèn)題，可采用“兩步法”相控隨機(jī)建模技術(shù)。相控建模時(shí)，先將巖相這一離散型變量改為連續(xù)型變量，以平面上不同微相內(nèi)部的巖相值域范圍作為約束條件建立相控巖相模型，然后建立該模型約束下的孔、滲屬性模型。此種建模方法不僅與所研究的地質(zhì)現(xiàn)象吻合，還能避免大多數(shù)連續(xù)變量模型對(duì)于平穩(wěn)性或均質(zhì)性的嚴(yán)格要求[8］。

2 研究區(qū)儲(chǔ)層特征

BLC油田產(chǎn)層為長(zhǎng)22油層，油層平均埋深為750 m，平均孔隙度為16%，平均含油飽和度為55%，平均滲透率為17.9 mD。油藏主體為CaCl2水型，部分為NaHCO3水型。

2.1 地層精細(xì)對(duì)比

BLC油田位于鄂爾多斯盆地陜北斜坡中東部[9］，含油層系屬長(zhǎng)22亞段。利用自然伽馬、自然電位、聲波時(shí)差的曲線形態(tài)分析該區(qū)213口井，通過(guò)剖面閉合對(duì)比，將目的層段長(zhǎng)22亞段劃分為5個(gè)沉積單元，即 5 個(gè)小層，分別為長(zhǎng) 221、長(zhǎng) 222、長(zhǎng) 223、長(zhǎng) 224、長(zhǎng)225小層。

2.2 儲(chǔ)層“四性”關(guān)系特征

2.2.1 巖性特征

根據(jù)巖心觀察，含油砂巖均為細(xì)砂巖，結(jié)合試采資料表明，BLC油田儲(chǔ)層巖性下限為細(xì)砂巖，即細(xì)砂巖以上儲(chǔ)層可作為有效儲(chǔ)層。

在2007年所取的87.46 m長(zhǎng)的巖心中，共取得油跡以上油砂41.13 m，其中油跡砂巖6.42 m，占油砂總數(shù)的16%；油斑砂巖28.3 m，占油砂總數(shù)的68%；油浸砂巖6.41 m，占油砂總數(shù)的16%。結(jié)合鄂爾多斯盆地同類油田對(duì)比分析，BLC油田儲(chǔ)層含油巖性下限為油斑細(xì)砂巖。

2.2.2 物性特征

該區(qū)儲(chǔ)層物性分析資料統(tǒng)計(jì)表明，油藏屬中孔、低滲儲(chǔ)層類型。孔隙度分布于4.0%～20.7%，集中分布于14%～20%，平均為16.1%，屬中孔儲(chǔ)層；滲透率分布于0.1～218.0 mD，集中分布于10～50 mD，平均為36.4 mD，屬低滲儲(chǔ)層。有效儲(chǔ)層物性下限為：孔隙度8%、滲透率0.43 mD。

儲(chǔ)層孔隙度與滲透率具有明顯的正相關(guān)性，其相關(guān)系數(shù)達(dá)到0.82，即隨著儲(chǔ)層孔隙度的增大，滲透率值也相應(yīng)增大（圖1）。

圖1 孔隙度與滲透率關(guān)系圖Fig.1 The relationship between porosity and permeability

2.2.3 孔隙度測(cè)井解釋模型

建立孔隙度測(cè)井解釋模型，可以看出聲波時(shí)差與孔隙度二者之間具有較高的正相關(guān)性（圖2），其測(cè)井解釋模型為式中：φ 為孔隙度，%；Δt為聲波時(shí)差，μs·m-1。

圖2 聲波時(shí)差與孔隙度關(guān)系圖Fig.2 The relationship between acoustic time and porosity

2.2.4 儲(chǔ)層含油性與電性關(guān)系

（1）有效儲(chǔ)層的確定

通過(guò)測(cè)井解釋或試油得出，有效儲(chǔ)層分為油層、油水同層、含油水層、水層，無(wú)效儲(chǔ)層為干層。根據(jù)試油規(guī)范，結(jié)合該區(qū)試采狀況，確定了層屬性的劃分標(biāo)準(zhǔn)：①油層，含水率＜40%；②油水同層，40%＜含水率＜80%；③含油水層，80%＜含水率＜95%；④水層，含水率＞95%。

試采資料分析表明，該區(qū)生產(chǎn)井以油水同層居多，其次是含油水層和油層，極少數(shù)為水層。

（2）含油性與電性關(guān)系

地層水礦化度的大小決定了油層電性響應(yīng)特征的高低。BLC油區(qū)地層水礦化度分布特征研究表明（圖3）：該油區(qū)存在2種地層水水型，一種為分布于斷層以北的構(gòu)造主體部位的CaCl2水型，礦化度為40 000～55 000 mg/L；另一種為分布于斷層以南的構(gòu)造東南端的NaHCO3水型，礦化度＜10 000 mg/L，2種水型的礦化度相差約5倍。

圖3 BLC油區(qū)地層水礦化度平面分布圖Fig.3 The planar distribution of reservoir water salinity in BLC oil region

根據(jù)射孔層段的電性特征以及試采層段的屬性特點(diǎn)，選取了歸一化處理后的深感應(yīng)電阻率Rild和地層電阻率Rt作為判別油水層的依據(jù)，按照地層水礦化度的高低分別建立測(cè)井解釋綜合圖，測(cè)井解釋標(biāo)準(zhǔn)見(jiàn)表1。

表1 不同地層水礦化度下的測(cè)井解釋標(biāo)準(zhǔn)Table 1 The standard of well logging interpretation underdifferent reservoir water salinity

圖4 高礦化度地區(qū)（左）和低礦化度地區(qū)（右）測(cè)井綜合解釋圖Fig.4 The comprehensive well logging interpretation in high salinity region （left） and low salinity region （right）

2.3 BLC油區(qū)沉積特征

延長(zhǎng)油田BLC油區(qū)沉積體系具有寬廣的三角洲平原相，為曲流河—三角洲平原亞相發(fā)育區(qū)，發(fā)育河道沉積[10］。BLC油區(qū)長(zhǎng)22亞段沉積模式主要為三角洲平原相沉積，曲流河河控平原亞相，沉積微相以河道砂微相為主，少量為堤岸微相[11-13］。

（1）河道砂微相。該微相沉積細(xì)砂巖，塊狀層理、無(wú)粒序構(gòu)造，顆粒磨圓為次棱—次圓，分選中等；單層厚度較大，一般在5～7 m，常為多期河道疊加形成；物性相對(duì)較好，平均孔隙度為14%～16%，是該油田主要儲(chǔ)集層。測(cè)井曲線特征為：自然電位、自然伽馬表現(xiàn)為箱狀結(jié)構(gòu)，層內(nèi)均質(zhì)程度較高，除非存在泥巖、鈣質(zhì)夾層；聲波時(shí)差曲線一般不發(fā)生漸變過(guò)程，而是垂向加積的過(guò)程。

（2）堤岸微相。該微相為水流過(guò)大時(shí)溢出河道的沉積產(chǎn)物，或在曲流河凸岸水體變淺、流速變緩時(shí)的產(chǎn)物，有時(shí)統(tǒng)稱為漫灘沉積，常發(fā)育交互水平層理。巖性主要是薄層細(xì)砂（厚度2～3 mm）與碳質(zhì)灰?guī)r，其次為質(zhì)純的灰色泥巖，沉積厚度一般為20～30 cm，夾于厚層砂巖之間，以隔夾層的形式出現(xiàn)，一般為一期河道沉積末期的產(chǎn)物。自然伽馬曲線反應(yīng)比較明顯，經(jīng)測(cè)試為非滲透層，水平滲透率一般小于0.2 mD，這對(duì)油氣的分布與聚集起到了一定的隔檔作用。

結(jié)合以上分析進(jìn)行巖相劃分，并繪制長(zhǎng)22儲(chǔ)層沉積微相平面圖（圖5）。

圖5 BLC油區(qū)長(zhǎng)22儲(chǔ)層沉積微相平面圖Fig.5 The sedimentary microfacies of Chang 22reservoir in BLC oil region

3 三維地質(zhì)模型的建立

3.1 三維構(gòu)造模型

精細(xì)的構(gòu)造模型是儲(chǔ)層描述的基礎(chǔ)[14］，它可為后續(xù)沉積微相模型和儲(chǔ)層物性參數(shù)模型提供邊界信息。構(gòu)造模型是根據(jù)斷層模型和層面模型建立起來(lái)的，斷層模型主要表現(xiàn)為單個(gè)斷裂面形態(tài)刻畫和多個(gè)斷層空間配置的擬合關(guān)系；層面模型反映的是地層界面的三維分布，疊合的層面模型即為地層格架模型[15-16］。

本文所建構(gòu)造模型采用了10 m×10 m×0.5 m的網(wǎng)格，建立BLC長(zhǎng)22地層構(gòu)造模型。從構(gòu)造模型來(lái)看（圖6），BLC油區(qū)構(gòu)造形態(tài)與區(qū)域構(gòu)造一致，為相對(duì)平緩向西傾斜的單斜構(gòu)造，油田范圍內(nèi)東西海拔高度差僅30 m。同時(shí)由于差異壓實(shí)作用，形成向西延伸多條長(zhǎng)度不一的小型鼻狀隆起，隆起幅度一般在4～6 m，隆起寬度一般在1～2個(gè)井距，延伸長(zhǎng)度為2～4個(gè)井距。該區(qū)構(gòu)造東南部的Z169井與Z4井之間發(fā)育一條小型高角度正斷層，斷面傾向東南，傾角約70°，斷距約40～50 m。

圖6 BLC油區(qū)長(zhǎng)22亞段構(gòu)造模型Fig.6 The tectonic model of Chang 22reservoir in BLC oil region

3.2 沉積相模型

以前文的沉積相分析和所繪制的沉積微相平面分布圖作為約束條件，采用“兩步法”相控隨機(jī)建模技術(shù)進(jìn)行巖相建模，建立BLC油區(qū)長(zhǎng)22油藏相控巖相三維模型（圖7）。該巖相模型描述了砂體的空間展布，河道與物源方向均為北東—南西向，沿河道方向砂體連續(xù)發(fā)育，縱向上單砂體相互疊置，自下而上砂體厚度逐漸增大、連片性逐漸增強(qiáng)，為建設(shè)性三角洲平原沉積。

圖7 BLC油區(qū)長(zhǎng)22亞段巖相模型Fig.7 The lithofacies model of Chang 22reservoir in BLC oil region

3.3 孔、滲屬性模型

孔隙度和滲透率，二者的分布規(guī)律受巖相的影響較大，且多呈正態(tài)分布。在建立孔、滲屬性模型時(shí)，以巖相為約束條件，對(duì)加載的屬性數(shù)據(jù)進(jìn)行地質(zhì)統(tǒng)計(jì)與數(shù)據(jù)分析，利用分析結(jié)果建立孔隙度和滲透率模型，能夠確保模型的真實(shí)可靠性。從圖8可以看出，孔隙度模型、滲透率模型是在沉積相邊界的約束下建立的，并且滲透率與孔隙度有很好的相關(guān)性。

結(jié)合以上建立的巖相與屬性模型可以得出，河道砂微相主體部位砂體厚度大、孔隙度和滲透率值較高，是該區(qū)有利的儲(chǔ)集砂體；平面上砂體連片發(fā)育，呈條帶狀分布，砂體連續(xù)性和連通性好。這與油田的生產(chǎn)認(rèn)識(shí)是一致的，能夠?yàn)橛吞锷a(chǎn)提供可靠的地質(zhì)依據(jù)。

圖8 BLC油區(qū)長(zhǎng)22儲(chǔ)層孔隙度模型（a）和滲透率模型（b）Fig.8 The porosity model（a） and permeability model（b） of Chang 22reservoir in BLC oil region

4 結(jié)論

（1）“兩步法”相控隨機(jī)建模技術(shù)，是先建立相控巖相模型，再以此為空間約束條件進(jìn)行孔、滲屬性建模，這樣不但可解決相控隨機(jī)建模的弊病，而且能更加準(zhǔn)確地反映實(shí)際地質(zhì)情況。

（2）BLC油區(qū)長(zhǎng)22亞段沉積模式主要為三角洲平原相沉積，曲流河河控平原亞相，沉積微相以河道砂微相為主，少量為堤岸微相。

（3）應(yīng)用“兩步法”相控隨機(jī)建模技術(shù)在延長(zhǎng)油田BLC油區(qū)建立的儲(chǔ)層地質(zhì)模型，其巖相模型與屬性模型在平面與縱向上的匹配性好，能夠?yàn)樵搮^(qū)油藏增儲(chǔ)上產(chǎn)和深度開(kāi)發(fā)調(diào)整提供地質(zhì)依據(jù)。

[1]宋子齊，唐長(zhǎng)久，林樹文，等.隨機(jī)建模在克拉瑪依油田礫巖儲(chǔ)層描述中的應(yīng)用[J］.特種油氣藏，2006，13（6）：25-28.

[2]宋子齊，伊軍鋒，龐振宇，等.三維儲(chǔ)層地質(zhì)建模與砂礫油層挖潛研究——以克拉瑪依油田七中區(qū)、七東區(qū)克拉瑪依組礫巖油藏為例[J］.巖性油氣藏，2007，19（4）：13-16.

[3]呂坐彬，趙春明，霍春亮，等.精細(xì)相控儲(chǔ)層地質(zhì)建模技術(shù)在老油田調(diào)整挖潛中的應(yīng)用——以綏中36-1油田為例[J］.巖性油氣藏，2010，22（3）：77-79.

[4]李少華，尹艷樹，張昌民.儲(chǔ)層隨機(jī)建模系列技術(shù)[M］.北京：石油工業(yè)出版社，2007：68-69.

[5]張熠，竇松江，賈玉梅，等.測(cè)井重構(gòu)約束反演在大港油田CH二區(qū)儲(chǔ)層建模中的應(yīng)用[J］.中國(guó)石油勘探，2010，15（6）：64-67.

[6]胡光義，古莉，王福利.礁灰?guī)r特征和塊狀體儲(chǔ)層建模方法探討[J］.巖性油氣藏，2007，19（2）：60-64.

[7]謝銳杰，張艷文，秦剛，等.相控建模在紅連油田精細(xì)油藏描述中的應(yīng)用[J］.巖性油氣藏，2011，23（1）：86-89.

[8]文華，孫娜，田華豐.基于灰色建模法的Logistic模型在油氣田開(kāi)發(fā)中的應(yīng)用[J］.巖性油氣藏，2009，21（3）：40-44.

[9]李士春，馮朝榮，殷世江.鄂爾多斯盆地南部中生界延長(zhǎng)組沉積體系與油氣富集[J］.巖性油氣藏，2010，22（2）：24-34.

[10]宋凱，呂劍文，杜金良，等.鄂爾多斯盆地中部上三疊統(tǒng)延長(zhǎng)組物源方向分析與三角洲沉積體系[J］.古地理學(xué)報(bào)，2002，4（3）：59-66.

[11]趙靖舟，武富禮，閆世可，等.陜北斜坡東部三疊系油氣富集規(guī)律研究[J］.石油學(xué)報(bào)，2006，27（5）：24-34.

[12]王華，鄭榮才，周祺，等.鄂爾多斯盆地長(zhǎng)北氣田山2段三角洲沉積體系和砂體展布特征[J］.巖性油氣藏，2008，20（2）：65-69.

[13]蘇進(jìn)昌，張嵐，馬新福.河流相儲(chǔ)層開(kāi)發(fā)初期地質(zhì)建模[J］.巖性油氣藏，2008，20（3）：114-117.

[14]尤少燕，時(shí)應(yīng)敏，喬玉雷.利用多尺度資料進(jìn)行油藏地質(zhì)建?！苑?24區(qū)塊為例[J］.油氣地質(zhì)與采收率，2005，12（3）：15-17.

[15]張永慶，代開(kāi)梅，陳舒薇.砂質(zhì)辮狀河儲(chǔ)層三維地質(zhì)建模研究[J］.大慶石油地質(zhì)與開(kāi)發(fā)，2002，21（5）：34-36.

[16]楊仁超.儲(chǔ)層地質(zhì)學(xué)研究新進(jìn)展[J］.特種油氣藏，2006，13（4）：1-5.

Application of“two-step” facies-constrained stochastic modeling to reservoir description in Yanchang Oilfield

WANG Hong， WANG Bei-lei， WANG Cheng-jun
（Research Institute of Shaanxi Yanchang Petroleum Co.Ltd.， Xi’an 710075， China）

The geological modelingfor lowpermeabilitysandstone reservoir in YanchangOilfield is

characterized by lack of data,poor accuracy degree and more uncertainty factors for reservoir description.Based on densewelllogginginterpretation,the“two-step”facies-constrainedstochasticmodelingmethodcombiningwith core analysis and production performance are used toestablish the facies model,and then carryout the porosityand the permeability modeling,in order todecrease the uncertainty for reservoir description.Taking BLC oil region in Yanchang Oilfield as an example,the“two-step”facies-constrained stochastic modeling method is applied to establish the lowpermeability sandstone reservoir model suitable for the actual production status in this area,which provides geological foundation for reasonable development decision and improvingdevelopment effect.

facies-constrained modeling；stochastic modeling；reservoir description；YanchangOilfield

TE348

A

1673-8926（2011）04-0110-05

2011-03-21；

2011-05-20

陜西延長(zhǎng)石油（集團(tuán)）有限責(zé)任公司研究院項(xiàng)目“特低滲透油藏提高采收率配套技術(shù)先導(dǎo)試驗(yàn)及推廣應(yīng)用”（編號(hào)：ycsy2010xjs-A-02）資助。

王宏，1982年生，男，碩士，主要從事石油開(kāi)發(fā)地質(zhì)、油藏描述、油藏?cái)?shù)值模擬與提高采收率的研究工作。地址：（710075）陜西省西安市高新區(qū)科技二路75號(hào)陜西延長(zhǎng)石油（集團(tuán)）有限責(zé)任公司研究院。E-mail：raskay@163.com

涂曉燕）