本申请公开了一种井控层析成像建模方法、装置、电子设备及介质。该方法可以包括：对于每口井，计算井数据的三维矩形窗体；在三维矩形窗体内，计算每口井的权系数；以权系数为井约束项的权系数，进行井控层析成像建模。本发明专利技术通过计算权系数的方法，平衡井约束项在方程求解中的影响，与常规权系数扫描法相比，可以快速寻找合适的权重系数来平衡井位置的地震道集和井信息的影响，适应性好、计算效率高，使层析运算快速收敛，提升反演模型可靠性。提升反演模型可靠性。提升反演模型可靠性。
【技术实现步骤摘要】
井控层析成像建模方法、装置、电子设备及介质
[0001]本专利技术涉及油气勘探开发领域，更具体地，涉及一种井控层析成像建模方法、装置、电子设备及介质。
技术介绍
[0002]对于层析成像速度建模方法，越来越多的研究表明，反演方程中加入井信息约束项，可以加快反演的收敛速度，提升反演模型可靠性。目前井控层析反演逐渐成为速度建模领域一项必备的技术，但是井信息约束项的权系数选取仍存在一定不足。
[0003]由于井信息相对于地震成像道集个数较少，因此在反演方程中需要对井信息加入权系数，来提升井信息的反演权重。如果权重过小，会导致井信息贡献很小，达不到约束反演的目的；如果权重太大，反演过程中会在井口位置产生速度极值，引起牛眼现象。目前常用的方法是权系数扫描法，选取多个不同权系数进行反演测试，通过测试效果来选取合适的权系数。这种方法耗费较多机时、降低计算效率，同时扫描方法无法建立空间变化的权系数。
[0004]因此，有必要开发一种井控层析成像建模方法、装置、电子设备及介质，针对井控层析成像方程中井约束项权系数难以求取的问题，自动计算权系数，平衡井约束项在方程求解中的影响。
[0005]公开于本专利技术
技术介绍
部分的信息仅仅旨在加深对本专利技术的一般
技术介绍
的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
技术实现思路
[0006]本专利技术提出了一种井控层析成像建模方法、装置、电子设备及介质，其能够通过计算权系数的方法，平衡井约束项在方程求解中的影响，与常规权系数扫描法相比，可以快速寻找合适的权重系数来平衡井位置的地震道集和井信息的影响，适应性好、计算效率高，使层析运算快速收敛，提升反演模型可靠性。
[0007]第一方面，本公开实施例提供了一种井控层析成像建模方法，包括：
[0008]对于每口井，计算井数据的三维矩形窗体；
[0009]在所述三维矩形窗体内，计算每口井的权系数；
[0010]以所述权系数为井约束项的权系数，进行井控层析成像建模。
[0011]优选地，计算井数据的三维矩形窗体包括：
[0012]确定窗口参数；
[0013]确定每口井实钻井点的空间位置的横、纵坐标以及最小、最大深度，根据所述窗口参数，计算井数据的三维矩形窗体。
[0014]优选地，所述三维矩形窗体为x∈(X
‑
WIN_SIZE,X+WIN_SIZE)、y∈(Y
‑
WIN_SIZE,Y+WIN_SIZE)、z∈(Zmin,Zmax)，其中，X为横坐标，Y为纵坐标，Zmin为最小深度，Zmax为最大深
度，WIN_SIZE为窗口参数。
[0015]优选地，在所述三维矩形窗体内，计算每口井的权系数包括：
[0016]计算该井所有井点处井震误差的平方和；
[0017]在三维矩形窗体内，计算地震道集走时残差的平方和；
[0018]根据井震误差的平方和与地震道集走时残差的平方和，计算所述权系数。
[0019]优选地，通过公式(1)计算所述权系数：
[0020]W＝ d_seismic/d_well
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0021]其中，W为权系数，d_seismic为地震道集走时残差的平方和，d_well为井震误差的平方和。
[0022]优选地，通过公式(2)进行井控层析成像建模：
[0023][0024]其中，S(m)为层析反演目标函数，为共成像点道集反演项，为井信息反演项，W
i
为第i个权系数。
[0025]作为本公开实施例的一种具体实现方式，
[0026]第二方面，本公开实施例还提供了一种井控层析成像建模装置，包括：
[0027]三维矩形窗体建立模块，对于每口井，计算井数据的三维矩形窗体；
[0028]权系数计算模块，在所述三维矩形窗体内，计算每口井的权系数；
[0029]建模模块，以所述权系数为井约束项的权系数，进行井控层析成像建模。
[0030]优选地，计算井数据的三维矩形窗体包括：
[0031]确定窗口参数；
[0032]确定每口井实钻井点的空间位置的横、纵坐标以及最小、最大深度，根据所述窗口参数，计算井数据的三维矩形窗体。
[0033]优选地，所述三维矩形窗体为x∈(X
‑
WIN_SIZE,X+WIN_SIZE)、y∈(Y
‑
WIN_SIZE,Y+WIN_SIZE)、z∈(Zmin,Zmax)，其中，X为横坐标，Y为纵坐标，Zmin为最小深度，Zmax为最大深度，WIN_SIZE为窗口参数。
[0034]优选地，在所述三维矩形窗体内，计算每口井的权系数包括：
[0035]计算该井所有井点处井震误差的平方和；
[0036]在三维矩形窗体内，计算地震道集走时残差的平方和；
[0037]根据井震误差的平方和与地震道集走时残差的平方和，计算所述权系数。
[0038]优选地，通过公式(1)计算所述权系数：
[0039]W＝d_seismic/d_well
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0040]其中，W为权系数，d_seismic为地震道集走时残差的平方和，d_well为井震误差的平方和。
[0041]优选地，通过公式(2)进行井控层析成像建模：
[0042][0043]其中，S(m)为层析反演目标函数，为共成像点道集反演项，为
井信息反演项，W
i
为第i个权系数。
[0044]第三方面，本公开实施例还提供了一种电子设备，该电子设备包括：
[0045]存储器，存储有可执行指令；
[0046]处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的井控层析成像建模方法。
[0047]第四方面，本公开实施例还提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的井控层析成像建模方法。
[0048]其有益效果在于：本专利技术能够快速计算井信息的最优权重系数，平衡井约束项在方程求解中的影响，与常规权系数扫描法相比，本专利技术可以根据每口井数据特点分别计算系数，实现权系数的快速求取、空间变化，具有适应性好、计算效率高的特点，使层析运算快速收敛，提升反演模型可靠性。
[0049]本专利技术的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本专利技术的特定原理。
附图说明
[0050]通过结合附图对本专利技术示例性实施例进行更详细的描述，本专利技术的上述以及其它目的、特征和优势将变得更加明显，其中，在本专利技术示例性实施例中，相同的参本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种井控层析成像建模方法，其特征在于，包括：对于每口井，计算井数据的三维矩形窗体；在所述三维矩形窗体内，计算每口井的权系数；以所述权系数为井约束项的权系数，进行井控层析成像建模。2.根据权利要求1所述的井控层析成像建模方法，其中，计算井数据的三维矩形窗体包括：确定窗口参数；确定每口井实钻井点的空间位置的横、纵坐标以及最小、最大深度，根据所述窗口参数，计算井数据的三维矩形窗体。3.根据权利要求2所述的井控层析成像建模方法，其中，所述三维矩形窗体为x∈(X
‑
WIN_SIZE,X+WIN_SIZE)、y∈(Y
‑
WIN_SIZE,Y+WIN_SIZE)、z∈(Zmin,Zmax)，其中，X为横坐标，Y为纵坐标，Zmin为最小深度，Zmax为最大深度，WIN_SIZE为窗口参数。4.根据权利要求1所述的井控层析成像建模方法，其中，在所述三维矩形窗体内，计算每口井的权系数包括：计算该井所有井点处井震误差的平方和；在三维矩形窗体内，计算地震道集走时残差的平方和；根据井震误差的平方和与地震道集走时残差的平方和，计算所述权系数。5.根据权利要求4所述的井控层析成像建模方法，其中，通过公式(1)计算所述权系数：W＝d_seismic/d_well
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)其中，W为权系数，d_seismic为地震道集走时残差的平方和，d_well为井震误差的平方和。6.根据权利要求1所述的井控...
【专利技术属性】
技术研发人员：齐鹏，胡玮，
申请(专利权)人：中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。