本申请公开了一种弹性参数确定方法、装置、计算机设备及存储介质，涉及油气田勘探领域，方法包括：对于勘探工区内每个地震道集，根据地震道集的位置，确定地震道集的初始弹性参数；根据勘探工区的测井资料和初始弹性参数，确定参数取值集合；将地震道集内目的层的一个波形周期的振幅值的总和，确定为降温的初始温度；基于初始弹性参数、参数取值集合、初始温度和降温函数，对地震道集进行参数迭代，得到地震道集对应的弹性参数。通过上述方案使得迭代的初始值与真实值更接近，大量减少了扰动个数，从而提高迭代效率，并且，由于扰动范围减小，使得降温函数与弹性参数的参数取值范围匹配，进一步降低了迭代次数，提高了迭代效率。提高了迭代效率。提高了迭代效率。
【技术实现步骤摘要】
弹性参数确定方法、装置、计算机设备及存储介质
[0001]本申请涉及油气田勘探领域，特别涉及一种弹性参数确定方法、装置、计算机设备及存储介质。
技术介绍
[0002]模拟退火算法(simulated annealing)也称蒙特卡罗(Monte Carlo)算法，因其独特的优势，在地震勘探的数据处理和地震反演中的应用越来越多。模拟退火算法是把全局寻优问题与物体在降温过程中结晶的热平衡问题进行类比得来的。在石油勘探中的应用也越来越多。例如，通过模拟退火算法，在石油勘探过程中，进行AVO(Amplitude variation with offset,振幅随偏移距的变化)弹性参数反演，以识别储层的岩性和含流体性质。
[0003]相关技术中，在给定的温度控制下，根据目标函数变化，修正退火的控制温度，随着温度的降低，随机从概率分布中随机确定地震道集中的一个弹性参数作为目标函数的解，通过遍历多个弹性参数，逐渐使目标函数向减小的方向搜寻，最终得到目标函数最小时对应的目标弹性参数，将该目标弹性参数确定为该地震道集的弹性参数。
[0004]虽然通过寻优算法不必遍历所有可能的解，使其效率大幅度提高，但在地震反演中由于基础数据量大，特别是近年来高密度、宽方位地震勘探方法的实施，叠前地震反演的数据量急剧增加，对地震反演的时效性影响更大，反演的效率仍然较低。
技术实现思路
[0005]本申请提供了一种弹性参数确定方法、装置、计算机设备及存储介质。所述技术方案如下：
[0006]根据本申请的一方面，提供了一种弹性参数确定方法，所述方法包括：
[0007]对于勘探工区内每个地震道集，根据所述地震道集的位置，确定所述地震道集的初始弹性参数；
[0008]根据所述勘探工区的测井资料和所述初始弹性参数，确定参数取值集合；
[0009]将所述地震道集内目的层的一个波形周期的振幅值的总和，确定为降温的初始温度；
[0010]基于所述初始弹性参数、所述参数取值集合、所述初始温度和降温函数，对所述地震道集进行参数迭代，得到所述地震道集对应的弹性参数。
[0011]在一些实施例中，所述根据所述地震道集的位置，确定所述地震道集的初始弹性参数，包括：
[0012]响应于所述地震道集的位置与测井所在的位置重合，获取所述地震道集的位置上的测井的第一测井资料；从所述第一测井资料中提取第一弹性参数，将所述第一弹性参数确定为所述初始弹性参数；或者，
[0013]响应于所述地震道集的位置与测井所在的位置不重合，获取所述地震道集的位置所在目标区域内的测井的第二测井资料，从所述第二测井资料中提取所述测井的第二弹性
参数，对所述第二弹性参数进行插值处理，得到插值处理结果，将所述插值处理结果确定为所述初始弹性参数。
[0014]在一些实施例中，所述对所述第二弹性参数进行插值处理，得到插值处理结果，包括：
[0015]响应于所述测井资料为二维地震测线对应的测井资料，采用线性内插方式进行插值处理，得到所述插值处理结果；
[0016]响应于所述测井资料为三维地震区块对应的测井资料，采用距离加权内插方式进行插值处理，得到所述插值处理结果。
[0017]在一些实施例中，所述根据所述勘探工区的测井资料和所述初始弹性参数，确定参数取值集合，包括：
[0018]根据所述勘探工区的测井资料的品质和储层的复杂程度，以初始弹性参数值为中心，向减小和增大方向，确定参数取值百分比；
[0019]基于所述初始弹性参数和所述参数取值百分比，依据最终弹性参数反演精度要求，确定弹性参数取值间隔，并得到取值个数；
[0020]基于所述取值个数，从所述参数取值范围中，确定所述取值个数的参数取值；
[0021]将所述取值个数的参数取值确定为所述参数取值集合。
[0022]在一些实施例中，所述将所述地震道集内目的层的一个波形周期的振幅值的总和，确定为降温的初始温度，包括：
[0023]基于如下公式，对所述地震道集内目的层的一个波形周期的振幅值进行求和，得到所述初始温度：
[0024][0025]其中，T0为初始温度，c为比例系数，NT为一个地震道集中的道数，t1为一个波形周期的起始时间，t2为一个波形周期的终止时间S
j
为第j个采样点的振幅值，abs( )为求绝对值函数。
[0026]在一些实施例中，所述基于所述初始弹性参数、所述参数取值集合、所述初始温度和降温函数，对所述地震道集进行参数迭代，得到所述地震道集对应的弹性参数，包括：
[0027]将所述初始弹性参数和所述初始温度作为迭代过程中算的初始参数，在迭代过程中，从所述参数取值集合中随机获取弹性参数作为目标函数的解；
[0028]响应于所述目标函数在所述弹性参数下的值不小于第一预设阈值，确定所述弹性参数对应的目标函数的增量；响应于所述目标函数的增量满足预设条件，将所述弹性参数确定为所述目标函数的当前解；其中，预设条件包括所述目标函数的增量小于第二预设阈值；或者，所述目标函数的增量不小于第二预设阈值，且所述增量与所述初始温度的比值的指数概率大于第三预设阈值；
[0029]响应于所述目标函数在所述弹性参数下的值小于所述第一预设阈值，结束迭代，将所述弹性参数确定为所述地震道集对应的弹性参数。
[0030]在一些实施例中，所述方法还包括：
[0031]响应于所述目标函数的增量不满足所述预设条件，重新从所述参数取值集合中选
取弹性参数，执行将所述弹性参数和所述初始温度作为迭代过程中算的初始参数，在迭代过程中，从所述参数取值集合中随机获取弹性参数作为目标函数的解的步骤。
[0032]根据本申请的另一方面，提供了一种弹性参数确定装置，所述装置包括：
[0033]第一确定模块，用于对于勘探工区内每个地震道集，根据所述地震道集的位置，确定所述地震道集的初始弹性参数；
[0034]第二确定模块，用于根据所述勘探工区的测井资料和所述初始弹性参数，确定参数取值集合；
[0035]第三确定模块，用于将所述地震道集内目的层的一个波形周期的振幅值的总和，确定为降温的初始温度；
[0036]迭代模块，用于基于所述初始弹性参数、所述参数取值集合、所述初始温度和降温函数，对所述地震道集进行参数迭代，得到所述地震道集对应的弹性参数。
[0037]在一些实施例中，所述第一确定模块，用于包括：
[0038]第一确定单元，用于响应于所述地震道集的位置与测井所在的位置重合，获取所述地震道集的位置上的测井的第一测井资料；从所述第一测井资料中提取第一弹性参数，将所述第一弹性参数确定为所述初始弹性参数；或者，
[0039]第二确定单元，用于响应于所述地震道集的位置与测井所在的位置不重合，获取所述地震道集的位置所在目标区域内的测井的第二测井资料，从所述第二测井资料中提取所述测井的第二弹性参数，对所述第二弹性参数进行插值处理，得到插本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种弹性参数确定方法，其特征在于，所述方法包括：对于勘探工区内每个地震道集，根据所述地震道集的位置，确定所述地震道集的初始弹性参数；根据所述勘探工区的测井资料和所述初始弹性参数，确定参数取值集合；将所述地震道集内目的层的一个波形周期的振幅值的总和，确定为降温的初始温度；基于所述初始弹性参数、所述参数取值集合、所述初始温度和降温函数，对所述地震道集进行参数迭代，得到所述地震道集对应的弹性参数。2.根据权利要求1所述的方法，其特征在于，所述根据所述地震道集的位置，确定所述地震道集的初始弹性参数，包括：响应于所述地震道集的位置与测井所在的位置重合，获取所述地震道集的位置上的测井的第一测井资料；从所述第一测井资料中提取第一弹性参数，将所述第一弹性参数确定为所述初始弹性参数；或者，响应于所述地震道集的位置与测井所在的位置不重合，获取所述地震道集的位置所在目标区域内的测井的第二测井资料，从所述第二测井资料中提取所述测井的第二弹性参数，对所述第二弹性参数进行插值处理，得到插值处理结果，将所述插值处理结果确定为所述初始弹性参数。3.根据权利要求2所述的方法，其特征在于，所述对所述第二弹性参数进行插值处理，得到插值处理结果，包括：响应于所述测井资料为二维地震测线对应的测井资料，采用线性内插方式进行插值处理，得到所述插值处理结果；响应于所述测井资料为三维地震区块对应的测井资料，采用距离加权内插方式进行插值处理，得到所述插值处理结果。4.根据权利要求1所述的方法，其特征在于，所述根据所述勘探工区的测井资料和所述初始弹性参数，确定参数取值集合，包括：根据所述勘探工区的测井资料的品质和储层的复杂程度，以初始弹性参数值为中心，向减小和增大方向确定参数取值范围的百分比；基于所述初始弹性参数和所述参数取值范围的百分比，依据最终弹性参数反演精度要求，确定弹性参数取值间隔，并得到取值个数；基于所述取值个数，从所述参数取值范围中，确定所述取值个数的参数取值；将所述取值个数的参数取值确定为所述参数取值集合。5.根据权利要求1所述的方法，其特征在于，所述将所述地震道集内目的层的一个波形周期的振幅值的总和，确定为降温的初始温度，包括：基于如下公式，对所述地震道集内目的层的一个波形周期的振幅值进行求和，得到所述初始温度：其中，T0为初始温度，c为比例系数，NT为一个地震道集中的道数，t1为一个波形周期...
【专利技术属性】
技术研发人员：张本健，郑超，王宇峰，尹宏，邓波，杨华，曹建，杨迅，裴森奇，胡欣，李容荣，赵漾，孙志昀，汤兴宇，唐明，刘柏，李洪林，罗静，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。