油氣儲運網(wǎng)

 找回密碼
 立即注冊

QQ登錄

QQ快速登錄

一步迅速開始

pdms培訓(xùn)&項目數(shù)字化360網(wǎng)
查看: 5694|回復(fù): 2
打印 上一主題 下一主題

[軟件介紹] 并行算法在射孔產(chǎn)能研究應(yīng)用

[復(fù)制鏈接]
跳轉(zhuǎn)到指定樓層
樓主
發(fā)表于 2012-7-9 14:30:28 | 只看該作者 回帖獎勵 |倒序瀏覽 |閱讀模式
本帖最后由 合肥辰工科技 于 2012-7-9 14:43 編輯

并行算法在射孔產(chǎn)能研究應(yīng)用
并行算法在射孔產(chǎn)能研究中的應(yīng)用. 天然氣工業(yè),1999;19 (6):63~65
摘 要 文章根據(jù)射孔完井的滲流特點,利用有限元方法建立了三維單相穩(wěn)定滲流的數(shù)值模型(詳見,數(shù)值模擬考慮了孔深
、孔徑、孔密、孔的相位、污染程度和厚度、壓實程度和厚度等射孔參數(shù)對油氣井產(chǎn)能的影響。有限元方程的求解是在曙光1000 大型并行機上利用區(qū)域分解法進行的,針對計算模型進行了并行算法加速比的研究,給出了不同計算結(jié)點的并行加速比,結(jié)果表明并行算法的計算速度與串行算法的計算速度明顯提高;通過各種射孔方案的計算,給出了各種射孔參數(shù)對油氣井產(chǎn)能影響的曲線,這些成果對于油氣田設(shè)計射孔方案具有一定的指導(dǎo)意義。
主題詞 射孔 射孔參數(shù) 并行計算機 有限元法 產(chǎn)能預(yù)測
目前,油氣田上射孔完井是較為廣泛的一種完井方式,射孔參數(shù)(孔深、孔徑、孔密、孔的相位、污染程度和厚度、壓實程度和厚度等) 對油氣井產(chǎn)能的影響是一個重要的研究課題。利用數(shù)值方法來研究射孔產(chǎn)能的機理起于60 年代〔1〕,較全面地考慮各種射孔參數(shù)對油氣井產(chǎn)能的影響始于80 年代〔2〕。由于射孔完井采油氣時形成的滲流場非常復(fù)雜,這就要求人們把滲流區(qū)域剖分得很細,才能滿足計算精度。由于受計算機內(nèi)存和計算速度的限制,人們在建立
計算模型時,總是要作一些假設(shè)或簡化。這次我們在曙光1000 上,利用目前普遍使用的區(qū)域分解法〔3〕來求解有限元方程組,取得了比較好的結(jié)果。通過各種方案的計算,給出了各種射孔參數(shù)對油氣井產(chǎn)能影響的曲線,這些成果對于設(shè)計射孔方案具有一定的指導(dǎo)意義。
計算模型
1. 數(shù)學(xué)模型
按照射孔完井滲流場的特點,數(shù)學(xué)模型采用三維單相穩(wěn)定滲流方程:
式中: p 為地層壓力,MPa ; K為地層滲透率, ;μ為地層原油氣粘度mPa ;Γ1 為已知壓力邊界;Γ2 為已知流量邊界(當q等于0 時,為不滲透邊界) ; nΓ2 的法線方向;Ω 為滲流區(qū)域。
2. 有限元方程的建立
在三維問題中最基本的單元是四面體單元,記為( e) ,它的四個頂點為有限元剖分的結(jié)點,其編號設(shè)為 i 、j k 、m ,根據(jù)里茨變分原理可得第 i個結(jié)點的有限元方程:
同理可寫出其他結(jié)點的有限元方程,就形成了有限元方程組。
3. 滲流區(qū)域的剖分
把三維滲流區(qū)域直接剖分成四面體單元,不僅難以繪出醒目的圖示,而且會使輸入的信息數(shù)量太大,因此在實際上都采用組合單元,最常用的是六面體單元,每個六面體組合單元又可剖分成五個四面體單元(有兩種不同的剖分形式)〔4〕。只要輸入組合單元的信息,可由計算機自動剖分成基本四面體單元進行計算(可按兩種剖分形式計算系數(shù)矩陣取每個系數(shù)的平均值) 。
六面體組合單元的剖分也是有計算機自動完成的,剖分的原則是: ①在平面上以井為中心,徑向剖分,夾角為10°。再以井為中心,做不同半徑的同心圓剖分,在射孔的頂端比較密,向外向內(nèi)變得比較疏(圖1 - a) ,共有30 層; ②在垂向上,以ΔZ 為步長進行剖分(圖1 - b) ; ③根據(jù)射孔的設(shè)計,在每個射孔處以孔為中心,用四個不同半徑的同心圓進行加密剖
分(圖1 - b) 。這樣就把整個滲流區(qū)域剖分成不同大小的六面體組合單元。
(a)平面上的剖分圖                   (b)射孔加密剖分圖
圖1  平面上的剖分和射孔加密剖分示意圖
有限元方程組的并行求解方法
有限元方程組的求解是在曙光1000 并行機上進行的。根據(jù)滲流區(qū)域剖分的特點,每一層同心圓上的剖分結(jié)點是相同的,比較適合區(qū)域分解法的并行算法。具體步驟如下。
有限元方程組形成一個N ×N 的大型稀疏系數(shù)矩陣如下:
由于系數(shù)矩陣很大,本文模型的剖分結(jié)點為36×30 ×20 = 21 600 ,再加上射孔處的加密結(jié)點,若按16 個射孔計算加密結(jié)點數(shù)為16 ×32 ×30 = 15 360 ,則系數(shù)矩陣是一個36 960 ×36 960 的矩陣。為了減少計算機內(nèi)存和計算結(jié)點之間的信息傳遞,在實際形成系數(shù)矩陣時,本文僅儲存非零元素,用兩個指示數(shù)組來確定元素的具體位置,用一個N ×1 的一維數(shù)組來儲存與每個結(jié)點有關(guān)的結(jié)點數(shù);用一個N ×22 的二維數(shù)組來儲存與每個結(jié)點有關(guān)的具體結(jié)點號,這兩個數(shù)組可在形成系數(shù)矩陣時同時形成。一般剖分成六面體組合單元每個結(jié)點周圍有26 個結(jié)點,但通過四面體建立基本單元后,實際上只有18個結(jié)點與其有關(guān),在射孔處有些結(jié)點周圍有33 個結(jié)點,實際有關(guān)的結(jié)點也僅有22 個結(jié)點,所以,系數(shù)矩陣變成了一個N ×22 的矩陣。
一維指示數(shù)組      二維指示數(shù)組      系數(shù)矩陣
                  
用超松弛迭代法求解有限元方程組,程序并行化采用部分重疊的區(qū)域分解法,具體做法如下:首先按30 層同心圓狀剖分不同層上的結(jié)點,平均地加載到不同計算結(jié)點上。比如加載到4 個計算結(jié)點上,1~8 層加載到0 號計算結(jié)點上;8~15 層加載到1 號計算結(jié)點上;15~22 層加載到2 號計算結(jié)點上;22~30 層加載到3 號計算結(jié)點上;每迭代一次,就將重疊層上的計算結(jié)果相互傳遞,取算術(shù)平均值作下一次迭代的初值,計算中僅需傳遞重疊層和相鄰層,如上
面重疊層(8 、15 、22 層) 和相鄰層(7 、9 、14 、16 、21 、23層) 。表1 是不同計算結(jié)點數(shù)計算時的并行加速比(串行程序計算時間與并行程序計算時間之比) 對比表,表明隨著計算結(jié)點的增多,計算速度明顯加快。
表1  并行計算加速比對比表
計算結(jié)點數(shù)
1
2
3
4
5
6
8
10
15
30
計算時間
(min)
計算加速比
78.9
1.00
51.6
1.53
35.9
2.20
29.1
2.71
25.7
3.07
22.4
3.52
19.2
4.11
15.9
4.96
13.3
5.93
11.0
7.17
射孔參數(shù)對油氣井產(chǎn)能的影響
這次研究所選取的數(shù)據(jù): 地層滲透率為0. 01 、油氣井半徑( Rw) 為0. 1 m ;根據(jù)前人的研究經(jīng)驗影響半徑(有限元計算的外邊界) 取30 倍的油氣井半徑,即S Rw = 30 Rw ;影響半徑處的地層壓力與射孔內(nèi)的壓力差為4 MPa 。在以上假定條件下,我們研究了射孔參數(shù)的變化對油氣井產(chǎn)能的影響。圖中的產(chǎn)率比是計算的射孔產(chǎn)能與裸眼井穩(wěn)定流產(chǎn)能的理論值之比。
圖2 是在沒有污染和壓實的情況下,射孔密度為每米16 個孔,孔徑為4 mm ,90°相位時射孔深度與孔徑的變化對產(chǎn)能的影響曲線,從圖中可看出,射孔深度和孔徑對產(chǎn)能的影響都比較大。
圖3 是在沒有污染和壓實的情況下,射孔密度為每米16 個孔,孔徑為4 mm ,90°相位時射孔相位不同對油氣井產(chǎn)能的影響曲線,圖中表明0°相位產(chǎn)能最低,螺旋式90°相位產(chǎn)能最大。
圖2  射孔深度和孔徑對產(chǎn)能影響曲線            圖3  射孔相位對產(chǎn)能的影響曲線
   
圖4   污染程度對產(chǎn)能的影響曲線              圖5   壓實程度對產(chǎn)能的影響曲線
圖4 是不考慮壓實,僅考慮污染的情況下,射孔密度為每米16 個孔,孔徑為4 mm ,90°相位,污染深度為0. 3 m 時,不同污染程度對產(chǎn)能的影響曲線,圖中Kd 代表污染區(qū)滲透率與原地層滲透率的比值。圖中可看出,在射孔深度小于污染深度時,污染程度對產(chǎn)能的影響較大;而射孔深度大于污染深度時,污染程度對產(chǎn)能的影響較小。
圖5 是不考慮污染,僅考慮壓實的情況下,射孔密度為每米16 個孔,孔徑為4 mm ,90°相位,壓實深度為0. 004 m 時,不同壓實程度對產(chǎn)能的影響曲線,圖中Kc 代表壓實區(qū)滲透率與原地層滲透率的比值。圖中可看出,壓實程度對產(chǎn)能的影響較大,這是因為壓實區(qū)包圍了整個射孔。
以上研究表明,利用并行算法求解地下滲流方程計算速度比串行算法的計算速度明顯提高,因此,并行算法應(yīng)該在射孔產(chǎn)能研究、油氣藏數(shù)值模擬等研究中被推廣應(yīng)用。本文給出的各種射孔參數(shù)對射孔產(chǎn)能的影響結(jié)果對油氣田設(shè)計射孔方案具有一定的指導(dǎo)意義。

變井筒存儲問題說明.doc

101 KB, 下載次數(shù): 1

售價: 3 金幣  [記錄]  [購買]

分形滲流分析研究.doc

594 KB, 下載次數(shù): 2

售價: 3 金幣  [記錄]  [購買]

并行算法在射孔產(chǎn)能研究應(yīng)用.doc

192 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

氣井的瞬時IPR 曲線.doc

154 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

閉式氣舉的井底壓力分析.doc

173 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

定產(chǎn)量煤層氣井底瞬時壓力計算.doc

339 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

間歇采油井的井底壓力計算.doc

163 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

考慮井筒熱效應(yīng)的氣體井底壓力計算.doc

419.5 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

考慮射孔的井底壓力計算.doc

194 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

扇形區(qū)域內(nèi)有井儲和表皮的井底壓力.doc

354.5 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

多層地層中的井筒及地層溫度解析解.doc

518 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

數(shù)值試井.doc

5.33 MB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

基于小信號提取及放大技術(shù)的早期試井項目研究開發(fā)報告.doc

1.4 MB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

早期補充材料.doc

289 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

早期分析理論.doc

430.5 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

早期試井分析軟件在大慶油田的應(yīng)用.doc

154 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

早期圖版適用條件.doc

26.5 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

多流量試井壓力偏差分析.doc

227.5 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

礫石充填完井多層油藏井底壓力評價.doc

260 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]

