Geological imaging method, apparatus, computer device, and storage medium

Abstract

The application relates to a geological imaging method, device, computer equipment and storage medium. The method comprises the following steps: in the process that a drilling machine moves to the deep part of a main borehole in a target detection area, acquiring a virtual seismic source signal at different depths in the main borehole detected by a main sensor arranged on the drilling machine, and acquiring an observation signal set obtained by detecting a propagated virtual seismic source signal by a secondary sensor arranged in a secondary borehole in the target detection area; performing denoising processing on an observation signal in the observation signal set according to the virtual seismic source signal at different depths, to obtain an ideal signal set; and performing geological imaging on the target detection area according to the ideal signal set based on a velocity model. The method can eliminate noise interference and improve the efficiency of geological detection in the target detection area.

Description

Technical Field

[0001] This application relates to the field of underground space exploration technology, and in particular to a geological imaging method, apparatus, computer equipment, and storage medium. Background Technology

[0002] With the rapid construction and development of densely populated cities, underground engineering projects, primarily consisting of urban subways, integrated utility tunnels, and underground roads, are gradually expanding towards larger-scale structures and more diverse functions. Numerous large-scale urban underground engineering projects have emerged, becoming an important tool for improving urban environmental quality and guiding the optimization of urban layout structure. However, during the construction of underground engineering projects, complex geological conditions such as isolated boulders, karst, and composite strata pose a serious threat to the safe excavation of tunnel boring machines.

[0003] Existing underground space exploration technologies primarily rely on direct borehole drilling, which involves drilling at equal intervals on the surface above underground projects such as urban shield tunnels to image the geology of the target area. However, traditional geological drilling methods struggle to eliminate interference from strong urban noise and are difficult to use for long-distance propagation and detection in shallow, strongly attenuated strata within urban areas, significantly limiting the efficiency of geological exploration. Summary of the Invention

[0004] Therefore, it is necessary to provide a geological imaging method, device, computer equipment, and storage medium to address the aforementioned technical problems, which can eliminate noise interference and improve the efficiency of geological exploration within the target detection area.

[0005] In a first aspect, this application provides a geological imaging method, including:

[0006] During the process of the drilling rig moving to the depth of the main borehole in the target detection area, the virtual seismic source signals at different depths in the main borehole detected by the main sensor deployed on the drilling rig, and the observation signal set obtained by the secondary sensor deployed in the secondary borehole in the target detection area detecting the propagating virtual seismic source signals are acquired.

[0007] Based on the virtual seismic source signals at different depths, the observation signals in the observation signal set are denoised to obtain an ideal signal set;

[0008] Based on the velocity model and according to the ideal signal set, geological imaging is performed on the target detection area.

[0009] In one embodiment, multiple secondary sensors are deployed in the secondary borehole, and different secondary sensors are deployed at different depths of the secondary borehole;

[0010] The virtual seismic source signal at each depth in the main borehole corresponds to an observation signal set; wherein, the observation signal set includes the observation signals obtained by each auxiliary sensor in the auxiliary borehole detecting the virtual seismic source signal at that depth.

[0011] In one embodiment, the observation signals in the observation signal set are denoised based on virtual seismic source signals at different depths to obtain an ideal signal set, including:

[0012] For each virtual seismic source signal at a given depth, cross-correlation is performed on the virtual seismic source signal at that depth and each observation signal in the observation signal set corresponding to the virtual seismic source signal at that depth to obtain the cross-correlation signal between the main sensor and each auxiliary sensor when the main borehole is at that depth.

[0013] Autocorrelation calculations are performed on each observation signal in the observation signal set corresponding to the virtual seismic source signal at that depth to obtain the autocorrelation signals of each sub-sensor;

[0014] Based on the cross-correlation signals between the main sensor and each auxiliary sensor at the depth of the main borehole, and the autocorrelation signals of each auxiliary sensor, the ideal signal set corresponding to the virtual seismic source signal at that depth is determined.

[0015] In one embodiment, based on the cross-correlation signals between the main sensor and each auxiliary sensor at that depth in the main borehole, and the autocorrelation signals of each auxiliary sensor, the ideal signal set corresponding to the virtual source signal at that depth is determined, including:

[0016] For each sensor, the ratio between the cross-correlation signal between the main sensor and the secondary sensor at that depth and the autocorrelation signal of the secondary sensor is used as an ideal signal in the set of ideal signals corresponding to the virtual source signal at that depth.

[0017] In one embodiment, before denoising the observed signals in the observed signal set based on virtual seismic source signals at different depths to obtain an ideal signal set, the method further includes:

[0018] Preprocessing is performed on the virtual source signal at each depth and the observation signal in the observation signal set corresponding to the virtual source signal at each depth; wherein the preprocessing includes at least one of the following: resampling, removal of instrument response, removal of mean, removal of linear trend, bandpass filtering, time-domain normalization, and spectral whitening.

[0019] In one embodiment, preprocessing is performed on the virtual source signal at each depth and the observation signal in the observation signal set corresponding to the virtual source signal at each depth, including:

[0020] According to the preset time window, the virtual seismic source signal at each depth and the observation signal in the observation signal set corresponding to the virtual seismic source signal at each depth are segmented and processed.

[0021] Preprocessing is performed on the virtual seismic source signal at each depth after segmentation and the segmented observation signal in the observation signal set corresponding to the virtual seismic source signal at each depth.

[0022] Secondly, this application also provides a geological imaging device, comprising:

