【技术实现步骤摘要】
本专利技术涉及含钾矿藏勘探，尤其涉及一种基于井震联合约束的深度域数据体获取方法、系统及存储介质。
技术介绍
1、当前，钾盐作为钾肥的主要原料是中国最为紧缺的战略性矿产资源之一,其广泛应用于农肥、化工、医药、纺织、印染等领域，特别是大量被用于制造复合肥。目前钾盐勘探采用的常规物探技术所使用和计算获得的成果数据，如地震数据、钾含量敏感因子反演数据等，大多是时间域的，为了能够进行准确的钻探井目标设计、地层厚度计算、资源量计算等工作，需要经过转换得到相对准确的地层参数深度域数据。
2、地震类数据(地震数据、反演数据等)从时间域到深度域的数据转换，目标是为了获得地质特征与钻探井位深度能够匹配对应关系的数据体，其精度主要取决于时深转换速度模型的合理性，而速度模型的准确与否，主要取决于该模型能否准确体现时间域地震类数据与深度域钻井信息的耦合关系。
3、现阶段，由时间域数据获取深度域数据有三种比较典型的方法：1、匀速速度场转换法；2、处理速度场转换法；3、反演速度场转换法。
4、方法1匀速速度场转换法，在一个地区，利用该区已有时间-深度统一速度尺开展时深转换，将时间域地震类数据转换为深度域数据。该类方法由于没有考虑地下速度的空间变化，显而易见会在井点位置出现深度域误差，得到的深度域地震数据与测井对应关系较差。
5、方法2处理速度场转换法，考虑到地层速度在空间上存在的变化，应用地震处理过程中得到的空间速度场，建立时间和深度的对应关系，以此为依据，实现空间速度变化的深度域数据转换。该类方法虽然考虑的速
6、方法3反演速度场转换法，应用地震资料和构造层位约束，可以反演得到符合空间变化的反演速度数据，以此为依据，开展深度域数据转换求取，可以得到相对精细的空间转换关系。但是，一方面目前反演方法得到的反演速度在井点处往往与井震对应速度存在误差，另一方面，反演速度的低频分量也与空间速度稳定变化信息存在差异，因此这种方法的步骤通常比较繁琐，需要不断的迭代修正深度转换误差，最终深度域数据在井点和空间任意点处会存在一定误差。
7、综上，现有技术中的三种转换法都存在精度不足的问题，特别是在井点部位存在一定的深度误差。
技术实现思路
1、为克服以上所述技术缺陷，本专利技术实施例提供一种基于井震联合约束的深度域数据体获取方法、系统及存储介质，能够获得与井点与构造层位精准符合的深度域数据体。
2、本专利技术一实施例提供了一种基于井震联合约束的深度域数据体获取方法，包括：
3、获取目标区域的测井数据标定结果，根据所述测井数据标定结果的时间-深度对应关系，修正测井层速度曲线，得到测井修正层速度数据；
4、通过时间域构造标准层约束，将井点对应层位深度进行全区插值，得到深度域插值层位数据；
5、获取所述目标区域中目标地层的构造解释层位，应用所述测井修正层速度数据，沿构造格架控制，进行测井速度空间插值，得到所述目标地层的地层层速度数据体；
6、在所述目标地层中选取基准面，根据所述深度域插值层位数据形成所述基准面的时间-深度基准对应关系，并以所述时间-深度基准对应关系为基准值进行速度积分运算，得到所述目标地层的时间域数据对应的深度域数据；
7、获取所述目标地层的时间域数据的采样间隔，根据所述时间域数据的采样间隔和所述地层层速度数据体计算得到深度域采样间隔；
8、根据所述深度域采样间隔对所述深度域数据进行重采样，得到深度域数据体。
9、作为上述方案的改进，所述通过时间域构造标准层约束，将井点对应层位深度进行全区插值，得到深度域插值层位数据，包括：
10、根据修正后的测井层速度曲线计算井点处目标地层的平均速度；
11、通过下式将多个井点关系在全区进行加权运算，得到所述目标地层的空变平均速度，并对所述目标地层的时间域层位数据进行深度域转换，获得所述目标地层深度域插值数据：
12、dxy＝∑(wwell*velwell)*txy
13、式中，dxy为所述目标地层在空间x，y位置处的深度域插值数据，wwell为权系数，velwell为所述井点处目标地层的平均速度，txy为空间x，y位置处的时间值。
14、作为上述方案的改进，所述在所述目标地层中选取基准面，根据所述深度域插值层位数据形成所述基准面的时间-深度基准对应关系，并以所述时间-深度基准对应关系为基准值进行速度积分运算，得到所述目标地层的时间域数据对应的深度域数据，包括：
15、选取所述目标地层的顶面或底面中的一个层面作为基准面；
16、将所述基准面的时间-深度基准对应关系(t0，d0)作为初始计算点，通过下式开展速度积分运算，得到所述目标地层的时间域数据对应的深度域数据：
17、
18、式中，(tk，dk)为所述目标地层对应位置任意点的时间-深度对应关系。
19、作为上述方案的改进，所述获取所述目标地层的时间域数据的采样间隔，根据所述时间域数据的采样间隔和所述地层层速度数据体计算得到深度域采样间隔，包括：
20、通过下式计算所述深度域采样间隔：
21、δd＝int(vmax*δt/2)
22、式中，δd为所述深度域采样间隔，vmax为所述地层层速度数据体的最大值，δt为所述目标地层的时间域数据的采样间隔。
23、作为上述方案的改进，所述进行测井速度空间插值，包括：
24、采用克里金插值、反距离插值、三角加权插值中的至少一种算法进行测井速度空间插值。
25、作为上述方案的改进，所述权系数wwell采用克里金插值、反距离插值、三角加权插值中的至少一种算法获得。
26、本专利技术另一实施例对应提供了一种基于井震联合约束的深度域数据体获取系统，包括：
27、速度修正模块，用于获取目标区域的测井数据标定结果，根据所述测井数据标定结果的时间-深度对应关系，修正测井层速度曲线，得到测井修正层速度数据；
28、深度插值模块，用于通过时间域构造标准层约束，将井点对应层位深度进行全区插值，得到深度域插值层位数据；
29、速度插值模块，用于获取所述目标区域中目标地层的构造解释层位，应用所述测井修正层速度数据，沿构造格架控制，进行测井速度空间插值，得到所述目标地层的地层层速度数据体；
30、深度域转换模块，用于在所述目标地层中选取基准面，根据所述深度域插值层位数据形成所述基准面的时间-深度基准对应关系，并以所述时间-深度基准对应关系为基准值进行速度积分运算，得到所述目标地层的时间域数据对应的深度域数据；
31、采样间隔计算模块，用于获取所述目标地层的时间域数据的采样间隔，根据所述时间域数据的采样间隔和所述地层层速度数据体计算得到深度本文档来自技高网...
【技术保护点】
1.一种基于井震联合约束的深度域数据体获取方法，其特征在于，包括：
2.如权利要求1所述的基于井震联合约束的深度域数据体获取方法，其特征在于，所述通过时间域构造标准层约束，将井点对应层位深度进行全区插值，得到深度域插值层位数据，包括：
3.如权利要求1所述的基于井震联合约束的深度域数据体获取方法，其特征在于，所述在所述目标地层中选取基准面，根据所述深度域插值层位数据形成所述基准面的时间-深度基准对应关系，并以所述时间-深度基准对应关系为基准值进行速度积分运算，得到所述目标地层的时间域数据对应的深度域数据，包括：
4.如权利要求1所述的基于井震联合约束的深度域数据体获取方法，其特征在于，所述获取所述目标地层的时间域数据的采样间隔，根据所述时间域数据的采样间隔和所述地层层速度数据体计算得到深度域采样间隔，包括：
5.如权利要求1所述的基于井震联合约束的深度域数据体获取方法，其特征在于，所述进行测井速度空间插值，包括：
6.如权利要求2所述的基于井震联合约束的深度域数据体获取方法，其特征在于，所述权系数Wwell采用克里金插值
7.一种基于井震联合约束的深度域数据体获取系统，其特征在于，包括：
8.一种基于井震联合约束的深度域数据体获取系统，其特征在于，所述通过时间域构造标准层约束，将井点对应层位深度进行全区插值，得到深度域插值层位数据，包括：
9.一种基于井震联合约束的深度域数据体获取系统，包括处理器、存储器以及存储在所述存储器中且被配置为由所述处理器执行的计算机程序，所述处理器执行所述计算机程序时实现如权利要求1至7中任意一项所述的基于井震联合约束的深度域数据体获取方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质包括存储的计算机程序，其中，在所述计算机程序运行时控制所述计算机可读存储介质所在设备执行如权利要求1至7中任意一项所述的基于井震联合约束的深度域数据体获取方法。
...【技术特征摘要】
1.一种基于井震联合约束的深度域数据体获取方法，其特征在于，包括：
2.如权利要求1所述的基于井震联合约束的深度域数据体获取方法，其特征在于，所述通过时间域构造标准层约束，将井点对应层位深度进行全区插值，得到深度域插值层位数据，包括：
3.如权利要求1所述的基于井震联合约束的深度域数据体获取方法，其特征在于，所述在所述目标地层中选取基准面，根据所述深度域插值层位数据形成所述基准面的时间-深度基准对应关系，并以所述时间-深度基准对应关系为基准值进行速度积分运算，得到所述目标地层的时间域数据对应的深度域数据，包括：
4.如权利要求1所述的基于井震联合约束的深度域数据体获取方法，其特征在于，所述获取所述目标地层的时间域数据的采样间隔，根据所述时间域数据的采样间隔和所述地层层速度数据体计算得到深度域采样间隔，包括：
5.如权利要求1所述的基于井震联合约束的深度域数据体获取方法，其特征在于，所述进行测井速度空间插值，包括：
6...
【专利技术属性】
技术研发人员：张永生，邢恩袁，桂宝玲，左璠璠，苏奎，牛新生，彭渊，商雯君，
申请(专利权)人：中国地质科学院矿产资源研究所，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。