使用隨鉆錄井資料計算地層靜溫.doc

168 KB, 下載次數(shù): 0

售價: 3 金幣  [記錄]  [購買]


油氣儲運網(wǎng) - 論壇版權(quán)1、本主題所有言論和圖片純屬會員個人意見,與本論壇立場無關(guān)
2、本站所有主題由該帖子作者發(fā)表,該帖子作者與油氣儲運網(wǎng)享有帖子相關(guān)版權(quán)
3、其他單位或個人使用、轉(zhuǎn)載或引用本文時必須同時征得該帖子作者和油氣儲運網(wǎng)的同意
4、帖子作者須承擔(dān)一切因本文發(fā)表而直接或間接導(dǎo)致的民事或刑事法律責(zé)任
5、本帖部分內(nèi)容轉(zhuǎn)載自其它媒體,但并不代表本站贊同其觀點和對其真實性負責(zé)
6、如本帖侵犯到任何版權(quán)問題,請立即告知本站,本站將及時予與刪除并致以最深的歉意
7、油氣儲運網(wǎng)管理員和版主有權(quán)不事先通知發(fā)貼者而刪除本文

沙發(fā)
發(fā)表于 2012-7-9 15:04:44 | 只看該作者
好多資料啊  好像是鉆井專業(yè)的
板凳
發(fā)表于 2012-7-13 13:10:57 | 只看該作者
好好好好!
您需要登錄后才可以回帖 登錄 | 立即注冊  

本版積分規(guī)則

掃碼訪問手機版

QQ|關(guān)于我們|網(wǎng)站地圖|油氣儲運網(wǎng) ( 魯ICP備11007657號-3 )

GMT+8, 2024-12-23 22:03 , Processed in 0.029468 second(s), 13 queries , Gzip On, MemCached On.

Powered by Discuz! X3.5

Copyright © 2001-2023, Tencent Cloud.

快速回復(fù) 返回列表