生產井開發指標的分層分方向定量計算方法

張繼成，任帥*，林立，趙帥，馮詩淼

（1.東北石油大學 石油工程學院，黑龍江 大慶 163318；2.大慶油田有限責任公司第三采油廠，黑龍江 大慶 163000）

生產井開發指標的分層分方向定量計算方法

張繼成1，任帥1*，林立2，趙帥1，馮詩淼2

（1.東北石油大學 石油工程學院，黑龍江 大慶 163318；2.大慶油田有限責任公司第三采油廠，黑龍江 大慶 163000）

針對注水開發的多層砂巖油藏分層動態分析難度大等問題，在常規井層開發指標計算基礎上，結合動、靜態劈分方法，綜合考慮滲透率、孔隙度、地層系數、含水飽和度、位置系數、措施系數及注水量系數，提出了一種既可將油、水井作為統一整體，又可對小層、方向流動分量開發指標進行定量計算的體現滲流力學本質的方法。用大慶油田N2-O1井組的產液剖面資料進行驗證。結果表明，所提方法的計算結果與測量結果吻合度較高，精度平均值達75.11%。用該方法計算開發指標，適用性強，能較真實地反映各小層、各方向的產液情況，對現場應用具有指導意義。

注采關系；開發效果；數學模型；分層注水；流動分量

為更有效地制訂分層調整措施和挖潛策略，提高油田的開發能力，須充分掌握油水井分層動用狀況，目前常見的計算井層開發指標的方法有地層系數法［1-4］、數值模擬法［5-7］、動態劈分方程法［8-10］、數學統計與回歸公式法［11-13］、物理實驗法［14-16］等。地層系數法，劈分方法簡單，參數單一，但因未考慮注采連通情況，只適合射孔層數少、儲層非均質性強的砂巖油藏；數值模擬法，需要大量開發動態數據進行擬合，適合開發時間長的油田；動態劈分方程法，是基于滲流力學進行劈分計算的一種方法，在計算過程中因涉及在實際油田開發中難以確定的多個參數，實用性不高；數學統計與回歸公式法，通常針對單個油田開發過程中的動靜態參數建立經驗公式，對其他油田不一定適用，不利于推廣應用；物理實驗法，成本高，且只能計算某時刻的開發指標，無法滿足實際需要。此外，以上方法只能計算至小層開發指標，無法定量計算油水井連通方向的開發指標。

本文在以往研究的基礎上，綜合考慮滲透率、孔隙度、地層系數、含水飽和度、位置系數、措施系數及注水量系數等，提出了一種可定量計算小層、方向流動分量開發指標的方法。

1 采油井流動分量系數

1.1 位置系數

先根據王治國等［17］提出的速率關聯度分析法確定井組的連通情況，再通過注采井之間的液量關聯度定量描述井間的連通情況，進而計算位置系數。某井某層某方向的位置系數

其中，n為周圍連通水井數，θ為某方向方位角。

如圖1所示，在以油井O為中心，水井A，B，C為方向的注采井連通關系中，求水井A方向的位置系數，需先計算∠1和∠2。

圖1 角度計算示意圖Fig.1 Schematic diagram of angle calculation

假設水井A坐標為（x1，y1），水井B坐標為（x2，y2），水井C坐標為（x3，y3），油井O坐標為（x0，y0），則AB連線長為

同理，可求得|BO|，|AO|，|CO|和|AC|，所以有

其中，M為AB的中點。已知|AM|，|AO|和∠1，可得

所以有

同理，可求得∠NOA，最終可得水井A方向的波及角為。

設中心油井連通的水井數為n，則水井A方向的位置系數為

1.2 平均地層系數

平均地層系數=（水井地層系數+油井地層系數）/2，計算式為

其中，Ki為水井i的有效滲透率（μm2）；Hi為水井i的砂巖厚度（m）；Kk為油井k的有效滲透率（μm2）；Hk為油井k的砂巖厚度（m）；為水井i與油井k連通方向的平均地層系數（μm2·m）。

