化探正演数值模拟方法、装置、计算机设备和存储介质与流程

1.本发明涉及油气勘探的数值模拟技术领域，特别涉及一种基于油藏模型的化探正演数值模拟方法、装置、计算机设备和存储介质。


背景技术：

2.油气化探技术作为油气勘探的一种手段，它的理论基础是建立在“地下油气藏中烃类气体以微弱但可检出的量近似垂直地向地表运移”的假设之上，多年来尽管取得一系列的成功案例，但地质界仍有许多学者对近地表油气化探异常形成机理存有质疑，亟需加强化探异常形成机理正演研究，其中数值模拟是开展正演研究的方法之一，可以从机理上描述烃类垂向微渗漏的过程及其近地表显示，夯实烃类垂向微渗漏的基础理论。
3.目前对于化探正演数值模拟研究，从油藏-直接盖层-上覆地层-地表基本上采用单一机制进行模拟。阮天健(1985)假设了6种盖层与气源组合状况，对扩散作用过程在地表形成的异常进行了数值模拟，模拟结果说明了小油田(窄小的气源)出现顶部异常，而大油田(大范围的气源)呈晕圈异常。肖伟、鲍征宇等(2003)在前人烃类物质微渗漏机制研究基础上，建立了烃类物质微渗漏的概念模型，并且推导了烃类物质垂向运移的动力学模型。在动力学模型的基础上，通过数值模拟认为烃类物质胶体气泡上升机制是较为合理的观点，并计算了不同情况下烃类物质垂向运移的平衡时间。黄志龙等(2007)认为轻烃呈胶束状态的微气泡形式发生准垂直向上的渗漏运移，建立了微渗漏散失的定量模型，并通过实际气藏定量估算了地史时期天然气通过盖层的微渗漏散失量。李萌(2009,2012),李志炜(2012)基于烃类垂向微渗漏基本概念和连续介质假设，构建理想单层和多层积木块地层介质模型；然后依据质量守恒定律，建立描述地层介质模型内部烃类垂向微渗漏过程的定量方程，即一类反应-对流-扩散偏微分方程。国外方面，对于扩散机制，krooss(1992)等进一步考虑了盖层中气相甲烷和水溶相甲烷的平衡，对前人的方程进一步修正，并计算了烃类在沉积柱中的扩散量并与实测结果进行对比，获得了较为吻合的结果。nelson和simmons对krooss(1992)的计算提出了批评，把扩散系数和孔隙度、渗透率、弯曲度结合，重新进行了数值计算，另外，jakel和klusman(1995),thomas和clouse(1990),都对扩散机制的数值模拟进行过研究[81]。ronald w.klusman(2005，2010)对浮力机制作用下油气藏上方饱和带和不饱和带中的烃类垂向运移过程进行了数值模拟，他认为微气泡的浮力上升机制是一种有活力的油气运移机制，为油气勘探中所观察到的地表地球化学现象提供支撑。对地表化探异常的形成，不同的学者在数值模拟时强调他们认为的主要机制，考虑相对单一。微渗漏在不同的层段，温压条件和地层水条件不同，其渗漏方式随着条件的变化而变化。同时盖层的封闭性能也决定了微渗漏的方式和规模。因此，从油藏-直接盖层-上覆地层-地表的整个地层介质系统中，必然存在各种微渗漏方式，以及影响因素。
[0004]
上述数值模拟研究，只强调了微渗漏的一种方式，地质条件下的影响因素考虑的也较少。因此，需要建立一种考虑各种微渗漏方式及其贡献，以及地质影响因素的数值模拟方法，来正演油气化探异常的形成机理。


技术实现要素：

