Depth domain seismic wavelet extraction method and device, storage medium and computer device
By applying time-depth conversion and convolution model to the seismic traces near the well, the problems of error and information loss in depth domain seismic wavelet extraction were solved, and more accurate depth domain seismic wavelet extraction was achieved.
Patent Information
- Application Number
- CN202011033989.8
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2020-09-27
- Publication Date
- 2025-12-05
- Estimated Expiration
- 2040-09-27
AI Technical Summary
In existing technologies, methods for extracting depth-domain seismic wavelets cannot directly utilize time-domain convolution models, leading to errors introduced by time shifts in the reflection coefficients and loss of detailed information, thus failing to accurately extract depth-domain seismic wavelets.
By converting the depth domain data to the time domain through the seismic traces near the well, the time domain seismic wavelet is extracted using a convolution model. The time interval is ensured to be uniform by interpolation. Then, the time domain data is converted to the depth domain to achieve accurate extraction of the depth domain seismic wavelet.
The method achieves accurate extraction of seismic wavelets in the depth domain, ensuring the correspondence between time and depth, reducing errors and information loss, and obtaining more accurate depth-domain seismic wavelet data.
Smart Images
Figure CN114428299B_ABST
Abstract
Description
Technical Field
[0001] This invention relates to the field of geophysical exploration technology in the petrochemical industry, and in particular to a method, apparatus, storage medium, and computer equipment for extracting depth-domain seismic wavelets. Background Technology
[0002] Currently, most conventional seismic data inversion methods are performed in the time domain, and the wavelets required for inversion are obtained based on time-domain convolution models. Convolution models are widely used in inversion processes due to their simplicity and applicability. In the time domain, the most common approach is to convolve the time-domain wavelet with the time-domain reflection coefficient to obtain a synthetic seismic record. The time-domain reflection coefficient is typically obtained from depth-domain well logging data after depth-time conversion. During the conversion of depth-domain well logging data to the time domain, stretching or compression of the logging curves in the time domain is unavoidable, inevitably causing a time shift in the reflection coefficients and introducing errors. Furthermore, the well logging curves after time-depth conversion may lose much important detailed information. Since the depth domain does not possess the condition of "linear time invariance," convolution models cannot be directly used. Therefore, how to accurately extract depth-domain seismic wavelets using convolution models is a crucial problem that urgently needs to be solved. Summary of the Invention
[0003] The main objective of this invention is to provide a method, apparatus, storage medium, and computer device for extracting depth domain seismic wavelets, so as to accurately extract depth domain seismic wavelets.
[0004] In a first aspect, this application provides a method for extracting depth-domain seismic wavelets, comprising the following steps: determining time-domain data of a well-side seismic trace based on depth-domain data of the well-side seismic trace; in the time domain, using a convolution model, determining the corresponding time-domain seismic wavelet based on the time-domain data of the well-side seismic trace; and determining the corresponding depth-domain seismic wavelet based on the time-domain seismic wavelet.
[0005] In one embodiment, the depth domain data of the well-side seismic trace includes multiple sampled data of the depth domain P-wave velocity of the well-side seismic trace; correspondingly, the time domain data of the well-side seismic trace includes the time domain P-wave velocity data of the well-side seismic trace.
[0006] In one embodiment, determining the time-domain data of the well-side seismic trace based on the depth-domain data of the well-side seismic trace includes: for each sampled data point of the depth-domain P-wave velocity of the well-side seismic trace, determining the sampling time of each sampled data point using the relationship between P-wave velocity, depth, and time; establishing a correspondence between P-wave velocity and sampling time in each sampled data point based on the sampling time of each sampled data point; and determining the correspondence between P-wave velocity and equally spaced time intervals based on the correspondence between P-wave velocity and sampling time in multiple sampled data points, thereby obtaining the time-domain P-wave velocity data of the well-side seismic trace.
[0007] In one embodiment, determining the correspondence between P-wave velocity and equal time intervals based on the correspondence between P-wave velocity and sampling time in multiple sampled data includes: determining the correspondence between P-wave velocity and equal time intervals using an interpolation method based on the correspondence between P-wave velocity and sampling time in multiple sampled data.
[0008] In one embodiment, the convolution model is:
[0009] s(t)=w(t)*r(t)
[0010] Where s(t) represents the time-domain data of the seismic trace near the well, r(t) represents the reflection coefficient sequence, and w(t) represents the seismic wavelet.
[0011] In one embodiment, the reflection coefficient sequence is determined by the following steps: determining the reflection coefficient of each formation in the target reservoir based on the density and P-wave velocity of each formation in the target reservoir, and determining the reflection coefficient sequence based on the reflection coefficient of each formation in the target reservoir.
[0012] In one embodiment, determining the corresponding depth-domain seismic wavelet based on the time-domain seismic wavelet includes: determining the depth reached by the P-wave at each time point using the relationship between P-wave velocity, depth, and time based on the time-domain seismic wavelet; establishing a correspondence between the P-wave velocity and the depth reached by the P-wave at each time point based on the depth reached by the P-wave at each time point; and determining a correspondence between the P-wave velocity and equally spaced depths based on the correspondence between the P-wave velocity and the depth reached by the P-wave at multiple time points, thereby obtaining the corresponding depth-domain seismic wavelet.
[0013] Secondly, this application provides a depth-domain seismic wavelet extraction device, comprising: a depth-time transformation module for determining the time-domain data of a well-side seismic trace based on the depth-domain data of the well-side seismic trace; a convolution operation module for determining the corresponding time-domain seismic wavelet based on the time-domain data of the well-side seismic trace using a convolution model in the time domain; and a time-depth transformation module for determining the corresponding depth-domain seismic wavelet based on the time-domain seismic wavelet.
[0014] Thirdly, this application provides a storage medium storing a computer program, characterized in that, when the computer program is executed by a processor, it implements the steps of the depth domain seismic wavelet extraction method described above.
[0015] Fourthly, this application provides a computer device, including a processor and a storage medium storing program code, wherein when the program code is executed by the processor, it implements the steps of the depth domain seismic wavelet extraction method described above.
[0016] This invention employs a time-depth conversion method for well-side seismic traces, converting depth-domain data from well-side seismic traces to the time domain. In the time domain, a convolution model is used to obtain the time-domain wavelet, and then the time-domain seismic wavelet is converted back to the depth domain, thereby obtaining a more accurate depth-domain time-invariant wavelet. Attached Figure Description
[0017] The accompanying drawings, which form part of this application, are used to provide a further understanding of the invention. The illustrative embodiments of the invention and their descriptions are used to explain the invention and do not constitute an undue limitation of the invention. In the drawings:
[0018] Figure 1 A flowchart of a method for extracting depth-domain seismic wavelets according to an exemplary embodiment of this application;
[0019] Figure 2 This is a flowchart of a method for extracting depth-domain seismic wavelets according to a specific embodiment of this application. Detailed Implementation
[0020] It should be noted that, unless otherwise specified, the embodiments and features described in this application can be combined with each other. The present invention will now be described in detail with reference to the accompanying drawings and embodiments.
[0021] Example 1
[0022] This embodiment provides a method for extracting seismic wavelets in the depth domain. Figure 1 This is a flowchart of a method for extracting depth-domain seismic wavelets according to an exemplary embodiment of this application. Figure 1 As shown, the method includes the following steps:
[0023] S100: Determine the time domain data of the seismic traces near the well based on the depth domain data of the seismic traces near the well.
[0024] S200: In the time domain, using the convolution model, the corresponding time domain seismic wavelet is determined based on the time domain data of the well-side seismic trace.
[0025] S300: Determine the corresponding depth-domain seismic wavelet based on the seismic wavelet in the time domain.
[0026] Example 2
[0027] This embodiment provides a method for extracting depth-domain seismic wavelets, including the following steps:
[0028] The first step is to determine the time domain data of the well-side seismic traces based on the depth domain data.
[0029] The depth domain data of the well-side seismic trace includes multiple sampled data of the depth domain P-wave velocity of the well-side seismic trace; correspondingly, the time domain data of the well-side seismic trace includes the time domain P-wave velocity data of the well-side seismic trace.
[0030] When the depth domain data of the well-side seismic trace includes multiple sampled data of the depth domain P-wave velocity of the well-side seismic trace, and the time domain data of the well-side seismic trace includes the time domain P-wave velocity data of the well-side seismic trace, determining the time domain data of the well-side seismic trace based on the depth domain data of the well-side seismic trace may include:
[0031] 1) For each sampled data of P-wave velocity in the depth domain of the seismic trace near the well, the sampling time of each sampled data is determined by using the relationship between P-wave velocity, depth and time, and the correspondence between P-wave velocity and sampling time in each sampled data is established based on the sampling time of each sampled data.
[0032] 2) Based on the correspondence between P-wave velocity and sampling time in multiple sampled data, determine the correspondence between P-wave velocity and equal time intervals, thereby obtaining the time-domain P-wave velocity data of the well-side seismic trace. Specifically, determining the correspondence between P-wave velocity and equal time intervals based on the correspondence between P-wave velocity and sampling time in multiple sampled data includes: using interpolation to determine the correspondence between P-wave velocity and equal time intervals based on the correspondence between P-wave velocity and sampling time in multiple sampled data.
[0033] The second step is to use a convolution model in the time domain to determine the corresponding seismic wavelet in the time domain based on the time domain data of the seismic traces near the well.
[0034] The convolution model can be:
[0035] s(t)=w(t)*r(t)
[0036] Wherein, s(t) represents the time-domain data of the well-side seismic trace, w(t) represents the seismic wavelet, and r(t) represents the reflection coefficient sequence. The reflection coefficient sequence can be determined by the following method: determining the reflection coefficient of each formation in the target reservoir based on the density and P-wave velocity of each formation in the target reservoir, and determining the reflection coefficient sequence based on the reflection coefficient of each formation in the target reservoir.
[0037] The third step is to determine the corresponding depth-domain seismic wavelet based on the time-domain seismic wavelet. Specifically, this may include:
[0038] 1) Based on the seismic wavelet in the time domain, the depth reached by the P-wave at each time point is determined by utilizing the relationship between P-wave velocity, depth and time. Based on the depth reached by the P-wave at each time point, the correspondence between the P-wave velocity and the depth reached by the P-wave at each time point is established.
[0039] 2) Based on the correspondence between P-wave velocity and depth reached by P-wave at multiple time points, determine the correspondence between P-wave velocity and equally spaced depths, thereby obtaining the corresponding seismic wavelet in the depth domain.
[0040] This invention employs a time-depth conversion method for well-side seismic traces. It converts the depth-domain data of well-side seismic traces into the time domain, uses a convolution model to obtain the time-domain wavelet, and then converts the time-domain seismic wavelet back into the depth domain, thus obtaining a more accurate depth-domain time-invariant wavelet. In this invention, after converting the well-side seismic traces from the depth domain to the time domain, interpolation is used to ensure uniform time intervals in the time-domain seismic data, which facilitates obtaining a more accurate time-domain seismic wavelet through convolution. Similarly, after converting the obtained time-domain seismic wavelet back to the depth domain, interpolation is used to ensure equal depth intervals between data points in the depth-domain seismic wavelet, ultimately yielding a more accurate depth-domain seismic wavelet.
[0041] Example 3
[0042] This embodiment provides a specific implementation of a method for extracting seismic wavelets in the depth domain, which includes the following steps:
[0043] The first step is to calculate the reflection coefficient based on the well logging data:
[0044] s(t)=w(t)*r(t) (1)
[0045] Where s(t) is the synthetic seismic record, r(t) is the reflection coefficient sequence, and w(t) is the seismic wavelet;
[0046] The reflection coefficient r(t) of an interface is obtained from the wave impedance of the upper and lower layers, and its expression is:
[0047]
[0048] In the formula, R is the reflection coefficient, ρ1 and ρ2 are the densities of the upper and lower layers of the interface, and v1 and v2 are the velocities of the upper and lower layers of the interface. The densities and velocities of the upper and lower layers of the interface can be obtained from well logging data.
[0049] The second step is to replace the synthetic seismic record s(t) in expression (1) with the well-side seismic trace s(m). Since the well-side seismic trace is depth domain seismic data, it needs to be converted into time domain seismic data in this step, that is, s(m) is converted into s(t). The specific method is as follows: based on the velocity information v(m) of the well-side trace, where m is the depth domain in meters, it is interpolated to data with a spacing of 0.1m. Then, the expression t = m / v is used to calculate the time information of each sampling point, thereby obtaining the well-side seismic trace s(t) in the time domain. The obtained data must not be arranged at equal time intervals, so it needs to be interpolated to the same sampling interval as the seismic data. If the seismic data is sampled in 1ms, the obtained data s(t) is interpolated to 1ms sampling, thus obtaining the seismic data in the time domain.
[0050] The third step is to obtain w(t) = s(t) * r(t) based on expression (1). -1 Thus, the seismic wavelet in the time domain can be obtained;
[0051] The fourth step is to reverse the steps of the depth-time conversion in the second step. At this time, w(t) is sampled at 1ms. The depth information w(m) of each point is calculated using the formula m = v*t. Then, it is interpolated to data with a spacing of 0.1m to finally obtain the seismic wavelet in the depth domain.
[0052] This invention employs a time-depth conversion method for well-side seismic traces, converting depth-domain data from well-side seismic traces into the time domain. In the time domain, a convolution model is used to obtain the time-domain wavelet, and then the time-domain seismic wavelet is converted back into the depth domain, thereby obtaining a more accurate depth-domain time-invariant wavelet.
[0053] Example 4
[0054] This embodiment provides a depth-domain seismic wavelet extraction device, comprising: a depth-time transformation module for determining the time-domain data of a well-side seismic trace based on the depth-domain data of the well-side seismic trace; a convolution operation module for determining the corresponding time-domain seismic wavelet based on the time-domain data of the well-side seismic trace using a convolution model in the time domain; and a time-depth transformation module for determining the corresponding depth-domain seismic wavelet based on the time-domain seismic wavelet.
[0055] In this embodiment, the depth domain seismic wavelet extraction device may further include a processor and a memory, wherein the processor is used to execute the following program modules stored in the memory: a depth-time transformation module, a convolution operation module, and a time-depth transformation module, to implement the steps of the depth domain seismic wavelet extraction method as described above:
[0056] The time domain data of the well-side seismic traces are determined based on the depth domain data of the well-side seismic traces.
[0057] In the time domain, the convolution model is used to determine the corresponding time domain seismic wavelet based on the time domain data of the well-side seismic trace;
[0058] Based on the seismic wavelet in the time domain, the corresponding seismic wavelet in the depth domain is determined.
[0059] Example 5
[0060] This embodiment provides a storage medium storing a computer program. When the computer program is executed by a processor, it implements the steps of the depth domain seismic wavelet extraction method described above:
[0061] The time domain data of the well-side seismic traces are determined based on the depth domain data of the well-side seismic traces.
[0062] In the time domain, the convolution model is used to determine the corresponding time domain seismic wavelet based on the time domain data of the well-side seismic trace;
[0063] Based on the seismic wavelet in the time domain, the corresponding seismic wavelet in the depth domain is determined.
[0064] In one embodiment, the depth domain data of the well-side seismic trace includes multiple sampled data of the depth domain P-wave velocity of the well-side seismic trace; correspondingly, the time domain data of the well-side seismic trace includes the time domain P-wave velocity data of the well-side seismic trace.
[0065] In one embodiment, determining the time-domain data of the well-side seismic trace based on the depth-domain data of the well-side seismic trace includes: for each sampled data point of the depth-domain P-wave velocity of the well-side seismic trace, determining the sampling time of each sampled data point using the relationship between P-wave velocity, depth, and time; establishing a correspondence between P-wave velocity and sampling time in each sampled data point based on the sampling time of each sampled data point; and determining the correspondence between P-wave velocity and equally spaced time intervals based on the correspondence between P-wave velocity and sampling time in multiple sampled data points, thereby obtaining the time-domain P-wave velocity data of the well-side seismic trace.
[0066] In one embodiment, determining the correspondence between P-wave velocity and equal time intervals based on the correspondence between P-wave velocity and sampling time in multiple sampled data includes: determining the correspondence between P-wave velocity and equal time intervals using an interpolation method based on the correspondence between P-wave velocity and sampling time in multiple sampled data.
[0067] In one embodiment, the convolution model is:
[0068] s(t)=w(t)*r(t)
[0069] Where s(t) represents the time-domain data of the seismic trace near the well, r(t) represents the reflection coefficient sequence, and w(t) represents the seismic wavelet.
[0070] In one embodiment, the reflection coefficient sequence is determined by the following steps: determining the reflection coefficient of each formation in the target reservoir based on the density and P-wave velocity of each formation in the target reservoir, and determining the reflection coefficient sequence based on the reflection coefficient of each formation in the target reservoir.
[0071] In one embodiment, determining the corresponding depth-domain seismic wavelet based on the time-domain seismic wavelet includes: determining the depth reached by the P-wave at each time point using the relationship between P-wave velocity, depth, and time based on the time-domain seismic wavelet; establishing a correspondence between the P-wave velocity and the depth reached by the P-wave at each time point based on the depth reached by the P-wave at each time point; and determining a correspondence between the P-wave velocity and equally spaced depths based on the correspondence between the P-wave velocity and the depth reached by the P-wave at multiple time points, thereby obtaining the corresponding depth-domain seismic wavelet.
[0072] Storage media, including permanent and non-permanent, removable and non-removable media, can be used to store information by any method or technology. Information can be computer-readable instructions, data structures, program modules, or other data. Examples of storage media include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, CD-ROM, digital versatile optical disc (DVD) or other optical storage, magnetic tape, magnetic magnetic disk storage or other magnetic storage devices, or any other non-transfer medium that can be used to store information accessible by computing devices.
[0073] Example 6
[0074] This embodiment provides a computer device, including a processor and a storage medium storing program code. When the program code is executed by the processor, it implements the steps of the depth domain seismic wavelet extraction method described above:
[0075] The time domain data of the well-side seismic traces are determined based on the depth domain data of the well-side seismic traces.
[0076] In the time domain, the convolution model is used to determine the corresponding time domain seismic wavelet based on the time domain data of the well-side seismic trace;
[0077] Based on the seismic wavelet in the time domain, the corresponding seismic wavelet in the depth domain is determined.
[0078] In one embodiment, the depth domain data of the well-side seismic trace includes multiple sampled data of the depth domain P-wave velocity of the well-side seismic trace; correspondingly, the time domain data of the well-side seismic trace includes the time domain P-wave velocity data of the well-side seismic trace.
[0079] In one embodiment, determining the time-domain data of the well-side seismic trace based on the depth-domain data of the well-side seismic trace includes: for each sampled data point of the depth-domain P-wave velocity of the well-side seismic trace, determining the sampling time of each sampled data point using the relationship between P-wave velocity, depth, and time; establishing a correspondence between P-wave velocity and sampling time in each sampled data point based on the sampling time of each sampled data point; and determining the correspondence between P-wave velocity and equally spaced time intervals based on the correspondence between P-wave velocity and sampling time in multiple sampled data points, thereby obtaining the time-domain P-wave velocity data of the well-side seismic trace.
[0080] In one embodiment, determining the correspondence between P-wave velocity and equal time intervals based on the correspondence between P-wave velocity and sampling time in multiple sampled data includes: determining the correspondence between P-wave velocity and equal time intervals using an interpolation method based on the correspondence between P-wave velocity and sampling time in multiple sampled data.
[0081] In one embodiment, the convolution model is:
[0082] s(t)=w(t)*r(t)
[0083] Where s(t) represents the time-domain data of the seismic trace near the well, r(t) represents the reflection coefficient sequence, and w(t) represents the seismic wavelet.
[0084] In one embodiment, the reflection coefficient sequence is determined by the following steps: determining the reflection coefficient of each formation in the target reservoir based on the density and P-wave velocity of each formation in the target reservoir, and determining the reflection coefficient sequence based on the reflection coefficient of each formation in the target reservoir.
[0085] In one embodiment, determining the corresponding depth-domain seismic wavelet based on the time-domain seismic wavelet includes: determining the depth reached by the P-wave at each time point using the relationship between P-wave velocity, depth, and time based on the time-domain seismic wavelet; establishing a correspondence between the P-wave velocity and the depth reached by the P-wave at each time point based on the depth reached by the P-wave at each time point; and determining a correspondence between the P-wave velocity and equally spaced depths based on the correspondence between the P-wave velocity and the depth reached by the P-wave at multiple time points, thereby obtaining the corresponding depth-domain seismic wavelet.
[0086] In one embodiment, a computer device includes one or more processors (CPUs), input / output interfaces, network interfaces, and memory.
[0087] Memory may include non-persistent storage in computer-readable media, such as random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash memory (flash FLASH RAM). Memory is an example of computer-readable media.
[0088] Any equivalent structural or procedural transformations made based on the description and drawings of this invention, or direct or indirect applications in other related technical fields, are similarly included within the scope of patent protection of this invention.
[0089] Through the above description of the embodiments, those skilled in the art can clearly understand that the methods of the above embodiments can be implemented by means of software plus necessary general-purpose hardware platforms. Of course, they can also be implemented by hardware, but in many cases the former is a better implementation method. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, can be embodied in the form of a software product. This computer software product is stored in a readable storage medium (such as ROM / RAM, magnetic disk, optical disk) as described above, and includes several instructions to cause a system device (which may be a mobile phone, computer, server, air conditioner, or network device, etc.) to execute the methods described in the various embodiments of the present invention.
[0090] It should be noted that the terminology used herein is for the purpose of describing particular implementations only and is not intended to limit the exemplary implementations according to this application. When the terms “comprising” and / or “including” are used in this specification, they indicate the presence of features, steps, operations, devices, components and / or combinations thereof.
[0091] It should be noted that the terms "first," "second," etc., used in the specification, claims, and drawings of this application are used to distinguish similar objects and are not used to describe a specific order or sequence. It should be understood that such terms can be used interchangeably where appropriate so that the embodiments of this application described herein can be implemented, for example, in sequences other than those illustrated or described herein.
[0092] It should be understood that the exemplary embodiments described herein can be implemented in many different forms and should not be construed as being limited to the embodiments set forth herein. These embodiments are provided so that the disclosure of this application is thorough and complete, and that the concept of these exemplary embodiments is fully conveyed to those skilled in the art, and should not be construed as limiting the invention.
Claims
1. A method for extracting a seismic wavelet in a depth domain, characterized in that, The method comprises the following steps: determining time domain data of the well seismic trace according to the depth domain data of the well seismic trace; in the time domain, determining the corresponding time domain seismic wavelet according to the time domain data of the well seismic trace by using a convolution model; determining the corresponding depth domain seismic wavelet according to the time domain seismic wavelet; the depth domain data of the well seismic trace comprises a plurality of sampling data of the depth domain P-wave velocity of the well seismic trace; and the time domain data of the well seismic trace comprises time domain P-wave velocity data of the well seismic trace; the determination of the time domain data of the well seismic trace according to the depth domain data of the well seismic trace comprises: for each sampling data of the depth domain P-wave velocity of the well seismic trace, determining the sampling time of each sampling data by using the relationship among the P-wave velocity, the depth and the time, and establishing the corresponding relationship between the P-wave velocity and the sampling time in each sampling data according to the sampling time of each sampling data; determining the corresponding relationship between the P-wave velocity and the equally spaced time by using an interpolation method according to the corresponding relationship between the P-wave velocity and the sampling time in the plurality of sampling data, so as to obtain the time domain P-wave velocity data of the well seismic trace; wherein the interpolation method is used to make the time interval of the seismic data of the well seismic trace in the time domain uniform; the determination of the corresponding depth domain seismic wavelet according to the time domain seismic wavelet comprises: determining the depth reached by the P-wave at each time point by using the relationship among the P-wave velocity, the depth and the time according to the time domain seismic wavelet, and establishing the corresponding relationship between the P-wave velocity and the depth reached by the P-wave at each time point according to the depth reached by the P-wave at each time point; determining the corresponding relationship between the P-wave velocity and the equally spaced depth by using an interpolation method according to the corresponding relationship between the P-wave velocity and the depth reached by the P-wave at a plurality of time points, so as to obtain the corresponding depth domain seismic wavelet; wherein the interpolation method is used to make the depth interval between each data point of the depth domain seismic wavelet equal.
2. The method of claim 1, wherein, the convolution model is: wherein s(t) represents the time domain data of the well seismic trace, r(t) represents a reflection coefficient sequence, and w(t) represents a seismic wavelet.
3. The method of claim 2, wherein, the reflection coefficient sequence is determined by the following steps: determining the reflection coefficient of each stratum in the target reservoir according to the density and the P-wave velocity of each stratum in the target reservoir, and determining the reflection coefficient sequence according to the reflection coefficient of each stratum in the target reservoir.
4. The method of claim 1, wherein, the determination of the corresponding depth domain seismic wavelet according to the time domain seismic wavelet comprises: determining the depth reached by the P-wave at each time point by using the relationship among the P-wave velocity, the depth and the time according to the time domain seismic wavelet, and establishing the corresponding relationship between the P-wave velocity and the depth reached by the P-wave at each time point according to the depth reached by the P-wave at each time point; determining the corresponding relationship between the P-wave velocity and the equally spaced depth by using an interpolation method according to the corresponding relationship between the P-wave velocity and the depth reached by the P-wave at a plurality of time points, so as to obtain the corresponding depth domain seismic wavelet.
5. An apparatus for extracting a seismic wavelet in a depth domain, characterized by comprising: comprise: a depth-time conversion module, configured to determine time domain data of the well seismic trace according to the depth domain data of the well seismic trace; a convolution operation module, configured to determine a corresponding time-domain seismic wavelet from time-domain data of a well seismic trace in a time domain by using a convolution model; a time-depth conversion module, configured to determine a corresponding depth-domain seismic wavelet from the time-domain seismic wavelet; the depth-domain data of the well seismic trace include a plurality of sampling data of depth-domain P-wave velocity of the well seismic trace; correspondingly, the time-domain data of the well seismic trace include time-domain P-wave velocity data of the well seismic trace; determining the time-domain data of the well seismic trace from the depth-domain data of the well seismic trace includes: for each sampling data of the depth-domain P-wave velocity of the well seismic trace, determining a sampling time of each sampling data by using a relationship among P-wave velocity, depth and time, and establishing a corresponding relationship between P-wave velocity and sampling time in each sampling data according to the sampling time of each sampling data; determining a corresponding relationship between P-wave velocity and equally-spaced time according to the corresponding relationship between P-wave velocity and sampling time in the plurality of sampling data, thereby obtaining the time-domain P-wave velocity data of the well seismic trace; wherein an interpolation method is used to make time intervals of the time-domain seismic data of the well seismic trace uniform; the determining the corresponding depth-domain seismic wavelet from the time-domain seismic wavelet includes: determining a depth reached by P-wave at each time point by using the relationship among P-wave velocity, depth and time according to the time-domain seismic wavelet, and establishing a corresponding relationship between P-wave velocity and the depth reached by P-wave at each time point according to the depth reached by P-wave at each time point; determining a corresponding relationship between P-wave velocity and equally-spaced depth by using an interpolation method according to the corresponding relationship between P-wave velocity and the depth reached by P-wave at the plurality of time points, thereby obtaining the corresponding depth-domain seismic wavelet; wherein the interpolation method is used to make depth intervals between each data point of the depth-domain seismic wavelet equal.
6. A storage medium storing a computer program, characterized by the computer program, when executed by a processor, implements the steps of the extraction method of the depth-domain seismic wavelet according to any one of claims 1-4. 7.A computer device, comprising a processor and a storage medium storing program code, the program code, when executed by the processor, implements the steps of the extraction method of the depth-domain seismic wavelet according to any one of claims 1-4.
Citation Information
Patent Citations
Depth-domain synthetic seismic record determining method and apparatus
CN107193040A