1.3 注水量系數

中心油井在各小層的產油量與產水量均可劈分到各方向，而各方向的產油量和產水量又受連通方向水井注水量的影響。因此，在進行各方向采出量劈分計算時，需引入注水量系數：

其中，Wik為水井i到油井k在j小層劈分的方向注水量（m3）。

1.4 連通狀況系數

在以油井為中心、水井為方向的連通關系中，連通狀況系數是衡量連通方向水井與中心油井砂體連通關系的一個指標，反映儲層的宏觀非均質性［18］，值越大，均質性越強。滲透率變異系數通常用于表征儲層滲透率的非均質性，值越小，均質性越強，反之，非均質性越強。所以，可用滲透率變異系數的倒數表示連通狀況系數：

其中，Z為水井i與油井k某小層連通方向的連通狀況系數；V為滲透率變異系數；Ki為第i口水井在該小層的滲透率（μm2）；hi為第i口水井在該小層的折算有效厚度（m）；為滲透率平均值（μm2）；l為連通方向小層滲透率變化量。

1.5 措施系數

根據該地區的現場情況，確定連通油井k的措施系數Mk，見表1。

表1 措施系數Table 1 Measure coefficient

1.6 井 距

其中，xi為水井i的橫坐標（m）；xk為油井k的橫坐標（m）；yi為水井i的縱坐標（m）；yk為油井k的縱坐標（m）；為對數化后的井距（m）。

1.7 相 滲

根據含水飽和度和大慶油田相滲關系曲線求相滲。計算相滲前，先計算不同時刻含油飽和度。

（1）含油飽和度。

根據油水井所在最優微相選擇相滲曲線，在相滲曲線上，水相相對滲透率零點所對應的含水飽和度為束縛水飽和度，1減去束縛水飽和度得到該方向的原始含油飽和度。

（2）各方向含油飽和度。

各方向原始地質儲量等于該方向所在區域體積與孔隙度和含油飽和度的乘積：

其中，N為方向域內的地質儲量（t）；H為小層折算厚度（m）；A為某四邊形區域面積（m2）；為加權平均后的孔隙度；為加權平均后的原始含油飽和度；為地面原油密度（g·cm-3）；Boi為原油體積系數。

對于剛開發的油田，其油相相對滲透率Kro=1，水相相對滲透率Krw=0，設水井與油井連通方向第1年的產油量為N，則第1年后，該方向剩余地質儲量為N1=N0-N，所以該方向剩余含油飽和度為

若受效油井數未發生變化，在計算第m年某方向波及區域的流動分量條件值時，該方向所有的油、水相相對滲透率均可取第m-1年剩余含油飽和度為o（m-1）時所對應的油、水相相對滲透率，得到第m年流動分量、剩余地質儲量和剩余含油飽和度。

若受效油井數發生變化，當出現新井F時，水井W的受效油井數增加，如圖2所示。其中WF方向波及區域為四邊形WrFs。四邊形WrFs的出現，導致四邊形WpEq和WqAm減小，即WE和WA方向波及區域減小。

（a）對比圖（b）出現新井圖

圖2 出現新井對比示意圖

Fig.2 New well comparison diagram

其中，三角形Wyq范圍內所有網格剩余含油飽和度與WE方向波及區域WpEr的剩余含油飽和度相同；三角形Wzq范圍內所有網格剩余含油飽和度與WA方向波及區域WsAm的剩余含油飽和度相同；五邊形Fszyr為原先未開發區域，其剩余含油飽和度仍為原始含油飽和度，見圖3。

圖3 出現新井含油飽和度計算示意Fig.3 The new well saturation calculation diagram appears

（3）相滲歸一化處理。

由于油、水井連通方向上存在不同微相，同一含油飽和度對應不同相滲，因此需對油水井連通方向不同微相的油、水相滲進行歸一化處理，使相滲更貼近油藏地質真實情況。

1.8 流動分量系數

