本發(fā)明涉及油氣勘探技術(shù)領(lǐng)域，尤其涉及一種估算開發(fā)井泄油面積的方法。

背景技術(shù)：

泄油面積是指每口采油井供油的面積，泄油半徑(供油半徑)則是指每口采油井供油面積相等的圓的半徑，有時也可以簡單地將兩口采油井之間的半徑當(dāng)作供油半徑?，F(xiàn)有技術(shù)中，計(jì)算開發(fā)區(qū)塊開發(fā)井泄油半徑的方法通常有兩種：一是根據(jù)壓力變化，利用流體力學(xué)原理計(jì)算壓力波及范圍，進(jìn)而估算開發(fā)井的泄油半徑；二是根據(jù)開發(fā)區(qū)塊面積和開發(fā)井的產(chǎn)量關(guān)系來計(jì)算泄油面積：就比較測試井產(chǎn)量與其鄰井產(chǎn)量而言，開發(fā)井泄油半徑的大小是測試井產(chǎn)量的函數(shù)，在油氣藏區(qū)域內(nèi)，假定地層的有效厚度為常數(shù)，那么在穩(wěn)定流動條件下，每口井的泄油面積與其產(chǎn)量成正比。因此可以根據(jù)開發(fā)井產(chǎn)量與時間的關(guān)系在典型井的產(chǎn)量遞減圖版上讀取相應(yīng)的參數(shù)值，以此計(jì)算單井泄油面積(m.j.fetkovich,etal.decline-curveanalysisusingtypecurves-casehistory，speformationevaluation,december1987)。然而在現(xiàn)有技術(shù)中，在確定單井泄油面積時，直接利用原始的生產(chǎn)數(shù)據(jù)在產(chǎn)量遞減圖版上投點(diǎn)，不能有效消除由非地質(zhì)因素引起的產(chǎn)量變化所帶來的影響，會直接影響參數(shù)讀取的準(zhǔn)確性。

技術(shù)實(shí)現(xiàn)要素：

為了解決上述技術(shù)問題，本發(fā)明提供了一種新的估算開發(fā)井泄油面積的方法。其包括以下步驟：

s101，對開發(fā)井的生產(chǎn)數(shù)據(jù)進(jìn)行預(yù)處理，所述預(yù)處理包括剔除異常數(shù)據(jù)；

s102，根據(jù)預(yù)處理后的生產(chǎn)數(shù)據(jù)建立時間-產(chǎn)量交匯圖，對其進(jìn)行模擬匹配以獲得產(chǎn)量遞減曲線；

s103，將模擬的產(chǎn)量遞減曲線投放到指定的油氣藏生產(chǎn)圖版中；

s104，通過比對模擬的產(chǎn)量遞減曲線在生產(chǎn)圖版中的投射位置來獲取計(jì)算該 開發(fā)井泄油半徑的相關(guān)參數(shù)，進(jìn)而確定該開發(fā)井的泄油半徑；

s105，根據(jù)泄油面積計(jì)算公式計(jì)算該開發(fā)井的泄油面積。

根據(jù)本發(fā)明的實(shí)施例，上述步驟s101中，所述異常數(shù)據(jù)主要指存在人為干預(yù)因素的生產(chǎn)數(shù)據(jù)。

根據(jù)本發(fā)明的實(shí)施例，上述步驟s102中，優(yōu)選采用蒙特卡洛模擬算法進(jìn)行模擬匹配以獲得產(chǎn)量遞減曲線。

根據(jù)本發(fā)明的實(shí)施例，上述步驟s103中，指定的油氣藏生產(chǎn)圖版是由一系列致密油氣藏生產(chǎn)數(shù)據(jù)所歸納出的生產(chǎn)圖版。

根據(jù)本發(fā)明的實(shí)施例，上述生產(chǎn)圖版為fetkovich圖版，其橫坐標(biāo)為相對生產(chǎn)時間對數(shù)坐標(biāo)，縱坐標(biāo)為產(chǎn)量對數(shù)坐標(biāo)，圖版中的曲線為不同的儲層半長與裂縫半長的比值所對應(yīng)的生產(chǎn)曲線。

