[發(fā)明專利]一種分段壓裂水平井的粗網(wǎng)格模擬方法及系統(tǒng)

有效

• 申請?zhí)枺?/td>: 201710432128.9
• 申請日: 2017-06-09
• 公開（公告）號: CN107169227B
• 公開（公告）日: 2019-04-05
• 發(fā)明（設(shè)計(jì)）人: 黃朝琴;梁夢吟;姚軍;黃濤
• 申請（專利權(quán)）人: 中國石油大學(xué)（華東）
• 主分類號: G06F17/50
• 分類號: G06F17/50;E21B43/26
• 代理公司: 北京高沃律師事務(wù)所 11569
• 代理人: 王戈

1.一種分段壓裂水平井的粗網(wǎng)格模擬方法，其特征在于，包括：

獲取油藏模型；

根據(jù)所述油藏模型確定多個裂縫信息；所述裂縫信息包括裂縫的長度、寬度和裂縫位置；

根據(jù)所述油藏模型采用邊界元法計(jì)算所述油藏模型的穩(wěn)態(tài)滲流方程的半解析解矩陣組；所述穩(wěn)態(tài)滲流方程為地層壓力達(dá)到穩(wěn)定，不再變化的狀態(tài)下的滲流方程；所述半解析解矩陣組是由壓力半解析解和流量半解析解構(gòu)成的矩陣組；

根據(jù)所述油藏模型建立笛卡爾網(wǎng)格；所述笛卡爾網(wǎng)格包括含有裂縫的網(wǎng)格、含有基質(zhì)的網(wǎng)格、含有井筒的網(wǎng)格；

根據(jù)所述半解析解矩陣組計(jì)算所述含有裂縫-基質(zhì)的網(wǎng)格的裂縫-基質(zhì)傳導(dǎo)率、含有井筒的網(wǎng)格的井指數(shù)以及相鄰網(wǎng)格之間的相鄰網(wǎng)格傳導(dǎo)率；所述裂縫-基質(zhì)傳導(dǎo)率表示為在某個含有裂縫的網(wǎng)格中，從基質(zhì)流入裂縫的流量與第一壓力差的比值；所述第一壓力差為裂縫所在網(wǎng)格和基質(zhì)所在網(wǎng)格的壓力差；所述井指數(shù)表示為流經(jīng)井筒的流量與第二壓力差的比值；所述第二壓力差為井所在網(wǎng)格壓力和井底的壓力差；所述相鄰網(wǎng)格傳導(dǎo)率表示為兩個相鄰網(wǎng)格之間流經(jīng)兩個相鄰網(wǎng)格交界面處的流量與第三壓力差的比值；所述第三壓力差為兩個相鄰網(wǎng)格之間的壓力差；

將所述裂縫-基質(zhì)傳導(dǎo)率、所述井指數(shù)和所述相鄰網(wǎng)格傳導(dǎo)率輸入到油藏?cái)?shù)值模擬器中，輸出得到油氣藏信息；所述油氣藏信息包括壓力、飽和度、產(chǎn)量。

2.根據(jù)權(quán)利要求1的所述粗網(wǎng)格模擬方法，其特征在于，所述根據(jù)所述半解析解矩陣組計(jì)算所述含有裂縫-基質(zhì)的網(wǎng)格的裂縫-基質(zhì)傳導(dǎo)率、含有井筒的網(wǎng)格的井指數(shù)以及相鄰網(wǎng)格之間的相鄰網(wǎng)格傳導(dǎo)率，具體包括：

獲取穩(wěn)態(tài)下含有基質(zhì)的網(wǎng)格的第一平均壓力；

獲取穩(wěn)態(tài)下所述含有裂縫的網(wǎng)格的第二平均壓力；

獲取穩(wěn)態(tài)下從所述含有基質(zhì)的網(wǎng)格流入所述含有裂縫的網(wǎng)格的第一流量；

根據(jù)所述第一平均壓力、所述第二平均壓力和所述第一流量計(jì)算所述裂縫-基質(zhì)傳導(dǎo)率。

3.根據(jù)權(quán)利要求1的所述粗網(wǎng)格模擬方法，其特征在于，所述根據(jù)所述半解析解矩陣組計(jì)算所述含有裂縫-基質(zhì)的網(wǎng)格的裂縫-基質(zhì)傳導(dǎo)率、含有井筒的網(wǎng)格的井指數(shù)以及相鄰網(wǎng)格之間的相鄰網(wǎng)格傳導(dǎo)率，具體包括：

獲取穩(wěn)態(tài)下的含有井筒網(wǎng)格的相鄰網(wǎng)格的第三平均壓力；

獲取穩(wěn)態(tài)下含有井筒的網(wǎng)格的第四平均壓力；

獲取所述含有井筒網(wǎng)格的相鄰網(wǎng)格流入所述含有井筒的網(wǎng)格的第二流量；

根據(jù)所述第三平均壓力、所述第四平均壓力以及所述第二流量計(jì)算所述井指數(shù)。

4.根據(jù)權(quán)利要求1的所述粗網(wǎng)格模擬方法，其特征在于，所述根據(jù)所述半解析解矩陣組計(jì)算所述含有裂縫-基質(zhì)的網(wǎng)格的裂縫-基質(zhì)傳導(dǎo)率、含有井筒的網(wǎng)格的井指數(shù)以及相鄰網(wǎng)格之間的相鄰網(wǎng)格傳導(dǎo)率，具體包括：

獲取任意兩個相鄰網(wǎng)格交界處的滲透率；所述兩個相鄰網(wǎng)格用第一網(wǎng)格和第二網(wǎng)格表示；

計(jì)算所述第一網(wǎng)格的中心點(diǎn)與所述交界處的第一傳導(dǎo)率；

計(jì)算所述第二網(wǎng)格的中心點(diǎn)與所述交界處的第二傳導(dǎo)率；

根據(jù)所述滲透率、所述第一傳導(dǎo)率和所述第二傳導(dǎo)率計(jì)算所述第一網(wǎng)格和第二網(wǎng)格之間的相鄰網(wǎng)格傳導(dǎo)率。

5.根據(jù)權(quán)利要求4的所述粗網(wǎng)格模擬方法，其特征在于，所述獲取任意兩個相鄰網(wǎng)格交界處的滲透率，具體包括：

獲取第一網(wǎng)格的第一滲透率；

獲取第二網(wǎng)格的第二滲透率；

根據(jù)所述第一傳導(dǎo)率計(jì)算所述第一網(wǎng)格到所述交界處的第一距離權(quán)函數(shù)；

根據(jù)所述第二傳導(dǎo)率計(jì)算所述第一網(wǎng)格到所述交界處的第二距離權(quán)函數(shù)；

根據(jù)所述第一滲透率、第二滲透率、第一距離權(quán)函數(shù)以及第二距離權(quán)函數(shù)計(jì)算所述兩個相鄰網(wǎng)格交界處的滲透率。