（1）單層單方向產水流動分量條件值：

（2）單層單方向產油流動分量條件值：

其中，Kro為油相相對滲透率。

（3）單層單方向產水流動分量系數。

將某小層各方向產水流動分量條件值相加，該小層某方向產水流動分量條件值與該小層產水流動分量條件值之和的比值就是該方向在該小層的產水流動分量系數：

其中，Cw為某方向產水流動分量系數；Fwik為某方向產水流動分量條件值。

（4）單層單方向產油流動分量系數。

將某小層各方向產油流動分量條件值相加，該小層某方向產油流動分量條件值與該小層產油流動分量條件值之和的比值就是該方向在該小層的產油流動分量系數：

其中，Co為某方向產油流動分量系數；Foik為某方向產油流動分量條件值。

（5）單層產水流動分量條件值。

將某小層各方向流動分量條件值相加，得該小層的產水流動分量條件值：

其中，Fwj為某小層產水流動分量條件值。

（6）單層產油流動分量條件值。

將某小層各方向產油流動分量條件值相加，得該小層的產油流動分量條件值：

其中，Foj為某小層產油流動分量條件值。

（7）單層產水流動分量系數。

將某井各小層產水流動分量條件值相加，某井某小層產水流動分量條件值與該井產水流動分量條件值總和的比值就是該小層在該井的產水流動分量系數：

其中，Swj為單層產水流動分量系數。

（8）單層產油流動分量系數。

將某井各小層產油流動分量條件值相加，某井某小層產油流動分量條件值與該井產油流動分量條件值總和的比值就是該小層在該井的產油流動分量系數：

其中，Soj為單層產油流動分量系數。

2 注水井流動分量系數

2.1 流動分量條件值

（1）單層單方向流動分量條件值：

（2）單層流動分量條件值：

其中，Yj為某小層流動分量條件值。

2.2 流動分量系數

（1）單層單方向流動分量系數：

其中，Cik為某方向流動分量系數。

（2）單層流動分量系數：

其中，Sj為單層流動分量系數。

3 井層方向指標

3.1 注采比

其中，Rip為方向注采比，Wik為方向月注水量（m3），Qwik為方向月產水量（m3），Qoik為方向月產油量（m3），Bo為原油體積系數。

3.2 含水率

其中，fw為含水率。

4 應用實例

4.1 井組概況

N2-O1油井于1998年7月投產，共射開19個小層，主要沉積類型為河道砂、主體砂、非主體砂和表外砂，小層平均孔隙度15.1%，滲透率0.112～0.390 μm2，平均滲透率0.161 μm2，截至2021年9月，該井累計產油量14 443 m3，累計產液量221 807 m3，月產液量為955 m3，含水率為95.49%。周圍連通的水井有N2-W1、N2-W2、N2-W3、N2-W4和N2-W5。N2-O1井組注采對應關系如圖4所示。

圖4 N2-O1井組注采關系示意Fig.4 Schematic of injection production relationship of well cluster N2-O1

4.2 注采指標

選取N2-O1井組的SIII3b小層作為研究對象。在SIII3b小層，N2-O1油井連通3個方向水井，注采關系如圖5所示。

圖5 N2-O1井組SIII3b小層注采關系示意Fig.5 Schematic of injection production relationship of SIII3b small layer of N2-O1 well cluster

整理N2-O1油井連通的N2-W4、N2-W3和N2-W2水井的物性參數，并計算SIII3b小層方向的流動分量系數，結果見表2～表4。

表2 N2-O1井組SIII3b小層物性參數Table 2 Physical property parameters of SIII3b sublayer of N2-O1 well cluster

表3 N2-O1油井SIII3b小層受效井各參數計算結果Table 3 Calculation results of parameters of affected wells in SIII3b sublayer of N2-O1 oil well

表4 N2-O1油井SIII3b小層受效井流動分量條件值計算結果Table 4 Calculation results of flow component condition parameters of effective well in SIII3b small layer of N2-O1 oil well

