[發明專利]化探正演數值模擬方法、裝置、計算機設備和存儲介質

在審

• 申請號: 202010883711.3
• 申請日: 2020-08-28
• 公開（公告）號: CN114117947A
• 公開（公告）日: 2022-03-01
• 發明（設計）人: 王國建;趙克斌;鮑征宇;楊俊;楊瑞琰;李武;盧麗
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發研究院
• 主分類號: G06F30/28
• 分類號: G06F30/28;G06F111/10;G06F113/08;G06F119/14
• 代理公司: 北京聿宏知識產權代理有限公司 11372
• 代理人: 吳大建;金淼

1.一種基于油藏模型的化探正演數值模擬方法，其特征在于，包括：

獲取油氣地質數據和地震數據，確定油藏地質模型；

獲取擴散數學模型、水溶相運移數學模型和浮力運移數學模型；

基于已確定的所述油藏地質模型，對所述擴散數學模型、所述水溶相運移數學模型和所述浮力運移數學模型耦合，建立油藏烴類微滲漏數學模型；

基于所述油藏烴類微滲漏數學模型進行分段模擬。

2.根據權利要求1所述的方法，其特征在于，所述基于已確定的所述油藏地質模型，對所述擴散數學模型、所述水溶相運移數學模型和所述浮力運移數學模型耦合，建立油藏烴類微滲漏數學模型的步驟包括：

獲取地層物理吸附量、碳酸鹽礦物化學吸附量和微生物降解量；

基于已確定的所述油藏地質模型，對所述擴散數學模型、所述水溶相運移數學模型和所述浮力運移數學模型耦合，結合所述地層物理吸附量、所述碳酸鹽礦物化學吸附量和所述微生物降解量，建立油藏烴類微滲漏數學模型。

3.根據權利要求2所述的方法，其特征在于，所述油藏烴類微滲漏數學模型為：

C＝C_k*α+C_s*β+C_f*γ-R₁-R₂-R₃

其中，C_k為擴散作用微滲漏的烴類濃度，C_s為水溶相運移烴類濃度，C_f為浮力運移烴類濃度，R₁為地層顆粒對烴類物質總吸附速率，R₂為碳酸巖對烴類物質總吸附速率，R₃為微生物對烴類物質的降解速率，α、β、γ為貢獻系數，且0≤α≤1，0≤β≤1，0≤γ≤1，α+β+γ＝1。

4.根據權利要求1所述的方法，其特征在于，所述擴散數學模型為：

其中，c是任一物質在直角坐標系中任一點上的濃度，D為三個方向上的擴散系數，單位為m²/s，其中，將C設為擴散作用微滲漏的烴類濃度C_k。

5.根據權利要求1所述的方法，其特征在于，所述水溶相運移數學模型為：

其中，q_ci表示烴類組分的通量，1/L；ε表示水的流速，L/T；D_T表示總分散系數；C_i表示水溶烴的濃度，表示r，1/L；

則水中烴類組分搬運公式為：

其中，將C設為水溶相運移烴類濃度C_s。

6.根據權利要求1所述的方法，其特征在于，所述獲取浮力運移數學模型的步驟為：

根據達西定律與每種相態的勢梯度成比例的關系，得到水和氣體在巖層的線速度計算式；

基于水和氣體在巖層的線速度計算式，計算得到烴類的量Q_f和濃度C_f，獲取所述浮力運移數學模型；

其中，水和氣體在巖層的線速度計算式分別為：

其中，V_g表示氣體速率；V_w表示水的速率；K表示巖層滲透率；Kg表示氣體相對滲透率；Kw表示水的相對滲透率；μg表示氣體粘滯度；μw表示水的粘滯度；Фg表示氣體勢，Фg＝P+ρgh；Фw表示水勢，Фw＝P+ρgh；

V＝Q/t，V＝Sh，Q為氣體通量，V為水的體積或者氣體的體積，S為橫截面積，h為高度，t為時間。

7.根據權利要求1所述的方法，其特征在于，所述基于所述油藏烴類微滲漏數學模型進行分段模擬步驟中的模擬內容包括烴類濃度、烴類滲漏速率和滲漏通量。