[0023] The signal acquisition model is used to acquire, during the process of the drilling rig moving to the depth of the main borehole in the target detection area, the virtual seismic source signals at different depths in the main borehole detected by the main sensor deployed on the drilling rig, and the set of observation signals obtained by the secondary sensors deployed in the secondary borehole in the target detection area detecting the propagating virtual seismic source signals.

[0024] The signal processing module is used to denoise the observation signals in the observation signal set based on the virtual seismic source signals at different depths to obtain an ideal signal set;

[0025] The geological imaging module is used to perform geological imaging of the target detection area based on a velocity model and an ideal signal set.

[0026] Thirdly, this application also provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to perform the following steps:

[0027] During the process of the drilling rig moving to the depth of the main borehole in the target detection area, the virtual seismic source signals at different depths in the main borehole detected by the main sensor deployed on the drilling rig, and the observation signal set obtained by the secondary sensor deployed in the secondary borehole in the target detection area detecting the propagating virtual seismic source signals are acquired.

[0028] Based on the virtual seismic source signals at different depths, the observation signals in the observation signal set are denoised to obtain an ideal signal set;

[0029] Based on the velocity model and according to the ideal signal set, geological imaging is performed on the target detection area.

[0030] Fourthly, this application also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, performs the following steps:

[0031] During the process of the drilling rig moving to the depth of the main borehole in the target detection area, the virtual seismic source signals at different depths in the main borehole detected by the main sensor deployed on the drilling rig, and the observation signal set obtained by the secondary sensor deployed in the secondary borehole in the target detection area detecting the propagating virtual seismic source signals are acquired.

[0032] Based on the virtual seismic source signals at different depths, the observation signals in the observation signal set are denoised to obtain an ideal signal set;

[0033] Based on the velocity model and according to the ideal signal set, geological imaging is performed on the target detection area.

[0034] Fifthly, this application also provides a computer program product, including a computer program that, when executed by a processor, performs the following steps:

[0035] During the process of the drilling rig moving to the depth of the main borehole in the target detection area, the virtual seismic source signals at different depths in the main borehole detected by the main sensor deployed on the drilling rig, and the observation signal set obtained by the secondary sensor deployed in the secondary borehole in the target detection area detecting the propagating virtual seismic source signals are acquired.

[0036] Based on the virtual seismic source signals at different depths, the observation signals in the observation signal set are denoised to obtain an ideal signal set;

[0037] Based on the velocity model and according to the ideal signal set, geological imaging is performed on the target detection area.

[0038] The aforementioned geological imaging method, apparatus, computer equipment, and storage medium, during the process of the drilling rig moving deeper into the main borehole within the target detection area, acquire virtual seismic source signals at different depths in the main borehole detected by the main sensor deployed on the drilling rig, and observation signal sets obtained by the secondary sensors deployed in the secondary boreholes within the target detection area detecting the propagating virtual seismic source signals. Based on the virtual seismic source signals at different depths, the observation signal sets are denoised to obtain an ideal signal set. Furthermore, based on a velocity model and the ideal signal set, geological imaging of the target detection area is performed. This scheme, through the main and secondary sensors, fully acquires the geological response wavefield information between the main and secondary boreholes; simultaneously, by acquiring virtual seismic source signals at different depths, it improves the detection accuracy of geological information in the target detection area; furthermore, this scheme, through denoising processing, eliminates strong noise interference and fully extracts the effective signals between the main and secondary boreholes, providing a data foundation and reliable means for effective geological imaging within the target detection area. Attached Figure Description

[0039] To more clearly illustrate the technical solutions in the embodiments or related technologies of this application, the accompanying drawings used in the description of the embodiments or related technologies will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of this application. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0040] Figure 1 This is a diagram illustrating the application environment of a geological imaging method in one embodiment;

[0041] Figure 2 This is a flowchart illustrating a geological imaging method in one embodiment;

[0042] Figure 3 This is a flowchart illustrating the process of determining an ideal signal set in one embodiment;

[0043] Figure 4 This is a schematic diagram of the structure of the observation system within the target detection area in one embodiment;

[0044] Figure 5 This is a flowchart illustrating a geological imaging method in another embodiment;

[0045] Figure 6 This is a structural block diagram of a geological imaging device in one embodiment;

[0046] Figure 7 This is a structural block diagram of a geological imaging device in another embodiment;

[0047] Figure 8 This is an internal structural diagram of a computer device in one embodiment. Detailed Implementation

[0048] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application.

[0049] The geological imaging method provided in this application embodiment can be applied to, for example... Figure 1In the application environment shown, sensor 101 communicates with main control device 102 via a network. Sensor 101 includes a main sensor 101-1 and a secondary sensor 101-2. Optionally, during the drilling rig's movement into the depth of the main borehole within the target detection area, the main control device 102 acquires virtual seismic source signals at different depths in the main borehole detected by the main sensor 101-1 deployed on the drilling rig, and observation signal sets obtained by the secondary sensor 101-2 deployed in the secondary borehole within the target detection area detecting the propagating virtual seismic source signals. Based on the virtual seismic source signals at different depths, the observation signal sets are denoised to obtain an ideal signal set. Further, based on a velocity model and the ideal signal set, the main control device 102 performs geological imaging of the target detection area. The main sensor 101-1 and secondary sensor 101-2 can be detectors, and the main control device 102 can be a server or terminal with powerful computing capabilities.

[0050] In one embodiment, such as Figure 2 As shown, a geological imaging method is provided, which can be applied to... Figure 1 Taking the main control device 102 as an example, the specific steps include:

[0051] S201, during the process of the drilling rig moving to the depth of the main borehole in the target detection area, the virtual seismic source signals at different depths in the main borehole detected by the main sensor deployed on the drilling rig, and the observation signal set obtained by the secondary sensors deployed in the secondary borehole in the target detection area detecting the propagating virtual seismic source signals.

[0052] The target detection area refers to the region where geological information needs to be detected; the observation signal set contains the observation signals. The observation signals refer to the signals obtained by the auxiliary sensors in the auxiliary borehole detecting the virtual seismic source signals.

[0053] Optionally, an observation system can be pre-deployed in the target detection area. This system includes a drilling rig with a main sensor, a main borehole, an auxiliary borehole, and auxiliary sensors deployed within the auxiliary borehole. Specifically, a suitable location is selected on the surface of the target detection area, and a vertical borehole is drilled using a drilling rig as an auxiliary borehole. A casing with an auxiliary sensor is placed in the auxiliary borehole; simultaneously, the main sensor is deployed on the drill string.

[0054] Furthermore, as the drilling rig moves deeper into the main borehole within the target detection area, noise signals are generated in the main borehole. These noise signals are used as virtual seismic sources, and the noise signals received by the sensors deployed on the drilling rig are used as virtual seismic source signals. The virtual seismic source signals propagate through the geology of the target detection area. Simultaneously, the secondary sensors deployed in the secondary boreholes within the target detection area detect the propagating virtual seismic source signals to obtain an observation signal set.

[0055] It should be noted that as the drilling depth in the main borehole changes, the auxiliary sensors in the secondary borehole simultaneously detect the virtual seismic source signal. For example, every time the drilling rig descends 20cm in the main borehole, the auxiliary sensors in the secondary borehole detect the virtual seismic source signal once.

[0056] S202, based on the virtual seismic source signals at different depths, performs noise reduction processing on the observation signals in the observation signal set to obtain an ideal signal set.

[0057] Among them, the ideal signal set includes the ideal signal; the ideal signal is the signal detected by the secondary sensor after only passing through the geological propagation of the target detection area, that is, the signal after the observation signal has been denoised.

[0058] Optionally, seismic interferometry can be used to denoise the observation signals in the observation signal set based on virtual seismic source signals at different depths. This removes interference signals from the observation signals in the observation signal set, resulting in a virtual seismic signal set, i.e., an ideal signal set.

[0059] S203, based on a velocity model, performs geological imaging of the target detection area according to an ideal signal set.

[0060] Optionally, a velocity model can be constructed, utilizing various principles of seismic wave propagation, such as Huygens' principle, Fermat's principle, and interchange principle. Based on an ideal signal set, the constructed velocity model can be forward modeled using linear interpolation ray tracing to model the travel time between the excitation point and the receiver point, thereby tracing the propagation path of the rays in the geology of the target detection area and then performing geological imaging of the target detection area.

[0061] In the aforementioned geological imaging method, as the drilling rig moves deeper into the main borehole within the target detection area, it acquires virtual seismic source signals at different depths within the main borehole detected by the main sensor deployed on the drilling rig, as well as observation signal sets obtained by the secondary sensors deployed in the secondary boreholes within the target detection area detecting the propagating virtual seismic source signals. Based on the virtual seismic source signals at different depths, the observation signals in the observation signal sets are denoised to obtain an ideal signal set. Furthermore, based on a velocity model and the ideal signal set, geological imaging of the target detection area is performed. This scheme, through the main and secondary sensors, fully acquires the geological response wavefield information between the main and secondary boreholes. Simultaneously, by acquiring virtual seismic source signals at different depths, it improves the detection accuracy of geological information in the target detection area. Furthermore, this scheme, through denoising processing, eliminates strong noise interference and fully extracts the effective signals between the main and secondary boreholes, providing a data foundation and reliable means for effective geological imaging within the target detection area.

[0062] It should be noted that, in order to improve the accuracy of geological information detection in the target detection area, multiple secondary sensors can be deployed in the secondary borehole, with different secondary sensors deployed at different depths in the secondary borehole; the virtual seismic source signal at each depth in the main borehole corresponds to an observation signal set, that is, each depth corresponds to a virtual seismic source signal and an observation signal set; wherein, the observation signal set includes the observation signals obtained by each secondary sensor in the secondary borehole detecting the virtual seismic source signal at that depth. Based on the above embodiments, in one embodiment, such as Figure 3 As shown, a method for determining an ideal signal set is provided, which specifically includes the following steps:

[0063] S301, for the virtual seismic source signal at each depth, perform cross-correlation calculation on the virtual seismic source signal at that depth and each observation signal in the observation signal set corresponding to the virtual seismic source signal at that depth to obtain the cross-correlation signal between the main sensor and each auxiliary sensor when the main borehole is at that depth.

[0064] Optionally, for each virtual seismic source signal at a given depth, a cross-correlation operation is performed between the virtual seismic source signal at that depth and each observation signal in the corresponding observation signal set. This can be achieved using the following formula (1):

[0065] (1)

[0066] in, This indicates that when the main borehole is at this depth, the main sensor and the... Cross-correlation signals between individual sensors; This represents the virtual seismic source signal at that depth; This represents the first observation signal in the set corresponding to the virtual seismic source signal at that depth. One observed signal.

[0067] S302, perform autocorrelation calculation on each observation signal in the observation signal set corresponding to the virtual seismic source signal at this depth to obtain the autocorrelation signal of each sub-sensor.

[0068] Optionally, autocorrelation calculation is performed on each observation signal in the observation signal set corresponding to the virtual seismic source signal at this depth, which can be achieved through the following formulas (2) and (3):

[0069] (2)

[0070] (3)