根據(jù)本發(fā)明的實(shí)施例，上述步驟s103中，具體投放過程如下：1)通過模擬的產(chǎn)量遞減曲線獲得模擬的產(chǎn)量數(shù)據(jù)；2)以模擬的產(chǎn)量數(shù)據(jù)為樣本數(shù)據(jù)繪制到指定的油氣藏生產(chǎn)圖版中，形成一條生產(chǎn)曲線。

根據(jù)本發(fā)明的實(shí)施例，上述步驟104中，計(jì)算開發(fā)井泄油半徑的相關(guān)參數(shù)為儲層半長與裂縫半長的比值，其中儲層半長等于泄油半徑。

根據(jù)本發(fā)明的實(shí)施例，上述步驟104中，需要先根據(jù)開發(fā)井的壓裂作業(yè)獲得裂縫半長，然后才能確定該開發(fā)井的儲層半長。

根據(jù)本發(fā)明的實(shí)施例，上述步驟105中，對于直井式的開發(fā)井，其泄油面積計(jì)算公式為a＝π*(xe-rw)²，對于水平井式的開發(fā)井，其泄油面積計(jì)算公式為a＝π*xe²+2*(xe-rw)*l，式中，a為泄油面積，xe為儲層半長，xf為裂縫半長，rw為井筒半徑，l為水平井段長度，π為圓周率。

與現(xiàn)有技術(shù)相比，本發(fā)明的一個或多個實(shí)施例可以具有如下優(yōu)點(diǎn)：

本發(fā)明根據(jù)開發(fā)井開發(fā)的實(shí)際數(shù)據(jù)，擬合開發(fā)井的產(chǎn)量遞減模式，模擬出新的井生產(chǎn)曲線，并利用fetkovich產(chǎn)量-時間圖版根據(jù)模擬的單井產(chǎn)量數(shù)據(jù)獲取相應(yīng)的參數(shù)，用來計(jì)算開發(fā)井的泄油半徑。該方法快捷便利，能夠在生產(chǎn)初期直觀地通過對產(chǎn)量遞減曲線和已有圖版進(jìn)行比較來得到計(jì)算泄油半徑的相關(guān)參數(shù)，從而計(jì)算開發(fā)井泄油半徑。

本發(fā)明的其它特征和優(yōu)點(diǎn)將在隨后的說明書中闡述，并且，部分地從說明書中變得顯而易見，或者通過實(shí)施本發(fā)明而了解。本發(fā)明的目的和其他優(yōu)點(diǎn)可通過 在說明書、權(quán)利要求書以及附圖中所特別指出的結(jié)構(gòu)來實(shí)現(xiàn)和獲得。

附圖說明

附圖用來提供對本發(fā)明的進(jìn)一步理解，并且構(gòu)成說明書的一部分，與本發(fā)明的實(shí)施例共同用于解釋本發(fā)明，并不構(gòu)成對本發(fā)明的限制。在附圖中：

圖1是本發(fā)明提供的估算開發(fā)井泄油面積的方法流程圖；

圖2示例性地顯示了根據(jù)開發(fā)井實(shí)際生產(chǎn)數(shù)據(jù)模擬的產(chǎn)量遞減模型曲線；

圖3示例性地顯示了根據(jù)產(chǎn)量遞減曲線在圖版中的投射位置確定相關(guān)參數(shù)的示意圖；

圖4是本發(fā)明實(shí)施例中w10井的時間-產(chǎn)量交匯圖；

圖5是本發(fā)明實(shí)施例中w10井的模擬的產(chǎn)量遞減曲線的示意圖；

圖6是本發(fā)明實(shí)施例中w10井的產(chǎn)量遞減曲線投放到圖版中的示意圖。

具體實(shí)施方式

為使本發(fā)明的目的、技術(shù)方案和優(yōu)點(diǎn)更加清楚，以下結(jié)合附圖對本發(fā)明作進(jìn)一步地詳細(xì)說明。本發(fā)明的主要核心思想是根據(jù)低滲透開發(fā)區(qū)塊開發(fā)井的實(shí)際生產(chǎn)數(shù)據(jù)分析產(chǎn)量遞減模式，并由此建立生產(chǎn)曲線，然后在此基礎(chǔ)上根據(jù)fetkovich典型生產(chǎn)曲線圖版讀取與儲層半長(也即泄油半徑)相關(guān)的參數(shù)，并根據(jù)相關(guān)公式估算泄油面積。