表5 N2-O1油井各小層開發指標計算結果Table 5 Calculation results of development indexes of each small layer of N2-O1 production well

當得到條件值中的各參數后，可計算小層各方向流動分量條件值，進而得到各小層流動分量條件值，最終得到各小層流動分量系數。將N2-O1井組各井資料代入相應計算式，與各井各小層流動分量系數相乘，得到各井各小層月產水量、月產油量、含水率等開發指標，見表5。

4.3 開發狀況評價

4.3.1 結果驗證

整理N2-O1油井各層段產液剖面測量數據，并與各層段流動分量計算結果進行對比，結果見圖6和表6。

圖6 產液剖面測量結果與計算結果對比Fig.6 Comparison between test and calculation results of liquid production profile of N2-O1 oil well

表6 N2-O1油井產液剖面測量值與計算結果對比Table 6 Comparison between test and calculation results of liquid production profile of N2-O1 oil well

可知，各層段產液占比測量結果與計算結果較為接近。層段1（SII5、SII6、SII13）與層段6（PII1a、PII9）的計算精度相對較低，其原因可能為實際測量時取的是單次測量結果，具有偶然性。整體看，各層段的產液占比計算值與測量值變化趨勢基本一致，吻合度較高，平均精度達75.11%。

4.3.2 開發現狀評價

由N2-O1油井各小層流動分量的計算結果和產液剖面實際產液量知，第4個層段SIII4b小層產液量最多，占全井產液量的30%左右，第2個層段SII15～SII16小層產液量最少，占全井產液量的3%左右。SII5小層的含水率最低，為71.28%，其次是SII15小層，含水率為78.58%。綜合分析得到，N2-O1油井的SII5、SII15和SIII6小層潛力較大，應加強注水開發。

5 結論

5.1在以往只考慮小層開發指標計算的基礎上，提出了方向開發指標計算方法，除考慮滲透率、厚度、位置系數、井距、含油飽和度、措施系數外，還加入了注水量系數，將油、水井作為整體考慮。

5.2用N2-O1井組的產液剖面現場資料驗證本文方法，各層段產液占比計算結果與測量結果吻合度較高，二者變化趨勢基本一致，平均精度達75.11%。對N2-O1井組進行綜合評價，認為該井組的SII5、SII15和SIII6小層潛力較大，應加強注水開發。

［1］李源流，楊兆平，潘多壽，等. 考慮相對滲透率變化的砂巖油藏產量劈分方法［J］. 西安石油大學學報（自然科學版），2020，35（3）：72-76， 85. DOI：10.3969/j.issn.1673-064X.2020.03.010

LI Y L， YANG Z P，PAN D S，et al. A Production splitting method of sandstone reservoir considering the change of relative permeability［J］. Journal of Xi′an Shiyou University （Natural Science Edition），2020，35（3）：72-76， 85. DOI：10.3969/j.issn.1673-064X. 2020.03.010

［2］夏靜，馮陽，陳新宇. 油井分層注水倍數計算方法［J］. 數學的實踐與認識，2017，47（3）：114-119.

XIA J， FENG Y，CHEN X Y. Computing method of separating layer water injected volume about oil wells［J］. Mathematics in Practice and Theory， 2017，47（3）：114-119.

［3］邸井東. 結合含水飽和度場分布的水淹級別評價體系［J］. 石油化工高等學校學報，2019，32（3）：51-57. DOI：10.3969/j.issn.1006-396X.2019.03.009

DI J D. Evaluation system of water flooded grade combined with distribution of water saturation field［J］. Journal of Petrochemical Universities，2019，32（3）：51-57. DOI：10.3969/j.issn.1006-396X.2019.03.009

［4］張繼成，王瀟悅. 考慮含水飽和度的產量劈分方法及應用［J］. 浙江大學學報（理學版），2015，42（5）：626-630. DOI：10.3785/j.issn.1008-9497.2015.05.020

