DAS VSP wave field simulation method and device, electronic equipment and medium
By converting conventional displacement wavefield simulation into strain/strain rate wavefield simulation, the technical gap in DAS VSP wavefield simulation is filled, achieving high-precision DAS VSP wavefield simulation, which is suitable for DAS VSP wavefield data processing.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-10-09
- Publication Date
- 2026-04-10
AI Technical Summary
Lacking mature DAS VSP wavefield simulation technology and software, traditional methods cannot effectively simulate the phase delay between the difference frequency signals at scattering points in optical fibers recorded by DAS.
By converting conventional displacement wave field simulation into strain/strain rate wave field along the fiber direction, DAS VSP wave field simulation is realized. Using parameters such as well trajectory, P-wave and S-wave velocity model, source wavelet, source excitation point position, receiver position and DAS demodulation equipment gauge length, combined with finite difference or finite element elastic wave simulation methods, displacement wave field data is calculated and converted into DAS VSP wave field data.
A simple and effective DAS VSP wavefield simulation was achieved, taking into account the influence of the gauge length of the DAS demodulation equipment, thus improving the simulation accuracy and efficiency.
Smart Images
Figure CN121831874A_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The present application relates to the field of DAS VSP data processing, and more particularly, to a DAS VSP wave field simulation method and device, electronic equipment and medium. BACKGROUND
[0002] Distributed acoustic sensing (DAS) is a new signal acquisition technology that characterizes seismic signals by demodulating the phase change of Rayleigh backscattering light signals. It has the characteristics of one-time full well coverage and high-density measurement in the well, and the construction efficiency and data consistency are greatly improved, so it has attracted widespread attention. With the rapid development of DAS acquisition technology in the well, high-density and full-well DAS VSP data acquisition can be realized. DAS-VSP technology can provide high-precision reservoir parameters, improve the fine description ability of the oil and gas reservoirs around the oil and gas well, and has the ability of dynamic monitoring of oil and gas reservoirs.
[0003] DAS data has great differences with traditional electronic geophone data due to its special sensing principle and acquisition method. DAS records the phase delay between the difference frequency signals of the scattering points in the optical fiber, which is essentially the strain / strain rate signal along the optical fiber direction. The signal detected by the traditional digital geophone is the velocity or acceleration signal. The conventional VSP wave field simulation is based on elastic wave simulation technology such as finite difference and finite element, and the simulated wave field is the displacement or velocity wave field at the geophone point; there is currently a lack of mature DAS VSP wave field simulation technology and software.
[0004] There is a need to develop a DAS VSP wave field simulation method.
[0005] The information disclosed in the background section of the present application is only intended to deepen the understanding of the general background of the present application, and should not be regarded as acknowledging or implying in any form that the information constitutes prior art known to those skilled in the art. SUMMARY
[0006] The present application provides a DAS VSP wave field simulation method, device, electronic equipment and medium, which can realize simple and effective DAS VSP wave field simulation by converting the displacement wave field simulation into strain / strain rate wave field along the optical fiber direction.
[0007] In a first aspect, the present application provides a DAS VSP wave field simulation method, comprising:
[0008] inputting initialization parameters of DAS VSP wave field simulation;
[0009] calculating displacement wave field data according to the initialization parameters;
[0010] According to the displacement wave field data, DAS VSP wave field data is calculated and output.
[0011] As a specific implementation manner of the embodiment of the present disclosure, the initialization parameters include: well trajectory, P-S wave velocity model, source wavelet, source excitation point position coordinates, receiver position coordinates, DAS demodulation device gauge and DAS VSP wave field type.
[0012] As a specific implementation manner of the embodiment of the present disclosure, the DAS VSP wave field type includes strain type and strain rate type.
[0013] As a specific implementation manner of the embodiment of the present disclosure, calculating displacement wave field data includes:
[0014] All discrete point position coordinates x along the well trajectory with fixed interval Δd are calculated i , and the displacement wave field data at each discrete point is calculated.
[0015] As a specific implementation manner of the embodiment of the present disclosure, the finite difference or finite element elastic wave simulation method is used to calculate the displacement wave field data.
[0016] As a specific implementation manner of the embodiment of the present disclosure, according to the displacement wave field data, calculating DAS VSP wave field data includes
[0017] According to the displacement wave field data, the wave field d(r i ,t) at all receivers is calculated.
[0018] If the DAS VSP wave field type is strain type, the DAS VSP wave field data is d(r i ,t); if the DAS VSP wave field type is strain rate type, the final DAS VSP wave field data is obtained by taking the time difference of the wave field d(r i ,t).
[0019] As a specific implementation manner of the embodiment of the present disclosure, the wave field at the receiver is:
[0020]
[0021] Wherein, L G is the DAS demodulation device gauge, and r i is any receiver position.
[0022] In a second aspect, the embodiment of the present disclosure also provides a DAS VSP wave field simulation device, which includes:
[0023] An input module inputs the initialization parameters of DAS VSP wave field simulation;
[0024] a calculation module configured to calculate displacement wavefield data according to the initialization parameters;
[0025] an output module configured to calculate and output DAS VSP wavefield data according to the displacement wavefield data.
[0026] As a specific implementation manner of the embodiment of the present disclosure, the initialization parameters include: well trajectory, P-S wave velocity model, source wavelet, source excitation point position coordinates, receiver position coordinates, DAS demodulation device gage length, and DAS VSP wavefield type.
[0027] As a specific implementation manner of the embodiment of the present disclosure, the DAS VSP wavefield type includes strain type and strain rate type.
[0028] As a specific implementation manner of the embodiment of the present disclosure, calculating displacement wavefield data includes:
[0029] calculating all discrete point position coordinates x along the well trajectory with a fixed interval of Δd i to obtain displacement wavefield data at each discrete point.
[0030] As a specific implementation manner of the embodiment of the present disclosure, the displacement wavefield data is calculated by a finite difference or finite element elastic wave simulation method.
[0031] As a specific implementation manner of the embodiment of the present disclosure, calculating DAS VSP wavefield data according to the displacement wavefield data includes
[0032] calculating wavefield d(r i ,t) at all receivers according to the displacement wavefield data.
[0033] If the DAS VSP wavefield type is strain type, the DAS VSP wavefield data is d(r i ,t); if the DAS VSP wavefield type is strain rate type, the DAS VSP wavefield data is obtained by taking time difference of wavefield d(r i ,t).
[0034] As a specific implementation manner of the embodiment of the present disclosure, the wavefield at the receiver is:
[0035]
[0036] wherein, L G is the DAS demodulation device gage length, and r i is any receiver position.
[0037] In a third aspect, the embodiment of the present disclosure further provides an electronic device, which comprises:
[0038] a memory storing executable instructions;
[0039] a processor running the executable instructions in the memory to implement the DAS VSP wave field simulation method.
[0040] In a fourth aspect, the embodiments of the present disclosure further provide a computer readable storage medium storing a computer program, which, when executed by a processor, implements the DAS VSP wave field simulation method.
[0041] The beneficial effects are that:
[0042] 1) The present application fully considers the influence of the DAS demodulation equipment gauge length;
[0043] 2) The present application is based on the conventional displacement wave field simulation, and is simple to implement.
[0044] The method and device of the present application have other characteristics and advantages, which will be apparent or will be described in detail in the accompanying drawings and subsequent detailed description incorporated herein, which together serve to explain the specific principles of the present application. BRIEF DESCRIPTION OF DRAWINGS
[0045] The above and other objects, features and advantages of the present application will become more apparent from the following detailed description when taken in conjunction with the accompanying drawings, in which like reference characters refer to like parts throughout the figures, and wherein:
[0046] Figure 1 A flow chart showing the steps of a DAS VSP wave field simulation method according to one embodiment of the present application is shown.
[0047] Figure 2a and Figure 2b Schematic diagrams showing the horizontal component and the vertical component of simulated VSP displacement wave field data according to one embodiment of the present application are shown, respectively.
[0048] Figure 3 A schematic diagram showing simulated DAS VSP strain wave field data according to one embodiment of the present application is shown.
[0049] Figure 4 A block diagram of a DAS VSP wave field simulation device according to one embodiment of the present application is shown.
[0050] BRIEF DESCRIPTION OF DRAWINGS
[0051] 201, input module; 202, calculation module; 203, output module. DETAILED DESCRIPTION
[0052] Preferred embodiments of the present application will be described in more detail below. Although the following describes preferred embodiments of the present application, it should be understood that the present application can be implemented in various forms and should not be limited by the embodiments set forth herein.
[0053] To facilitate understanding of the solutions and effects of the embodiments of the present application, six specific application examples are given below. Those skilled in the art should understand that the examples are only for facilitating understanding of the present application, and any specific details thereof are not intended to limit the present application in any way.
[0054] Example 1
[0055] Figure 1 A flowchart showing steps of a DAS VSP wavefield simulation method according to an embodiment of the present application is shown.
[0056] As shown in Figure 1 , the DAS VSP wavefield simulation method comprises:
[0057] Step 101, input initialization parameters of the DAS VSP wavefield simulation;
[0058] Step 102, calculate displacement wavefield data according to the initialization parameters;
[0059] Step 103, calculate DAS VSP wavefield data according to the displacement wavefield data and output.
[0060] In one example, the initialization parameters include: well trajectory, P-SV velocity model, source wavelet, source excitation point position coordinates, receiver position coordinates, DAS demodulation device gage length, and DAS VSP wavefield type.
[0061] In one example, the DAS VSP wavefield type includes strain type and strain rate type.
[0062] In one example, calculating displacement wavefield data comprises:
[0063] calculating displacement wavefield data at each discrete point. i
[0064] In one example, the finite difference or finite element elastic wave simulation method calculates displacement wavefield data.
[0065] In one example, calculating DAS VSP wavefield data according to displacement wavefield data comprises
[0066] calculating wavefield d(r i t);
[0067] If the DAS VSP wave field type is strain type, the DAS VSP wave field data is d(r i t) is differentiated with respect to time to obtain the final DAS VSP wave field data. i t) is differentiated with respect to time to obtain the final DAS VSP wave field data.
[0068] In one example, the wave field at the receiver point is:
[0069]
[0070] wherein L G is the DAS demodulation device pitch, and r i is any receiver point position.
[0071] Specifically, the initialization parameters of the DAS VSP wave field simulation are input, including: well trajectory, P and S wave velocity model, source wavelet, source excitation point position coordinates, receiver point position coordinates, DAS demodulation device pitch, and DAS VSP wave field type, etc., wherein the DAS wave field type includes two types of strain or strain rate.
[0072] Displacement wave field data is calculated; all discrete point position coordinates x i along the well trajectory with a fixed interval of Δd are calculated, wherein the fixed interval Δd can be artificially selected and is generally less than or equal to the receiver point coordinate interval; the displacement wave field U(x i , t) at each discrete point is calculated, wherein the displacement wave field U is a vector including a horizontal component and a depth component; x i is the discrete point position coordinate, and t is time. The displacement wave field can be calculated by using conventional elastic wave simulation methods such as finite difference and finite element.
[0073] DAS VSP wave field data is calculated: the simulated displacement wave field U(x i , t) is input; the wave field d(r i , t) at all receiver points is calculated by using the displacement wave field U(x i , t):
[0074]
[0075] wherein L G is the DAS demodulation device pitch, r i is any receiver point position, and the discrete point x j in the formula needs to satisfy the condition that the distance from the receiver point position r i is less than 1 / 2 pitch, i.e., |x j -r i | < LG / 2.
[0076] If the DAS wave field type is a strain type, the DAS VSP wave field data is d(r i ,t) directly; if it is a strain rate type, the wave field d(r i ,t) needs to be differentiated with respect to time to obtain the final DAS VSP wave field data.
[0077] Example 2
[0078] The application also provides a DAS VSP wave field simulation device, comprising:
[0079] An input module for inputting initialization parameters of DAS VSP wave field simulation;
[0080] A calculation module for calculating displacement wave field data according to the initialization parameters;
[0081] An output module for calculating and outputting DAS VSP wave field data according to the displacement wave field data.
[0082] In one example, the initialization parameters include: well trajectory, P-S wave velocity model, source wavelet, source excitation point position coordinates, receiver position coordinates, DAS demodulation device gauge length, and DAS VSP wave field type.
[0083] In one example, the DAS VSP wave field type includes a strain type and a strain rate type.
[0084] In one example, calculating displacement wave field data includes:
[0085] Calculating all discrete point position coordinates x i along the well trajectory with a fixed interval of Δd, and calculating displacement wave field data at each discrete point.
[0086] In one example, the displacement wave field data is calculated by a finite difference or finite element elastic wave simulation method.
[0087] In one example, calculating DAS VSP wave field data according to the displacement wave field data includes
[0088] Calculating wave field d(r i ,t) at all receivers according to the displacement wave field data;
[0089] If the DAS VSP wave field type is a strain type, the DAS VSP wave field data is d(r i ,t); if the DAS VSP wave field type is a strain rate type, the wave field d(r i ,t) is differentiated with respect to time to obtain the final DAS VSP wave field data.
[0090] In one example, the wave field at the receiver point is:
[0091]
[0092] Among them, L G For the gauge length of the DAS demodulation device, r i For any receiver location.
[0093] Specifically, the initialization parameters for the DAS VSP wavefield simulation are input, including: well trajectory, P-wave and S-wave velocity model, source wavelet, source excitation point location coordinates, receiver location coordinates, DAS demodulation equipment gauge length, and DAS VSP wavefield type, among which the DAS wavefield type includes two types: strain or strain rate.
[0094] Calculate displacement wavefield data; calculate the position coordinates x of all discrete points along the well trajectory at a fixed interval Δd. i The fixed spacing Δd can be selected manually, and is generally less than or equal to the coordinate spacing of the receiver points; the displacement wave field U(x) at each discrete point is calculated. i ,t), where the displacement wave field U is a vector, including horizontal and depth components; x i Let be the coordinates of the discrete point, and t be the time. The displacement wave field can be calculated using conventional elastic wave simulation methods such as finite difference and finite element methods.
[0095] The DAS VSP wavefield data was calculated as follows: Input simulated displacement wavefield U(x) i ,t); using displacement wave field U(x i The wave field d(r) at all receiver points was calculated using the method t). i ,t):
[0096]
[0097] Among them, L G For the gauge length of the DAS demodulation device, r i For any receiver location, the discrete point x in the formula j The distance r from the detector point must be satisfied. i Less than 1 / 2 gauge length, i.e., |x j -r i | <L G / 2.
[0098] If the DAS wavefield type is strain type, then the DAS VSP wavefield data is directly d(r i If it is a strain rate type, then a wave field d(r) is required. i The final DAS VSP wavefield data is obtained by differentiating the time (t) with respect to time.
[0099] Example 3
[0100] This embodiment uses the method to simulate DAS VSP wave field. In this embodiment, a simple layered stratum velocity model is used to simulate the full well section DAS VSP wave field of a horizontal well to verify the method.
[0101] Figure 2a and Figure 2b Figures respectively show schematic diagrams of horizontal and vertical components of simulated VSP displacement wave field data according to an embodiment of the present application.
[0102] Figure 3 Figure shows a schematic diagram of simulated DAS VSP strain wave field data according to an embodiment of the present application.
[0103] As Figure 2a , Figure 2b , Figure 3 can be seen, the wave field characteristics of DAS VSP, such as polarity reversal, are effectively simulated. Model testing shows the effectiveness of the simulation method of the present application.
[0104] Example 4
[0105] Figure 4 Figure shows a block diagram of a DAS VSP wave field simulation device according to an embodiment of the present application.
[0106] As Figure 4 shown, the DAS VSP wave field simulation device comprises:
[0107] an input module 201 for inputting initialization parameters of DAS VSP wave field simulation;
[0108] a calculation module 202 for calculating displacement wave field data according to the initialization parameters;
[0109] an output module 203 for calculating DAS VSP wave field data according to the displacement wave field data and outputting.
[0110] In one example, the initialization parameters comprise: well trajectory, P-S wave velocity model, source wavelet, source excitation point position coordinates, receiver position coordinates, DAS demodulation device gauge length, and DAS VSP wave field type.
[0111] In one example, the DAS VSP wave field type comprises strain type and strain rate type.
[0112] In one example, calculating displacement wave field data comprises:
[0113] Compute all discrete point position coordinates x along the well trajectory with fixed interval of Δd i , and the displacement wavefield data at each discrete point is calculated.
[0114] In one example, the finite-difference or finite-element elastic wave simulation method calculates the displacement wavefield data.
[0115] In one example, the computing DAS VSP wavefield data according to the displacement wavefield data includes
[0116] According to the displacement wavefield data, the wavefield d(r i ,t) at all geophones is calculated.
[0117] If the DAS VSP wavefield type is strain type, the DAS VSP wavefield data is d(r i ,t); if the DAS VSP wavefield type is strain rate type, the DAS VSP wavefield data is obtained by taking the time difference of the wavefield d(r i ,t).
[0118] In one example, the wavefield at the geophone is:
[0119]
[0120] where L G is the DAS demodulation device gauge length, and r i is any geophone position.
[0121] Example 5
[0122] The embodiment provides an electronic device, which includes a memory storing executable instructions, and a processor running the executable instructions in the memory to implement the DAS VSP wavefield simulation method.
[0123] The electronic device according to the embodiment of the present disclosure includes a memory and a processor.
[0124] The memory is configured to store non-transitory computer-readable instructions. Specifically, the memory can include one or more computer program products, which can include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM), cache memory, and / or the like. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, and / or the like.
[0125] The processor can be a central processing unit (CPU) or other form of processing unit that has data processing and / or instruction execution capabilities, and can control other components in the electronic device to perform desired functions. In one embodiment of the present disclosure, the processor is configured to execute the computer-readable instructions stored in the memory.
[0126] Those skilled in the art will understand that, in order to solve the technical problem of how to obtain a good user experience effect, the embodiments can also include well-known structures such as a communication bus, an interface, and the like, which should also be included in the protection scope of the present disclosure.
[0127] Detailed descriptions of the embodiments are described above with reference to the corresponding descriptions of the previous embodiments, and will not be repeated here.
[0128] Example 6
[0129] The embodiments provide a computer readable storage medium, which stores a computer program. The computer program is executed by a processor to implement the DAS VSP wave field simulation method.
[0130] The computer readable storage medium according to the embodiments of the present disclosure has non-transitory computer readable instructions stored thereon. When the non-transitory computer readable instructions are run by a processor, all or part of the steps of the method of the embodiments of the present disclosure are executed.
[0131] The computer readable storage medium described above includes, but is not limited to, an optical storage medium (for example, CD-ROM and DVD), a magneto-optical storage medium (for example, MO), a magnetic storage medium (for example, magnetic tape or a mobile hard disk), a medium with a built-in rewritable non-volatile memory (for example, a memory card), and a medium with a built-in ROM (for example, a ROM cartridge).
[0132] Those skilled in the art will understand that the above description of the embodiments of the present disclosure is only for the purpose of exemplarily illustrating the beneficial effects of the embodiments of the present disclosure, and is not intended to limit the embodiments of the present disclosure to any examples given.
[0133] The embodiments of the present disclosure have been described above, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes are obvious to those skilled in the art without departing from the scope and spirit of the described embodiments.
Claims
1. A DAS VSP wavefield simulation method, characterized in that, include: Input the initialization parameters for the DAS VSP wavefield simulation; Calculate the displacement wave field data based on the initialization parameters; Based on the displacement wave field data, calculate and output the DAS VSP wave field data.
2. The DAS VSP wavefield simulation method according to claim 1, wherein, The initialization parameters include: well trajectory, P-wave and S-wave velocity model, source wavelet, source excitation point location coordinates, detector point location coordinates, DAS demodulation equipment gauge length, and DAS VSP wave field type.
3. The DAS VSP wavefield simulation method according to claim 2, wherein, DAS VSP wave field types include strain type and strain rate type.
4. The DAS VSP wavefield simulation method according to claim 1, wherein, The displacement wavefield data to be calculated includes: Calculate the coordinates x of all discrete points along the well trajectory at a fixed interval Δd. i The displacement wave field data at each discrete point were calculated.
5. The DAS VSP wavefield simulation method according to claim 4, wherein, The displacement wave field data are calculated using finite difference or finite element elastic wave simulation methods.
6. The DAS VSP wavefield simulation method according to claim 3, wherein, Based on the displacement wavefield data, the calculation of DAS VSP wavefield data includes... Based on the displacement wavefield data, calculate the wavefield d(r) at all receiver points. i ,t); If the DAS VSP wavefield type is strain type, then the DAS VSP wavefield data is d(r i If the DAS VSP wave field type is strain rate type, then for the wave field d(r) i The final DAS VSP wavefield data is obtained by differentiating the time (t) with respect to time.
7. The DAS VSP wavefield simulation method according to claim 6, wherein, The wave field at the receiver point is: Among them, L G For the gauge length of the DAS demodulation device, r i For any receiver location.
8. A DAS VSP wavefield simulation device, characterized in that, include: Input module: Input the initialization parameters for the DAS VSP wavefield simulation; The calculation module calculates the displacement wave field data based on the initialization parameters; The output module calculates and outputs DAS VSP wavefield data based on the displacement wavefield data.
9. An electronic device, characterized in that, The electronic device includes: Memory, which stores executable instructions; A processor that executes the executable instructions in the memory to implement the DAS VSP wavefield simulation method according to any one of claims 1-7.
10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, implements the DAS VSP wavefield simulation method according to any one of claims 1-7.