Wave field separation method, device and equipment of DAS data, storage medium and computer program
By using channel grouping and shaping median filtering, the difficulties in wavefield separation caused by the small difference between P-wave and S-wave velocities and high-density acquisition in DAS-VSP data were solved, achieving high-precision wavefield separation.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- CHINA PETROLEUM & CHEMICAL CORP
- Filing Date
- 2024-12-27
- Publication Date
- 2026-06-30
AI Technical Summary
In the processing of vertical seismic profile (VSP) data in wells, conventional wavefield separation methods such as median filtering and mean filtering are ineffective when the difference between P-wave and S-wave velocities is small. Furthermore, high-density acquisition results in small inter-channel time differences, making it difficult to effectively separate the overlapping P-wave and S-wave fields.
A method combining channel grouping and shaping median filtering is adopted. By setting the channel interval, channel groups are formed and median filtering and shaping are performed on each channel of the new channel group to restore the original channel group and achieve wavefield separation.
High-precision DAS-VSP data wavefield separation was achieved, clearly separating the uplink and downlink longitudinal and transverse wavefields, solving the separation difficulties of conventional methods under conditions of small velocity difference and high-density acquisition.
Smart Images
Figure CN122307649A_ABST
Abstract
Description
Technical Field
[0001] This disclosure pertains to the application field of fiber optic sensing technology, and specifically relates to a wavefield separation method, conversion device, equipment, storage medium, and computer program for distributed fiber optic sensing (DAS) data. Background Technology
[0002] Distributed fiber acoustic sensing (DAS) technology, as an emerging and transformative technology, has developed rapidly in fields such as oil and gas exploration and development and engineering monitoring, but there are still some problems in DAS data processing.
[0003] In the acquisition of Vertical Seismic Profile (VSP) data, DAS technology uses single-component data. Within the same component, uplink and downlink P-waves, S-waves, and other interfering waves are superimposed. Conventional vector wavefield separation methods based on three-component reception cannot be used; only scalar wavefield separation methods are applicable. However, these scalar wavefield separation methods also have limitations in DAS data. For example, median filtering is ineffective when the P-wave and S-wave velocity differences are small; mean filtering, on the other hand, exhibits severe boundary smoothing effects. Furthermore, due to the high-density acquisition of DAS-VSP data, the gather spacing is typically on the order of meters. Such small gather spacing results in very small inter-trace time differences, hindering the use of conventional wavefield separation methods such as median filtering. Summary of the Invention
[0004] Based on the above, and addressing the issue that conventional vector wavefield separation methods based on three-component reception and scalar wavefield separation methods such as median filtering are unsuitable for vertical seismic profile (VSP) data acquired by distributed optical fiber acoustic sensing (DAS) technology, which involves single-component reception, high sampling density, and complex wavefield, the purpose of this disclosure is to propose a wavefield separation method and device for distributed optical fiber sensing data. This wavefield separation method uses a channel-level combination of shaping median filtering to perform wavefield separation of DAS-VSP data, thereby achieving high-precision wavefield separation of DAS-VSP data.
[0005] In a first aspect, this disclosure provides a wavefield separation method for distributed optical fiber sensing data, including:
[0006] S1. Obtain DAS-VSP data, wherein the DAS-VSP data contains aliased wavefields that need to be separated.
[0007] S2. Divide the sampling channels of the DAS-VSP data into several new channel sets;
[0008] S3. The new track set is shaped by applying median filtering to each track to obtain the shaped new track set;
[0009] S4. Perform gather recovery on the shaped new gather to obtain DAS-VSP data after wavefield separation.
[0010] S5. Output the DAS-VSP data after wavefield separation.
[0011] In this disclosure, the DAS-VSP data refers to well vertical seismic profile data acquired by distributed fiber optic sensing.
[0012] In some embodiments, the aliased wave field includes an uplink aliased wave field and a crosslink aliased wave field.
[0013] In some embodiments, step S2, which involves grouping the sampling channels of the DAS-VSP data into several new channel sets, includes:
[0014] S21. Set a channel interval d, where the original interval of the sampling channel < d ≤ 20;
[0015] S22. For the sampling channels of the DAS-VSP data, starting from the first channel with the original number, another channel is extracted every d channels until the channel interval d is no longer met, at which point the extraction stops, and the extracted sampling channels are combined into a new channel set according to the extraction order.
[0016] S23. For the remaining sampling channels after step S22, start sampling from the first channel, and sample one channel every d channels until the channel interval d is no longer met, then stop sampling and form another new channel set by sampling the above-mentioned sampling channels in the order of sampling.
[0017] S24. Repeat step S23 until all sampling channels have been extracted to form a new gather, and there are no remaining sampling channels.
[0018] In this disclosure, the channel grouping in step S2 does not change the original numbering of each sampling channel.
[0019] In some embodiments, the size of the channel spacing d is determined based on the final wavefield separation effect test.
[0020] In some embodiments, step S3, which involves shaping the new gather channel by channel using median filtering to obtain the shaped new gather, includes:
[0021] S31. Pick up the original travel time of each sampling channel in the new channel set and compare them to obtain a minimum travel time;
[0022] S32. Based on the minimum travel time, the travel times of all sampling channels are aligned to obtain a new set of channels after alignment;
[0023] S33. For the new channel set after alignment, calculate the midpoint of the reshape at all sampling times for all sampling channels;
[0024] S34. Restore the travel time of the sampling channel to the original travel time to obtain the new set after anti-alignment as the new set after shaping.
[0025] In some embodiments, in step S32, the alignment means converting the travel times of all acquisition channels from the original travel times to the minimum travel times.
[0026] In some embodiments, step S33, calculating the shaped midpoint over all sampling times for all sampling channels in the aligned new gather, includes:
[0027] For each sampling channel in the newly aligned channel set, perform the following steps:
[0028] S331. Set multiple sampling times. At any sampling time t, take the current sampling channel c. With the sample point (t,c) as the midpoint, take a time window w. Calculate the median value med and the average value mean of the sample points within the time window w. Calculate the new value of the midpoint as the integer midpoint (t',c') according to the following formula:
[0029] The midpoint of the shape (t', c') = med + (mean - med) × α(I)
[0030] Where α is the shaping factor, 0 < α < 1.
[0031] In some embodiments, the time window w is an odd number greater than or equal to 10 and less than or equal to 100. The size of the shaping factor α and the time window w is determined by the effect of the wavefield-separated DAS-VSP data finally output in step S5.
[0032] In some embodiments, step S4, which involves recovering the shaped new gather to obtain wavefield-separated DAS-VSP data, includes:
[0033] S41. Restore all sampling channels in the reshaped new gather according to their original numbers to obtain the target gather and the DAS-VSP data after wavefield separation.
[0034] Secondly, this disclosure provides a wavefield separation device for distributed optical fiber data, comprising:
[0035] The acquisition module is used to acquire DAS-VSP data, which contains aliased wavefields that need to be separated.
[0036] The channel set grouping module is used to divide the sampling channels of the DAS-VSP data into several new channel sets;
[0037] The shaping module is used to shape the new track set track by track using median filtering to obtain the shaped new track set.
[0038] The gather group recovery module is used to recover the shaped new gathers to obtain DAS-VSP data after wavefield separation.
[0039] The output module is used to output the DAS-VSP data after wavefield separation.
[0040] In some embodiments, the aliased wave field includes an uplink aliased wave field and a crosslink aliased wave field.
[0041] In some embodiments, the lane grouping module includes:
[0042] The setting unit is used to set a channel interval d, where the original interval of the sampling channel is <d≤20;
[0043] The first extraction unit extracts samples from the first channel of the original number for the sampling channels of the DAS-VSP data. It then extracts one channel every d channels until the channel interval d is no longer met, at which point the extraction stops. The extracted sampling channels are then combined into a new channel set according to the extraction order.
[0044] The second extraction unit extracts from the remaining sampling channels after the first extraction unit extracts, starting from the first channel. It then extracts one channel every d channels until the channel interval d is no longer met, at which point it stops extracting and combines the extracted sampling channels into another new channel set in the order of extraction.
[0045] Repeat the extraction unit, repeating the operation of the second extraction unit, until all sampling channels have been extracted to form a new gather, and there are no remaining sampling channels.
[0046] In some embodiments, the size of the channel spacing d is determined based on the final wavefield separation effect test.
[0047] In some embodiments, the shaping module includes:
[0048] The second picking unit is used to pick up the original travel time of each sampling channel in the new channel set and compare it to obtain a minimum travel time;
[0049] The alignment unit is used to align the travel times of all sampling channels based on the minimum travel time to obtain a new channel set after alignment.
[0050] The shaping unit is used to calculate the shaping midpoint of all sampling channels at all sampling times for the new, aligned channel set.
[0051] The de-aligning unit is used to restore the travel time of the sampling channel to the original travel time, and obtain a new track set after de-aligning as the new track set after shaping.
[0052] In some embodiments, the shaping unit is configured to perform the following steps for each of the sampling channels in the aligned new channel set:
[0053] S331. Set multiple sampling times. At any sampling time t, take the current sampling channel c. With the sample point (t,c) as the midpoint, take a time window w. Calculate the median value med and the average value mean of the sample points within the time window w. Calculate the new value of the midpoint as the integer midpoint (t',c') according to the following formula:
[0054] The midpoint of the shape (t', c') = med + (mean - med) × α(I)
[0055] Where α is the shaping factor, 0 < α < 1.
[0056] In some embodiments, the time window w is an odd number greater than or equal to 10 and less than or equal to 100. The size of the shaping factor α and the time window w is determined by the effect of the wavefield-separated DAS-VSP data finally output by the output module.
[0057] In some embodiments, the gather group recovery module is used to perform the following steps:
[0058] S41. Restore all sampling channels in the reshaped new gather to the original gather according to their original numbers.
[0059] Thirdly, this disclosure provides a computer device including a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the above-described method.
[0060] Fourthly, this disclosure provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps of the above-described method.
[0061] Fifthly, this disclosure provides a computer program product, including a computer program that, when executed by a processor, implements the steps of the above-described method.
[0062] The beneficial effects of this disclosure are as follows:
[0063] To address the issue that conventional vector wavefield separation methods based on three-component reception and scalar wavefield separation methods such as median filtering are unsuitable for vertical seismic profile (VSP) data acquired by distributed acoustic sensing (DAS) technology, which involves single-component reception, high sampling density, and complex wavefields, this disclosure presents a wavefield separation method based on channel grouping combined with shaping median filtering to process DAS-VSP data, achieving effective separation of the wavefield in DAS-VSP data. Attached Figure Description
[0064] The present disclosure will be described in more detail below based on embodiments and with reference to the accompanying drawings:
[0065] Figure 1 This is a schematic diagram of an application scenario of the wave field separation method for distributed optical fiber sensing data provided in the embodiments of this disclosure.
[0066] Figure 2 This is a flowchart of a second embodiment of the wavefield separation method for distributed optical fiber sensing data provided in accordance with the present disclosure.
[0067] Figure 3 This is the original DAS-VSP data before wavefield separation.
[0068] Figure 4a To Figure 3 The original DAS-VSP data uses the downlink longitudinal wave after wavefield separation as disclosed in this publication.
[0069] Figure 4b To Figure 3 The original DAS-VSP data uses the uplink P-wave after wavefield separation as disclosed in this publication.
[0070] Figure 4c To Figure 3 The original DAS-VSP data uses the downlink shear wave after wavefield separation as disclosed in this publication.
[0071] Figure 4d To Figure 3 The original DAS-VSP data uses the up-going shear wave after wavefield separation as disclosed in this publication.
[0072] from Figures 4a to 4d As can be seen, each wave field is completely separated, and the wave field phase axis is clear, indicating that this disclosure has achieved effective separation of the wave field of DAS-VSP data.
[0073] Figure 5 This is a simplified structural diagram of a wave field separation device for distributed optical fiber sensing data provided according to an embodiment of the present disclosure.
[0074] Figure 6 This is a schematic diagram of an electronic device provided according to an embodiment of the present disclosure.
[0075] In the accompanying drawings, the same parts are referred to by the same reference numerals, and the drawings are not drawn to scale. Detailed Implementation
[0076] To enable those skilled in the art to better understand the technical solutions of this disclosure, and to fully understand and implement the process of how this disclosure applies technical means to solve technical problems and achieve corresponding technical effects, the technical solutions in the embodiments of this disclosure will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this disclosure, not all embodiments. The embodiments of this disclosure and the various features within them can be combined with each other without conflict, and the resulting technical solutions are all within the protection scope of this disclosure. All other embodiments obtained by those skilled in the art based on the embodiments of this disclosure without creative effort should fall within the protection scope of this disclosure.
[0077] It should be noted that the terms "first," "second," etc., in the specification, claims, and accompanying drawings of this disclosure are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that the embodiments of this disclosure described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.
[0078] It should be noted that the steps shown in the flowchart in the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical order is shown in the flowchart, in some cases the steps shown or described may be executed in a different order than that shown here.
[0079] Example 1
[0080] Figure 1This is a schematic diagram illustrating an application scenario of an embodiment of this disclosure. The application scenario may include a server 4, a network 5, and a terminal device, which can be hardware or software. When the terminal device is hardware, it can be various electronic devices with an LED display screen that support communication with the server 1, including but not limited to smartphones 1, tablets 3, laptops 2, and desktop computers; when the terminal device is software, it can be installed in the aforementioned electronic devices. The terminal device can be implemented as multiple software programs or software modules, or as a single software program or software module; this disclosure does not limit this. Furthermore, various applications can be installed on the terminal device, such as data processing applications, instant messaging tools, social platform software, search applications, shopping applications, etc.
[0081] Server 4 can be a server that provides various services, such as a backend server that receives requests sent by terminal devices with which it has established communication connections. This backend server can receive and analyze the requests sent by the terminal devices and generate processing results. Server 4 can be a single server, a server cluster consisting of several servers, or a cloud computing service center. This disclosure embodiment does not limit this.
[0082] It should be noted that server 4 can be either hardware or software. When server 1 is hardware, it can be any electronic device that provides various services to the terminal device. When server 1 is software, it can implement multiple software programs or software modules that provide various services to the terminal device, or it can implement a single software program or software module that provides various services to the terminal device. This disclosure does not impose any limitations on this aspect.
[0083] Network 5 can be a wired network using coaxial cable, twisted pair, and fiber optic connection, or it can be a wireless network that enables interconnection of various communication devices without wiring, such as Bluetooth, Near Field Communication (NFC), and Infrared. This disclosure does not limit the scope of the embodiments.
[0084] Users can establish a communication connection with server 4 via network 5 through terminal devices to receive or send information, etc.
[0085] Specifically, firstly, server 4 can acquire the raw data from the distributed fiber optic sensing. Secondly, server 4 can perform a first transformation process, converting the raw data from a phase response into a strain response to obtain first data. Thirdly, server 4 can perform a second transformation process, converting the first data from a strain response into a displacement response to obtain second data. Fourthly, server 4 can perform a third transformation process, converting the second data from a displacement response into a velocity response to obtain third data. Finally, server 4 can output the third data.
[0086] It should be noted that the specific types, quantities, and combinations of server 4, network 5, and terminal devices can be adjusted according to the actual needs of the application scenario, and this disclosure embodiment does not impose any restrictions on this.
[0087] Example 2
[0088] Continue to refer to Figure 2 , Figure 2 A schematic flowchart of a wavefield separation method for distributed optical fiber sensing data provided in an embodiment of this disclosure is shown. This method can be... Figure 1 It is executed by electronic devices within the system. For example... Figure 2 As shown, the wavefield separation method for distributed optical fiber sensing data includes:
[0089] S1. Obtain DAS-VSP data, wherein the DAS-VSP data contains aliased wavefields that need to be separated.
[0090] S2. Divide the sampling channels of the DAS-VSP data into several new channel sets;
[0091] S3. The new track set is shaped by applying median filtering to each track to obtain the shaped new track set;
[0092] S4. Perform gather recovery on the shaped new gather to obtain DAS-VSP data after wavefield separation.
[0093] S5, Output DAS-VSP data after wavefield separation.
[0094] In some embodiments, the entity performing the conversion method (such as...) Figure 1 The electronic device shown can connect to the target device via wired or wireless connection, and then acquire DAS-VSP data. DAS data refers to vertical seismic profile data acquired from wells using distributed fiber optic sensing. DAS-VSP data contains aliased wavefields that require wavefield separation, for example... Figure 3 The data in the middle.
[0095] In some embodiments, the aliased wave field includes an uplink aliased wave field and a crosslink aliased wave field.
[0096] It should be noted that the aforementioned wireless connection methods may include, but are not limited to, 3G / 4G / 5G connections, WiFi connections, Bluetooth connections, WiMAX connections, Zigbee connections, UWB (ultra wideband) connections, and other currently known or future wireless connection methods.
[0097] In this embodiment, the channel grouping in step S2 does not change the original number of each sampling channel.
[0098] In some embodiments, step S2, which involves grouping the sampling channels of the DAS-VSP data into several new channel sets, includes:
[0099] S21. Set a channel interval d, where the original interval of the sampling channel < d ≤ 20;
[0100] S22. For the sampling channels of the DAS-VSP data, starting from the first channel with the original number, another channel is extracted every d channels until the channel interval d is no longer met, at which point the extraction stops, and the extracted sampling channels are combined into a new channel set according to the extraction order.
[0101] S23. For the remaining sampling channels after step S22, start sampling from the first channel, and sample one channel every d channels until the channel interval d is no longer met, then stop sampling and form another new channel set by sampling the above-mentioned sampling channels in the order of sampling.
[0102] S24. Repeat step S23 until all sampling channels have been extracted to form a new gather, and there are no remaining sampling channels.
[0103] In some embodiments, the size of the channel spacing d is determined based on the final wavefield separation effect test.
[0104] In some embodiments, step S3, which involves shaping the new gather channel by channel using median filtering to obtain the shaped new gather, includes:
[0105] S31. Pick up the original travel time of each sampling channel in the new channel set and compare them to obtain a minimum travel time;
[0106] S32. Based on the minimum travel time, the travel times of all sampling channels are aligned to obtain a new set of channels after alignment;
[0107] S33. For the new channel set after alignment, calculate the midpoint of the reshape at all sampling times for all sampling channels;
[0108] S34. Restore the travel time of the sampling channel to the original travel time to obtain the new set after anti-alignment as the new set after shaping.
[0109] In this embodiment, in step S32, the alignment means converting the travel time of all acquisition channels from the original travel time to the minimum travel time.
[0110] In some embodiments, step S33, calculating the shaped midpoint over all sampling times for all sampling channels in the aligned new gather, includes:
[0111] For each sampling channel in the newly aligned channel set, perform the following steps:
[0112] S331. Set multiple sampling times. At any sampling time t, take the current sampling channel c. With the sample point (t,c) as the midpoint, take a time window w. Calculate the median value med and the average value mean of the sample points within the time window w. Calculate the new value of the midpoint as the integer midpoint (t',c') according to the following formula:
[0113] The midpoint of the shape (t', c') = med + (mean - med) × α(I)
[0114] Where α is the shaping factor, 0 < α < 1.
[0115] In some embodiments, the time window w is an odd number greater than or equal to 10 and less than or equal to 100. The size of the shaping factor α and the time window w is determined by the effect of the wavefield-separated DAS-VSP data finally output in step S5.
[0116] In some embodiments, step S4, which involves recovering the shaped new gather to obtain wavefield-separated DAS-VSP data, includes:
[0117] S41. Restore all sampling channels in the reshaped new gather according to their original numbers to obtain the target gather and the DAS-VSP data after wavefield separation.
[0118] The beneficial effects of one of the above embodiments of this disclosure include at least the following: effective separation of the wavefield of DAS-VSP data is achieved by processing DAS-VSP data based on channel grouping combined with shaping median filtering.
[0119] The present disclosure is illustrated below through a specific embodiment:
[0120] Figure 3 Raw DAS-VSP data was collected for a certain work area before wavefield separation. Figure 3The original DAS-VSP data is obtained after processing according to steps S2 to S5 in the above embodiment. Figures 4a to 4d The separated wavefield data shown are as follows:
[0121] Figure 4a To Figure 3 The original DAS-VSP data uses the downlink P-wave after wavefield separation as disclosed in this publication;
[0122] Figure 4b To Figure 3 The original DAS-VSP data uses the uplink P-wave after wavefield separation as disclosed in this publication;
[0123] Figure 4c To Figure 3 The original DAS-VSP data uses the downlink shear wave after wavefield separation as disclosed in this publication;
[0124] Figure 4d To Figure 3 The original DAS-VSP data uses the up-going shear wave after wavefield separation as disclosed in this publication.
[0125] from Figures 4a to 4d It can be seen that each wave field is completely separated, and the phase axis of the wave field is clear, indicating that the present invention has achieved effective separation of the wave field of DAS-VSP data.
[0126] Example 3
[0127] This embodiment describes some further implementation steps of the wavefield separation method for distributed optical fiber sensor data according to this disclosure. The wavefield separation method for distributed optical fiber sensor data includes:
[0128] S100. Obtain DAS-VSP data. The DAS-VSP data contains aliased wave fields that need to be separated, such as uplink and downlink aliased wave fields and crosslink and crosslink aliased wave fields.
[0129] S200. The sampling channels of the DAS-VSP data are grouped into several new channel sets, including:
[0130] S210. Set a channel interval d, where the original interval of the sampling channel < d ≤ 20;
[0131] S220. For the sampling channels of the DAS-VSP data, starting from the first channel with the original number, another channel is extracted every d channels until the channel interval d is no longer met, at which point the extraction stops, and the extracted sampling channels are combined into a new channel set according to the extraction order.
[0132] S230. For the remaining sampling channels after extraction in step S22, start extraction from the first channel, and extract one channel every d channels until the channel interval d is no longer met, then stop extraction and form another new channel set by extracting each of the above sampling channels in the order of extraction.
[0133] S240. Repeat step S23 until all sampling channels have been extracted to form a new gather, and there are no remaining sampling channels.
[0134] S300. The new track set obtained in step S200 is shaped by median filtering for each track to obtain a shaped new track set, including:
[0135] S310. Pick up the original travel time of each sampling channel in the new channel set and compare them to obtain a minimum travel time;
[0136] S320. Based on the minimum travel time, the travel times of all sampling channels are aligned to obtain a new set of channels after alignment.
[0137] S330. For the aligned new gather, calculate the shaped midpoint of all sampling channels at all sampling times, including:
[0138] For each sampling channel in the newly aligned channel set, perform the following steps:
[0139] S331. Set multiple sampling times. At any sampling time t, take the current sampling channel c. With the sample point (t,c) as the midpoint, take a time window w. Calculate the median value med and the average value mean of the sample points within the time window w. Calculate the new value of the midpoint as the integer midpoint (t',c') according to the following formula:
[0140] The midpoint of the shape (t', c') = med + (mean - med) × α(I)
[0141] Where α is the shaping factor, 0 < α < 1;
[0142] S340. Restore the travel time of the sampling channel to the original travel time to obtain the new channel set after de-alignment as the new channel set after shaping;
[0143] S400. Perform gather recovery on the shaped new gather to obtain DAS-VSP data after wavefield separation, including:
[0144] S410. Restore all sampling channels in the reshaped new channel set to the original channel set according to their original numbers;
[0145] S500 outputs DAS-VSP data after wavefield separation.
[0146] In some embodiments, the specific implementation of steps S100 to S500 and the resulting technical effects can be referred to Figure 2 The steps in the corresponding embodiments will not be repeated here.
[0147] All of the above-mentioned optional technical solutions can be combined in any way to form the optional embodiments of this application, and will not be described in detail here.
[0148] Example 4
[0149] This embodiment describes some structural embodiments of the frequency response correction device for fractional fiber optic sensing data according to this disclosure. It can be used to perform the method embodiments of this disclosure. For details not disclosed in the device embodiments of this disclosure, please refer to the method embodiments of this disclosure.
[0150] refer to Figure 5 In some embodiments, the wavefield separation device 6 for distributed optical fiber data includes:
[0151] Acquisition module 7 is used to acquire DAS-VSP data, which contains aliased wavefields that need to be separated.
[0152] Channel set grouping module 8 is used to divide the sampling channels of the DAS-VSP data into several new channel sets;
[0153] Shaping module 9 is used to shape the new track set by applying median filtering to each track to obtain the shaped new track set;
[0154] The gather group recovery module 10 is used to recover the shaped new gathers to obtain DAS-VSP data after wavefield separation.
[0155] Output module 11 is used to output DAS-VSP data after wavefield separation.
[0156] In some embodiments, the aliased wave field includes an uplink aliased wave field and a crosslink aliased wave field.
[0157] In some embodiments, the lane grouping module includes:
[0158] The setting unit is used to set a channel interval d, where the original interval of the sampling channel is <d≤20;
[0159] The first extraction unit extracts samples from the first channel of the original number for the sampling channels of the DAS-VSP data. It then extracts one channel every d channels until the channel interval d is no longer met, at which point the extraction stops. The extracted sampling channels are then combined into a new channel set according to the extraction order.
[0160] The second extraction unit extracts from the remaining sampling channels after the first extraction unit extracts, starting from the first channel. It then extracts one channel every d channels until the channel interval d is no longer met, at which point it stops extracting and combines the extracted sampling channels into another new channel set in the order of extraction.
[0161] Repeat the extraction unit, repeating the operation of the second extraction unit, until all sampling channels have been extracted to form a new gather, and there are no remaining sampling channels.
[0162] In this embodiment, the channel grouping does not change the original numbering of each sampling channel.
[0163] In some embodiments, the size of the channel spacing d is determined based on the final wavefield separation effect test.
[0164] In some embodiments, the shaping module includes:
[0165] The second picking unit is used to pick up the original travel time of each sampling channel in the new channel set and compare it to obtain a minimum travel time;
[0166] The alignment unit is used to align the travel times of all sampling channels based on the minimum travel time to obtain a new channel set after alignment.
[0167] The shaping unit is used to calculate the shaping midpoint of all sampling channels at all sampling times for the new, aligned channel set.
[0168] The de-aligning unit is used to restore the travel time of the sampling channel to the original travel time, and obtain a new track set after de-aligning as the new track set after shaping.
[0169] In some embodiments, the shaping unit is configured to perform the following steps for each of the sampling channels in the aligned new channel set:
[0170] S331. Set multiple sampling times. At any sampling time t, take the current sampling channel c. With the sample point (t,c) as the midpoint, take a time window w. Calculate the median value med and the average value mean of the sample points within the time window w. Calculate the new value of the midpoint as the integer midpoint (t',c') according to the following formula:
[0171] The midpoint of the shape (t', c') = med + (mean - med) × α(I)
[0172] Where α is the shaping factor, 0 < α < 1.
[0173] In some embodiments, in step S331, the time window w is an odd number greater than or equal to 10 and less than or equal to 100. The size of the shaping factor α and the time window w is determined by the effect of the wavefield-separated DAS-VSP data finally output by the output module.
[0174] In some embodiments, the gather group recovery module is specifically used to perform the following steps:
[0175] S41. Restore all sampling channels in the reshaped new gather to the original gather according to their original numbers.
[0176] It is understandable that the modules described in the wavefield separation device 6 for distributed optical fiber data are similar to those in the reference device. Figure 2 The steps in the described method correspond to each other. Therefore, the operations, features, and beneficial effects described above for the method also apply to the wavefield separation device 6 and the modules contained therein, and will not be repeated here.
[0177] Example 5
[0178] like Figure 6 As shown, the electronic device 12 may include a processing unit (e.g., a central processing unit, a graphics processor, etc.) 12-1, which can perform various appropriate actions and processes according to a program stored in a read-only memory (ROM) 12-2 or a program loaded from a storage device 12-8 into a random access memory (RAM) 12-3. The RAM 12-3 also stores various programs and data required for the operation of the electronic device 12. The processing unit 12-1, ROM 12-2, and RAM 12-3 are interconnected via a bus 12-4. An input / output (I / O) interface 12-5 is also connected to the bus 12-4.
[0179] Typically, the following devices can be connected to I / O interface 12-5: input devices 12-6 including, for example, touchscreens, touchpads, keyboards, mice, cameras, microphones, accelerometers, gyroscopes, etc.; output devices 12-7 including, for example, liquid crystal displays (LCDs), speakers, vibrators, etc.; storage devices 12-8 including, for example, magnetic tapes, hard disks, etc.; and communication devices 12-9. Communication device 12-9 allows electronic device 12 to communicate wirelessly or wiredly with other devices to exchange data. Although Figure 6 An electronic device 12 with various devices is shown, but it should be understood that it is not required to implement or have all of the devices shown. More or fewer devices may be implemented or have instead. Figure 6 Each box shown can represent a device or multiple devices as needed.
[0180] Specifically, according to this embodiment, the process described in the above-mentioned flowchart can be implemented as a computer software program. For example, this embodiment includes a computer program product comprising a computer program or instructions carried on a computer-readable medium, the computer program or instructions containing program code for performing the methods shown in the flowchart. In this embodiment, the computer program can be downloaded and installed from a network via communication device 12-9, or installed from storage device 12-8, or installed from ROM 12-2. When the computer program is executed by processing device 12-1, the wavefield separation method for distributed optical fiber sensing data described above can be performed.
[0181] It should be noted that the computer-readable medium described above in this embodiment can be a computer-readable signal medium, a computer-readable storage medium, or any combination of the two. A computer-readable storage medium can be, for example,—but not limited to—an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination thereof. More specific examples of a computer-readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fiber, and a portable compact disk read-only memory. Optical storage devices, magnetic storage devices, or any suitable combination thereof. In this embodiment, a computer-readable storage medium can be any tangible medium containing or storing a program that can be used by or in connection with an instruction execution system, apparatus, or device. In this embodiment, a computer-readable signal medium can include a data signal propagated in baseband or as part of a carrier wave, carrying computer-readable program code. Such propagated data signals can take various forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination thereof. A computer-readable signal medium can also be any computer-readable medium other than a computer-readable storage medium, which can send, propagate, or transmit a program for use by or in connection with an instruction execution system, apparatus, or device. The program code contained on the computer-readable medium can be transmitted using any suitable medium, including but not limited to: wires, optical cables, RF (radio frequency), etc., or any suitable combination thereof.
[0182] In some implementations, clients and servers can communicate using any currently known or future-developed network protocol such as HTTP (Hypertext Transfer Protocol) and can interconnect with digital data communication (e.g., communication networks) of any form or medium. Examples of communication networks include local area networks (“LANs”), wide area networks (“WANs”), the Internet (e.g., the Internet of Things), and peer-to-peer networks (e.g., ad hoc peer-to-peer networks), as well as any currently known or future-developed networks.
[0183] The aforementioned computer-readable medium may be included in the aforementioned device; or it may exist independently and not assembled into the electronic device. The aforementioned computer-readable medium carries one or more programs that, when executed by the electronic device, cause the electronic device to perform the aforementioned wavefield separation method for distributed optical fiber sensing data.
[0184] Computer program code for performing the operations of this embodiment can be written in one or more programming languages or a combination thereof. These programming languages include object-oriented programming languages—such as Java, Smalltalk, and C++—and conventional procedural programming languages—such as the "C" language or similar programming languages. The program code can be executed entirely on the user's computer, partially on the user's computer, as a standalone software package, partially on the user's computer and partially on a remote computer, or entirely on a remote computer or server. In cases involving remote computers, the remote computer can be connected to the user's computer via any type of network—including a local area network (LAN) or a wide area network (WAN)—or can be connected to an external computer (e.g., via the Internet using an Internet service provider).
[0185] The flowcharts and block diagrams in the accompanying drawings illustrate the architecture, functionality, and operation of possible implementations of systems, methods, and computer program products according to various embodiments of this disclosure. In this regard, each block in a flowchart or block diagram may represent a module, segment, or portion of code containing one or more executable instructions for implementing a specified logical function. It should also be noted that in some alternative implementations, the functions indicated in the blocks may occur in a different order than those indicated in the drawings. For example, two consecutively indicated blocks may actually be executed substantially in parallel, and they may sometimes be executed in reverse order, depending on the functions involved. It should also be noted that each block in the block diagrams and / or flowcharts, and combinations of blocks in the block diagrams and / or flowcharts, can be implemented using a dedicated hardware-based system that performs the specified function or operation, or using a combination of dedicated hardware and computer instructions.
[0186] The modules described in this embodiment can be implemented in software or hardware. The described modules can also be located in a processor; for example, they can be described as:
[0187] The module comprises an acquisition module, a generation module, and a calculation module. For example, the acquisition module can also be described as "a module for acquiring seismic wave data collected by distributed optical fibers for a target area."
[0188] The functions described above in this document can be performed, at least in part, by one or more hardware logic components. For example, exemplary types of hardware logic components that can be used, without limitation, include: Field Programmable Gate Arrays (FPGAs), Application-Specific Integrated Circuits (ASICs), Application Standard Products (ASSPs), System-on-Chip (SoCs), Complex Programmable Logic Devices (CPLDs), and so on.
[0189] The above description is merely a selection of preferred embodiments of this disclosure and an explanation of the technical principles employed. Those skilled in the art should understand that the scope of the invention involved in the embodiments of this disclosure is not limited to technical solutions formed by specific combinations of the above-described technical features, but should also cover other technical solutions formed by arbitrary combinations of the above-described technical features or their equivalents without departing from the above-described inventive concept. For example, technical solutions formed by substituting the above-described features with (but not limited to) technical features with similar functions disclosed in the embodiments of this disclosure.
Claims
1. A wavefield separation method for distributed optical fiber sensor data, characterized in that, include: Acquire DAS-VSP data, which contains aliased wavefields that need to be separated. The sampling channels of the DAS-VSP data are divided into several new channel sets; The new track gather is shaped by applying median filtering to each track to obtain the shaped new track gather. The reshaped new gather is used to recover the gather and obtain the DAS-VSP data after wavefield separation. Output the DAS-VSP data after wavefield separation.
2. The wave field separation method according to claim 1, characterized in that, The step of dividing the sampling channels of the DAS-VSP data into several new channel sets includes: Set a channel interval d, where the original interval of the sampling channels is <d≤20; For the sampling channels of the DAS-VSP data, sampling starts from the first channel with the original number. Then, sampling is performed every channel interval d until the channel interval d is no longer met. Sampling stops then, and the sampling channels are combined into a new channel set according to the order of sampling. For the remaining sampling channels after the previous step, start sampling from the first channel, and sample one channel every d channels until the channel interval d is no longer met. Stop sampling and combine the above-mentioned sampling channels into another new channel set according to the order of sampling. Repeat the previous step until all sampling channels have been extracted to form a new gather, and there are no remaining sampling channels.
3. The wavefield separation method according to claim 1, characterized in that, The step of shaping the new gather by applying median filtering to each channel to obtain the shaped new gather includes: The original travel time of each sampling channel in the new channel set is picked up and compared to obtain a minimum travel time; Based on the minimum travel time, the travel times of all sampling channels are aligned to obtain a new set of channels after alignment; For the new set after alignment, calculate the midpoint of the shaped curve at all sampling times for all sampling channels; The travel time of the sampling channel is restored to the original travel time, and the new set after anti-alignment is used as the new set after shaping.
4. The wavefield separation method according to claim 3, characterized in that, The step of calculating the shaped midpoint of all sampling channels at all sampling times for the aligned new gather includes: For each sampling channel in the newly aligned channel set, perform the following steps: Multiple sampling times are set. At any sampling time t, the current sampling channel c is taken. Taking the sample point (t,c) as the midpoint, a time window w is taken. The median value med and the average value mean of the sample points within the time window w are calculated. The new value of the midpoint is then calculated according to the following formula as the integer midpoint (t',c'): The midpoint of the shape (t', c') = med + (mean - med) × α(I) Where α is the shaping factor, 0 < α < 1.
5. The wave field separation method according to claim 4, characterized in that, The time window w is an odd number greater than or equal to 10 and less than or equal to 100.
6. The wavefield separation method according to any one of claims 1-5, characterized in that, The step of recovering the wavefield-separated DAS-VSP data from the shaped new gather includes: All sampling channels in the reshaped new gather are restored according to their original numbers to obtain the target gather and the DAS-VSP data after wavefield separation.
7. A wavefield separation device for distributed optical fiber data, characterized in that, include: The acquisition module is used to acquire DAS-VSP data; The channel set grouping module is used to divide the sampling channels of the DAS-VSP data into several new channel sets; The shaping module is used to shape the new track set track by track using median filtering to obtain the shaped new track set. The gather group recovery module is used to recover the shaped new gathers to obtain DAS-VSP data after wavefield separation. The output module is used to output the DAS-VSP data after wavefield separation.
8. A computer device, comprising a memory, a processor, and a computer program stored in the memory, characterized in that, The processor executes the computer program to implement the steps of the method according to any one of claims 1 to 6.
9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by a processor, the computer program implements the steps of the method according to any one of claims 1 to 6.
10. A computer program product, comprising a computer program, characterized in that, When executed by a processor, the computer program implements the steps of the method according to any one of claims 1 to 6.