一種低滲氣藏產(chǎn)能計算方法
【專利摘要】本發(fā)明公開了一種低滲氣藏產(chǎn)能計算方法，該方法包括：將油藏對象投影到二維平面上，獲得油藏對象的二維油藏模型；將二維油藏模型中的所有油藏對象表示成有序點集；對有序點集進行德勞內(nèi)Delaunay三角剖分，獲得Delaunay三角形網(wǎng)；基于Delaunay三角形網(wǎng)，生成Voronoi網(wǎng)格圖；根據(jù)Voronoi網(wǎng)格圖計算低滲氣藏產(chǎn)能。通過本發(fā)明的方案，能夠精細模擬壓裂井的產(chǎn)能變化，計算出地層達到真實產(chǎn)能的時間和大小，節(jié)省了測試時間和費用。
【專利說明】
一種低滲氣藏產(chǎn)能計算方法

【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及低滲氣井的壓裂測試技術(shù)，尤其涉及一種低滲氣藏產(chǎn)能計算方法。

【背景技術(shù)】
[0002]目前的業(yè)務發(fā)展需要油藏能夠準確的預測和評價壓裂產(chǎn)能，正確的指導壓裂工藝，減小無效壓裂的次數(shù)，提高壓裂服務的水平。
[0003]但目前還沒有一種有效的方法和軟件能夠準確的描述裂縫形態(tài)，同時精確地預測壓裂產(chǎn)能。主要的研宄都是在常規(guī)的裂縫形態(tài)描述方法基礎(chǔ)上進行改進和完善。近期大量的實踐表明，這些現(xiàn)有技術(shù)存在網(wǎng)格極差變化大，網(wǎng)格模型在進行數(shù)值模擬計算時容易產(chǎn)生不收斂的問題。


【發(fā)明內(nèi)容】