ZHANG J C， WANG X Y. The method of oil production splitting based on water saturation and its application［J］. Journal of Zhejiang University （Science Edition），2015，42（5）：626-630. DOI：10. 3785/j.issn.1008-9497.2015.05.020

［5］陳建華，晏慶輝，駱逸婷，等. 基于歷史生產數據的多層合采井產量劈分新方法［J］.中國海上油氣，2022，34（1）：110-116. DOI：10.11935/j.issn.1673-1506.2022.01.013

CHEN J H， YAN Q H，LUO Y T，et al. A historical production data based method for production splitting of multi-layer commingled gas wells［J］. China Offshore Oil and Gas， 2022， 34（1）： 110-116. DOI：10.11935/j.issn.1673-1506.2022.01.013

［6］LIU X，WANG P，LIANG X，et al. Research progress and surface/interfacial regulation methods for electrophotocatalytic hydrogen production from water splitting［J］. Materials Today Energy， 2020，18：100524. DOI：10.1016/j.mtener.2020.100524

［7］林孟雄，成育紅，張林，等. 蘇里格氣田蘇東區塊多層產量劈分新方法［J］. 新疆地質，2019，37（3）：419-421. DOI：10.3969/j.issn.1000-8845.2019.03.024

LIN M X， CHENG Y H，ZHANG L，et al. A new method of production dividing for commingling production in the Sudong area［J］. Xinjiang Geology， 2019，37（3）：419-421. DOI：10.3969/j.issn.1000-8845.2019.03.024

［8］HEBA H F，SAYED Y S，NAGEH K A，et al. Emerging nanoporous anodized stainless steel for hydrogen production from solar water splitting［J］. Journal of Cleaner Production， 2020，274：122826. DOI：10.1016/j.jclepro.2020.122826

［9］周立國. 低滲砂巖油藏優勢注采方向量化表征技術及應用［J］. 特種油氣藏，2021，28（6）：98-104. DOI：10.3969/j.issn.1006-6535.2021.06.013

ZHOU L G. Quantitative characterization technology and application of dominant injection-production direction in low-permeability sandstone reservoirs［J］. Special Oil and Gas Reservoirs， 2021，28（6）：98-104. DOI：10.3969/j.issn.1006-6535.2021.06.013

［10］RODOLPHE D，ROBERT L P，ALANA G，et al. Reinventing residual reserves in the sea： Are we favouring ease of establishment over need for protection？［J］. Aquatic Conservation： Marine and Freshwater Ecosystems，2015，25（4）：480-504. DOI：10.1002/aqc.2445

［11］楊兆平，岳世俊，鄭長龍，等. 薄互層砂巖油藏多因素綜合約束的產量劈分方法［J］. 巖性油氣藏，2018，30（6）：117-124. DOI：10.12108/yxyqc.20180614

YANG Z P， YUE S J，ZHENG C L，et al. Production split method restricted synthetically by multi-factors in thin interbed sandstone reservoirs［J］. Lithologic Reservoirs， 2018，30（6）：117-124. DOI：10.12108/yxyqc.20180614

［12］王立，喻高明，傅宣豪，等. 基于反距離加權插值法的產量劈分新方法［J］. 斷塊油氣田，2018，25（5）：617-621. DOI：10.6056/dkyqt201805015

WANG L， YU G M，FU X H，et al. New method for production cleavage by inverse distance weighted interpolation［J］. Fault-block Oil and Gas Field， 2018，25（5）：617-621. DOI：10.6056/dkyqt201805015

［13］張繼成，單五一，盧光夫，等. 一種預測注采關系調整效果的數學模型［J］. 數學的實踐與認識，2020，50（4）：162-170.

ZHANG J C， SHAN W Y，LU G F，et al. A mathematical model to predict the effects of adjustment between injection and production［J］. Mathematics in Practice And Theory， 2020，50（4）：162-170.

［14］趙輝，劉鄧，江厚順，等. 基于數據物理驅動的油藏調剖動態預測新方法［J］. 長江大學學報（自然科學版），2021，18（1）：60-69. DOI：10.3969/j.issn.1673-1409.2021.01.007

