[发明专利]化探正演数值模拟方法、装置、计算机设备和存储介质

说明书

本发明提供一种化探正演数值模拟方法、装置、计算机设备和存储介质，所述方法包括获取油气地质数据和地震数据，确定油藏地质模型；获取扩散数学模型、水溶相运移数学模型和浮力运移数学模型；基于已确定的所述油藏地质模型，对所述扩散数学模型、所述水溶相运移数学模型和所述浮力运移数学模型耦合，建立油藏烃类微渗漏数学模型；基于所述油藏烃类微渗漏数学模型进行分段模拟。采用数值模拟方法较好地反映实际地质条件下微渗漏三种方式的转变及其微渗漏浓度变化，考虑了多种地质因素的影响，能实现一维、二维、三维地层空间的数值模拟及其可视化效果，正演油气化探形成机理，夯实了烃类垂向微渗漏理论。

技术领域

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

背景技术

油气化探技术作为油气勘探的一种手段，它的理论基础是建立在“地下油气藏中烃类气体以微弱但可检出的量近似垂直地向地表运移”的假设之上，多年来尽管取得一系列的成功案例，但地质界仍有许多学者对近地表油气化探异常形成机理存有质疑，亟需加强化探异常形成机理正演研究，其中数值模拟是开展正演研究的方法之一，可以从机理上描述烃类垂向微渗漏的过程及其近地表显示，夯实烃类垂向微渗漏的基础理论。

目前对于化探正演数值模拟研究，从油藏-直接盖层-上覆地层-地表基本上采用单一机制进行模拟。阮天健(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)对浮力机制作用下油气藏上方饱和带和不饱和带中的烃类垂向运移过程进行了数值模拟，他认为微气泡的浮力上升机制是一种有活力的油气运移机制，为油气勘探中所观察到的地表地球化学现象提供支撑。对地表化探异常的形成，不同的学者在数值模拟时强调他们认为的主要机制，考虑相对单一。微渗漏在不同的层段，温压条件和地层水条件不同，其渗漏方式随着条件的变化而变化。同时盖层的封闭性能也决定了微渗漏的方式和规模。因此，从油藏-直接盖层-上覆地层-地表的整个地层介质系统中，必然存在各种微渗漏方式，以及影响因素。

上述数值模拟研究，只强调了微渗漏的一种方式，地质条件下的影响因素考虑的也较少。因此，需要建立一种考虑各种微渗漏方式及其贡献，以及地质影响因素的数值模拟方法，来正演油气化探异常的形成机理。

发明内容

基于此，有必要针对上述技术问题，提供一种化探正演数值模拟方法、装置、计算机设备和存储介质。

一种化探正演数值模拟方法，包括：

获取油气地质数据和地震数据，确定油藏地质模型；

获取扩散数学模型、水溶相运移数学模型和浮力运移数学模型；