Method and apparatus for determining geological change information, device, and storage medium
By employing a forward and reverse double-difference inversion strategy, benchmark and monitoring data were inverted sequentially, solving the problem of low accuracy in geological change information and achieving higher accuracy.
Patent Information
- Authority / Receiving Office
- WO · WO
- Patent Type
- Applications
- Current Assignee / Owner
- CHINA NAT PETROLEUM CORP
- Filing Date
- 2025-08-21
- Publication Date
- 2026-06-04
AI Technical Summary
In existing technologies, when constructing high-resolution velocity models of geological structures using full-waveform inversion methods, spurious components in the inversion results of baseline data and monitoring data are added together rather than canceled out, resulting in low accuracy of geological change information.
A forward and reverse double-difference inversion strategy was adopted. First, the baseline data was inverted, and then the monitoring data was inverted to obtain the velocity time-shift changes in the forward and reverse double-difference inversion, respectively. Then, the average value of the two was used as the change information of the target geology.
The accuracy of geological change information has been improved. The double-difference inversion strategy has eliminated the influence of spurious components and enhanced the accuracy of geological change information.
Smart Images
Figure CN2025116249_04062026_PF_FP_ABST
Abstract
Description
Methods, devices, equipment, and storage media for determining geological change information
[0001] This application claims priority to Chinese Patent Application No. 202411728913.5, filed on November 28, 2024, entitled “Method, Apparatus, Equipment and Storage Medium for Determining Geological Change Information”, the entire contents of which are incorporated herein by reference. Technical Field
[0002] This application relates to the field of computer technology, and specifically to a method, apparatus, device, and storage medium for determining geological change information. Background Technology
[0003] In some scenarios, high-resolution velocity models of geological structures can be constructed using the information contained in seismic data. For example, in oil reservoir development and carbon dioxide sequestration monitoring, velocity models can be used to obtain information about changes in geological structures.
[0004] In related technologies, four-dimensional seismic data can be processed using the Full-waveform inversion (FWI) method to construct high-resolution subsurface models. Parallel inversion strategies can be employed to invert baseline data and monitoring data separately, and then the difference between the inversion results of the baseline data and monitoring data can be calculated to obtain the temporal changes of the subsurface structure.
[0005] However, when subtracting the inversion results from the baseline data and the monitoring data, the spurious components in the two inversion results may be added together rather than canceled out, resulting in lower accuracy. Summary of the Invention
[0006] This application provides a method, apparatus, equipment, and storage medium for determining geological change information, in order to solve the problem of low accuracy of geological change information.
[0007] Firstly, this application provides a method for determining geological change information, including:
[0008] Obtain a data processing request, wherein the data processing request is used to request the acquisition of geological change information of the target geology;
[0009] According to the data processing request, the geological data of the target geology is obtained. The geological data includes the baseline data of the target geology at a historical time, the monitoring data at the most recent time, the first earthquake data at the historical time, and the second earthquake data at the most recent time.
[0010] The reference data is inverted to obtain a first seismic wavelet. Based on the first seismic wavelet and the geological data, the reference data is inverted, and the monitoring data is inverted to obtain the first geological change information of the target geology.
[0011] The monitoring data is inverted to obtain a second seismic wavelet. Based on the second seismic wavelet and the geological data, the monitoring data is inverted, and the reference data is inverted to obtain the second geological change information of the target geology.
[0012] Based on the first geological change information and the second geological change information, the target geological change information of the target geology is determined.
[0013] In one possible design, based on the first seismic wavelet and the geological data, the reference data is inverted, and the monitoring data is inverted to obtain the first geological change information of the target geology, including:
[0014] Based on the baseline data, the first seismic wavelet, and the first seismic data, a first baseline model for the target geology is determined;
[0015] Based on the baseline data, the monitoring data, the second seismic data, and the first baseline model, a first monitoring model for the target geology is determined.
[0016] The first geological change information is obtained by performing difference processing on the first monitoring model and the first benchmark model.
[0017] In one possible design, a first reference model of the target geology is determined based on the reference data, the first seismic wavelet, and the first seismic data, including:
[0018] The first seismic data is subjected to grid tomography to obtain an initial velocity model;
[0019] The reference data is bandpass filtered to obtain sub-reference data of at least one frequency band, and the first seismic wavelet is bandpass filtered to obtain sub-seismic wavelets of at least one frequency band.
[0020] The first reference model is determined based on the initial velocity model, the sub-reference data of the at least one frequency band, and the sub-seismic wavelet of the at least one frequency band.
[0021] In one possible design, the at least one frequency band includes a first frequency band, a second frequency band, and a third frequency band; the first reference model is determined based on the initial velocity model, the sub-reference data of the at least one frequency band, and the sub-seismic wavelet of the at least one frequency band, including:
[0022] The initial velocity model, the sub-reference data of the first frequency band, and the sub-seismic wavelet of the first frequency band are inverted to obtain the low-frequency model;
[0023] The low-frequency model, the sub-reference data of the second frequency band, and the sub-seismic wavelet of the second frequency band are inverted to obtain the mid-frequency model;
[0024] The first reference model is obtained by inverting the mid-frequency model, the sub-reference data of the third frequency band, and the sub-seismic wavelet of the third frequency band.
[0025] In one possible design, a first monitoring model for the target geology is determined based on the baseline data, the monitoring data, the second seismic data, and the first baseline model, including:
[0026] Based on the first benchmark model, the monitoring data, and the benchmark data, the pseudo-observation data and the seismic inversion wavelet are determined;
[0027] The pseudo-observation data is bandpass filtered to obtain sub-pseudo-observation data of at least one frequency band, and the seismic inversion wavelet is bandpass filtered to obtain sub-seismic inversion wavelet of at least one frequency band.
[0028] The first monitoring model is determined based on the first benchmark model, the sub-pseudo-observation data of the at least one frequency band, and the sub-seismic inversion wavelet of the at least one frequency band.
[0029] In one possible design, based on the first reference model, the monitoring data, and the reference data, the pseudo-observation data and the seismic inversion wavelet are determined, including:
[0030] Retrieve the preset seismic wavelet from the preset database;
[0031] The first reference model and the preset seismic wavelet are subjected to forward modeling to obtain simulation data;
[0032] The monitoring data and the benchmark data are differentially processed to obtain differential data;
[0033] The simulated data and the difference data are summed to obtain the pseudo-observation data;
[0034] The pseudo-observation data is inverted to obtain the seismic inversion wavelet.
[0035] In one possible design, the monitoring data is inverted based on the second seismic wavelet and the geological data, and the reference data is also inverted to obtain second geological change information of the target geology, including:
[0036] Based on the monitoring data, the second seismic wavelet, and the second seismic data, a second monitoring model for the target geology is determined;
[0037] Based on the baseline data, the monitoring data, the first seismic data, and the second monitoring model, a second baseline model for the target geology is determined.
[0038] The second geological change information is obtained by performing difference processing on the second monitoring model and the second benchmark model.
[0039] Secondly, this application provides a device for determining geological change information, comprising: a first acquisition module, a second acquisition module, a first inversion processing module, a second inversion processing module, and a determination module, wherein,
[0040] The first acquisition module is used to acquire a data processing request, wherein the data processing request is used to request the acquisition of geological change information of the target geology;
[0041] The second acquisition module is used to acquire geological data of the target geology according to the data processing request. The geological data includes baseline data of the target geology at a historical time, monitoring data at a recent time, first earthquake data at the historical time, and second earthquake data at the recent time.
[0042] The first inversion processing module is used to perform inversion processing on the reference data to obtain a first seismic wavelet, perform inversion processing on the reference data based on the first seismic wavelet and the geological data, and perform inversion processing on the monitoring data to obtain the first geological change information of the target geology;
[0043] The second inversion processing module is used to perform inversion processing on the monitoring data to obtain a second seismic wavelet, and to perform inversion processing on the monitoring data and the reference data based on the second seismic wavelet and the geological data to obtain the second geological change information of the target geology.
[0044] The determining module is used to determine the target geological change information of the target geology based on the first geological change information and the second geological change information.
[0045] In one possible design, the first inversion processing module is specifically used for:
[0046] Based on the baseline data, the first seismic wavelet, and the first seismic data, a first baseline model for the target geology is determined;
[0047] Based on the baseline data, the monitoring data, the second seismic data, and the first baseline model, a first monitoring model for the target geology is determined.
[0048] The first geological change information is obtained by performing difference processing on the first monitoring model and the first benchmark model.
[0049] In one possible design, the first inversion processing module is specifically used for:
[0050] The first seismic data is subjected to grid tomography to obtain an initial velocity model;
[0051] The reference data is bandpass filtered to obtain sub-reference data of at least one frequency band, and the first seismic wavelet is bandpass filtered to obtain sub-seismic wavelets of at least one frequency band.
[0052] The first reference model is determined based on the initial velocity model, the sub-reference data of the at least one frequency band, and the sub-seismic wavelet of the at least one frequency band.
[0053] In one possible design, the first inversion processing module is specifically used for:
[0054] The initial velocity model, the sub-reference data of the first frequency band, and the sub-seismic wavelet of the first frequency band are inverted to obtain the low-frequency model;
[0055] The low-frequency model, the sub-reference data of the second frequency band, and the sub-seismic wavelet of the second frequency band are inverted to obtain the mid-frequency model;
[0056] The first reference model is obtained by inverting the mid-frequency model, the sub-reference data of the third frequency band, and the sub-seismic wavelet of the third frequency band.
[0057] In one possible design, the first inversion processing module is specifically used for:
[0058] Based on the first benchmark model, the monitoring data, and the benchmark data, the pseudo-observation data and the seismic inversion wavelet are determined;
[0059] The pseudo-observation data is bandpass filtered to obtain sub-pseudo-observation data of at least one frequency band, and the seismic inversion wavelet is bandpass filtered to obtain sub-seismic inversion wavelet of at least one frequency band.
[0060] The first monitoring model is determined based on the first benchmark model, the sub-pseudo-observation data of the at least one frequency band, and the sub-seismic inversion wavelet of the at least one frequency band.
[0061] In one possible design, the first inversion processing module is specifically used for:
[0062] Retrieve the preset seismic wavelet from the preset database;
[0063] The first reference model and the preset seismic wavelet are subjected to forward modeling to obtain simulation data;
[0064] The monitoring data and the benchmark data are differentially processed to obtain differential data;
[0065] The simulated data and the difference data are summed to obtain the pseudo-observation data;
[0066] The pseudo-observation data is inverted to obtain the seismic inversion wavelet.
[0067] In one possible design, the second inversion processing module is specifically used for:
[0068] Based on the monitoring data, the second seismic wavelet, and the second seismic data, a second monitoring model for the target geology is determined;
[0069] Based on the baseline data, the monitoring data, the first seismic data, and the second monitoring model, a second baseline model for the target geology is determined.
[0070] The second geological change information is obtained by performing difference processing on the second monitoring model and the second benchmark model.
[0071] Thirdly, embodiments of this application provide an electronic device, including: at least one processor and a memory; the memory stores computer-executable instructions; the at least one processor executes the computer-executable instructions stored in the memory, causing the at least one processor to perform the method for determining geological change information as described in the first aspect above and various possible designs of the first aspect.
[0072] Fourthly, embodiments of this application provide a computer-readable storage medium storing computer-executable instructions. When a processor executes the computer-executable instructions, it implements the method for determining geological change information as described in the first aspect above and various possible designs of the first aspect.
[0073] Fifthly, embodiments of this application provide a computer program product, including a computer program that, when executed by a processor, implements the method for determining geological change information as described in the first aspect above and various possible designs of the first aspect.
[0074] The method, apparatus, equipment, and storage medium for determining geological change information provided in this application can, when it is necessary to determine the change information of a target geological location from a historical time to the most recent time, determine the baseline data of the historical time and the monitoring data of the most recent time based on the geological data of the target geological location. By first inverting the baseline data and then inverting the monitoring data, the forward double-difference velocity-time change is obtained. Then, by first inverting the monitoring data and then inverting the baseline data, the reverse double-difference velocity-time change is obtained. Finally, the average value of the two velocity-time changes is determined as the change information of the target geological location, thereby improving the accuracy of geological change information. Attached Figure Description
[0075] The accompanying drawings, which are incorporated in and form part of this specification, illustrate embodiments consistent with this application and, together with the description, serve to explain the principles of this application.
[0076] Figure 1 is a schematic diagram of the system architecture provided in an embodiment of this application;
[0077] Figure 2 is a flowchart illustrating a method for determining geological change information provided in an embodiment of this application;
[0078] Figure 3 is a schematic diagram of the process for determining the first geological change information provided in an embodiment of this application;
[0079] Figure 4 is a schematic diagram of the process of determining the first reference model provided in the embodiments of this application;
[0080] Figure 5 is a schematic diagram of another process for determining first geological change information provided in an embodiment of this application;
[0081] Figure 6 is a schematic diagram of the process for determining the second geological change information provided in an embodiment of this application;
[0082] Figure 7 is a schematic diagram of the determination process of the second monitoring model provided in the embodiments of this application;
[0083] Figure 8 is a schematic diagram of another process for determining second geological change information provided in an embodiment of this application;
[0084] Figure 9 is a schematic diagram of a two-dimensional model of the baseline data;
[0085] Figure 10 is a schematic diagram of the two-dimensional model of the initial velocity model;
[0086] Figure 11 is a schematic diagram of the two-dimensional model of the first reference model;
[0087] Figure 12 is a schematic diagram of a two-dimensional model of the first geological change information;
[0088] Figure 13 is a schematic diagram of a two-dimensional model of the second geological change information;
[0089] Figure 14 is a schematic diagram of a two-dimensional model of the target geological change information;
[0090] Figure 15 is a schematic diagram of a two-dimensional model of actual geological changes of the target.
[0091] Figure 16 is a schematic diagram of a device for determining geological change information provided in an embodiment of this application;
[0092] Figure 17 is a schematic diagram of the structure of the electronic device provided in the embodiment of this application.
[0093] The accompanying drawings illustrate specific embodiments of this application, which will be described in more detail below. These drawings and descriptions are not intended to limit the scope of the concept in any way, but rather to illustrate the concept of this application to those skilled in the art through reference to particular embodiments. Detailed Implementation
[0094] Exemplary embodiments will now be described in detail, examples of which are illustrated in the accompanying drawings. When the following description relates to the drawings, unless otherwise indicated, the same numbers in different drawings denote the same or similar elements. The embodiments described in the following exemplary embodiments do not represent all embodiments consistent with this application. Rather, they are merely examples of apparatuses and methods consistent with some aspects of this application as detailed in the appended claims.
[0095] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, stored data, displayed data, etc.) involved in one or more embodiments of this specification are all information and data authorized by the user or fully authorized by all parties. Furthermore, the collection, use and processing of related data must comply with relevant laws, regulations and standards, and corresponding operation entry points are provided for users to choose to authorize or refuse.
[0096] It should be noted that in the embodiments of this application, certain software, components, models and other existing solutions in the industry may be mentioned. These should be regarded as exemplary and are only intended to illustrate the feasibility of implementing the technical solution of this application. However, it does not mean that the applicant has used or necessarily used the solution.
[0097] For ease of understanding, the system architecture applicable to the embodiments of this application will be described below with reference to Figure 1.
[0098] Figure 1 is a schematic diagram of the system architecture provided in an embodiment of this application. Referring to Figure 1, it includes a terminal device 101 and a storage device 102. The terminal device 101 can be a user's local device, and it may contain data processing software. The user can use the data processing software to determine geological change information. The storage device 102 may contain a preset database for storing geological data, preset seismic wavelets, and other data information. The terminal device 101 can retrieve the data stored in the storage device 102. The storage device 102 can be a USB flash drive, hard drive, or cloud storage device, etc.
[0099] In related technologies, FWI (Focused Inversion Theory) can be used to process four-dimensional seismic data to construct high-resolution subsurface models. Parallel inversion strategies can be employed to invert baseline and monitoring data separately, and then the difference between the inversion results of the baseline and monitoring data is calculated to obtain the temporal changes in subsurface structures. However, when the inversion results of the baseline and monitoring data are subtracted, spurious components in the two inversion results may be additive rather than canceled out, leading to lower accuracy.
[0100] To address the aforementioned technical issues, in this embodiment of the application, when it is necessary to determine the change information of the target geology from a historical time to the most recent time, the baseline data of the historical time and the monitoring data of the most recent time can be determined based on the geological data of the target geology. By first inverting the baseline data and then inverting the monitoring data, the forward double-difference velocity time-shift change is obtained. Then, by first inverting the monitoring data and then inverting the baseline data, the reverse double-difference velocity time-shift change is obtained. Finally, the average value of the two velocity time-shift changes is determined as the change information of the target geology, thereby improving the accuracy of the geological change information.
[0101] The technical solution of this application and how the technical solution of this application solves the above-mentioned technical problems are described in detail below with specific embodiments. These specific embodiments can be combined with each other, and the same or similar concepts or processes may not be described again in some embodiments. The embodiments of this application will now be described with reference to the accompanying drawings.
[0102] Figure 2 is a flowchart illustrating a method for determining geological change information according to an embodiment of this application. Referring to Figure 2, the method may include the following steps:
[0103] S201, Obtain data processing request.
[0104] The execution entity in this application embodiment can be a terminal device or a device for determining geological change information installed in the terminal device. The device for determining geological change information can be implemented by software or by a combination of software and hardware.
[0105] A data processing request is used to request geological change information of a target geological site. The terminal device may have a built-in data processing system with a visual interface. Users can select the target geological site, historical time, and most recent time on the visual interface, and generate a data processing request by clicking the "Process" control.
[0106] Among them, target geology can refer to the geological area where the user needs to obtain information on the changes in the target geology from historical time to the most recent time. For example, the target address can be a geological area with oil and gas reservoirs or a geological area containing carbon dioxide.
[0107] Geological change information of the target geology can refer to the temporal shift of the target geology from a historical moment to a recent moment. For example, a historical moment can refer to a certain survey moment ten years ago, and a recent moment can refer to a certain survey moment this year. Geological change information can refer to the temporal shift between the velocity of the target geology ten years ago and the velocity of the target geology this year.
[0108] S202. Obtain geological data of the target geology according to the data processing request.
[0109] The geological data includes baseline data of the target geology at historical times, monitoring data at the most recent time, first earthquake data at historical times, and second earthquake data at the most recent time.
[0110] Seismic data for the target geology can be obtained through on-site surveys. The seismic data can be stored in a preset database on the terminal device, which can be located in the terminal device's storage space.
[0111] Reference data refers to the magnitude of the signal received by the receiver at a target geological area at a historical moment, originating from the shot point. Here, the receiver can refer to the instrument that receives the signal, such as a seismic detector, and the shot point can refer to the location of the seismic wave source in seismic exploration. Reference data can be represented by d1(x,y,t), where x can refer to the receiver's coordinates along the longitudinal survey line, y can refer to the receiver's coordinates along a non-longitudinal survey line, and t can refer to the travel time coordinates, which indicate the time at which the receiver received the measurement signal.
[0112] Monitoring data refers to the magnitude of the signal received by the receiver from the shot point in the target geological area at the most recent moment. Monitoring data can be represented by d2(x,y,t), where x represents the receiver's coordinates along the longitudinal survey line, y represents the receiver's coordinates along a non-longitudinal survey line, and t represents the travel time coordinates.
[0113] The first type of earthquake data can refer to earthquake velocity data at historical moments, such as the P-wave velocity. The second type of earthquake data can refer to earthquake velocity data at the most recent moment.
[0114] S203. Perform inversion processing on the baseline data to obtain the first seismic wavelet. Based on the first seismic wavelet and geological data, perform inversion processing on the baseline data and the monitoring data to obtain the first geological change information of the target geology.
[0115] Inversion processing refers to the process of inferring the underground geological structure and properties from observation data. The terminal equipment can be equipped with an inversion processing algorithm module. Through inversion processing, relevant data on the geological structure can be extracted from the geological data. For example, by inverting the reference data, the first seismic wavelet can be obtained.
[0116] The first seismic wavelet can refer to the seismic wave characteristics extracted from the reference data after inversion processing. The first seismic wavelet can be represented by w1(x,y,t), where x can refer to the coordinates of the receiver along the longitudinal line, y can refer to the coordinates of the receiver along the non-longitudinal line, and t can refer to the travel time coordinates.
[0117] The primary geological change information can refer to the temporal variation of subsurface velocity in a target geological area, determined through a forward double-difference inversion strategy. The forward double-difference inversion strategy involves first inverting the baseline data and then inverting the monitoring data.
[0118] The first geological change information can be determined as follows: Based on the baseline data, the first seismic wavelet, and the first seismic data, determine the first baseline model of the target geology; based on the baseline data, monitoring data, the second seismic data, and the first baseline model, determine the first monitoring model of the target geology; perform difference processing on the first monitoring model and the first baseline model to obtain the first geological change information.
[0119] The first reference model can refer to the velocity model obtained after inverting the reference data, the first seismic wavelet, and the first seismic data. The first monitoring model can refer to the velocity model obtained after inverting the reference data, monitoring data, the second seismic data, and the first reference model. Differential processing refers to subtracting the first monitoring model from the first reference model.
[0120] S204. Perform inversion processing on the monitoring data to obtain the second seismic wavelet. Based on the second seismic wavelet and geological data, perform inversion processing on the monitoring data and the reference data to obtain the second geological change information of the target geology.
[0121] The second seismic wavelet refers to the seismic wave characteristics extracted from the monitoring data after inversion processing. The second seismic wavelet can be represented by w2(x,y,t), where x represents the coordinates of the receiver along the longitudinal survey line, y represents the coordinates of the receiver along a non-longitudinal survey line, and t represents the travel time coordinates.
[0122] The second geological change information can refer to the time-shift changes in the subsurface velocity of the target geology, determined through a reverse double-difference inversion strategy. The reverse double-difference inversion strategy refers to a method that first inverts the monitoring data and then inverts the baseline data.
[0123] The second geological change information can be determined as follows: Based on monitoring data, the second seismic wavelet, and the second seismic data, determine the second monitoring model of the target geology; based on the baseline data, monitoring data, the first seismic data, and the second monitoring model, determine the second baseline model of the target geology; perform difference processing on the second monitoring model and the second baseline model to obtain the second geological change information.
[0124] The second monitoring model can refer to the velocity model obtained by inverting monitoring data, the second seismic wavelet, and the second seismic data. The second reference model can refer to the velocity model obtained by inverting reference data, monitoring data, the first seismic data, and the second monitoring model. Difference processing refers to subtracting the second monitoring model from the second reference model.
[0125] S205. Based on the first geological change information and the second geological change information, determine the target geological change information of the target geology.
[0126] The target geological change information can be determined as follows: obtain the first geological change information and the second geological change information, calculate the average value of the first geological change information and the second geological change information, and determine the average value as the target geological change information.
[0127] For example, assuming the target geological change information is Δm, and the first geological change information is Δm ＝> The second geological change information is Δm <＝ The first monitoring model is m′ m The first baseline model is m b The second monitoring model is m m The second baseline model is m′ b Then the target geological change information can be determined as follows:
[0128] In this embodiment, when it is necessary to determine the change information of a target geology from a historical time to the most recent time, the geological data of the target geology can be determined according to the data processing request. This includes the baseline data of the target geology at the historical time, the monitoring data at the most recent time, the first seismic data at the historical time, and the second seismic data at the most recent time. Inversion processing of the baseline data yields a first seismic wavelet. Based on the first seismic wavelet and the geological data, inversion processing of the baseline data yields a first baseline model. Inversion processing of the monitoring data yields a first monitoring model. The difference between the first monitoring model and the first baseline model yields the first geological change information. Inversion processing of the monitoring data yields a second seismic wavelet. Based on the second seismic wavelet and the geological data, inversion processing of the monitoring data yields a second monitoring model. Inversion processing of the baseline data yields a second baseline model. The difference between the second monitoring model and the second baseline model yields the second geological change information. The average of the first geological change information and the second geological change information yields the target geological change information, thus improving the accuracy of the geological change information.
[0129] Based on any of the above embodiments, the process of determining the first geological change information (S203 in the embodiment of Figure 2) will be described in detail below with reference to Figure 3.
[0130] Figure 3 is a schematic diagram illustrating the process of determining the first geological change information provided in an embodiment of this application. Referring to Figure 3, the method may include:
[0131] S301. Obtain geological data of the target geology according to the data processing request.
[0132] The geological data may include baseline data of the target geology at historical times, monitoring data at the most recent time, first earthquake data at historical times, and second earthquake data at the most recent time.
[0133] It should be noted that the execution process of S301 above can be found in S202, and will not be repeated here.
[0134] S302. Perform grid tomography on the first earthquake data to obtain the initial velocity model.
[0135] First earthquake data can refer to earthquake velocity data at a historical moment, such as the velocity of earthquake P-waves.
[0136] Grid tomography refers to velocity modeling processing technology in seismic exploration. The terminal equipment can be equipped with a grid tomography algorithm module. Through grid tomography, the first seismic data can be converted into discrete grid data, and the grid data can be updated using iterative optimization methods to obtain an initial velocity model.
[0137] S303. Perform bandpass filtering on the reference data to obtain sub-reference data for at least one frequency band, and perform bandpass filtering on the first seismic wavelet to obtain sub-seismic wavelet for at least one frequency band.
[0138] At least one frequency band may include a first frequency band, a second frequency band, and a third frequency band. For example, the first frequency band may be a low-frequency band, the second frequency band may be a mid-frequency band, and the third frequency band may be a high-frequency band.
[0139] Sub-reference data can be determined as follows: determine the frequency ranges of the first, second, and third frequency bands, and perform bandpass filtering on the reference data according to the frequency ranges of the first, second, and third frequency bands to obtain the sub-reference data of the first, second, and third frequency bands, respectively.
[0140] The sub-seismic wavelet can be determined as follows: determine the frequency ranges of the first, second, and third frequency bands, and perform bandpass filtering on the first seismic wavelet according to the frequency ranges of the first, second, and third frequency bands to obtain the sub-seismic wavelet of the first frequency band, the sub-seismic wavelet of the second frequency band, and the sub-seismic wavelet of the third frequency band, respectively.
[0141] S304. Determine the first reference model based on the initial velocity model, sub-reference data of at least one frequency band, and sub-seismic wavelet of at least one frequency band.
[0142] The first reference model can be determined as follows: the initial velocity model, the sub-reference data of the first frequency band, and the sub-seismic wavelet of the first frequency band are inverted to obtain the low-frequency model; the low-frequency model, the sub-reference data of the second frequency band, and the sub-seismic wavelet of the second frequency band are inverted to obtain the mid-frequency model; the mid-frequency model, the sub-reference data of the third frequency band, and the sub-seismic wavelet of the third frequency band are inverted to obtain the first reference model.
[0143] The process of determining the first benchmark model will be explained below with reference to Figure 4 and a specific example.
[0144] Figure 4 is a schematic diagram of the process for determining the first reference model provided in an embodiment of this application. Referring to Figure 4, the first reference model m can be obtained by inverting the initial velocity model, sub-reference data of at least one frequency band, and sub-seismic wavelet of at least one frequency band. b .
[0145] S305. Based on the first benchmark model, monitoring data, and benchmark data, determine the pseudo-observation data and seismic inversion wavelet.
[0146] Pseudo-observation data can be determined as follows: obtain a preset seismic wavelet from a preset database; perform forward modeling on the first reference model and the preset seismic wavelet to obtain simulated data; perform differential processing on the monitoring data and the reference data to obtain differential data; and sum the simulated data and the differential data to obtain pseudo-observation data.
[0147] Among them, the preset seismic wavelet can refer to the seismic wavelet under idealized conditions. The preset seismic wavelet can be set by the user in advance and stored in the preset database of the terminal device.
[0148] Forward modeling refers to generating theoretical observation data by simulating the propagation of seismic waves based on a known geological model. Terminal devices can be equipped with a forward modeling module. For example, forward modeling can be performed on a first reference model and a preset seismic wavelet to obtain simulated data.
[0149] Differential processing refers to subtracting the monitoring data from the baseline data to obtain differential data.
[0150] Summation processing can refer to summing simulated data and differenced data to obtain pseudo-observation data.
[0151] For example, suppose the simulation data is u baseline The monitoring data is d monitor The simulated data is d baseline Then it can be determined that the difference data can be represented as d monitor -d baseline , pseudo-observation data d composite It can be represented as:
[0152] d composite =d monitor -d baseline +u baseline
[0153] Among them, differential processing and summation processing require interpolating the monitoring data and the benchmark data to the same position before processing.
[0154] The seismic inversion wavelet can be determined as follows: obtain pseudo-observation data, perform inversion processing on the pseudo-observation data, and thus obtain the seismic inversion wavelet.
[0155] S306. Perform bandpass filtering on the pseudo-observation data to obtain sub-pseudo-observation data of at least one frequency band, and perform bandpass filtering on the seismic inversion wavelet to obtain sub-seismic inversion wavelet of at least one frequency band.
[0156] Sub-pseudo-observation data can be determined as follows: determine the frequency ranges of the first, second, and third frequency bands, and perform bandpass filtering on the pseudo-observation data according to the frequency ranges of the first, second, and third frequency bands to obtain the sub-pseudo-observation data of the first, second, and third frequency bands, respectively.
[0157] The sub-seismic inversion wavelet can be determined as follows: determine the frequency ranges of the first, second, and third frequency bands, and perform bandpass filtering on the seismic inversion wavelet according to the frequency ranges of the first, second, and third frequency bands to obtain the sub-seismic inversion wavelet of the first frequency band, the sub-seismic inversion wavelet of the second frequency band, and the sub-seismic inversion wavelet of the third frequency band, respectively.
[0158] S307. Determine the first monitoring model based on the first benchmark model, the pseudo-observation data of at least one frequency band, and the seismic inversion wavelet of at least one frequency band.
[0159] The first monitoring model can be determined as follows: Inversion processing is performed on the first baseline model, the pseudo-observation data of the first frequency band, and the sub-seismic inversion wavelet of the first frequency band to obtain the low-frequency monitoring model; inversion processing is performed on the low-frequency monitoring model, the pseudo-observation data of the second frequency band, and the sub-seismic inversion wavelet of the second frequency band to obtain the mid-frequency monitoring model; inversion processing is performed on the mid-frequency monitoring model, the pseudo-observation data of the third frequency band, and the sub-seismic inversion wavelet of the third frequency band to obtain the first monitoring model.
[0160] S308. Perform difference processing on the first monitoring model and the first benchmark model to obtain the first geological change information.
[0161] Differential processing can refer to subtracting the first monitoring model from the first benchmark model to obtain the first geological change information.
[0162] The process of determining the first geological change information will be explained below with reference to Figure 5 and through specific examples.
[0163] Figure 5 is a schematic diagram of another process for determining the first geological change information provided in an embodiment of this application. Referring to Figure 5, by performing forward modeling on the first reference model and the preset seismic wavelet, the simulation data u can be obtained. baseline , for monitoring data d monitor and benchmark data d baseline By performing difference processing, we can obtain the difference data d. monitor -d baseline For the simulated data u baseline By summing the summation of the difference data, we can obtain the pseudo-observation data d. compositeBy inverting the first benchmark model, the pseudo-observation data of at least one frequency band, and the seismic inversion wavelet of at least one frequency band, the first monitoring model m′ can be obtained. m For the first benchmark model m b and the first monitoring model m′ m By performing interpolation, the first geological change information Δm can be obtained. ＝> .
[0164] In this embodiment of the application, when it is necessary to determine the change information of the target geology from a historical moment to the most recent moment, the first geological change information can be determined first. Based on the data processing request, the following steps can be taken: Geological data of the target geology can be determined; grid tomography processing can be performed on the first seismic data to obtain an initial velocity model; bandpass filtering can be performed on the reference data to obtain sub-reference data for at least one frequency band, and bandpass filtering can be performed on the first seismic wavelet to obtain sub-seismic wavelets for at least one frequency band; a first reference model can be determined based on the initial velocity model, the sub-reference data for at least one frequency band, and the sub-seismic wavelets for at least one frequency band; pseudo-observation data and seismic inversion wavelets can be determined based on the first reference model, monitoring data, and reference data; bandpass filtering can be performed on the pseudo-observation data to obtain sub-pseudo-observation data for at least one frequency band, and bandpass filtering can be performed on the seismic inversion wavelet to obtain sub-seismic inversion wavelets for at least one frequency band; a first monitoring model can be determined based on the first reference model, the sub-pseudo-observation data for at least one frequency band, and the sub-seismic inversion wavelets for at least one frequency band; difference processing can be performed on the first monitoring model and the first reference model to obtain the first geological change information, thus improving the accuracy of the geological change information.
[0165] Based on any of the above embodiments, the process of determining the second geological change information will be described in detail below with reference to Figure 6.
[0166] Figure 6 is a schematic diagram illustrating the process of determining the second geological change information provided in an embodiment of this application. Referring to Figure 6, the method may include:
[0167] S601. Based on the data processing request, determine the geological data of the target geology.
[0168] S602. Perform grid tomography on the second seismic data to obtain the second initial velocity model.
[0169] S603. Perform bandpass filtering on the monitoring data to obtain sub-monitoring data of at least one frequency band, and perform bandpass filtering on the second seismic wavelet to obtain the second sub-seismic wavelet of at least one frequency band.
[0170] S604. Determine the second monitoring model based on the second velocity model, sub-monitoring data of at least one frequency band, and the second sub-seismic wavelet of at least one frequency band.
[0171] The process of determining the second monitoring model will be explained below with reference to Figure 7 and a specific example.
[0172] Figure 7 is a schematic diagram of the determination process of the second monitoring model provided in the embodiment of this application. Referring to Figure 7, the second monitoring model can be obtained by performing inversion processing on the second initial velocity model, the second sub-reference data of at least one frequency band, and the second sub-seismic wavelet of at least one frequency band.
[0173] S605. Based on the second monitoring model, monitoring data, and benchmark data, determine the second pseudo-observation data and the second seismic inversion wavelet.
[0174] S606. Bandpass filtering is applied to the second pseudo-observation data to obtain at least one frequency band of second sub-pseudo-observation data, and bandpass filtering is applied to the second seismic inversion wavelet to obtain at least one frequency band of second sub-seismic inversion wavelet.
[0175] S607. Determine the second benchmark model based on the second monitoring model, the second sub-pseudo-observation data of at least one frequency band, and the second sub-seismic inversion wavelet of at least one frequency band.
[0176] S608. Perform differential processing on the second monitoring model and the second benchmark model to obtain the second geological change information.
[0177] The process of determining the second geological change information will be explained below with reference to Figure 8 and through specific examples.
[0178] Figure 8 is a schematic diagram of another process for determining the second geological change information provided in an embodiment of this application. Referring to Figure 8, by performing forward modeling on the second monitoring model and the preset seismic wavelet, the simulated data u can be obtained. monitor , for monitoring data d monitor and benchmark data d baseline By performing difference processing, the second difference data d can be obtained. baseline -d monitor For the simulated data u monitor By summing the second difference data, we can obtain the second pseudo-observation data d. composite By inverting the second monitoring model, the pseudo-observation data of at least one frequency band, and the seismic inversion wavelet of at least one frequency band, a second benchmark model can be obtained. By performing difference processing on the second benchmark model and the second monitoring model, the second geological change information can be obtained.
[0179] The method for determining second geological change information provided in this application embodiment can execute the technical solution shown in the above method embodiment. Its implementation principle and beneficial effects are similar, and will not be described again here.
[0180] The data processing effects of the embodiments of this application will be explained below with reference to Figures 9-15 through specific examples.
[0181] Figure 9 is a schematic diagram of the two-dimensional model of the baseline data, Figure 10 is a schematic diagram of the two-dimensional model of the initial velocity model, Figure 11 is a schematic diagram of the two-dimensional model of the first baseline model, Figure 12 is a schematic diagram of the two-dimensional model of the first geological change information, Figure 13 is a schematic diagram of the two-dimensional model of the second geological change information, Figure 14 is a schematic diagram of the two-dimensional model of the target geological change information, and Figure 15 is a schematic diagram of the two-dimensional model of the actual target geological change information. Please refer to Figures 9-15. The first geological change information obtained based on the forward double difference inversion (Figure 12) and the second geological change information obtained based on the reverse double difference inversion (Figure 13) have a large error compared with the actual target geological change information. However, after averaging the first geological change information and the second geological change information, the target geological change information (Figure 14) is obtained, which is very close to the actual target geological change information (Figure 15). It can be seen that after data processing according to the method provided in the embodiments of this application, a more accurate time-shift change of the underground structure can be obtained.
[0182] Figure 16 is a schematic diagram of a geological change information determination device provided in an embodiment of this application. Referring to Figure 16, the geological change information determination device 10 includes: a first acquisition module 11, a second acquisition module 12, a first inversion processing module 13, a second inversion processing module 14, and a determination module 15, wherein...
[0183] The first acquisition module 11 is used to acquire a data processing request, wherein the data processing request is used to request the acquisition of geological change information of the target geology;
[0184] The second acquisition module 12 is used to acquire geological data of the target geology according to the data processing request. The geological data includes the baseline data of the target geology at a historical time, the monitoring data at the most recent time, the first earthquake data at the historical time, and the second earthquake data at the most recent time.
[0185] The first inversion processing module 13 is used to perform inversion processing on the reference data to obtain a first seismic wavelet, perform inversion processing on the reference data based on the first seismic wavelet and the geological data, and perform inversion processing on the monitoring data to obtain the first geological change information of the target geology;
[0186] The second inversion processing module 14 is used to perform inversion processing on the monitoring data to obtain a second seismic wavelet, perform inversion processing on the monitoring data based on the second seismic wavelet and the geological data, and perform inversion processing on the reference data to obtain the second geological change information of the target geology;
[0187] The determining module 15 is used to determine the target geological change information of the target geology based on the first geological change information and the second geological change information.
[0188] The geological change information determination device provided in this application embodiment can execute the technical solution shown in the above method embodiment. Its implementation principle and beneficial effects are similar, and will not be described again here.
[0189] In one possible design, the first inversion processing module 13 is specifically used for:
[0190] Based on the baseline data, the first seismic wavelet, and the first seismic data, a first baseline model for the target geology is determined;
[0191] Based on the baseline data, the monitoring data, the second seismic data, and the first baseline model, a first monitoring model for the target geology is determined.
[0192] The first geological change information is obtained by performing difference processing on the first monitoring model and the first benchmark model.
[0193] In one possible design, the first inversion processing module 13 is specifically used for:
[0194] The first seismic data is subjected to grid tomography to obtain an initial velocity model;
[0195] The reference data is bandpass filtered to obtain sub-reference data of at least one frequency band, and the first seismic wavelet is bandpass filtered to obtain sub-seismic wavelets of at least one frequency band.
[0196] The first reference model is determined based on the initial velocity model, the sub-reference data of the at least one frequency band, and the sub-seismic wavelet of the at least one frequency band.
[0197] In one possible design, the first inversion processing module 13 is specifically used for:
[0198] The initial velocity model, the sub-reference data of the first frequency band, and the sub-seismic wavelet of the first frequency band are inverted to obtain the low-frequency model;
[0199] The low-frequency model, the sub-reference data of the second frequency band, and the sub-seismic wavelet of the second frequency band are inverted to obtain the mid-frequency model;
[0200] The first reference model is obtained by inverting the mid-frequency model, the sub-reference data of the third frequency band, and the sub-seismic wavelet of the third frequency band.
[0201] In one possible design, the first inversion processing module 13 is specifically used for:
[0202] Based on the first benchmark model, the monitoring data, and the benchmark data, the pseudo-observation data and the seismic inversion wavelet are determined;
[0203] The pseudo-observation data is bandpass filtered to obtain sub-pseudo-observation data of at least one frequency band, and the seismic inversion wavelet is bandpass filtered to obtain sub-seismic inversion wavelet of at least one frequency band.
[0204] The first monitoring model is determined based on the first benchmark model, the sub-pseudo-observation data of the at least one frequency band, and the sub-seismic inversion wavelet of the at least one frequency band.
[0205] In one possible design, the first inversion processing module 13 is specifically used for:
[0206] Retrieve the preset seismic wavelet from the preset database;
[0207] The first reference model and the preset seismic wavelet are subjected to forward modeling to obtain simulation data;
[0208] The monitoring data and the benchmark data are differentially processed to obtain differential data;
[0209] The simulated data and the difference data are summed to obtain the pseudo-observation data;
[0210] The pseudo-observation data is inverted to obtain the seismic inversion wavelet.
[0211] In one possible design, the second inversion processing module 14 is specifically used for:
[0212] Based on the monitoring data, the second seismic wavelet, and the second seismic data, a second monitoring model for the target geology is determined;
[0213] Based on the baseline data, the monitoring data, the first seismic data, and the second monitoring model, a second baseline model for the target geology is determined.
[0214] The second geological change information is obtained by performing difference processing on the second monitoring model and the second benchmark model.
[0215] The geological change information determination device provided in this application embodiment can execute the technical solution shown in the above method embodiment. Its implementation principle and beneficial effects are similar, and will not be described again here.
[0216] Figure 17 is a schematic diagram of the structure of an electronic device provided in an embodiment of this application. As shown in Figure 17, the electronic device 20 may include: a transceiver 21, a processor 22, and a memory 23.
[0217] Processor 22 executes computer execution instructions stored in memory, causing processor 22 to perform the scheme in the above embodiments. Processor 22 can be a general-purpose processor, including a central processing unit (CPU), a network processor (NP), etc.; it can also be a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components.
[0218] The memory 23 is connected to the processor 22 via the system bus and completes communication between them. The memory 23 is used to store computer program instructions.
[0219] Transceiver 21 can be used to obtain the task to be run and the configuration information of the task to be run.
[0220] The system bus can be a Peripheral Component Interconnect (PCI) bus or an Extended Industry Standard Architecture (EISA) bus, etc. The system bus can be divided into address bus, data bus, control bus, etc. For ease of representation, only one thick line is used in the diagram, but this does not indicate that there is only one bus or one type of bus. Transceivers are used to enable communication between database access devices and other computers (e.g., clients, read-write libraries, and read-only libraries). Memory may include random access memory (RAM) and may also include non-volatile memory.
[0221] The electronic device provided in this application embodiment can be the terminal device described in the above embodiments.
[0222] This application also provides a chip for executing instructions, which is used to execute the technical solution of the method for determining geological change information in the above embodiments.
[0223] This application also provides a computer-readable storage medium storing computer instructions. When the computer instructions are executed on a computer, the computer performs the technical solution of the method for determining geological change information described in the above embodiments.
[0224] This application also provides a computer program product, which includes a computer program stored in a computer-readable storage medium. At least one processor can read the computer program from the computer-readable storage medium. When the at least one processor executes the computer program, it can implement the technical solution of the method for determining geological change information in the above embodiments.
[0225] In the several embodiments provided in this application, it should be understood that the disclosed devices and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative; for instance, the division of modules is only a logical functional division, and in actual implementation, there may be other division methods. For example, multiple modules may be combined or integrated into another system, or some features may be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection shown or discussed may be indirect coupling or communication connection through some interfaces, devices, or modules, and may be electrical, mechanical, or other forms.
[0226] The modules described as separate components may or may not be physically separate. The components shown as modules may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the modules can be selected to implement the solution of this embodiment according to actual needs.
[0227] Furthermore, the functional modules in the various embodiments of this application can be integrated into one processing unit, or each module can exist physically separately, or two or more modules can be integrated into one unit. The unit composed of the above modules can be implemented in hardware or in the form of hardware plus software functional units.
[0228] The integrated modules described above, implemented as software functional modules, can be stored in a computer-readable storage medium. These software functional modules, stored in a storage medium, include several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) or processor to execute some steps of the methods of the various embodiments of this application.
[0229] It should be understood that the aforementioned processor can be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), etc. A general-purpose processor can be a microprocessor or any conventional processor. The steps of the method disclosed in this invention can be directly manifested as being executed by a hardware processor, or executed by a combination of hardware and software modules within the processor.
[0230] The memory may include high-speed RAM, and may also include non-volatile storage (NVM), such as at least one disk storage device, and may also be a USB flash drive, external hard drive, read-only memory, disk or optical disc, etc.
[0231] The bus can be an Industry Standard Architecture (ISA) bus, a Peripheral Component Interconnect (PCI) bus, or an Extended Industry Standard Architecture (EISA) bus, etc. Buses can be categorized as address buses, data buses, control buses, etc. For ease of illustration, the buses shown in the accompanying drawings are not limited to a single bus or a single type of bus.
[0232] The aforementioned storage medium can be implemented from any type of volatile or non-volatile storage device or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic storage, flash memory, magnetic disk, or optical disk. The storage medium can be any available medium accessible to general-purpose or special-purpose computers.
[0233] An exemplary storage medium is coupled to a processor, enabling the processor to read information from and write information to the storage medium. Alternatively, the storage medium can be an integral part of the processor. The processor and storage medium can reside in an Application Specific Integrated Circuit (ASIC). Alternatively, the processor and storage medium can exist as discrete components in an electronic control unit or main control device.
[0234] Those skilled in the art will understand that all or part of the steps of the above-described method embodiments can be implemented by hardware related to program instructions. The aforementioned program can be stored in a computer-readable storage medium. When executed, the program performs the steps of the above-described method embodiments; and the aforementioned storage medium includes various media capable of storing program code, such as ROM, RAM, magnetic disks, or optical disks.
[0235] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of this application, and are not intended to limit them. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some or all of the technical features therein. Such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the scope of the technical solutions of the embodiments of this application.
Claims
1. A method for determining geological change information, characterized in that, include: Obtain a data processing request, wherein the data processing request is used to request the acquisition of geological change information of the target geology; According to the data processing request, the geological data of the target geology is obtained. The geological data includes the baseline data of the target geology at a historical time, the monitoring data at the most recent time, the first earthquake data at the historical time, and the second earthquake data at the most recent time. The reference data is inverted to obtain a first seismic wavelet. Based on the first seismic wavelet and the geological data, the reference data is inverted, and the monitoring data is inverted to obtain the first geological change information of the target geology. The monitoring data is inverted to obtain a second seismic wavelet. Based on the second seismic wavelet and the geological data, the monitoring data is inverted, and the reference data is inverted to obtain the second geological change information of the target geology. Based on the first geological change information and the second geological change information, the target geological change information of the target geology is determined.
2. The method according to claim 1, characterized in that, Based on the first seismic wavelet and the geological data, the reference data is inverted, and the monitoring data is inverted to obtain the first geological change information of the target geology, including: Based on the baseline data, the first seismic wavelet, and the first seismic data, a first baseline model for the target geology is determined; Based on the baseline data, the monitoring data, the second seismic data, and the first baseline model, a first monitoring model for the target geology is determined. The first geological change information is obtained by performing difference processing on the first monitoring model and the first benchmark model.
3. The method according to claim 2, characterized in that, Based on the baseline data, the first seismic wavelet, and the first seismic data, a first baseline model for the target geology is determined, including: The first seismic data is subjected to grid tomography to obtain an initial velocity model; The reference data is bandpass filtered to obtain sub-reference data of at least one frequency band, and the first seismic wavelet is bandpass filtered to obtain sub-seismic wavelets of at least one frequency band. The first reference model is determined based on the initial velocity model, the sub-reference data of the at least one frequency band, and the sub-seismic wavelet of the at least one frequency band.
4. The method according to claim 3, characterized in that, The at least one frequency band includes a first frequency band, a second frequency band, and a third frequency band; the first reference model is determined based on the initial velocity model, the sub-reference data of the at least one frequency band, and the sub-seismic wavelet of the at least one frequency band, including: The initial velocity model, the sub-reference data of the first frequency band, and the sub-seismic wavelet of the first frequency band are inverted to obtain the low-frequency model; The low-frequency model, the sub-reference data of the second frequency band, and the sub-seismic wavelet of the second frequency band are inverted to obtain the mid-frequency model; The first reference model is obtained by inverting the mid-frequency model, the sub-reference data of the third frequency band, and the sub-seismic wavelet of the third frequency band.
5. The method according to any one of claims 2-4, characterized in that, Based on the baseline data, the monitoring data, the second seismic data, and the first baseline model, a first monitoring model for the target geology is determined, including: Based on the first benchmark model, the monitoring data, and the benchmark data, the pseudo-observation data and the seismic inversion wavelet are determined; The pseudo-observation data is bandpass filtered to obtain sub-pseudo-observation data of at least one frequency band, and the seismic inversion wavelet is bandpass filtered to obtain sub-seismic inversion wavelet of at least one frequency band. The first monitoring model is determined based on the first benchmark model, the sub-pseudo-observation data of the at least one frequency band, and the sub-seismic inversion wavelet of the at least one frequency band.
6. The method according to claim 5, characterized in that, Based on the first benchmark model, the monitoring data, and the benchmark data, the pseudo-observation data and seismic inversion wavelet are determined, including: Retrieve the preset seismic wavelet from the preset database; The first reference model and the preset seismic wavelet are subjected to forward modeling to obtain simulation data; The monitoring data and the benchmark data are differentially processed to obtain differential data; The simulated data and the difference data are summed to obtain the pseudo-observation data; The pseudo-observation data is inverted to obtain the seismic inversion wavelet.
7. The method according to any one of claims 1-6, characterized in that, Based on the second seismic wavelet and the geological data, the monitoring data is inverted, and the reference data is inverted to obtain the second geological change information of the target geology, including: Based on the monitoring data, the second seismic wavelet, and the second seismic data, a second monitoring model for the target geology is determined; Based on the baseline data, the monitoring data, the first seismic data, and the second monitoring model, a second baseline model for the target geology is determined. The second geological change information is obtained by performing difference processing on the second monitoring model and the second benchmark model.
8. A device for determining geological change information, characterized in that, include: The system comprises a first acquisition module, a second acquisition module, a first inversion processing module, a second inversion processing module, and a determination module, wherein... The first acquisition module is used to acquire a data processing request, wherein the data processing request is used to request the acquisition of geological change information of the target geology; The second acquisition module is used to acquire geological data of the target geology according to the data processing request. The geological data includes baseline data of the target geology at a historical time, monitoring data at a recent time, first earthquake data at the historical time, and second earthquake data at the recent time. The first inversion processing module is used to perform inversion processing on the reference data to obtain a first seismic wavelet, perform inversion processing on the reference data based on the first seismic wavelet and the geological data, and perform inversion processing on the monitoring data to obtain the first geological change information of the target geology; The second inversion processing module is used to perform inversion processing on the monitoring data to obtain a second seismic wavelet, and to perform inversion processing on the monitoring data and the reference data based on the second seismic wavelet and the geological data to obtain the second geological change information of the target geology. The determining module is used to determine the target geological change information of the target geology based on the first geological change information and the second geological change information.
9. An electronic device, characterized in that, include: A processor, and a memory communicatively connected to the processor; The memory stores computer-executed instructions; The processor executes computer execution instructions stored in the memory to implement the method as described in any one of claims 1 to 7.
10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer-executable instructions, which, when executed by a processor, are used to implement the method as described in any one of claims 1 to 7.
11. A computer program product, characterized in that, Includes a computer program that, when executed by a processor, implements the method of any one of claims 1 to 7.