如1所示，本發(fā)明提供的估算開發(fā)井泄油面積的方法可以包括五大步驟。

步驟s101：對開發(fā)井的生產(chǎn)數(shù)據(jù)進(jìn)行整理。

一般情況下，開發(fā)井的生產(chǎn)數(shù)據(jù)具有一定的規(guī)律，隨著地層壓力以及地下其它地質(zhì)條件的變化，自然產(chǎn)量隨著生產(chǎn)時間的推移呈逐漸下降的趨勢。但是，由于在生產(chǎn)中還存在人為干預(yù)(比如一些增產(chǎn)或減產(chǎn)措施)的影響，產(chǎn)量的變化可能不是單一的遞減，有時局部產(chǎn)量還可能會有上升，甚至劇烈下降，因此需要先對開發(fā)井的生產(chǎn)數(shù)據(jù)進(jìn)行整理，剔除某些異常數(shù)據(jù)，特別是剔除存在人為干預(yù)因素的生產(chǎn)數(shù)據(jù)。

步驟s102：根據(jù)整理后的生產(chǎn)數(shù)據(jù)建立產(chǎn)量遞減模型。

經(jīng)過整理后的生產(chǎn)數(shù)據(jù)，以相對生產(chǎn)時間(即生產(chǎn)時長)為橫坐標(biāo)、產(chǎn)量為縱坐標(biāo)建立時間-產(chǎn)量交匯圖，并根據(jù)時間-產(chǎn)量的關(guān)系選擇合適的產(chǎn)量遞減模型 進(jìn)行模擬匹配，獲得如圖2所示的產(chǎn)量遞減曲線(也稱生產(chǎn)曲線)。

步驟s103：將模擬的產(chǎn)量遞減曲線投放到指定的油氣藏生產(chǎn)圖版中。

具體地，在步驟2的基礎(chǔ)上，根據(jù)模擬的產(chǎn)量遞減曲線得到模擬的生產(chǎn)數(shù)據(jù)，以此為樣本數(shù)據(jù)將其繪制在指定的油氣藏生產(chǎn)圖版中。在此，指定的油氣藏生產(chǎn)圖版通常是由一系列致密油氣藏生產(chǎn)數(shù)據(jù)所歸納出的生產(chǎn)圖版，例如圖3所示的fetkovich圖版。

fetkovich圖版是fetkovich根據(jù)一系列致密油氣藏生產(chǎn)數(shù)據(jù)歸納總結(jié)的典型生產(chǎn)圖版，該圖版在低滲透油氣藏開發(fā)中具有廣泛的適用性。圖版中的橫坐標(biāo)為相對生產(chǎn)時間(即生產(chǎn)時長)，縱坐標(biāo)為產(chǎn)量，均以10為底的對數(shù)坐標(biāo)繪制。圖版中的曲線表示不同的儲層半長xe與裂縫半長xf的比值所對應(yīng)的生產(chǎn)曲線(也稱為xe/xf線)，這些生產(chǎn)曲線之間等分。圖3中以實(shí)線顯示的曲線是已有的根據(jù)一系列致密油氣藏生產(chǎn)數(shù)據(jù)歸納總結(jié)的生產(chǎn)曲線；以點(diǎn)劃線顯示的曲線則是由開發(fā)井模擬的產(chǎn)量數(shù)據(jù)在該圖版中的投影所形成的一條生產(chǎn)曲線，也即開發(fā)井模擬的產(chǎn)量遞減曲線在圖版中的投影。根據(jù)開發(fā)井模擬的產(chǎn)量遞減曲線在圖版中的投射位置就能夠確定開發(fā)井模擬的產(chǎn)量遞減曲線所對應(yīng)的xe/xf的比值。

步驟4：通過比對模擬的產(chǎn)量遞減曲線在生產(chǎn)圖版中的投射位置來獲取計(jì)算該開發(fā)井泄油半徑的相關(guān)參數(shù)，進(jìn)而確定該開發(fā)井的泄油半徑。