[0071] in, Indicates the first The autocorrelation signal of each secondary sensor; For regularization terms; The regularization factor can be determined based on the percentage of the average energy of the virtual source signal spectrum at that depth.

[0072] S303. Based on the cross-correlation signals between the main sensor and each auxiliary sensor at this depth, and the autocorrelation signals of each auxiliary sensor, determine the ideal signal set corresponding to the virtual seismic source signal at this depth.

[0073] Optionally, by performing corresponding calculations on the cross-correlation signals between the main sensor and each auxiliary sensor at this depth, as well as the autocorrelation signals of each auxiliary sensor, the ideal signal set corresponding to the virtual seismic source signal at this depth can be obtained.

[0074] Specifically, to improve the accuracy of the ideal signal in the ideal signal set, for each sensor, the ratio between the cross-correlation signal between the main sensor and the auxiliary sensor at that depth and the autocorrelation signal of the auxiliary sensor can be used as an ideal signal in the ideal signal set corresponding to the virtual source signal at that depth. Specifically, the ideal signal can be determined using the following formula (4):

[0075] (4)

[0076] in, This represents the first ideal signal in the set corresponding to the virtual seismic source signal at that depth. An ideal signal. This provides a method for accurately determining an ideal signal from a set of ideal signals corresponding to a virtual seismic source signal at each depth.

[0077] In this embodiment, through autocorrelation and cross-correlation operations, the signals and noise common to each observation signal in the observation signal set corresponding to the virtual seismic source signal are preserved, while the incoherent noise components that exist independently outside each observation signal in the observation signal set corresponding to the virtual seismic source signal are attenuated.

[0078] It should be noted that, in order to improve the signal quality of the virtual source signals at each depth and the corresponding observation signal sets at each depth, and to remove some interference, before denoising the observation signals at each depth and obtaining the ideal signal set based on the virtual source signals at different depths, preprocessing is required for both the virtual source signals at each depth and the corresponding observation signal sets at each depth. This preprocessing includes at least one of the following: resampling, removal of instrument response, removal of mean, removal of linear trends, bandpass filtering, time-domain normalization, and spectral whitening. Understandably, by introducing preprocessing, strong vibrations generated by the operation of the drilling rig motor and other components can be removed, signal quality can be controlled, and the observation signals, primarily carrying geological information wavefields, can be preserved.

[0079] Furthermore, to improve the anti-interference capability of the signals during preprocessing, the virtual source signals at each depth and the observation signals in the corresponding observation signal sets at each depth can be segmented beforehand. Optionally, the segmentation process is as follows: according to a preset time window, the virtual source signals at each depth and the observation signals in the corresponding observation signal sets at each depth are segmented; preprocessing is then performed on the segmented virtual source signals at each depth and the segmented observation signals in the corresponding observation signal sets at each depth. The size of the time window can be set according to the propagation time of the virtual source signals; for example, the time window size can be set to 3-4 times the propagation time of the virtual source signals. It is understood that by introducing the segmentation process, the anti-interference capability of the virtual source signals at each depth and the observation signals in the corresponding observation signal sets at each depth is improved.

[0080] In one embodiment, such as Figure 4 As shown, taking the target detection area as an urban environment as the application scenario, a schematic diagram of the observation system within the target detection area is provided. Among them, Figure 4 The 102 module is the main control device; the 402 module is the sensor data transmission and synchronization system, used to transmit the virtual seismic source signal detected by the main sensor and the observation signals detected by each auxiliary sensor; furthermore, Figure 4 Multiple secondary sensors were deployed in the middle, namely Figure 4 401 in the diagram is used to detect the virtual seismic source signal propagating in the main borehole. Optionally, the secondary sensor represented by 401 can... Figure 1The secondary sensor shown in 101-2 is presented. 403 is the drill bit of the drilling rig in the main borehole, and the main sensor (not shown in the figure) is deployed on the drill bit's shaft; 405 is the surface overburden layer; 406 is the deep underground rock and soil structure; 407 is the surface drilling rig control system, used to control the movement of the drill bit in the deep part of the main borehole; 408 is unfavorable geological bodies such as boulders.

[0081] Optionally, an observation system can be deployed in advance in the target detection area. Specifically, a suitable location can be selected on the ground surface of the target detection area to drill a vertical low-level borehole as a secondary borehole, and a casing with multiple secondary sensors can be placed in the secondary borehole. At the same time, the main sensor can be deployed on the drill string.

[0082] Furthermore, during the movement of the drill bit deep within the main borehole, the surface drilling rig control system generates noise. The main sensor in the main borehole detects this noise as a virtual seismic source signal. This virtual seismic source signal propagates through the geology, and the auxiliary sensors in the secondary boreholes detect the propagating virtual seismic source signal to obtain an observation signal set. The main control equipment, through a sensor transmission and synchronization system, acquires virtual seismic source signals at different depths, as well as the observation signal sets obtained from the detection of the propagating virtual seismic source signals by the auxiliary sensors deployed in the secondary boreholes within the target detection area. Based on the virtual seismic source signals at different depths and the corresponding observation signal sets, geological imaging of the target detection area can be performed.

[0083] Figure 5 This is a flowchart illustrating a geological imaging method in another embodiment. Based on the above embodiments, this embodiment provides an optional example of a geological imaging method. (Combined with...) Figure 5 The specific implementation process is as follows:

[0084] S501, during the process of the drilling rig moving to the depth of the main borehole in the target detection area, acquire the virtual seismic source signals at different depths in the main borehole detected by the main sensor deployed on the drilling rig, and the set of observation signals obtained by each auxiliary sensor deployed in the auxiliary borehole in the target detection area detecting each propagating virtual seismic source signal.