[0004]為了解決上述問題，本發(fā)明提出了一種低滲氣藏產(chǎn)能計算方法，能夠精細模擬壓裂井的產(chǎn)能變化，計算出地層達到真實產(chǎn)能的時間和大小，節(jié)省了測試時間和費用。
[0005]為了達到上述目的，本發(fā)明提出了一種低滲氣藏產(chǎn)能計算方法，該方法包括:
[0006]將油藏對象投影到二維平面上，獲得油藏對象的二維油藏模型。
[0007]將二維油藏模型中的所有油藏對象表示成有序點集。
[0008]對有序點集進行德勞內(nèi)Delaunay三角剖分，獲得Delaunay三角形網(wǎng)。
[0009]基于Delaunay三角形網(wǎng)，生成Voronoi網(wǎng)格圖。
[0010]根據(jù)Voronoi網(wǎng)格圖計算低滲氣藏產(chǎn)能。
[0011]優(yōu)選地，油藏對象包括:油藏邊界、層面、斷層、井和水力壓裂裂縫。
[0012]優(yōu)選地，將油藏對象投影到二維平面上，獲得油藏對象的二維油藏模型是指:
[0013]將油藏對象投影到二維平面上，邊界轉(zhuǎn)化為包絡(luò)線，斷層轉(zhuǎn)化為多邊形，井轉(zhuǎn)化為點，斜度大于預定斜度閾值的井或者水平井轉(zhuǎn)化為多邊形，水力壓裂裂縫轉(zhuǎn)化為線段。
[0014]優(yōu)選地，將二維油藏模型中的所有油藏對象表示成有序點集是指:將二維油藏模型中的包絡(luò)線、多邊形、線段以及點轉(zhuǎn)化為二維平面上不同的有序點集。
[0015]優(yōu)選地，該方法還包括:當油藏對象非均勻分布和/或油藏對象的有序點集中的散點個數(shù)超過預定的散點閾值時，增加二維油藏模型中的油藏對象的限定條件，將油藏對象的有序點集進彳丁細分。
[0016]其中，增加二維油藏模型中的油藏對象的限定條件，將油藏對象的有序點集進行細分包括:增加每個和/或個別油藏對象的局部邊界區(qū)域，在局部邊界區(qū)域內(nèi)部對該油藏對象的有序點集進行細分。
[0017]優(yōu)選地，對有序點集進行德勞內(nèi)DeIaunay三角剖分，獲得DeIaunay三角形網(wǎng)的方法包括以下步驟:
[0018]61、遍歷所述有序點集的所有散點，得出有序點集的包容盒，獲得作為有序點集的凸殼的初始三角形并放入Delaunay三角形鏈表。
[0019]62、將有序點集中的一個散點插入Delaunay三角形鏈表，在Delaunay三角形鏈表中找出其外接圓包含插入的散點的三角形，刪除影響三角形的公共邊，將插入的散點與影響三角形的全部頂點連接起來，完成一個散點在Delaunay三角形鏈表中的插入。
[0020]63、根據(jù)Delaunay三角形優(yōu)化準則對局部新形成的三角形進行優(yōu)化，將優(yōu)化后的三角形放入Delaunay三角形鏈表。
[0021]64、判斷有序點集中是否還有未插入的散點，如果有則進行步驟62 ;如果沒有則結(jié)束。
[0022]優(yōu)選地，當引入地面的地性線和特征線時，對有序點集進行德勞內(nèi)Delaunay三角剖分，獲得Delaunay三角形網(wǎng)的方法包括以下步驟:
[0023]71、根據(jù)已有的地性線和特征線，形成控制邊鏈表。
[0024]72、以控制邊鏈表中的一個線段為基邊，從有序點集中找出同基邊的兩端點的距離和最小的點，以該點為頂點，以基邊為邊，向外擴展一個三角形并放入三角形鏈表中。
[0025]73、判斷控制邊鏈表中是否還有未擴展的線段，如果有則進行步驟72 ;如果沒有則進行步驟74。
[0026]74、依次將新形成的三角形的邊作為基邊，形成新的控制邊鏈表，按照步驟72，對新的控制邊鏈表中的所有的線段進行擴展。
[0027]75、判斷新的控制邊鏈表中是否還有未擴展的新形成的三角形，如果有則進行步驟74;如果沒有則結(jié)束。
[0028]優(yōu)選地，該方法還包括:當油藏對象非均勻分布和/或油藏對象的有序點集中的散點個數(shù)超過預定的散點閾值時，在對點集進行Delaunay三角剖分時，如果翻邊后獲得的三角形仍然不能滿足Delaunay三角形的基本準則，則在該點集內(nèi)部增加新的點作為輔助點，或增加新的點集作為輔助點集繼續(xù)進行Delaunay三角剖分。
[0029]其中，在該點集內(nèi)部增加新的點作為輔助點，或增加新的點集作為輔助點集采用隨機增量構(gòu)造法，隨機增量構(gòu)造法包括以下步驟:
[0030]構(gòu)造相對于原始點集P大于一定閾值的點集P-1PjPy使用該點集把原始點集P包圍起來，并以隨機的次序依次把原始點集P中的各個點插入到點集P-A2P-的三角剖分中，在插入的過程中始終保證得到的三角剖分是Delaunay三角剖分，在原始點集P的所有的點都插入以后，獲得{P-!, P_2，P_3} U P的Delaunay三角剖分，將P^1, P_2，P_3以及與其相關(guān)聯(lián)的各邊刪除，獲得原始點集P的Delaunay三角剖分。
[0031]優(yōu)選地，
[0032]Voronoi網(wǎng)格圖是指:由一組連接兩相鄰點的直線的垂直平分線組成的連續(xù)多邊形所組成的網(wǎng)格；
[0033]基于Delaunay三角形網(wǎng)，生成Voronoi網(wǎng)格圖是指:基于Delaunay三角形網(wǎng)中的Delaunay三角形，獲得Delaunay三角形中的每條邊的垂直平分線，連接Delaunay三角形網(wǎng)中的全部垂直平分線，并使每個連接點與Delaunay三角形的外接圓的圓心重合，連接全部連接點，獲得Voronoi網(wǎng)格圖。
[0034]優(yōu)選地，
[0035]二維油藏模型為二維PEBI網(wǎng)格；Voronoi網(wǎng)格是三維PEBI網(wǎng)格。
[0036]與現(xiàn)有技術(shù)相比，本發(fā)明包括:將油藏對象投影到二維平面上，獲得油藏對象的二維油藏模型；將二維油藏模型中的所有油藏對象表示成有序點集；對有序點集進行德勞內(nèi)Delaunay三角剖分，獲得Delaunay三角形網(wǎng)；基于Delaunay三角形網(wǎng)，生成Voronoi網(wǎng)格圖；根據(jù)Voronoi網(wǎng)格圖計算低滲氣藏產(chǎn)能。通過本發(fā)明的方案，能夠精細模擬壓裂井的產(chǎn)能變化，計算出地層達到真實產(chǎn)能的時間和大小，節(jié)省了測試時間和費用。