ZHAO H， LIU D，JIANG H S，et al. A new method of reservoir profile control dynamic prediction based on the data-driven model［J］. Journal of Yangtze University （Natural Science Edition），2021，18（1）：60-69. DOI：10.3969/j.issn.1673-1409.2021.01.007

［15］HARIMURTI D. Indonesian milestone in production-sharing contract in perspective of government take， contractor take， cost recovery and production target［C］// Spe/Iatmi Asia Pacific Oil amp; Gas Conference amp; Exhibition. Jakarta：Silverchair，2017. DOI：10.2118/187008-MS

［16］謝書劍，周玉輝，丁艷雪，等. 水驅油藏分層注水效率評價新方法［J］. 中國科技論文，2021，16（5）：475-481. DOI：10.3969/j.issn.2095-2783.2021.05.004

XIE S J， ZHOU Y H，DING Y X，et al. A new method to evaluate the injection efficiency of different layers in waterflooding reservoir［J］. China Sciencepaper， 2021，16（5）：475-481. DOI：10. 3969/j.issn.2095-2783.2021.05.004

［17］王治國，張繼成，宋考平，等. 應用速率關聯度分析井間連通情況［J］. 特種油氣藏，2011，18（3）：91-93， 140. DOI：10.3969/j.issn.1006-6535.2011.03.024

WANG Z G， ZHANG J C，SONG K P，et al. Analyze interwell communication by rate correlation［J］. Special Oil and Gas Reservoirs， 2011，18（3）：91-93， 140. DOI：10.3969/j.issn.1006-6535.2011.03.024

［18］婁瑞，劉宗堡，張云峰，等. 齊家地區高Ⅲ、高Ⅳ組砂巖輸導層靜態連通評價［J］. 斷塊油氣田，2019，26（3）：285-289. DOI：10.6056/dkyqt201903003

LOU R， LIU Z B，ZHANG Y F，et al. Static connectivity evaluation of sandstone carrier bed in GⅢ-GⅣ groups in Qijia area［J］. Fault Block Oil and Gas Field， 2019，26 （3）：285-289. DOI：10.6056/dkyqt201903003

Quantitative calculation method of development indexes for layered and directional of production wells

ZHANG Jicheng1, REN Shuai1, LIN Li2, ZHAO Shuai1, FENG Shimiao2

（1.School of Petroleum Engineering，Northeast Petroleum University，Daqing163318，Heilongjiang Province，China；2. No.3Oil Production Plant，Daqing Oilfield Co.，Ltd.，Daqing163000，Heilongjiang Province，China）

On account of the difficulty of Layered Dynamic Analysis of multi-layer sandstone reservoir developed by water flooding, based on the calculation approach of conventional well layer development indexes, combined with dynamic and static splitting methods, and comprehensively considering the parameters such as permeability, porosity, formation coefficient, water saturation, position coefficient, measure coefficient and water injection coefficient, this paper propose a new method which can not only deal with oil and water wells as a unified whole, but also for small layers. The mathematical method for quantitative calculation of development index of directional flow component reflects the essence of seepage mechanics. The proposed method is verified by the liquid production profile data of N2-O1 well group in Daqing Oilfield, and shows high coincidence with the actual production, with average accuracy 75.11%. It has strong applicability to the calculation of development indicators and plays an important guiding role in field application.

injection-production relationship; development effect; mathematical model; layered water injection; flow component

P 548

A

1008?9497（2022）05?598?08

2021?11?03.

黑龍江省自然科學基金資助項目（E201407）.

張繼成（1972—），ORCID：https：//orcid.org/0000-0002-6791-7767，男，博士，教授，主要從事石油與天然氣工程研究，E-mail：zhangjc777@163.com.

通信作者，ORCID：https：//orcid.org/ 0000-0002-5274-1674，E-mail：769277251@qq.com.

10.3785/j.issn.1008-9497.2022.05.011