[0085] In this system, multiple secondary sensors are deployed in the secondary borehole, and different secondary sensors are deployed at different depths in the secondary borehole; the virtual seismic source signal at each depth in the main borehole corresponds to an observation signal set; the observation signal set includes the observation signals obtained by each secondary sensor in the secondary borehole detecting the virtual seismic source signal at that depth.

[0086] S502, according to the preset time window, performs segmented processing on the virtual seismic source signal at each depth and the observation signal in the corresponding observation signal set at each depth.

[0087] S503 preprocesses the virtual source signal at each depth after segmentation and the segmented observation signal in the corresponding observation signal set at each depth.

[0088] The preprocessing includes at least one of the following: resampling, removal of instrument response, removal of mean, removal of linear trend, bandpass filtering, time-domain normalization, and spectral whitening.

[0089] S504 performs cross-correlation calculations on the virtual seismic source signal and the corresponding observation signal at each depth after preprocessing, and obtains the cross-correlation signal between the main sensor and each auxiliary sensor at each depth of the main borehole.

[0090] S505 performs autocorrelation calculation on each observation signal after centralized preprocessing of the observation signals corresponding to each depth, to obtain the autocorrelation signal of each sub-sensor corresponding to each depth.

[0091] S506. Based on the cross-correlation signals between the main sensor and each auxiliary sensor at each depth of the main borehole, and the autocorrelation signals of each auxiliary sensor at each depth, determine the ideal signal set corresponding to the virtual source signal at each depth.

[0092] S507, based on a velocity model, performs geological imaging of the target detection area according to various ideal signal sets.

[0093] The specific processes of S501-S507 described above can be found in the description of the above method embodiments. Their implementation principles and technical effects are similar, and will not be repeated here.

[0094] It should be understood that although the steps in the flowcharts of the embodiments described above are shown sequentially according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the flowcharts of the embodiments described above may include multiple steps or multiple stages. These steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these steps or stages is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the steps or stages of other steps.

[0095] Based on the same inventive concept, this application also provides a geological imaging apparatus for implementing the geological imaging method described above. The solution provided by this apparatus is similar to the implementation scheme described in the above method; therefore, the specific limitations in one or more geological imaging apparatus embodiments provided below can be found in the limitations of the geological imaging method described above, and will not be repeated here.

[0096] In one embodiment, such as Figure 6 As shown, a geological imaging device 1 is provided, comprising: a signal acquisition module 10, a signal processing module 20, and a geological imaging module 30, wherein:

[0097] Signal acquisition model 10 is used to acquire, during the process of the drilling rig moving to the depth of the main borehole in the target detection area, the virtual seismic source signals at different depths in the main borehole detected by the main sensor deployed on the drilling rig, and the observation signal set obtained by the secondary sensors deployed in the secondary borehole in the target detection area detecting the propagating virtual seismic source signals.

[0098] The signal denoising module 20 is used to denoise the observation signals in the observation signal set based on the virtual seismic source signals at different depths, so as to obtain an ideal signal set.

[0099] The geological imaging module 30 is used to perform geological imaging of the target detection area based on a velocity model and an ideal signal set.

[0100] The aforementioned geological imaging device, during the drilling rig's movement into the depths of the main borehole within the target detection area, acquires virtual seismic source signals at different depths within the main borehole detected by the main sensor deployed on the drilling rig, as well as observation signal sets obtained by the secondary sensors deployed in the auxiliary boreholes within the target detection area detecting the propagating virtual seismic source signals. Based on the virtual seismic source signals at different depths, the observation signal sets are denoised to obtain an ideal signal set. Furthermore, based on a velocity model and the ideal signal set, geological imaging of the target detection area is performed. This scheme, through the main and secondary sensors, fully acquires the geological response wavefield information between the main and auxiliary boreholes; simultaneously, by acquiring virtual seismic source signals at different depths, it improves the detection accuracy of geological information in the target detection area; furthermore, this scheme, through denoising processing, eliminates strong noise interference and fully extracts the effective signals between the main and auxiliary boreholes, providing a data foundation and reliable means for effective geological imaging within the target detection area.

[0101] In one embodiment, multiple secondary sensors are deployed in the secondary borehole, and different secondary sensors are deployed at different depths in the secondary borehole; the virtual seismic source signal at each depth in the main borehole corresponds to an observation signal set; wherein, the observation signal set includes the observation signals obtained by each secondary sensor in the secondary borehole detecting the virtual seismic source signal at that depth.

[0102] In one embodiment, in Figure 6 On the basis of, such as Figure 7 As shown, the signal denoising module 20 includes:

[0103] The first arithmetic unit 21 is used to perform cross-correlation calculation on the virtual seismic source signal at each depth and each observation signal in the observation signal set corresponding to the virtual seismic source signal at that depth, so as to obtain the cross-correlation signal between the main sensor and each auxiliary sensor when the main borehole is at that depth.

[0104] The second processing unit 22 is used to perform autocorrelation calculation on each observation signal in the observation signal set corresponding to the virtual seismic source signal at this depth, so as to obtain the autocorrelation signal of each sub-sensor.

[0105] The third calculation unit 23 is used to determine the ideal signal set corresponding to the virtual source signal at that depth based on the cross-correlation signal between the main sensor and each auxiliary sensor when the main borehole is at that depth, as well as the autocorrelation signal of each auxiliary sensor.

[0106] In one embodiment, the third arithmetic unit 23 is specifically used for:

[0107] For each sensor, the ratio between the cross-correlation signal between the main sensor and the secondary sensor at that depth and the autocorrelation signal of the secondary sensor is used as an ideal signal in the set of ideal signals corresponding to the virtual source signal at that depth.

[0108] In one embodiment, the geological imaging device 1 further includes:

[0109] The preprocessing module is used to preprocess the virtual source signal at each depth and the observation signal in the observation signal set corresponding to the virtual source signal at each depth, respectively; wherein, the preprocessing includes at least one of the following: resampling, removal of instrument response, removal of mean, removal of linear trend, bandpass filtering, time domain normalization and spectral whitening.

[0110] In one embodiment, the preprocessing module is specifically used for:

[0111] According to the preset time window, the virtual source signal at each depth and the observation signal in the observation signal set corresponding to the virtual source signal at each depth are segmented and processed; the segmented virtual source signal at each depth and the segmented observation signal in the observation signal set corresponding to the virtual source signal at each depth are preprocessed respectively.

[0112] Each module in the aforementioned geological imaging device can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in or independent of the processor in a computer device, or stored in the memory of a computer device as software, so that the processor can call and execute the operations corresponding to each module.

[0113] In one embodiment, a computer device is provided, which may be a server, and its internal structure diagram may be as follows: Figure 8 As shown, this computer device includes a processor, memory, input / output (I / O) interfaces, and a communication interface. The processor, memory, and I / O interfaces are connected via a system bus, and the communication interface is also connected to the system bus via the I / O interfaces. The processor provides computational and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system, computer programs, and a database. The internal memory provides the environment for the operating system and computer programs stored in the non-volatile storage media. The database stores virtual seismic source signals at different depths and sets of observation signals obtained by detecting propagating virtual seismic source signals. The I / O interfaces are used for exchanging information between the processor and external devices. The communication interface is used for communication with external terminals via a network connection. When the computer program is executed by the processor, it implements a geological imaging method.

[0114] Those skilled in the art will understand that Figure 8 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0115] In one embodiment, a computer device is provided, including a memory and a processor, wherein the memory stores a computer program, and the processor executes the computer program to perform the following steps:

[0116] During the process of the drilling rig moving to the depth of the main borehole in the target detection area, the virtual seismic source signals at different depths in the main borehole detected by the main sensor deployed on the drilling rig, as well as the observation signal set obtained by the secondary sensors deployed in the secondary borehole in the target detection area detecting the propagating virtual seismic source signals;

[0117] Based on the virtual seismic source signals at different depths, the observed signals in the observation signal set are denoised to obtain an ideal signal set;

[0118] Based on the velocity model and the ideal signal set, geological imaging is performed on the target detection area.

[0119] In one embodiment, when the processor executes a computer program, multiple secondary sensors are deployed in the secondary borehole, and different secondary sensors are deployed at different depths in the secondary borehole; the virtual seismic source signal at each depth in the main borehole corresponds to an observation signal set; wherein, the observation signal set includes the observation signal obtained by each secondary sensor in the secondary borehole detecting the virtual seismic source signal at that depth.

[0120] In one embodiment, when the processor executes a computer program to denoise the observed signals in the observed signal set based on virtual seismic source signals at different depths to obtain an ideal signal set, it also performs the following steps:

[0121] For each virtual seismic source signal at a given depth, cross-correlation is performed on the virtual seismic source signal at that depth and each observation signal in the corresponding observation signal set to obtain the cross-correlation signal between the main sensor and each auxiliary sensor at that depth in the main borehole. Autocorrelation is then performed on each observation signal in the corresponding observation signal set to obtain the autocorrelation signal of each auxiliary sensor. Based on the cross-correlation signal between the main sensor and each auxiliary sensor at that depth in the main borehole, and the autocorrelation signal of each auxiliary sensor, the ideal signal set corresponding to the virtual seismic source signal at that depth is determined.

[0122] In one embodiment, when the processor executes a computer program to determine the ideal signal set corresponding to the virtual source signal at that depth based on the cross-correlation signals between the main sensor and each auxiliary sensor at that depth in the main borehole, as well as the autocorrelation signals of each auxiliary sensor, the processor also performs the following steps:

[0123] For each sensor, the ratio between the cross-correlation signal between the main sensor and the secondary sensor at that depth and the autocorrelation signal of the secondary sensor is used as an ideal signal in the set of ideal signals corresponding to the virtual source signal at that depth.

[0124] In one embodiment, before the processor executes a computer program to denoise the observed signals in the observed signal set based on virtual seismic source signals at different depths to obtain an ideal signal set, the following steps are also performed:

[0125] Preprocessing is performed on the virtual source signal at each depth and the observation signal in the observation signal set corresponding to the virtual source signal at each depth. The preprocessing includes at least one of the following: resampling, removal of instrument response, removal of mean, removal of linear trend, bandpass filtering, time-domain normalization, and spectral whitening.

[0126] In one embodiment, when the processor executes a computer program to preprocess the virtual seismic source signal at each depth and the observation signal in the observation signal set corresponding to the virtual seismic source signal at each depth, it also performs the following steps:

[0127] According to the preset time window, the virtual source signal at each depth and the observation signal in the observation signal set corresponding to the virtual source signal at each depth are segmented and processed; the segmented virtual source signal at each depth and the segmented observation signal in the observation signal set corresponding to the virtual source signal at each depth are preprocessed respectively.

[0128] In one embodiment, a computer-readable storage medium is provided having a computer program stored thereon, the computer program performing the following steps when executed by a processor:

