CN-121995444-A - Chromatography inversion iterative computation quality control method, device, equipment, storage medium and product

CN121995444ACN 121995444 ACN121995444 ACN 121995444ACN-121995444-A

Abstract

The application provides a chromatography inversion iterative computation quality control method, a device, equipment, a storage medium and a product, relates to the technical field of seismic exploration, and aims to solve the problem that the related chromatography inversion technology is difficult to intuitively analyze the specific cause of the error result. The method comprises the steps of performing forward modeling on an initial velocity model, calculating second first arrival time of each receiving channel, calculating difference values between the second arrival time of each receiving channel and the corresponding picked first arrival time to obtain errors corresponding to each receiving channel, performing root mean square error statistics according to the errors corresponding to all receiving channels of each shot point to obtain root mean square errors corresponding to each shot point, performing multi-round iterative updating on the initial velocity model based on a tomographic inversion algorithm, visually displaying seismic wave inversion errors of all shot points after each round iterative calculation to obtain a multi-round error distribution map, calculating seismic wave inversion error reduction percentages of all shot points after multiple rounds of iterative calculation, and analyzing a tomographic inversion effect.

Inventors

ZHANG HENG
WANG ZHAOLEI
WANG GUANGYIN
LIU HONG
XIN KEFENG
JIN DEGANG

Assignees

中国石油天然气集团有限公司
中国石油集团东方地球物理勘探有限责任公司
中油油气勘探软件国家工程研究中心有限公司

Dates

Publication Date: 20260508
Application Date: 20241101

Claims (10)

1. A chromatography-inversion iterative computation quality control method, characterized in that the method comprises: Acquiring a shot set record of a target area, wherein the shot set record comprises seismic wave information of a plurality of receiving channels corresponding to each shot point, and picking up the first arrival time of the seismic waves of each receiving channel from the shot set record; establishing an initial velocity model according to prior information, and performing forward modeling on the initial velocity model based on a ray tracing method to obtain second first arrival time of seismic waves of each receiving channel; calculating the difference value between the second first arrival time of the seismic waves of each receiving channel and the corresponding first arrival time to obtain the corresponding error of each receiving channel; Carrying out root mean square error statistics according to the errors corresponding to all the receiving channels of each shot point to obtain root mean square errors corresponding to each shot point; According to the root mean square error corresponding to each shot point, performing multi-round iterative updating on the initial velocity model based on a tomographic inversion algorithm, and performing visual display on the seismic wave inversion errors of all shot points after each round of iterative calculation to obtain a multi-round error distribution diagram; And respectively carrying out weight reduction calculation on the seismic wave inversion errors of the seismic waves of each shot point in different iteration rounds based on the multi-round error distribution diagram to obtain the seismic wave inversion error reduction percentage of each shot point after the multi-round iteration calculation, and analyzing the tomographic inversion effect based on the seismic wave inversion error reduction percentage.
2. The method for controlling quality of tomographic inversion iterative computation according to claim 1, wherein said performing root mean square error statistics according to the errors corresponding to all the receiving channels of each shot to obtain root mean square error corresponding to each shot comprises: the root mean square error corresponding to each shot point is calculated according to the following formula: , Wherein X rms represents the root mean square error corresponding to each shot point, X i represents the error corresponding to the ith receiving channel, and n represents the number of receiving channels corresponding to the shot point.
3. The method of claim 1, wherein the performing multiple iterative updates on the initial velocity model based on the tomographic inversion algorithm according to the root mean square error corresponding to each shot point comprises: and adjusting a speed value in the initial speed model according to the seismic wave inversion error to obtain an initial speed model after multiple rounds of iterative updating.
4. The method of claim 3, wherein the analyzing the tomographic inversion effect based on the percentage reduction of the seismic inversion error comprises: and under the condition that the reduction percentages of the seismic wave inversion errors of all shot points are smaller than a preset value, taking the initial velocity model after iterative updating of the round as a target velocity model.
5. The method of claim 4, wherein analyzing the tomographic inversion effect based on the percentage reduction of the seismic inversion error comprises: And under the condition that the reduction percentage of the seismic wave inversion error of a certain shot point is not smaller than a preset value, performing quality control adjustment on the certain shot point based on the multi-round error distribution map.
6. The method of claim 5, wherein performing the quality control adjustment on the shot point based on the multi-round error distribution map comprises: Analyzing the waveform characteristics of the certain shot point, wherein the waveform characteristics comprise amplitude, frequency and phase; and under the condition that abnormal amplitude, frequency deficiency or phase discontinuity exists in the waveform characteristics, corresponding adjustment is carried out, wherein the adjustment comprises the adjustment of source parameters, receiving equipment setting, sampling rate, filter setting and data alignment.
7. A tomographic inversion iterative computation quality control apparatus, the apparatus comprising: The acquisition module is used for acquiring a shot set record of the target area, wherein the shot set record comprises seismic wave information of a plurality of receiving channels corresponding to each shot point, and the first arrival time of the seismic waves of each receiving channel is picked up from the shot set record; The simulation module is used for establishing an initial velocity model according to the priori information, and performing forward modeling on the initial velocity model based on a ray tracing method to obtain second first arrival time of the seismic waves of each receiving channel; the first calculation module is used for calculating the difference value between the second first arrival time of the seismic waves of each receiving channel and the corresponding first arrival time to obtain the corresponding error of each receiving channel; The second calculation module is used for carrying out root mean square error statistics according to the errors corresponding to all the receiving channels of each shot point to obtain root mean square errors corresponding to each shot point; The inversion display module is used for carrying out multi-round iterative updating on the initial speed model based on a tomographic inversion algorithm according to the root mean square error corresponding to each shot point, and carrying out visual display on the seismic wave inversion errors of all the shot points after each round of iterative calculation to obtain a multi-round error distribution diagram; And the analysis module is used for respectively carrying out weight reduction calculation on the seismic wave inversion errors of the seismic waves of each shot point in different iteration rounds based on the multi-round error distribution diagram to obtain the seismic wave inversion error reduction percentage of each shot point after the multi-round iteration calculation, and analyzing the tomographic inversion effect based on the seismic wave inversion error reduction percentage.
8. An electronic device, comprising: A processor, a memory, and a computer program stored on the memory and executable on the processor, the processor implementing the tomographic inversion iterative computational quality control method of any one of claims 1 to 6 when the program is executed.
9. A computer readable storage medium, characterized in that instructions in the storage medium, when executed by a processor of an electronic device, enable the electronic device to perform the tomographic inversion iterative computational quality control method of any one of claims 1 to 6.
10. A computer program product comprising a computer program which, when executed by a processor, implements a tomographic inversion iterative computational quality control method as claimed in any one of claims 1 to 6.