【專利附圖】

【附圖說明】
[0037]下面對本發(fā)明實施例中的附圖進行說明，實施例中的附圖是用于對本發(fā)明的進一步理解，與說明書一起用于解釋本發(fā)明，并不構(gòu)成對本發(fā)明保護范圍的限制。
[0038]圖1為本發(fā)明的低滲氣藏產(chǎn)能計算方法流程圖；
[0039]圖2為本發(fā)明的油藏對象的二維油藏模型示意圖；
[0040]圖3為本發(fā)明的增加二維油藏模型中的油藏對象的限定條件示意圖；
[0041]圖4為常規(guī)的三角剖分技術(shù)不意圖；
[0042]圖5為常規(guī)的Delaunay三角剖分示意圖；
[0043]圖6為常規(guī)的Delaunay三角剖分中非法邊的翻轉(zhuǎn)示意圖；
[0044]圖7為常規(guī)的互為對偶的Voronoi圖與Delaunay圖示意圖。

【具體實施方式】
[0045]為了便于本領(lǐng)域技術(shù)人員的理解，下面結(jié)合附圖對本發(fā)明作進一步的描述，并不能用來限制本發(fā)明的保護范圍。
[0046]目前的業(yè)務發(fā)展需要油藏能夠準確的預測和評價壓裂產(chǎn)能，正確的指導壓裂工藝，減小無效壓裂的次數(shù)，提高壓裂服務的水平。但現(xiàn)有的壓裂產(chǎn)能計算中存在以下問題:
[0047]1、現(xiàn)有技術(shù)對裂縫網(wǎng)格的描述都是直線或者長條狀的網(wǎng)格，與主體大網(wǎng)格極差較大，不利于計算。
[0048]2、現(xiàn)有技術(shù)在裂縫網(wǎng)格創(chuàng)建中三角化和生成PEBI網(wǎng)格時存在問題。
[0049]3、現(xiàn)有技術(shù)生成的裂縫網(wǎng)格在進行數(shù)值模擬計算時存在計算極易不收斂的問題。
[0050]針對上述問題，本發(fā)明主要提供一種毫米級別的非結(jié)構(gòu)化網(wǎng)格進行低滲氣井的產(chǎn)能計算方法，包括模型的描述和網(wǎng)格生成方法。主要應用于低滲氣井的壓裂精細模擬和產(chǎn)能計算。在壓裂時由于儲層容易受到外來入井液的傷害，導致壓裂測試初期的產(chǎn)能低，此時的產(chǎn)量并不是地層的真實產(chǎn)能，隨著壓裂液的逐漸返出，壓裂井產(chǎn)能才能漸漸恢復到儲層的真實產(chǎn)能。而要達到測試出真實產(chǎn)能的目的則需要延長測試時間。如TWT-A5井壓裂后180天才到達地層真實產(chǎn)能。對于海上探井而言，延長測試時間則意味著需要指大量的鉆井船費用，所以東海前期壓裂探井的測試時間都很短，一般在一周以內(nèi)，導致部分壓裂井測試產(chǎn)能低于預期，影響了事業(yè)部的壓裂服務水平。通過采用該技術(shù)，精細模擬壓裂井的產(chǎn)能變化，計算出地層達到真實產(chǎn)能的時間，以及產(chǎn)能大小。節(jié)省測試時間和費用，提高壓裂服務水平，降低壓裂風險。
[0051]具體地，本發(fā)明提出了一種低滲氣藏產(chǎn)能計算方法，如圖1所示，該方法包括:
[0052]S101、將油藏對象投影到二維平面上，獲得油藏對象的二維油藏模型。
[0053]優(yōu)選地，油藏對象包括:油藏邊界、層面、斷層、井和水力壓裂裂縫。
[0054]優(yōu)選地，將油藏對象投影到二維平面上，獲得油藏對象的二維油藏模型是指:
[0055]將油藏對象投影到二維平面上，邊界轉(zhuǎn)化為包絡(luò)線，斷層轉(zhuǎn)化為多邊形，井轉(zhuǎn)化為點，斜度大于預定斜度閾值的井或者水平井轉(zhuǎn)化為多邊形，水力壓裂裂縫轉(zhuǎn)化為線段。
[0056]油藏模型由油藏邊界、層面、斷層、井這些對象組成，將這些對象投影到二維平面上，邊界可以轉(zhuǎn)化為包絡(luò)線，斷層可以轉(zhuǎn)化為多邊形、井可以轉(zhuǎn)化為點，大斜度井或者水平井可以轉(zhuǎn)化成多邊形(有序點集)，水力壓裂裂縫可以轉(zhuǎn)化為線段。由于多邊形和線段都可以由有序點集組成，因此，歸納起來，油藏模型中的全部對象可以轉(zhuǎn)化為平面上不同的有序點集。如圖2所示。
[0057]S102、將二維油藏模型中的所有油藏對象表示成有序點集。
[0058]優(yōu)選地，將二維油藏模型中的所有油藏對象表示成有序點集是指:將二維油藏模型中的包絡(luò)線、多邊形、線段以及點轉(zhuǎn)化為二維平面上不同的有序點集。
[0059]當油藏對象以上述模型進行描述后，所有的油藏對象在平面上都表示成了有序點集，這樣就為后續(xù)的網(wǎng)格化提供了數(shù)據(jù)。
[0060]優(yōu)選地，該方法還包括:當油藏對象非均勻分布且比較復雜和/或油藏對象的有序點集中的散點個數(shù)超過預定的散點閾值時，增加二維油藏模型中的油藏對象的限定條件，將油藏對象的有序點集進行細分，以便進行局部的網(wǎng)格化處理。如圖3所示。
[0061]其中，增加二維油藏模型中的油藏對象的限定條件，將油藏對象的有序點集進行細分包括:增加每個和/或個別油藏對象的局部邊界區(qū)域，在局部邊界區(qū)域內(nèi)部對該油藏對象的有序點集進行細分。
[0062]S103、對有序點集進行德勞內(nèi)Delaunay三角剖分，獲得Delaunay三角形網(wǎng)。
[0063]散點集的三角剖分，對數(shù)值分析以及圖形學來說，都是極為重要的一項預處理技術(shù)，該技術(shù)的示意圖如圖4所示。
[0064]在實際中運用Delaunay三角剖分，它是一種特殊的三角剖分。它的特點是生成的三角形都是銳角三角形。這樣一般可以保證任意一個三角形的外接圓的圓心在三角形內(nèi)。
[0065]圖5中所示為平面上4個點的兩種三角剖分，顯然左圖的剖分得到兩個狹長的鈍角三角形，而右圖的剖分則接近正三角形。通?？梢詫⒏魅瞧史种兴淖钚A角作為指標來衡量它們的優(yōu)劣，如果最小夾角相等，就比較它們的次小夾角；以此類推，對于任一給定的點集P，存在的三角剖分的總數(shù)必然是有限的，因此必然存在一個使最小夾角達到最大的三角剖分，這個三角剖分就稱為Delaunay三角剖分，圖5中左圖的三角劃分由六個角度，如果所有角度中有超過最大角度(例如90度的夾角)，則刪除PR，然后連接P1Pk再重復判斷，這個操作通常稱為翻邊；翻邊后如此時存在的所有夾角集合眾所有的夾角都滿足要求，則保存該邊，繼續(xù)下一個判斷，直到所有的邊都被遍歷；這個三角剖分就是Delaunay三角劃分，因此可以得到從任意三角剖分到Delaunay三角剖分的算法。
[0066]如圖6所示，當P,落入三角形P PjPk時，只需要把P ,分別與三個頂點連線即可得到三角剖分。這個剖分中出現(xiàn)非法邊，假設(shè)非法的只需對這三條邊遞歸進行邊反轉(zhuǎn)操作即可得到合理的Delaunay三角剖分。
[0067]優(yōu)選地，對有序點集進行德勞內(nèi)DeIaunay三角剖分，獲得DeIaunay三角形網(wǎng)的方法包括以下步驟:
[0068]61、遍歷所述有序點集的所有散點，得出有序點集的包容盒，獲得作為有序點集的凸殼的初始三角形并放入Delaunay三角形鏈表。
[0069]62、將有序點集中的一個散點插入Delaunay三角形鏈表，在Delaunay三角形鏈表中找出其外接圓包含插入的散點的三角形，刪除影響三角形的公共邊，將插入的散點與影響三角形的全部頂點連接起來，完成一個散點在Delaunay三角形鏈表中的插入。
[0070]63、根據(jù)Delaunay三角形優(yōu)化準則對局部新形成的三角形進行優(yōu)化，將優(yōu)化后的三角形放入Delaunay三角形鏈表。
[0071]64、判斷有序點集中是否還有未插入的散點，如果有則進行步驟62 ;如果沒有則結(jié)束。
[0072]上述基于散點的構(gòu)網(wǎng)算法理論嚴密、唯一性好，網(wǎng)格滿足空圓特性，較為理想。由其逐點插入的構(gòu)網(wǎng)過程可知，在完成構(gòu)網(wǎng)后，增加新點時，無需對所有的點進行重新構(gòu)網(wǎng)，只需對新點的影響三角形范圍進行局部聯(lián)網(wǎng)，且局部聯(lián)網(wǎng)的方法簡單易行。同樣，點的刪除、移動也可快速動態(tài)地進行。但在實際應用當中，這種構(gòu)網(wǎng)算法不易引入地面的地性線和特征線，當點集較大時構(gòu)網(wǎng)速度也較慢，如果點集范圍是非凸區(qū)域或者存在內(nèi)環(huán)，則會產(chǎn)生非法三角形。
[0073]為了克服基于散點構(gòu)網(wǎng)算法的上述缺點，特別是為了提高算法效率，可以對網(wǎng)格中三角形的空圓特性稍加放松，亦即采用基于邊的構(gòu)網(wǎng)方法。
[0074]優(yōu)選地，當引入地面的地性線和特征線時，對有序點集進行德勞內(nèi)Delaunay三角剖分，獲得Delaunay三角形網(wǎng)的方法包括以下步驟:
[0075]71、根據(jù)已有的地性線和特征線，形成控制邊鏈表。
[0076]72、以控制邊鏈表中的一個線段為基邊，從有序點集中找出同基邊的兩端點的距離和最小的點，以該點為頂點，以基邊為邊，向外擴展一個三角形并放入三角形鏈表中。
[0077]73、判斷控制邊鏈表中是否還有未擴展的線段，如果有則進行步驟72 ;如果沒有則進行步驟74。
[0078]74、依次將新形成的三角形的邊作為基邊，形成新的控制邊鏈表，按照步驟72，對新的控制邊鏈表中的所有的線段進行擴展。
[0079]75、判斷新的控制邊鏈表中是否還有未擴展的新形成的三角形，如果有則進行步驟74;如果沒有則結(jié)束。
[0080]研宄表明，對η個實數(shù)進行排序的問題，可以在線性時間內(nèi)歸約為Voronoi圖的構(gòu)造問題。
[0081]在實際的應用中，由于情況可能非常復雜，對三角形的要求較高，翻邊后常常也不滿足要求，就需要在該集合內(nèi)部增加新的點作為輔助點，以輔助點生成一個局部點集再繼續(xù)三角剖分。
[0082]優(yōu)選地，該發(fā)明方法還包括:當油藏對象非均勻分布和/或油藏對象的有序點集中的散點個數(shù)超過預定的散點閾值時，在對點集進行Delaunay三角剖分時，如果翻邊后獲得的三角形仍然不能滿足Delaunay三角形的基本準則，則在該點集內(nèi)部增加新的點作為輔助點，或增加新的點集作為輔助點集繼續(xù)進行Delaunay三角剖分。
[0083]其中，在該點集內(nèi)部增加新的點作為輔助點，或增加新的點集作為輔助點集采用隨機增量構(gòu)造法，隨機增量構(gòu)造法包括以下步驟:
[0084]構(gòu)造相對于原始點集P大于一定閾值的點集P-1PjPy使用該點集把原始點集P包圍起來，并以隨機的次序依次把原始點集P中的各個點插入到點集P-A2P-的三角剖分中，在插入的過程中始終保證得到的三角剖分是Delaunay三角剖分，在原始點集P的所有的點都插入以后，獲得{P-!, P_2，P_3} U P的Delaunay三角剖分，將P^1, P_2，P_3以及與其相關(guān)聯(lián)的各邊刪除，獲得原始點集P的Delaunay三角剖分。
[0085]S104、基于Delaunay三角形網(wǎng)，生成Voronoi網(wǎng)格圖。
[0086]Voronoi圖，又叫泰森多邊形或Dirichlet圖，它是由一組由連接兩鄰點直線的垂直平分線組成的連續(xù)多邊形組成。N個在平面上有區(qū)別的點，按照最鄰近原則劃分平面；每個點與它的最近鄰區(qū)域相關(guān)聯(lián)。Delaunay三角形是由與相鄰Voronoi多邊形共享一條邊的相關(guān)點連接而成的三角形。Delaunay三角形的外接圓圓心是與三角形相關(guān)的Voronoi多邊形的一個頂點。Voronoi三角形是Delaunay圖的對偶圖，如圖7所示，為互為對偶的Voronoi 圖與 Delaunay 圖。
[0087]優(yōu)選地，Voronoi網(wǎng)格圖是指:由一組連接兩相鄰點的直線的垂直平分線組成的連續(xù)多邊形所組成的網(wǎng)格。
[0088]基于Delaunay三角形網(wǎng)，生成Voronoi網(wǎng)格圖是指:基于Delaunay三角形網(wǎng)中的Delaunay三角形，獲得Delaunay三角形中的每條邊的垂直平分線，連接Delaunay三角形網(wǎng)中的全部垂直平分線，并使每個連接點與Delaunay三角形的外接圓的圓心重合，連接全部連接點，獲得Voronoi網(wǎng)格圖。
[0089]S105、根據(jù)Voronoi網(wǎng)格圖計算低滲氣藏產(chǎn)能。
[0090]優(yōu)選地，本發(fā)明中的二維油藏模型為二維PEBI網(wǎng)格；Voronoi網(wǎng)格是三維PEBI網(wǎng)格。
[0091]綜上所述，本發(fā)明是將油藏對象加上限定條件和檢驗條件，得到平面二維PEBI網(wǎng)格，再由二維PEBI網(wǎng)格獲得三維PEBI網(wǎng)格，利用獲得的三維PEBI網(wǎng)格對低滲氣藏產(chǎn)能進行計算。
[0092]生成的PEBI網(wǎng)格取決于油藏區(qū)域的點集內(nèi)點的分布。按照模型描述方法，油藏中的構(gòu)造在平面上被抽象為點、多邊形和圓等圖形元素。這些元素在點集處理模塊中會被處理為平面上的有序點。平面模型預處理模塊就是完成這部分的工作。除了元素的有序點集化以外，還需要對點集之間的影響進行處理，例如要將一些距離較近可能存在干擾的點進行判斷和處理，避免后期三角化的時候影響。
[0093]利用平面網(wǎng)格化方法得到的二維網(wǎng)格可以使用一定的算法擴展為三維網(wǎng)格，經(jīng)過研宄得出，在平面網(wǎng)格化的基礎(chǔ)上，利用其網(wǎng)格化框架，再將其網(wǎng)格化結(jié)果在垂直方向上擴展。即可得到三維網(wǎng)格。按照這種方法將平面PEBI網(wǎng)格擴展到三維網(wǎng)格，得到的是近似的三維PEBI，即平面上連接面是正交網(wǎng)格，但是縱向連接的網(wǎng)格連接面并不是正交的。不過，這樣的網(wǎng)格已經(jīng)比通常的矩形或長方體網(wǎng)格具有更好的正交性，計算的平面上由于網(wǎng)格尺寸帶來的誤差影響已經(jīng)大大減小。
[0094]本發(fā)明的技術(shù)方案具有以下技術(shù)效果:
[0095]1、本方法能夠得到毫米級別的精細網(wǎng)格，提高了常規(guī)非結(jié)構(gòu)化網(wǎng)格的精度。
[0096]2、本發(fā)明增加了輔助點與輔助線，并優(yōu)化了增加的輔助點、線與固定圖形的三角化方法。
[0097]3、本發(fā)明建立的模型將毫米級別的網(wǎng)格與50米級別的網(wǎng)格耦合在一起，克服了計算不收斂等數(shù)值計算問題。
[0098]需要說明的是，以上所述的實施例僅是為了便于本領(lǐng)域的技術(shù)人員理解而已，并不用于限制本發(fā)明的保護范圍，在不脫離本發(fā)明的發(fā)明構(gòu)思的前提下，本領(lǐng)域技術(shù)人員對本發(fā)明所做出的任何顯而易見的替換和改進等均在本發(fā)明的保護范圍之內(nèi)。
【權(quán)利要求】
1.一種低滲氣藏產(chǎn)能計算方法，其特征在于，所述方法包括: 將油藏對象投影到二維平面上，獲得所述油藏對象的二維油藏模型； 將所述二維油藏模型中的所有所述油藏對象表示成有序點集； 對所述有序點集進行德勞內(nèi)Delaunay三角剖分，獲得Delaunay三角形網(wǎng)； 基于所述Delaunay三角形網(wǎng)，生成Voronoi網(wǎng)格圖； 根據(jù)所述Voronoi網(wǎng)格圖計算低滲氣藏產(chǎn)能。
2.如權(quán)利要求1所述的方法，其特征在于，所述油藏對象包括:油藏邊界、層面、斷層、井和水力壓裂裂縫。
3.如權(quán)利要求1或2所述的方法，其特征在于，所述將油藏對象投影到二維平面上，獲得所述油藏對象的二維油藏模型是指: 將所述油藏對象投影到二維平面上，所述邊界轉(zhuǎn)化為包絡(luò)線，所述斷層轉(zhuǎn)化為多邊形，所述井轉(zhuǎn)化為點，斜度大于預定斜度閾值的井或者水平井轉(zhuǎn)化為多邊形，所述水力壓裂裂縫轉(zhuǎn)化為線段。
4.如權(quán)利要求3所述的方法，其特征在于，所述將所述二維油藏模型中的所有所述油藏對象表示成有序點集是指:將所述二維油藏模型中的所述包絡(luò)線、所述多邊形、所述線段以及所述點轉(zhuǎn)化為二維平面上不同的有序點集。
5.如權(quán)利要求1所述的方法，其特征在于，所述方法還包括:當所述油藏對象非均勻分布和/或所述油藏對象的所述有序點集中的散點個數(shù)超過預定的散點閾值時，增加所述二維油藏模型中的所述油藏對象的限定條件，將所述油藏對象的所述有序點集進行細分； 其中，所述增加所述二維油藏模型中的所述油藏對象的限定條件，將所述油藏對象的所述有序點集進行細分包括:增加每個和/或個別所述油藏對象的局部邊界區(qū)域，在所述局部邊界區(qū)域內(nèi)部對該油藏對象的所述有序點集進行細分。
6.如權(quán)利要求1所述的方法，其特征在于，所述對所述有序點集進行德勞內(nèi)Delaunay三角剖分，獲得Delaunay三角形網(wǎng)的方法包括以下步驟: 61、遍歷所述有序點集的所有散點，得出所述有序點集的包容盒，獲得作為所述有序點集的凸殼的初始三角形并放入Delaunay三角形鏈表； 62、將所述有序點集中的一個散點插入所述Delaunay三角形鏈表，在所述Delaunay三角形鏈表中找出其外接圓包含插入的所述散點的三角形，刪除影響所述三角形的公共邊，將插入的所述散點與影響所述三角形的全部頂點連接起來，完成一個所述散點在Delaunay三角形鏈表中的插入； 63、根據(jù)Delaunay三角形優(yōu)化準則對局部新形成的三角形進行優(yōu)化，將優(yōu)化后的所述三角形放入所述Delaunay三角形鏈表； 64、判斷所述有序點集中是否還有未插入的散點，如果有則進行步驟62;如果沒有則結(jié)束。
7.如權(quán)利要求6所述的方法，其特征在于，當引入地面的地性線和特征線時，所述對所述有序點集進行德勞內(nèi)Delaunay三角剖分，獲得Delaunay三角形網(wǎng)的方法包括以下步驟: 71、根據(jù)已有的所述地性線和所述特征線，形成控制邊鏈表； 72、以所述控制邊鏈表中的一個線段為基邊，從所述有序點集中找出同所述基邊的兩端點的距離和最小的點，以該點為頂點，以所述基邊為邊，向外擴展一個三角形并放入三角形鏈表中； 73、判斷所述控制邊鏈表中是否還有未擴展的線段，如果有則進行步驟72;如果沒有則進行步驟74 ； 74、依次將新形成的三角形的邊作為基邊，形成新的控制邊鏈表，按照步驟72，對所述新的控制邊鏈表中的所有的線段進行擴展； 75、判斷所述新的控制邊鏈表中是否還有未擴展的所述新形成的三角形，如果有則進行步驟74;如果沒有則結(jié)束。
8.如權(quán)利要求1所述的方法，其特征在于，所述方法還包括:當所述油藏對象非均勻分布和/或所述油藏對象的所述有序點集中的散點個數(shù)超過預定的散點閾值時，在對所述點集進行所述Delaunay三角剖分時，如果翻邊后獲得的三角形仍然不能滿足Delaunay三角形的基本準則，則在該點集內(nèi)部增加新的點作為輔助點，或增加新的點集作為輔助點集繼續(xù)進行所述Delaunay三角剖分； 其中，在所述該點集內(nèi)部增加新的點作為輔助點，或增加新的點集作為輔助點集采用隨機增量構(gòu)造法，所述隨機增量構(gòu)造法包括以下步驟: 構(gòu)造相對于原始點集P大于一定閾值的點集PJjPy使用該點集把所述原始點集P包圍起來，并以隨機的次序依次把所述原始點集P中的各個點插入到所述點集P_iP_2P_3的三角剖分中，在插入的過程中始終保證得到的三角剖分是Delaunay三角剖分，在所述原始點集P的所有的點都插入以后，獲得{P-!,，P_2，P_3} U P的Delaunay三角剖分，將P+ P_2，P_3以及與其相關(guān)聯(lián)的各邊刪除，獲得所述原始點集P的Delaunay三角剖分。
9.如權(quán)利要求1所述的方法，其特征在于， 所述VOTonoi網(wǎng)格圖是指:由一組連接兩相鄰點的直線的垂直平分線組成的連續(xù)多邊形所組成的網(wǎng)格； 所述基于所述Delaunay三角形網(wǎng)，生成Voronoi網(wǎng)格圖是指:基于所述Delaunay三角形網(wǎng)中的Delaunay三角形，獲得所述Delaunay三角形中的每條邊的垂直平分線，連接所述Delaunay三角形網(wǎng)中的全部所述垂直平分線，并使每個連接點與所述Delaunay三角形的外接圓的圓心重合，連接全部所述連接點，獲得所述VOTonoi網(wǎng)格圖。
10.如權(quán)利要求9所述的方法，其特征在于， 所述二維油藏模型為二維PEBI網(wǎng)格；所述Voronoi網(wǎng)格是三維PEBI網(wǎng)格。
【文檔編號】G06F19/00GK104504298SQ201510038773
【公開日】2015年4月8日 申請日期:2015年1月26日 優(yōu)先權(quán)日:2015年1月26日
【發(fā)明者】劉子雄, 王杏尊, 吳英, 李敬松, 黃子俊, 秦俊杰, 曾鳴 申請人:中國海洋石油總公司, 中海油田服務股份有限公司