[0129] During the process of the drilling rig moving to the depth of the main borehole in the target detection area, the virtual seismic source signals at different depths in the main borehole detected by the main sensor deployed on the drilling rig, as well as the observation signal set obtained by the secondary sensors deployed in the secondary borehole in the target detection area detecting the propagating virtual seismic source signals;

[0130] Based on the virtual seismic source signals at different depths, the observed signals in the observation signal set are denoised to obtain an ideal signal set;

[0131] Based on the velocity model and the ideal signal set, geological imaging is performed on the target detection area.

[0132] In one embodiment, when the computer program is executed by the processor, multiple secondary sensors are deployed in the secondary borehole, and different secondary sensors are deployed at different depths in the secondary borehole; the virtual seismic source signal at each depth in the main borehole corresponds to an observation signal set; wherein, the observation signal set includes the observation signal obtained by each secondary sensor in the secondary borehole detecting the virtual seismic source signal at that depth.

[0133] In one embodiment, the computer program performs noise reduction processing on the observed signals in the observed signal set based on virtual seismic source signals at different depths to obtain an ideal signal set. When the processor executes the ideal signal set, it also performs the following steps:

[0134] For each virtual seismic source signal at a given depth, cross-correlation is performed on the virtual seismic source signal at that depth and each observation signal in the corresponding observation signal set to obtain the cross-correlation signal between the main sensor and each auxiliary sensor at that depth in the main borehole. Autocorrelation is then performed on each observation signal in the corresponding observation signal set to obtain the autocorrelation signal of each auxiliary sensor. Based on the cross-correlation signal between the main sensor and each auxiliary sensor at that depth in the main borehole, and the autocorrelation signal of each auxiliary sensor, the ideal signal set corresponding to the virtual seismic source signal at that depth is determined.

[0135] In one embodiment, when the computer program determines the ideal signal set corresponding to the virtual source signal at that depth based on the cross-correlation signals between the main sensor and each auxiliary sensor at that depth, and the autocorrelation signals of each auxiliary sensor, the process further implements the following steps:

[0136] For each sensor, the ratio between the cross-correlation signal between the main sensor and the secondary sensor at that depth and the autocorrelation signal of the secondary sensor is used as an ideal signal in the set of ideal signals corresponding to the virtual source signal at that depth.

[0137] In one embodiment, before the computer program performs noise reduction processing on the observed signals in the observed signal set based on virtual seismic source signals at different depths to obtain an ideal signal set, the following steps are also performed:

[0138] Preprocessing is performed on the virtual source signal at each depth and the observation signal in the observation signal set corresponding to the virtual source signal at each depth. The preprocessing includes at least one of the following: resampling, removal of instrument response, removal of mean, removal of linear trend, bandpass filtering, time-domain normalization, and spectral whitening.

[0139] In one embodiment, when the computer program preprocesses the virtual seismic source signal at each depth and the observation signal in the observation signal set corresponding to the virtual seismic source signal at each depth, it also performs the following steps when executed by the processor:

[0140] According to the preset time window, the virtual source signal at each depth and the observation signal in the observation signal set corresponding to the virtual source signal at each depth are segmented and processed; the segmented virtual source signal at each depth and the segmented observation signal in the observation signal set corresponding to the virtual source signal at each depth are preprocessed respectively.

[0141] In one embodiment, a computer program product is provided, including a computer program that, when executed by a processor, performs the following steps:

[0142] During the process of the drilling rig moving to the depth of the main borehole in the target detection area, the virtual seismic source signals at different depths in the main borehole detected by the main sensor deployed on the drilling rig, as well as the observation signal set obtained by the secondary sensors deployed in the secondary borehole in the target detection area detecting the propagating virtual seismic source signals;

[0143] Based on the virtual seismic source signals at different depths, the observed signals in the observation signal set are denoised to obtain an ideal signal set;

[0144] Based on the velocity model and the ideal signal set, geological imaging is performed on the target detection area.

[0145] In one embodiment, when the computer program is executed by the processor, multiple secondary sensors are deployed in the secondary borehole, and different secondary sensors are deployed at different depths in the secondary borehole; the virtual seismic source signal at each depth in the main borehole corresponds to an observation signal set; wherein, the observation signal set includes the observation signal obtained by each secondary sensor in the secondary borehole detecting the virtual seismic source signal at that depth.

[0146] In one embodiment, the computer program performs noise reduction processing on the observed signals in the observed signal set based on virtual seismic source signals at different depths to obtain an ideal signal set. When the processor executes the ideal signal set, it also performs the following steps:

[0147] For each virtual seismic source signal at a given depth, cross-correlation is performed on the virtual seismic source signal at that depth and each observation signal in the corresponding observation signal set to obtain the cross-correlation signal between the main sensor and each auxiliary sensor at that depth in the main borehole. Autocorrelation is then performed on each observation signal in the corresponding observation signal set to obtain the autocorrelation signal of each auxiliary sensor. Based on the cross-correlation signal between the main sensor and each auxiliary sensor at that depth in the main borehole, and the autocorrelation signal of each auxiliary sensor, the ideal signal set corresponding to the virtual seismic source signal at that depth is determined.

[0148] In one embodiment, when the computer program determines the ideal signal set corresponding to the virtual source signal at that depth based on the cross-correlation signals between the main sensor and each auxiliary sensor at that depth, and the autocorrelation signals of each auxiliary sensor, the process further implements the following steps:

[0149] For each sensor, the ratio between the cross-correlation signal between the main sensor and the secondary sensor at that depth and the autocorrelation signal of the secondary sensor is used as an ideal signal in the set of ideal signals corresponding to the virtual source signal at that depth.