[0005]
基于此，有必要针对上述技术问题，提供一种化探正演数值模拟方法、装置、计算机设备和存储介质。
[0006]
一种化探正演数值模拟方法，包括：
[0007]
获取油气地质数据和地震数据，确定油藏地质模型；
[0008]
获取扩散数学模型、水溶相运移数学模型和浮力运移数学模型；
[0009]
基于已确定的所述油藏地质模型，对所述扩散数学模型、所述水溶相运移数学模型和所述浮力运移数学模型耦合，建立油藏烃类微渗漏数学模型；
[0010]
基于所述油藏烃类微渗漏数学模型进行分段模拟。
[0011]
在一个实施例中，所述基于已确定的所述油藏地质模型，对所述扩散数学模型、所述水溶相运移数学模型和所述浮力运移数学模型耦合，建立油藏烃类微渗漏数学模型的步骤包括：
[0012]
获取地层物理吸附量、碳酸盐矿物化学吸附量和微生物降解量；
[0013]
基于已确定的所述油藏地质模型，对所述扩散数学模型、所述水溶相运移数学模型和所述浮力运移数学模型耦合，结合所述地层物理吸附量、所述碳酸盐矿物化学吸附量和所述微生物降解量，建立油藏烃类微渗漏数学模型。
[0014]
在一个实施例中，所述油藏烃类微渗漏数学模型为：
[0015]
c＝ck*α+cs*β+cf*γ-r
1-r
2-r3[0016]
其中，ck为扩散作用微渗漏的烃类浓度，cs为水溶相运移烃类浓度，cf为浮力运移烃类浓度，r1为地层颗粒对烃类物质总吸附速率，r2为碳酸岩对烃类物质总吸附速率，r3为微生物对烃类物质的降解速率，α、β、γ为贡献系数，且0≤α≤1，0≤β≤1，0≤γ≤1，α+β+γ＝1。
[0017]
在一个实施例中，所述扩散数学模型为：
[0018][0019]
其中，c是任一物质在直角坐标系中任一点上的浓度，d为三个方向上的扩散系数，单位为m2/s，其中，将c设为扩散作用微渗漏的烃类浓度ck。
[0020]
在一个实施例中，所述水溶相运移数学模型为：
[0021][0022]
其中，q
ci
表示烃类组分的通量，1/l；ε表示水的流速，l/t；d
t
表示总分散系数；ci表示水溶烃的浓度，表示1/l；
[0023]
则水中烃类组分搬运公式为：
[0024][0025]
，其中，将c设为水溶相运移烃类浓度cs。
[0026]
在一个实施例中，所述获取浮力运移数学模型的步骤为：
[0027]
根据达西定律与每种相态的势梯度成比例的关系，得到水和气体在岩层的线速度计算式；
[0028]
基于水和气体在岩层的线速度计算式，计算得到烃类的量qf和浓度cf，获取所述浮力运移数学模型；
[0029]
其中，水和气体在岩层的线速度计算式分别为：
[0030][0031][0032]
其中，vg表示气体速率；vw表示水的速率；k表示岩层渗透率；kg表示气体相对渗透率；kw表示水的相对渗透率；μg表示气体粘滞度；μw表示水的粘滞度；фg表示气体势，фg＝p+ρgh；фw表示水势，фw＝p+ρgh；
[0033]
v＝q/t，v＝sh，q为气体通量，v为水的体积或者气体的体积，s为横截面积，h为高度，t为时间。
[0034]
在一个实施例中，所述基于所述油藏烃类微渗漏数学模型进行分段模拟步骤中的模拟内容包括烃类浓度、烃类渗漏速率和渗漏通量。
[0035]
一种基于油藏模型的化探正演数值模拟装置，包括：
[0036]
油藏地质模型获取模块，用于获取油气地质数据和地震数据，确定油藏地质模型；
[0037]
数学模型获取模块，用于获取扩散数学模型、水溶相运移数学模型和浮力运移数学模型；
[0038]
模型耦合模块，用于基于已确定的所述油藏地质模型，对所述扩散数学模型、所述水溶相运移数学模型和所述浮力运移数学模型耦合，建立油藏烃类微渗漏数学模型；
[0039]
分段模拟模块，用于基于所述油藏烃类微渗漏数学模型进行分段模拟。
[0040]
一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现以下步骤：
[0041]
获取油气地质数据和地震数据，确定油藏地质模型；
[0042]
获取扩散数学模型、水溶相运移数学模型和浮力运移数学模型；
[0043]
基于已确定的所述油藏地质模型，对所述扩散数学模型、所述水溶相运移数学模型和所述浮力运移数学模型耦合，建立油藏烃类微渗漏数学模型；
[0044]
基于所述油藏烃类微渗漏数学模型进行分段模拟。
[0045]
一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
[0046]
获取油气地质数据和地震数据，确定油藏地质模型；
[0047]
获取扩散数学模型、水溶相运移数学模型和浮力运移数学模型；
[0048]
基于已确定的所述油藏地质模型，对所述扩散数学模型、所述水溶相运移数学模型和所述浮力运移数学模型耦合，建立油藏烃类微渗漏数学模型；
[0049]
基于所述油藏烃类微渗漏数学模型进行分段模拟。
[0068]
其中，ck为扩散作用微渗漏的烃类浓度，cs为水溶相运移烃类浓度，cf为浮力运移烃类浓度，r1为地层颗粒对烃类物质总吸附速率，r2为碳酸岩对烃类物质总吸附速率，r3为微生物对烃类物质的降解速率，α、β、γ为贡献系数，且0≤α≤1，0≤β≤1，0≤γ≤1，α+β+γ＝1。
[0069]
在一个实施例中，所述扩散数学模型为：
[0070][0071]
其中，c是任一物质在直角坐标系中任一点上的浓度，单位mg/l，d为三个方向上的扩散系数，单位为m2/s，其中，将c设为扩散作用微渗漏的烃类浓度ck。
[0072]
在一个实施例中，所述水溶相运移数学模型为：
[0073][0074]
上述计算式(2)为采用水溶相单个烃类分子因对流和分散而形成的通量得到(bear，1979)，其中，q
ci
表示烃类组分的通量，1/l；ε表示水的流速，l/t；d
t
表示总分散系数；ci表示水溶烃的浓度，表示1/l；
[0075]
在缺乏化学反应与吸附作用(或吸附点已饱和)的情况下，水溶烃的质量平衡方程可表述如下：
[0076][0077]
将等式(2)代入等式(3)，则得到水中烃类组分搬运公式，其中，将c设为水溶相运移烃类浓度cs。
[0078]
在一个实施例中，所述获取浮力运移数学模型的步骤为：根据达西定律与每种相态的势梯度成比例的关系，得到水和气体在岩层的线速度计算式；基于水和气体在岩层的线速度计算式，计算得到烃类的量qf和浓度cf，获取所述浮力运移数学模型。
[0079]
应该理解的是，不论是对于气体还是对于水，达西定律与每种相态的势梯度成比例，可以表达如下(据saeed，1991)：
[0080]
水和气体在岩层的线速度计算式分别为：
[0081][0082][0083]
其中，vg表示气体速率；vw表示水的速率；k表示岩层渗透率；kg表示气体相对渗透率；kw表示水的相对渗透率；μg表示气体粘滞度；μw表示水的粘滞度；фg表示气体势，фg＝p+ρgh；фw表示水势，фw＝p+ρgh；
[0084][0085]
v＝q/t，v＝sh，q为气体通量，v为水的体积或者气体的体积，s为横截面积，h为高度，t为时间。等式(4)和(5)分别表示水和气体的线速度。如气体通量q＝v/t，对于圆柱体而言，体积v＝sh，代入至上述通量等式，h/t即为水和气体的线速度，等于q/s。由于浓度cf与烃类的量qf有关，qf与线速度有关，因此，可得到浓度cf也是线速度有关。
[0086]
气体的相对渗透率取决于水的饱和度。水的饱和度降低，气体的饱和度就会增加，气体相对渗透率就会加大。然后根据油气组分、气和水的饱和度就可以求出烃类的量qf和浓度cf。
[0087]
除了上述烃类微渗漏的三种数学模型，还考虑了以下因素对微渗漏烃类浓度的影响及其数学模型：
[0088]
(一)地层的物理吸附量：
[0089]
采用指数吸附模式，经过数据拟合得到地层吸附量v(cm3/g)与地层的压力p(mpa)经验的公式：
[0090]
v＝2.24p
1/1
·
30
ꢀꢀꢀꢀꢀꢀ
(6)
[0091]
应该理解的是，上式(6)只是表示总的吸附量，而没有考虑吸附率的问题。为此地层颗粒对某种烃类物质总吸附速率r1与介质的孔隙度、体比表面积值、以及岩性有关。可假设为与该烃类物质在地层中的游离烃浓度c成正比、而与吸附在地层颗粒表面的烃浓度cr成反比，比例系数为－kr(称为单位时间单位长度(体积)是土壤对该烃类物质的吸附率，1/mt，负号是因为烃类物质的浓度变小)，当吸附浓度接近吸附量v时中止吸附过程，若为地层介质的孔隙度，于是：
[0092]
r1＝-krψc(v-cr)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(7)
[0093]
(二)碳酸盐矿物对烃类组分的化学吸附作用：
[0094]
烃类垂向微渗漏的过程中，一部分烃类溶解于水，水岩相互作用导致烃类进入碳酸盐矿物晶格，以化学吸附态存在于地层中。当温度压力相同时，砂、泥岩的吸附量具有一定的比例，在得到某种条件下砂岩背景后，即可求得相同条件下泥岩酸解烃背景。为获得砂、泥岩间相同条件下吸附烃的比例，采用freundlich公式对文献数据拟合得(王振平等，1996)：：
[0095]q1s
＝4.86644
×
10-3
p
0.938423
，q
1m
＝1.38999
×
10-2
p
0.940838
[0096]q3s
＝0.213835p
0.838807
，q
3m
＝0.400188p
0.842807
[0097]q4s
＝0.87435p
0.74826
，q
4m
＝1.36484p
0.738051
ꢀꢀ
(8)
[0098]
上列式(8)中q
1s
，q
1m
，q
3s
，q
3m
，q
4s
，q
4m
，，分别为砂岩、泥岩对甲烷、丙烷及丁烷的吸附量(cm3/100g)，为压力(
×
133.3pa)，各式中的数值为回归系数。
[0099]
据上式，可以得到泥岩、砂岩对同种烃类吸附量的比值：
[0100]
[0101]
式(9)中a1，a3，a4分别为相同条件下泥岩与砂岩对甲烷、丙烷、丁烷吸附量之比。
[0102]
为此碳酸岩对某种烃类物质总吸附速率r2可假设与该类物质在地层中的游离烃浓度c与吸附在碳酸演岩中的cc的差成正比－kc(称为单位时间单位长度(体积)是碳酸盐岩对该烃类物质的吸附率，1/mt，负号是因为烃类物质的浓度变小)，为地层介质的孔隙度，于是：
[0103]
r2＝-ψkcc(q-cc)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(10)
[0104]
(三)地表微生物的作用：
[0105]
地表氧化带的微生物(烃类氧化菌)活动较为强烈，垂向迁移的烃类通过时浓度会大大降低；同时微生物的活动受气候的影响较大。因此在发明中，简单地假设微生物对某个烃类物质的降解速率r3为该烃类物质的一次函数，即：
[0106]
r3＝-ψk1c
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(11)
[0107]
负号是因为烃类物质的浓度变小，k1为单位时间单位长度(体积)里微生物对该烃类物质的消耗率，1/mt。
[0108]
因此，将三种机制耦合后，并考虑地层物理吸附量、碳酸盐矿物化学吸附量、微生物降解量，耦合后的数学模型任一点微渗漏烃类浓度c可表示为：
[0109]
c＝ck*α+cs*β+cf*γ-r1-r2-r3
[0110]
其中α(0≤α≤1)、β(0≤β≤1)、γ(0≤γ≤1)为贡献系数，且α+β+γ＝1
[0111]
步骤140，基于所述油藏烃类微渗漏数学模型进行分段模拟。
[0112]
具体地，分段模拟就是由于油气藏上覆地层封盖条件和通道发育不同，烃类在垂向上微渗漏会以不同地层条件中以不同运移方式为主，致密(孔隙度小、渗透率小)、无裂隙发育岩层以扩散方式为主，比如，采用扩散运移方程；孔隙度大、渗透率大的岩层以水溶方式为主，比如，采用水溶运移方程；裂隙发育的岩层以微气泡方式为主，比如，采用浮力运移方程。
[0113]
基于上述数学模型编制计算程序，不仅可对烃类垂向微渗漏单种机制进行数值模拟，还可对三种机制进行耦合数值模拟。模拟方法中将扩散方式、水溶方式、微气泡方式的数学模型耦合成一个模型，不同地层中主要烃类微渗漏方式可能不同，耦合方程中贡献系数的选取可能就不同。采取分段模拟：一、可排除主要微渗漏方式不同带来的影响，可以随意截取一个地层进行模拟，将下一层的最终状态作为紧邻上一地层的初始状态进行模拟；二、软件能将任一实测点或者模拟点(插值点)作为初始条件，进行模拟计算。
[0114]
本技术还能够实现一维、二维和三维模拟功能。数值模拟方法可进行一维、二维和三维地层空间模拟。模拟内容包括烃类浓度、烃类渗漏速率和渗漏通量。具体的模拟，根据不同的运移方式的数学模型进行计算，不同时间、距离的浓度大小，然后将这些浓度数值变化用图表达出来。
[0115]
此外，本技术还具有模拟结果可视化功能，数值模拟方法不仅给出数值结果，还能将模拟结果图像化：
①
一维模拟，任意时间点地层中或者地表测线土壤中烃类组分浓度或比值分布图、地层中某一点浓度或比值随时间变化图；
②
二维模拟，任意时间点垂向剖面或者平面上烃类组分浓度或比值分布图、地层中某一点浓度或比值随时间变化图；
③
三维模拟，任意线或者切面或者体的浓度图或比值图，以等值面形成显示、空间地层中某一点浓度或比值随时间变化图。三维等值面显示的优点是能够看出三维空间中的优势通道展布；二
维平面数据的三维显示，只能显示优势通道在平面上的展布。
[0116]
在一个实施例中，所述基于所述油藏烃类微渗漏数学模型进行分段模拟步骤中的模拟内容包括烃类浓度、烃类渗漏速率和渗漏通量。
[0117]
实施例二
[0118]
本实施例中，请结合图4，提供一种基于油藏模型的化探正演数值模拟方法，包括：
[0119]
步骤一、确定已知油藏地质模型。
[0120]
由油气地质、地震资料，明确已知油藏埋深、油气组分、盖层厚度、含油面积、油藏剖面、多口井全井地质录井岩性剖面，统计各层砂泥岩厚度等。
[0121]
步骤二、建立油藏烃类微渗漏数学模型，包括扩散数学模型、水溶相运移数学模型、浮力运移数学模型，以及三种机制耦合的数学模型。
[0122]
本实施例中，各数学模型可同时加载地层的物理吸附量、碳酸盐矿物化学吸附量、地表微生物消耗量等。根据数学模型，编制相应计算机计算程序，对油藏油气组份通过盖层及上覆地层、第四系沉积层微渗漏的浓度进行计算。
[0123]
扩散作用服从费克定律，在直角坐标系中，其表达式为：
[0124][0125]
其中，c是任一物质在直角坐标系中任一点上的浓度，单位mg/l，d为三个方向上的扩散系数，单位m2/s。
[0126]
将扩散作用微渗漏的烃类浓度设为ck。
[0127]
水溶相单个烃类分子因对流和分散而形成的通量可采用下式计算(bear，1979)：
[0128][0129]
式中，q
ci
表示烃类组分的通量，1/l；ε表示水的流速，l/t；dt表示总分散系数；ci表示水溶烃的浓度。表示1/l。
[0130]
在缺乏化学反应与吸附作用(或吸附点已饱和)的情况下，水溶烃的质量平衡方程可表述如下：
[0131]
将等式(2)代入等式(3)，可得如下水中烃类组分搬运公式：
[0132][0133]
将水溶相运移烃类浓度设为cs。
[0134]
浮力运移，不论是对于气体还是对于水，达西定律与每种相态的势梯度成比例，可以表达如下(据saeed，1991)：
[0135]
[0136]
相对渗透率；kw表示水的相对渗透率；μg表示气体粘滞度；μw表示水的粘滞度；фg表示气体势，фg＝p+ρgh；фw表示水势，фw＝p+ρgh；
[0137][0138]
气体的相对渗透率取决于水的饱和度。水的饱和度降低，气体的饱和度就会增加，气体相对渗透率就会加大。然后根据油气组分、气和水的饱和度就可以求出烃类的量qf和浓度cf。
[0139]
除了上述烃类微渗漏的三种数学模型，还考虑了以下因素对微渗漏烃类浓度的影响及其数学模型：
[0140]
(一)地层的物理吸附量：
[0141]
采用指数吸附模式，经过数据拟合得到地层吸附量v(cm3/g)与地层的压力p(mpa)经验的公式：
[0142]
v＝2.24p
1/1
·
30
ꢀꢀ
(6)
[0143]
应该注意的是，上式只是表示总的吸附量，而没有考虑吸附率的问题。为此地层颗粒对某种烃类物质总吸附速率r1与介质的孔隙度、体比表面积值、以及岩性有关。可假设为与该烃类物质在地层中的游离烃浓度c成正比、而与吸附在地层颗粒表面的烃浓度cr成反比，比例系数为－kr(称为单位时间单位长度(体积)是土壤对该烃类物质的吸附率，1/mt，负号是因为烃类物质的浓度变小)，当吸附浓度接近吸附量v时中止吸附过程，若为地层介质的孔隙度，于是：
[0144]
r1＝-krψc(v-cr)
ꢀꢀ
(7)
[0145]
(二)碳酸盐矿物对烃类组分的化学吸附作用：
[0146]
烃类垂向微渗漏的过程中，一部分烃类溶解于水，水岩相互作用导致烃类进入碳酸盐矿物晶格，以化学吸附态存在于地层中。当温度压力相同时，砂、泥岩的吸附量具有一定的比例，在得到某种条件下砂岩背景后，即可求得相同条件下泥岩酸解烃背景。为获得砂、泥岩间相同条件下吸附烃的比例，采用freundlich公式对文献数据拟合得(王振平等，1996)：：
[0147]q1s
＝4.86644
×
10-3
p
0.938423q1m
＝1.38999
×
10-2
p
0.940838
[0148]q3s
＝0.213835p
0.838807q3m
＝0.400188p
0.842807
ꢀꢀ
(8)
[0149]q4s
＝0.87435p
0.74826q4m
＝1.36484p
0.738051
[0150]
上列各式中q
1s
，q
1m
，q
3s
，q
3m
，q
4s
，q
4m
，分别为砂岩、泥岩对甲烷、丙烷及丁烷的吸附量(cm3/100g)，p为压力(
×
133.3pa)，各式中的数值为回归系数。
[0151]
据上式，可以得到泥岩、砂岩对同种烃类吸附量的比值：
[0152][0153]
式中a1，a3，a4分别为相同条件下泥岩与砂岩对甲烷、丙烷、丁烷吸附量之比。
[0154]
为此碳酸岩对某种烃类物质总吸附速率r2可假设与该类物质在地层中的游离烃
[0176]
其中，ck为扩散作用微渗漏的烃类浓度，cs为水溶相运移烃类浓度，cf为浮力运移烃类浓度，r1为地层颗粒对烃类物质总吸附速率，r2为碳酸岩对烃类物质总吸附速率，r3为微生物对烃类物质的降解速率，α、β、γ为贡献系数，且0≤α≤1，0≤β≤1，0≤γ≤1，α+β+γ＝1。
[0177]
在一个实施例中，所述扩散数学模型为：
[0178][0179]
其中，c是任一物质在直角坐标系中任一点上的浓度，d为三个方向上的扩散系数，单位为m2/s，其中，将c设为扩散作用微渗漏的烃类浓度ck。
[0180]
在一个实施例中，所述水溶相运移数学模型为：
[0181][0182]
其中，q
ci
表示烃类组分的通量，1/l；ε表示水的流速，l/t；d
t
表示总分散系数；ci表示水溶烃的浓度，表示1/l；
[0183]
则水中烃类组分搬运公式为：
[0184][0185]
，其中，将c设为水溶相运移烃类浓度cs。
[0186]
在一个实施例中，所述获取浮力运移数学模型的步骤为：
[0187]
根据达西定律与每种相态的势梯度成比例的关系，得到水和气体在岩层的线速度计算式；
[0188]
基于水和气体在岩层的线速度计算式，计算得到烃类的量qf和浓度cf，获取所述浮力运移数学模型；
[0189]
其中，水和气体在岩层的线速度计算式分别为：
[0190][0191][0192]
其中，vg表示气体速率；vw表示水的速率；k表示岩层渗透率；kg表示气体相对渗透率；kw表示水的相对渗透率；μg表示气体粘滞度；μw表示水的粘滞度；фg表示气体势，фg＝p+ρgh；фw表示水势，фw＝p+ρgh；
[0193][0194]
v＝q/t，v＝sh，q为气体通量，v为水的体积或者气体的体积，s为横截面积，h为高度，t为时间。
[0195]
在一个实施例中，所述基于所述油藏烃类微渗漏数学模型进行分段模拟步骤中的模拟内容包括烃类浓度、烃类渗漏速率和渗漏通量。
[0196]
关于化探正演数值模拟装置的具体限定可以参见上文中对于化探正演数值模拟方法的限定，在此不再赘述。上述化探正演数值模拟装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。
[0197]
实施例四
[0198]
提供了计算机设备，其内部结构图可以如图3所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与其他计算机设备连接通信。该计算机程序被处理器执行时以实现一种化探正演数值模拟方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
[0199]
本领域技术人员可以理解，图3中示出的结构，仅仅是与本技术方案相关的部分结构的框图，并不构成对本技术方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
[0200]
实施例五
[0201]
提供了一种计算机设备，包括存储器和处理器，该存储器存储有计算机程序，该处理器执行计算机程序时实现以下步骤：
[0202]
获取油气地质数据和地震数据，确定油藏地质模型；
[0203]
获取扩散数学模型、水溶相运移数学模型和浮力运移数学模型；
[0204]
基于已确定的所述油藏地质模型，对所述扩散数学模型、所述水溶相运移数学模型和所述浮力运移数学模型耦合，建立油藏烃类微渗漏数学模型；
[0205]
基于所述油藏烃类微渗漏数学模型进行分段模拟。
[0206]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
[0207]
获取地层物理吸附量、碳酸盐矿物化学吸附量和微生物降解量；
[0208]
基于已确定的所述油藏地质模型，对所述扩散数学模型、所述水溶相运移数学模型和所述浮力运移数学模型耦合，结合所述地层物理吸附量、所述碳酸盐矿物化学吸附量和所述微生物降解量，建立油藏烃类微渗漏数学模型。
[0209]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：获取油藏烃类微渗漏数学模型，所述油藏烃类微渗漏数学模型为：
[0210]
c＝ck*α+cs*β+cf*γ-r
1-r
2-r3[0211]
其中，ck为扩散作用微渗漏的烃类浓度，cs为水溶相运移烃类浓度，cf为浮力运移烃类浓度，r1为地层颗粒对烃类物质总吸附速率，r2为碳酸岩对烃类物质总吸附速率，r3为微生物对烃类物质的降解速率，α、β、γ为贡献系数，且0≤α≤1，0≤β≤1，0≤γ≤1，α+β+γ
＝1。
[0212]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：获取扩散数学模型，所述扩散数学模型为：
[0213][0214]
其中，c是任一物质在直角坐标系中任一点上的浓度，d为三个方向上的扩散系数，单位为m2/s，其中，将c设为扩散作用微渗漏的烃类浓度ck。
[0215]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：获取水溶相运移数学模型，所述水溶相运移数学模型为：
[0216][0217]
其中，q
ci
表示烃类组分的通量，1/l；ε表示水的流速，l/t；d
t
表示总分散系数；ci表示水溶烃的浓度，表示1/l；
[0218]
则水中烃类组分搬运公式为：
[0219][0220]
，其中，将c设为水溶相运移烃类浓度cs。
[0221]
在一个实施例中，处理器执行计算机程序时还实现以下步骤：
[0222]
根据达西定律与每种相态的势梯度成比例的关系，得到水和气体在岩层的线速度计算式；
[0223]
基于水和气体在岩层的线速度计算式，计算得到烃类的量qf和浓度cf，获取所述浮力运移数学模型；
[0224]
其中，水和气体在岩层的线速度计算式分别为：
[0225][0226][0227]
其中，vg表示气体速率；vw表示水的速率；k表示岩层渗透率；kg表示气体相对渗透率；kw表示水的相对渗透率；μg表示气体粘滞度；μw表示水的粘滞度；фg表示气体势，фg＝p+ρgh；фw表示水势，фw＝p+ρgh；
[0228][0229]
v＝q/t，v＝sh，q为气体通量，v为水的体积或者气体的体积，s为横截面积，h为高度，t为时间。
[0230]
实施例六
[0231]
提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：
[0232]
获取油气地质数据和地震数据，确定油藏地质模型；
[0233]
获取扩散数学模型、水溶相运移数学模型和浮力运移数学模型；
[0234]
基于已确定的所述油藏地质模型，对所述扩散数学模型、所述水溶相运移数学模型和所述浮力运移数学模型耦合，建立油藏烃类微渗漏数学模型；
[0235]
基于所述油藏烃类微渗漏数学模型进行分段模拟。
[0236]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
[0237]
获取地层物理吸附量、碳酸盐矿物化学吸附量和微生物降解量；
[0238]
基于已确定的所述油藏地质模型，对所述扩散数学模型、所述水溶相运移数学模型和所述浮力运移数学模型耦合，结合所述地层物理吸附量、所述碳酸盐矿物化学吸附量和所述微生物降解量，建立油藏烃类微渗漏数学模型。
[0239]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：获取油藏烃类微渗漏数学模型，所述油藏烃类微渗漏数学模型为：
[0240]
c＝ck*α+cs*β+cf*γ-r
1-r
2-r3[0241]
其中，ck为扩散作用微渗漏的烃类浓度，cs为水溶相运移烃类浓度，cf为浮力运移烃类浓度，r1为地层颗粒对烃类物质总吸附速率，r2为碳酸岩对烃类物质总吸附速率，r3为微生物对烃类物质的降解速率，α、β、γ为贡献系数，且0≤α≤1，0≤β≤1，0≤γ≤1，α+β+γ＝1。
[0242]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：获取扩散数学模型，所述扩散数学模型为：
[0243][0244]
其中，c是任一物质在直角坐标系中任一点上的浓度，d为三个方向上的扩散系数，单位为m2/s，其中，将c设为扩散作用微渗漏的烃类浓度ck。
[0245]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：获取水溶相运移数学模型，所述水溶相运移数学模型为：
[0246][0247]
其中，q
ci
表示烃类组分的通量，1/l；ε表示水的流速，l/t；d
t
表示总分散系数；ci表示水溶烃的浓度，表示1/l；
[0248]
则水中烃类组分搬运公式为：
[0249][0250]
，其中，将c设为水溶相运移烃类浓度cs。
[0251]
在一个实施例中，计算机程序被处理器执行时还实现以下步骤：
[0252]
根据达西定律与每种相态的势梯度成比例的关系，得到水和气体在岩层的线速度
计算式；
[0253]
基于水和气体在岩层的线速度计算式，计算得到烃类的量qf和浓度cf，获取所述浮力运移数学模型；
[0254]
其中，水和气体在岩层的线速度计算式分别为：
[0255][0256][0257]
其中，vg表示气体速率；vw表示水的速率；k表示岩层渗透率；kg表示气体相对渗透率；kw表示水的相对渗透率；μg表示气体粘滞度；μw表示水的粘滞度；фg表示气体势，фg＝p+ρgh；фw表示水势，фw＝p+ρgh；
[0258][0259]
v＝q/t，v＝sh，q为气体通量，v为水的体积或者气体的体积，s为横截面积，h为高度，t为时间。
[0260]
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本技术所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(rom)、可编程rom(prom)、电可编程rom(eprom)、电可擦除可编程rom(eeprom)或闪存。易失性存储器可包括随机存取存储器(ram)或者外部高速缓冲存储器。作为说明而非局限，ram以多种形式可得，诸如静态ram(sram)、动态ram(dram)、同步dram(sdram)、双数据率sdram(ddrsdram)、增强型sdram(esdram)、同步链路(synchlink)dram(sldram)、存储器总线(rambus)直接ram(rdram)、直接存储器总线动态ram(drdram)、以及存储器总线动态ram(rdram)等。
[0261]
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。
[0262]
以上所述实施例仅表达了本技术的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本技术构思的前提下，还可以做出若干变形和改进，这些都属于本技术的保护范围。因此，本技术专利的保护范围应以所附权利要求为准。