Description

Chromatography inversion iterative computation quality control method, device, equipment, storage medium and product Technical Field The embodiment of the application relates to the technical field of seismic exploration, in particular to a chromatography inversion iterative computation quality control method, a device, equipment, a storage medium and a product. Background With the increase of exploration force, the exploration target is deep, the ground surface of a work area is complex, the lateral change of the shallow speed is severe, the traditional near-surface modeling and static correction technology cannot completely adapt to the processing requirements of the complex work area, and the near-surface speed cannot be finely depicted. The first-arrival chromatography inversion technology is a mainstream modeling technical means because of strong adaptability and more advanced principle. However, in practical application, the quality control means has limitations, such as root mean square error during fitting travel and picking travel, the iterative calculation result of each round is a value, the average effect is high, and the details of the calculation result are difficult to comprehensively reflect. Disclosure of Invention The embodiment of the application provides a chromatography inversion iterative computation quality control method, a device, equipment, a storage medium and a product, and aims to solve the problems that the computation process of a related first-arrival chromatography inversion technology is not visual, and the specific cause causing the error result is difficult to analyze. An embodiment of the present application provides a chromatography inversion iterative computation quality control method, which includes: Acquiring a shot set record of a target area, wherein the shot set record comprises seismic wave information of a plurality of receiving channels corresponding to each shot point, and picking up the first arrival time of the seismic waves of each receiving channel from the shot set record; establishing an initial velocity model according to prior information, and performing forward modeling on the initial velocity model based on a ray tracing method to obtain second first arrival time of seismic waves of each receiving channel; calculating the difference value between the second first arrival time of the seismic waves of each receiving channel and the corresponding first arrival time to obtain the corresponding error of each receiving channel; Carrying out root mean square error statistics according to the errors corresponding to all the receiving channels of each shot point to obtain root mean square errors corresponding to each shot point; According to the root mean square error corresponding to each shot point, performing multi-round iterative updating on the initial velocity model based on a tomographic inversion algorithm, and performing visual display on the seismic wave inversion errors of all shot points after each round of iterative calculation to obtain a multi-round error distribution diagram; And respectively carrying out weight reduction calculation on the seismic wave inversion errors of the seismic waves of each shot point in different iteration rounds based on the multi-round error distribution diagram to obtain the seismic wave inversion error reduction percentage of each shot point after the multi-round iteration calculation, and analyzing the tomographic inversion effect based on the seismic wave inversion error reduction percentage. In an optional implementation manner, the performing root mean square error statistics according to the errors corresponding to all the receiving channels of each shot point to obtain root mean square errors corresponding to each shot point includes: the root mean square error corresponding to each shot point is calculated according to the following formula: , Wherein X rms represents the root mean square error corresponding to each shot point, X i represents the error corresponding to the ith receiving channel, and n represents the number of receiving channels corresponding to the shot point. In an optional implementation manner, the performing multiple iterative updating on the initial velocity model based on a tomographic inversion algorithm according to the root mean square error corresponding to each shot point includes: and adjusting a speed value in the initial speed model according to the seismic wave inversion error to obtain an initial speed model after multiple rounds of iterative updating. In an alternative embodiment, said analyzing tomographic inversion effects based on said percent reduction in seismic inversion error comprises: and under the condition that the reduction percentages of the seismic wave inversion errors of all shot points are smaller than a preset value, taking the initial velocity model after iterative updating of the round as a target velocity model. In an alternative embodiment, said analy