[0150] In one embodiment, before the computer program performs noise reduction processing on the observed signals in the observed signal set based on virtual seismic source signals at different depths to obtain an ideal signal set, the following steps are also performed:

[0151] Preprocessing is performed on the virtual source signal at each depth and the observation signal in the observation signal set corresponding to the virtual source signal at each depth. The preprocessing includes at least one of the following: resampling, removal of instrument response, removal of mean, removal of linear trend, bandpass filtering, time-domain normalization, and spectral whitening.

[0152] In one embodiment, when the computer program preprocesses the virtual seismic source signal at each depth and the observation signal in the observation signal set corresponding to the virtual seismic source signal at each depth, it also performs the following steps when executed by the processor:

[0153] According to the preset time window, the virtual source signal at each depth and the observation signal in the observation signal set corresponding to the virtual source signal at each depth are segmented and processed; the segmented virtual source signal at each depth and the segmented observation signal in the observation signal set corresponding to the virtual source signal at each depth are preprocessed respectively.

[0154] It should be noted that the data involved in this application (including but not limited to data used for analysis, data stored, data displayed, etc.) are all data that have been fully authorized by all parties, and the collection, use and processing of the relevant data must comply with relevant regulations.

[0155] Those skilled in the art will understand that all or part of the processes in the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium. When executed, the computer program can include the processes of the embodiments described above. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM). The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, etc., and are not limited to these.

[0156] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0157] The embodiments described above are merely illustrative of several implementation methods of this application, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of this patent application. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this application should be determined by the appended claims.

Claims

1. A geological imaging method, characterized in that, The method includes: During the process of the drilling rig moving to the depth of the main borehole in the target detection area, the virtual seismic source signals at different depths in the main borehole detected by the main sensor deployed on the drilling rig, and the observation signal set obtained by the secondary sensor deployed in the secondary borehole in the target detection area detecting the propagating virtual seismic source signals are acquired. Based on the virtual seismic source signals at different depths, the observation signals in the observation signal set are denoised to obtain an ideal signal set; Based on the velocity model and according to the ideal signal set, geological imaging is performed on the target detection area.

2. The method according to claim 1, characterized in that, Multiple secondary sensors are deployed in the secondary borehole, and different secondary sensors are deployed at different depths in the secondary borehole; The virtual seismic source signal at each depth in the main borehole corresponds to an observation signal set; wherein, the observation signal set includes the observation signals obtained by each auxiliary sensor in the auxiliary borehole detecting the virtual seismic source signal at that depth.

3. The method according to claim 2, characterized in that, The step of denoising the observed signals in the observed signal set based on virtual seismic source signals at different depths to obtain an ideal signal set includes: For each virtual seismic source signal at a given depth, cross-correlation is performed on the virtual seismic source signal at that depth and each observation signal in the observation signal set corresponding to the virtual seismic source signal at that depth to obtain the cross-correlation signal between the main sensor and each auxiliary sensor when the main borehole is at that depth. Autocorrelation calculations are performed on each observation signal in the observation signal set corresponding to the virtual seismic source signal at that depth to obtain the autocorrelation signals of each sub-sensor; Based on the cross-correlation signals between the main sensor and each auxiliary sensor at the depth of the main borehole, and the autocorrelation signals of each auxiliary sensor, the ideal signal set corresponding to the virtual seismic source signal at that depth is determined.

4. The method according to claim 3, characterized in that, The step of determining the ideal signal set corresponding to the virtual seismic source signal at that depth based on the cross-correlation signals between the main sensor and each auxiliary sensor at that depth in the main borehole, and the autocorrelation signals of each auxiliary sensor, includes: For each sensor, the ratio between the cross-correlation signal between the main sensor and the secondary sensor at that depth and the autocorrelation signal of the secondary sensor is used as an ideal signal in the set of ideal signals corresponding to the virtual source signal at that depth.

5. The method according to claim 2, characterized in that, Before denoising the observed signals in the observed signal set based on virtual seismic source signals at different depths to obtain an ideal signal set, the method further includes: Preprocessing is performed on the virtual source signal at each depth and the observation signal in the observation signal set corresponding to the virtual source signal at each depth; wherein the preprocessing includes at least one of the following: resampling, removal of instrument response, removal of mean, removal of linear trend, bandpass filtering, time-domain normalization, and spectral whitening.

6. The method according to claim 5, characterized in that, The preprocessing of the virtual seismic source signal at each depth and the observation signal set corresponding to the virtual seismic source signal at each depth includes: According to the preset time window, the virtual seismic source signal at each depth and the observation signal in the observation signal set corresponding to the virtual seismic source signal at each depth are segmented and processed. Preprocessing is performed on the virtual seismic source signal at each depth after segmentation and the segmented observation signal in the observation signal set corresponding to the virtual seismic source signal at each depth.

7. A geological imaging device, characterized in that, The device includes: The signal acquisition model is used to acquire, during the process of the drilling rig moving to the depth of the main borehole in the target detection area, the virtual seismic source signals at different depths in the main borehole detected by the main sensor deployed on the drilling rig, and the set of observation signals obtained by the secondary sensors deployed in the secondary borehole in the target detection area detecting the propagating virtual seismic source signals. The signal processing module is used to denoise the observation signals in the observation signal set based on the virtual seismic source signals at different depths to obtain an ideal signal set; The geological imaging module is used to perform geological imaging of the target detection area based on a velocity model and an ideal signal set.

8. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements 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 the computer program is executed by a processor, it 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 the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 6.

Images