常州大學(xué)史文洋獲國家專利權(quán)
買專利賣專利找龍圖騰,真高效! 查專利查商標(biāo)用IPTOP,全免費(fèi)!專利年費(fèi)監(jiān)控用IP管家,真方便!
龍圖騰網(wǎng)獲悉常州大學(xué)申請的專利一種多分支斷溶體油藏縫洞參數(shù)的確定方法及系統(tǒng)獲國家發(fā)明授權(quán)專利權(quán),本發(fā)明授權(quán)專利權(quán)由國家知識(shí)產(chǎn)權(quán)局授予,授權(quán)公告號(hào)為:CN115559714B 。
龍圖騰網(wǎng)通過國家知識(shí)產(chǎn)權(quán)局官網(wǎng)在2025-08-15發(fā)布的發(fā)明授權(quán)授權(quán)公告中獲悉:該發(fā)明授權(quán)的專利申請?zhí)?專利號(hào)為:202211062078.7,技術(shù)領(lǐng)域涉及:E21B49/00;該發(fā)明授權(quán)一種多分支斷溶體油藏縫洞參數(shù)的確定方法及系統(tǒng)是由史文洋;蔣征霖;陶磊;高敏;劉永川;朱慶杰;白佳佳;張娜;劉現(xiàn)坤設(shè)計(jì)研發(fā)完成,并于2022-08-31向國家知識(shí)產(chǎn)權(quán)局提交的專利申請。
本一種多分支斷溶體油藏縫洞參數(shù)的確定方法及系統(tǒng)在說明書摘要公布了:本發(fā)明公開了一種多分支斷溶體油藏縫洞參數(shù)的確定方法及系統(tǒng),所述方法包括:獲取關(guān)井后測試井井底壓力隨時(shí)間變化的實(shí)測數(shù)據(jù);根據(jù)多分支斷溶體油藏縫洞的基礎(chǔ)地質(zhì)資料,確定預(yù)先構(gòu)建的多分支斷溶體合采油藏井底壓力數(shù)學(xué)模型的初始參數(shù),將初始參數(shù)作為擬合參數(shù)的初始值,求解得到多分支斷溶體合采油藏井底壓力解;計(jì)算實(shí)測數(shù)據(jù)的壓力與求解得到的壓力解之間的差值;若差值小于預(yù)設(shè)的擬合誤差,則將初始參數(shù)作為最終擬合參數(shù);若差值大于等于預(yù)設(shè)的擬合誤差,修正擬合參數(shù),直到差值小于預(yù)設(shè)的擬合誤差,得到最終擬合參數(shù)。本發(fā)明僅通過一次壓力數(shù)據(jù)獲取就能夠確定多分支斷溶油藏縫洞空間分布。
本發(fā)明授權(quán)一種多分支斷溶體油藏縫洞參數(shù)的確定方法及系統(tǒng)在權(quán)利要求書中公布了:1.一種多分支斷溶體油藏縫洞參數(shù)的確定方法,其特征在于,包括: 獲取關(guān)井后測試井井底壓力隨時(shí)間變化的實(shí)測數(shù)據(jù); 根據(jù)多分支斷溶體油藏縫洞的基礎(chǔ)地質(zhì)資料,確定預(yù)先構(gòu)建的多分支斷溶體合采油藏井底壓力數(shù)學(xué)模型的初始參數(shù),將初始參數(shù)作為擬合參數(shù)的初始值,求解得到多分支斷溶體合采油藏井底壓力解; 其中,所述預(yù)先構(gòu)建的多分支斷溶體合采油藏井底壓力數(shù)學(xué)模型,是根據(jù)水平井筒的水平管流、裂縫區(qū)域線性滲流、溶洞區(qū)域儲(chǔ)集流的特征,考慮縫洞初始條件和邊界條件構(gòu)建的; 其中,根據(jù)水平井筒的水平管流的特征:水平井穩(wěn)定滿管層流的流量與驅(qū)動(dòng)壓差之間關(guān)系通過以下數(shù)學(xué)模型表征: 式1中,△pw為井筒長度對(duì)應(yīng)的壓力差,Pa;qw為井筒區(qū)域流量,m3s;μ為粘度,Pa·s;d為井筒長度,m;rw為井筒半徑,m; 其中,根據(jù)裂縫區(qū)域線性滲流的特征,描述流體滲流的數(shù)學(xué)微分方程為: 式2中,k為裂縫區(qū)域滲透率,m2;μ為粘度,Pa·s;pf為裂縫區(qū)域壓力,Pa;h為深度,m;ρ為流體密度kgm3;g為重力加速度,9.8ms2;cl為流體壓縮系數(shù),Pa-1;φ為孔隙度,%;ct為儲(chǔ)層綜合壓縮系數(shù),Pa-1;t為時(shí)間,s; 儲(chǔ)層裂縫區(qū)域的流量由壓力變化得到: 式3中,qf為裂縫區(qū)域流量,m3s;Af為裂縫區(qū)域滲流面積,m2;k為裂縫區(qū)域滲透率,m2;pf為裂縫區(qū)域壓力,Pa;t為時(shí)間,s; 其中,根據(jù)溶洞區(qū)域儲(chǔ)集流的特征,描述溶洞流體的凈流量與溶洞壓力變化之間關(guān)系的數(shù)學(xué)方程為: 式4中,qv為溶洞區(qū)域流量,m3s;Cv為溶洞儲(chǔ)集系數(shù),m3Pa;pv為溶洞區(qū)域壓力,Pa;t為時(shí)間,s; 其中,縫洞初始條件為開井生產(chǎn)前儲(chǔ)層各處壓力相同且等于儲(chǔ)層原始?jí)毫Γ瑒t初始?jí)毫Φ臄?shù)學(xué)模型表達(dá)為: ptfh,0=pv0=pbfh,0=pi5 式5中,ptf為溶洞上部裂縫區(qū)域壓力,Pa;pv為溶洞區(qū)域壓力,Pa;pbf為溶洞下部裂縫區(qū)域壓力,Pa;pi為儲(chǔ)層原始?jí)毫Γ琍a;h為深度,m; 其中,縫洞邊界條件包括: 內(nèi)邊界條件: 井筒產(chǎn)量等于所有溶洞上部裂縫區(qū)域流量之和: 井底壓力等于近井分支儲(chǔ)層壓力與近井水平段壓力差之和: pwt=Δpw1+ptf1t7 式6、7中,Q為井底流量,m3s;qtf為溶洞上部裂縫區(qū)域流量,m3s;qtfj為j分支溶洞上部裂縫區(qū)域流量,m3s;pw為井底壓力,Pa;△pw1為近井段井筒長度對(duì)應(yīng)的壓力差,Pa;t為時(shí)間,s;ptf1為第1分支裂縫在井筒處的壓力,Pa; 交界面連接條件: 流體從溶洞下部裂縫系統(tǒng),經(jīng)過溶洞區(qū)域,流向溶洞上部裂縫區(qū);在溶洞與上部、下部裂縫區(qū)域的界面處,流量相同: 在溶洞與上部、下部裂縫區(qū)域的界面處,壓力相同: ptfhtv,t=pvt=pbfhbv,t9 式8、9中,Cv為溶洞儲(chǔ)集系數(shù),m3Pa;htv為溶洞頂面深度,m;hbv為溶洞底面深度,m; 外邊界條件: 儲(chǔ)層外邊界為封閉邊界,無流體流動(dòng): 式中:hfb為溶洞下部裂縫區(qū)域底度,m; 其中,所述求解得到多分支斷溶體合采油藏井底壓力解,包括: 將數(shù)學(xué)模型中的參數(shù)進(jìn)行無量綱化,得到:無量綱壓力k1為第1分支儲(chǔ)層滲透率,m2,Af1為第1分支儲(chǔ)層投影面積,m2;p為壓力,Pa;無量綱時(shí)間φ為孔隙度,%;j為j分支;無量綱流量q為流量,m3d;無量綱儲(chǔ)集系數(shù)C為儲(chǔ)集系數(shù),m3Pa;無量綱深度無量綱重力系數(shù)無量綱分支儲(chǔ)層距離無量綱流動(dòng)系數(shù)比Afj為第j分支儲(chǔ)層投影面積,m2;無量綱壓力傳導(dǎo)系數(shù)比 經(jīng)過無量綱化和拉普拉斯空間變換后,任意j分支溶洞上部裂縫、下部裂縫壓力為: 式11中,為拉普拉斯空間下的第j分支上部裂縫無量綱壓力,hD為垂直方向上無量綱深度,無量綱;為拉普拉斯空間下的第j分支下部裂縫無量綱壓力,無量綱; 任意j分支溶洞上部裂縫流動(dòng)方程為: 式12中,htvDj為第j分支溶洞頂面無量綱深度,無量綱; 任意j分支溶洞流動(dòng)方程為: 式13中,CvDj為第j分支溶洞區(qū)域無量綱儲(chǔ)集系數(shù),無量綱;為拉普拉斯空間下第j分支溶洞區(qū)域無量綱壓力,無量綱; 任意j段井筒流動(dòng)方程為: 井筒流量條件為: 井筒壓力條件為: 式16中,為拉普拉斯空間下第j分支上部裂縫區(qū)域無量綱壓力,無量綱;dD1為第1分支儲(chǔ)層到井底的無量綱距離,無量綱;為拉普拉斯空間下井底無量綱壓力,無量綱; 任意j分支下部裂縫-溶洞-上部裂縫,流量條件為: 式17中,CvDj為第j分支溶洞區(qū)域無量綱儲(chǔ)集系數(shù),無量綱;任意j分支下部裂縫-溶洞-上部裂縫,壓力條件為: 式18中,hbvDj為第j分支溶洞底面無量綱深度,無量綱; 儲(chǔ)層外邊界處為封邊邊界,則流量為0: 式19中,hbfDj為第j分支下部裂縫底面無量綱深度,無量綱; 方程11對(duì)應(yīng)解的形式為: 式20中,c為待求系數(shù),r為共軛特征根,其表達(dá)形式為: 式21中,為第j分支上部裂縫區(qū)域流動(dòng)方程對(duì)應(yīng)特征根,無量綱;為第j分支下部裂縫區(qū)域流動(dòng)方程對(duì)應(yīng)特征根,無量綱; 將方程20帶入內(nèi)邊界條件15-16、交界面連接條件17-18、外邊界條件19,得到關(guān)于待求系數(shù)c的線性方程組: 式22中,ctj為第j分支上部裂縫區(qū)域流動(dòng)方程中待求系數(shù),無量綱,下標(biāo)j={1,…,n};cbj為第j分支下部裂縫區(qū)域流動(dòng)方程中待求系數(shù),無量綱,下標(biāo)j={1,…,n}; 矩陣D為4n×4n個(gè)元素,除了以下元素,其余均為0; 第一行元素為: 第2~n行元素為: 式22-2中,λ為矩陣D行序號(hào),無量綱,取值為λ={2,…,n};dDj為第j分支儲(chǔ)層到井底的無量綱距離,無量綱;第n+1~2n行元素為: 第2n+1~3n行元素為: 第3n+1~4n行元素為: 對(duì)于雙分支模型,D具體為: 對(duì)于單分支,D具體為: 解方程組22得到所有待求系數(shù)ct1,則根據(jù)公式16可得井底壓力解為 式23中,為拉普拉斯空間下無量綱井底壓力,無量綱;為第1分支上部裂縫區(qū)域滲流方程中待求系數(shù),無量綱;為第1分支上部裂縫區(qū)域滲流方程中待求共軛系數(shù),無量綱;dD1為第1分支儲(chǔ)層到井底的無量綱距離,無量綱; 通過Stehfest數(shù)值積分算法對(duì)Laplace空間無量綱井底壓力23進(jìn)行反演,得到實(shí)空間下無因次井底壓力隨時(shí)間變化的求解數(shù)據(jù)pwD; 計(jì)算實(shí)測數(shù)據(jù)的壓力與求解得到的壓力解之間的差值;若差值小于預(yù)設(shè)的擬合誤差,則將初始參數(shù)作為最終擬合參數(shù);若差值大于等于預(yù)設(shè)的擬合誤差,則修正擬合參數(shù),直到差值小于預(yù)設(shè)的擬合誤差,將最后修正的擬合參數(shù)作為最終擬合參數(shù);所述最終擬合參數(shù)即為多分支斷溶體油藏縫洞的參數(shù)。
如需購買、轉(zhuǎn)讓、實(shí)施、許可或投資類似專利技術(shù),可聯(lián)系本專利的申請人或?qū)@麢?quán)人常州大學(xué),其通訊地址為:213164 江蘇省常州市武進(jìn)區(qū)滆湖中路21號(hào);或者聯(lián)系龍圖騰網(wǎng)官方客服,聯(lián)系龍圖騰網(wǎng)可撥打電話0551-65771310或微信搜索“龍圖騰網(wǎng)”。
1、本報(bào)告根據(jù)公開、合法渠道獲得相關(guān)數(shù)據(jù)和信息,力求客觀、公正,但并不保證數(shù)據(jù)的最終完整性和準(zhǔn)確性。
2、報(bào)告中的分析和結(jié)論僅反映本公司于發(fā)布本報(bào)告當(dāng)日的職業(yè)理解,僅供參考使用,不能作為本公司承擔(dān)任何法律責(zé)任的依據(jù)或者憑證。
- 江蘇大學(xué)王晗昊獲國家專利權(quán)
- 波士頓科學(xué)國際有限公司曼·明·恩古耶恩獲國家專利權(quán)
- 深圳盈特創(chuàng)智能科技有限公司瞿晗昕獲國家專利權(quán)
- 上海鴻研物流技術(shù)有限公司請求不公布姓名獲國家專利權(quán)
- 巴斯夫公司C·塔蘭塔獲國家專利權(quán)
- 廣州市華濾環(huán)保設(shè)備有限公司曾志松獲國家專利權(quán)
- 合肥泰禾光電科技股份有限公司許大紅獲國家專利權(quán)
- 浙江富瑞浦衛(wèi)浴科技有限公司宣品安獲國家專利權(quán)
- 伊利諾斯工具制品有限公司馮吉林獲國家專利權(quán)
- 博西華電器(江蘇)有限公司薄濤獲國家專利權(quán)


熱門推薦
- 科尼全球公司N·勞卡寧獲國家專利權(quán)
- 中國石油化工股份有限公司張?jiān)诗@國家專利權(quán)
- 成都市瓏熙科技有限公司肖仁旺獲國家專利權(quán)
- 中興通訊股份有限公司邱徵虹獲國家專利權(quán)
- 四川瑞可達(dá)連接系統(tǒng)有限公司黃應(yīng)德獲國家專利權(quán)
- 博西華電器(江蘇)有限公司張玲獲國家專利權(quán)
- 深圳市斯永威科技有限公司歐陽俊偉獲國家專利權(quán)
- 江蘇風(fēng)和醫(yī)療器材股份有限公司孫寶峰獲國家專利權(quán)
- 易泰博商業(yè)設(shè)備(中國)有限公司丁萬春獲國家專利權(quán)
- 青島海爾智能技術(shù)研發(fā)有限公司楊少強(qiáng)獲國家專利權(quán)