如上所述，fetkovich圖版中已經(jīng)存在一些根據(jù)一系列致密油氣藏生產(chǎn)數(shù)據(jù)歸納總結(jié)的生產(chǎn)曲線，以及這些生產(chǎn)曲線各自所具有的xe/xf的比值。因此，根據(jù)開發(fā)井模擬的產(chǎn)量遞減曲線在fetkovich圖版中的投射位置能夠確定該產(chǎn)量遞減曲線所對應(yīng)的xe/xf的大小。然后，進(jìn)一步地結(jié)合開發(fā)井的裂縫半長xf(xf參數(shù)可以根據(jù)開發(fā)井的壓裂作業(yè)獲得)就可以計(jì)算出開發(fā)井的儲層半長xe，也即泄油半徑。

s105，根據(jù)泄油面積計(jì)算公式計(jì)算該開發(fā)井的泄油面積。

對于直井式的開發(fā)井，其泄油面積計(jì)算公式為a＝π*(xe-rw)²；

對于水平井式的開發(fā)井，其泄油面積計(jì)算公式為a＝π*xe²+2*(xe-rw)*l；

式中，a為泄油面積，xe為儲層半長，xf為裂縫半長，rw為井筒半徑，l為水平井段長度，π為圓周率。

下面以某頁巖氣生產(chǎn)區(qū)塊為例進(jìn)行說明。該區(qū)塊含氣頁巖厚度平均為93.8m，平均孔隙度為3.5％，平均滲透率為6.35md。該區(qū)塊內(nèi)現(xiàn)有開發(fā)井10口，其中生 產(chǎn)時間最長的為w2井，累計(jì)生產(chǎn)136天，生產(chǎn)時間最短為w5井，累計(jì)生產(chǎn)85天，累計(jì)產(chǎn)量最大為571901.69m3，最小為40767.84m3。其中w10井為水平井，該井投產(chǎn)126天；水平段為裸眼完井，井徑118mm，水平井段長300m，該井壓裂后裂縫半長xf為23.3m。

表1某頁巖氣生產(chǎn)區(qū)塊開發(fā)井產(chǎn)量數(shù)據(jù)

以上述方法計(jì)算w10井的泄油面積。

首先根據(jù)步驟s101，首先對w10井的產(chǎn)量數(shù)據(jù)進(jìn)行預(yù)處理，消除由人為因素所引起的異常產(chǎn)量值，繪制產(chǎn)量-時間交匯圖(如圖4所示)。

然后對圖4中的w10井的生產(chǎn)數(shù)據(jù)進(jìn)行產(chǎn)量遞減模擬。例如，可以通過蒙特卡洛模擬算法對圖4中的生產(chǎn)數(shù)據(jù)進(jìn)行擬合，獲得如圖5所示的產(chǎn)量遞減曲線。

然后根據(jù)模擬的產(chǎn)量遞減曲線得到模擬的產(chǎn)量數(shù)據(jù)，以它們?yōu)闃颖緮?shù)據(jù)投放到如圖6所示的fetkovich圖版中，也即將產(chǎn)量遞減曲線投放到帶有xe/xf線的fetkovich圖版中。接著，比對產(chǎn)量遞減曲線與fetkovich圖版中的xe/xf線的位置，以確定該產(chǎn)量遞減曲線所對應(yīng)的xe/xf的值。從圖6的圖版中可以讀出，w10井的模擬的產(chǎn)量數(shù)據(jù)所對應(yīng)的xe/xf的值約為5.2，因此該井泄油半徑為xe＝xf×5.2＝23.3×5.2＝121.16m。最后，該開發(fā)井為水平井，因此泄油面積a＝π*xe²+2*(xe-rw)*l＝3.14×121.16²+2×(121.16-0.118/2)×300＝118755m²。

以上所述，僅為本發(fā)明的具體實(shí)施案例，本發(fā)明的保護(hù)范圍并不局限于此，任何熟悉本技術(shù)的技術(shù)人員在本發(fā)明所述的技術(shù)規(guī)范內(nèi)，對本發(fā)明的修改或替換，都應(yīng)在本發(fā)明的保護(hù)范圍之內(nèi)。