【技术实现步骤摘要】
本专利技术属于石油地球物理勘探，更具体地，涉及一种快速优选岩石物理建模参数的方法、装置、电子设备及介质。
技术介绍
1、岩石物理模型能够反映岩石微观物性参数(孔隙度、渗透率、含水饱和度等)与宏观弹性参数(纵波速度、横波速度等)之间的物理联系。通过岩石物理建模，可以使用测井解释结果进行地区纵、横波速度的预测，达到预测无横波测井地区的横波速度的目的。在使用岩石物理模型实施横波速度预测时，需要输入的参数有：基质的密度、模量(或者速度)、孔隙纵横比等。由于岩石的沉积环境和成岩作用不同，同一种矿物的模量在不同地区差异较大，如粘土矿物高岭石的平均体积模量为1.5gpa，而墨西哥湾粘土矿物的体积模量可以达到25gpa，因此，基质模量的准确选取往往是非常模糊的，进而导致横波速度预测偏差。相较于其他参数，矿物模量对预测横波速度的影响占很大权重。因此，针对岩石基质矿物模量的预测就成为了亟待解决的问题。
2、岩石物理建模中骨架组分较多，每种组分的参数随地区和沉积环境的变化而不同，骨架组分参数难以有效的获取，致使岩石物理建模技术和测井横波预测技术的应用受到限制。
技术实现思路
1、本专利技术的目的是针对岩石物理建模中骨架组分参数难以确定的问题，本专利技术提供一种快速优选岩石物理建模参数的方法、装置、电子设备及介质，通过采用数学最优化的方式对岩石物理建模中骨架组分的参数进行最优化选取，为基于岩石物理模型的测井横波预测提供最优骨架组分参数。
2、为了实现上述目的，本专利技术的第一方面提
3、(1)获取目标区测井数据，岩石物理模型及建模所需参数；
4、(2)基于所述目标区测井数据和岩石物理模型，利用纵波速度、横波速度联合约束，建立骨架组分参数的最优化目标函数；
5、(3)基于所述目标区测井数据，获取关于纵波速度和横波速度的骨架组分参数的数值分布范围；
6、(4)对所述关于纵波速度和横波速度的骨架组分参数的数值分布范围进行分割，得到骨架组分参数的数值序列；
7、(5)将所述骨架组分参数的数值序列利用枚举法代入所述岩石物理模型，基于所述骨架组分参数的最优化目标函数进行优化计算，得到岩石物理模型的最优化的骨架参数。
8、根据本专利技术，优选地，所述目标区测井数据包括储层物性参数曲线和储层弹性参数曲线；
9、所述储层物性参数曲线优选包括：孔隙度测井曲线、泥质含量测井曲线和含水饱和度测井曲线。
10、所述储层弹性参数曲线优选包括：密度曲线、纵波速度曲线和横波速度曲线。
11、根据本专利技术，优选地，所述建模所需参数包括：密度、纵波速度、横波速度和孔隙纵横比。
12、本专利技术中，岩石物理模型为用户给定的岩石物理模型。用户需给定岩石物理模型及其它建模所需参数，如简单组分的密度、纵横波速度、孔隙纵横比等参数。
13、根据本专利技术，优选地，所述骨架组分参数的最优化目标函数如公式(1)所示：
14、f(m1,m2,…,mn)＝min(||vpmod-vpwell||+2*||vsmod-vswell||) 公式(1)
15、其中，mi，i＝1,...,n是以向量形式表示第i个骨架组分参数，vpmod和vsmod分别表示所述岩石物理模型正演计算的纵波速度和横波速度，vpwell和vswell分别表示所述目标区测井数据中的纵波速度和横波速度，||·||表示二范数。
16、本专利技术中，横波速度项的权重系数为2，该设置既可以提高横波速度预测的精度，同时也能保证纵波速度预测精度。
17、根据本专利技术，优选地，步骤(3)中，所述关于纵波速度和横波速度的骨架组分参数为mi＝(vpi,vsi)；
18、所述关于纵波速度和横波速度的骨架组分参数的数值分布范围为mi∈[milow,miup]；
19、其中，milow＝(vpilow,vsilow)，miup＝(vpiup,vsiup)，vpilow表示基于目标区测井数据的骨架组分的纵波速度最小值，vsilow表示基于目标区测井数据的骨架组分的横波速度最小值，vpiup表示基于目标区测井数据的骨架组分的纵波速度最大值，vsiup表示基于目标区测井数据的骨架组分的横波速度最大值。
20、本专利技术中，关于纵波速度和横波速度的骨架组分参数组分的纵波速度和横波速度，即mi＝(vpi,vsi)；关于纵波速度和横波速度的骨架组分参数的数值分布范围为mi∈[milow,miup]从测井数据中获取；vpilow、vsilow、vpiup和vsiup均为已知参数。
21、根据本专利技术，优选地，步骤(4)中，按照1％的间隔对所述关于纵波速度和横波速度的骨架组分参数的数值分布范围进行分割，得到骨架组分参数的数值序列。
22、根据本专利技术，优选地，步骤(4)中，所述骨架组分参数的数值序列包括基于纵波速度的骨架组分参数的数值序列和基于横波速度的骨架组分参数的数值序列；
23、δvpi＝(vpiup-vpilow)/100，…vsi＝(vsiup-vsilow)/100，
24、vpik＝vpilow+k*δvpi，k＝1,2,3,…,100，
25、vsik＝vsilow+k*δvsi，k＝1,2,3,…,100，
26、其中，vpik表示基于纵波速度的骨架组分参数的数值序列，vsik表示基于横波速度的骨架组分参数的数值序列，k表示间隔序号。
27、本专利技术中，步骤(5)中，将骨架组分参数的数值序列vpik和vsik利用枚举法代入用户自定义的岩石物理模型，得到vpmod和vsmod，结合已知参数vpwell和vswell，根据最优化目标函数(公式(1))进行优化计算，得到最优解，即岩石物理模型的最优化的骨架参数。
28、本专利技术中，关于最优解的应用：使用步骤(5)得到的岩石物理模型的最优化的骨架参数来实施目标区无横波测井地区的横波速度预测以及定量解释量版制作等工作。
29、本专利技术的第二方面提供一种快速优选岩石物理建模参数的装置，该装置包括：
30、数据获取模块，获取目标区测井数据，岩石物理模型及建模所需参数；
31、目标函数建立模块，基于所述目标区测井数据和岩石物理模型，利用纵波速度、横波速度联合约束，建立骨架组分参数的最优化目标函数；
32、骨架参数分布范围获取模块，基于所述目标区测井数据，获取关于纵波速度和横波速度的骨架组分参数的数值分布范围；
33、骨架参数分割模块，对所述关于纵波速度和横波速度的骨架组分参数的数值分布范围进行分割，得到骨架组分参数的数值序列；
34、优化计算模块，将所述骨架组分参数的数值序列利用枚举法代入所述岩石物理模型，基于所述骨架组分参数的最优化目标函数进行优化计算，得到岩石物理模型的最优化的骨架参数。
<本文档来自技高网...【技术保护点】
1.一种快速优选岩石物理建模参数的方法，其特征在于，该方法包括：
2.根据权利要求1所述的方法，其中，所述目标区测井数据包括储层物性参数曲线和储层弹性参数曲线；
3.根据权利要求1所述的方法，其中，所述建模所需参数包括：密度、纵波速度、横波速度和孔隙纵横比。
4.根据权利要求1所述的方法，其中，所述骨架组分参数的最优化目标函数如公式(1)所示：
5.根据权利要求1所述的方法，其中，步骤(3)中，所述关于纵波速度和横波速度的骨架组分参数为mi＝(vpi,vsi)；
6.根据权利要求1所述的方法，其中，步骤(4)中，按照1％的间隔对所述关于纵波速度和横波速度的骨架组分参数的数值分布范围进行分割，得到骨架组分参数的数值序列。
7.根据权利要求5所述的方法，其中，步骤(4)中，所述骨架组分参数的数值序列包括基于纵波速度的骨架组分参数的数值序列和基于横波速度的骨架组分参数的数值序列；
8.一种快速优选岩石物理建模参数的装置，其特征在于，该装置包括：
9.一种电子设备，其特征在于，该电子设备包括：
10.一种计算机可读存储介质，其特征在于，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现权利要求1-7中任一项所述的快速优选岩石物理建模参数的方法。
...【技术特征摘要】
1.一种快速优选岩石物理建模参数的方法，其特征在于，该方法包括：
2.根据权利要求1所述的方法，其中，所述目标区测井数据包括储层物性参数曲线和储层弹性参数曲线；
3.根据权利要求1所述的方法，其中，所述建模所需参数包括：密度、纵波速度、横波速度和孔隙纵横比。
4.根据权利要求1所述的方法，其中，所述骨架组分参数的最优化目标函数如公式(1)所示：
5.根据权利要求1所述的方法，其中，步骤(3)中，所述关于纵波速度和横波速度的骨架组分参数为mi＝(vpi,vsi)；
6.根据权利要求1所述的方法，其中，步骤(4)中，按照1％的间...
【专利技术属性】
技术研发人员：李敏，王洋，刘卫华，
申请(专利权)人：中国石油化工股份有限公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。