Horizontal component orientation method of shallow well geophone based on stable source signal

Abstract

The present invention belongs to the field of microseismic geophone orientation, and specifically relates to a method for orienting the horizontal component of shallow well seismic geophones based on stable source signals, aiming to solve the problems of inaccurate correction angles and low detection quality of geophones. The present invention includes: calculating the relative orientation of each interference source signal relative to each geophone; calculating the waveform matching coefficient of each representative record, screening the representative record based on the waveform matching coefficient, and obtaining each stability record detected by each geophone after screening; the interference source signal corresponding to each stability record is recorded as a stable source signal; superimposing all stability records obtained by each geophone to obtain a superimposed channel record with enhanced signal-to-noise ratio of each geophone; calculating the calibration angle of each geophone relative to each stable source signal, and calculating the calibration angle of the horizontal component of each geophone based on each calibration angle. The present invention improves the calibration angle of the geophone and improves the detection quality.

Description

Background Art

[0002] Industrial production activities, such as unconventional oil and gas reservoir fracturing, coal mining, underground gas storage injection and production, carbon dioxide storage and reservoir water storage, etc., will cause changes in the stress and strain state of the formation, resulting in rock fractures or pre-existing faults or fractures that become unstable and slide, causing earthquakes. In order to understand the state of production activities and avoid geological disasters, it is necessary to use microseismic monitoring technology to conduct real-time dynamic monitoring and evaluation of industrial production activities. Three-component sensors are often used to collect seismic monitoring data. Since the horizontal component azimuth of the detector is random during the process of lowering the detector into the well, the azimuth of the detector needs to be determined by the active source event in the well (perforation event, ball-casting sleeve or detonating cord signal, etc.) or the active source event on the ground (controllable source, heavy hammer or explosive source) at a known position. The existing technology usually analyzes the polarization of the direct P wave of the active source event based on the position coordinates of the active source event relative to the detector, and calculates the horizontal azimuth of the detector by rotating the detector that receives the seismic wave signal. Methods for obtaining the polarization angle of multi-component records include polarization analysis methods, linear fitting methods, and waveform cross-correlation methods. The results calculated by these methods will be affected by the signal-to-noise ratio of active source events, resulting in certain errors in the polarization angle, which in turn affects the orientation of the horizontal component of the detector and the azimuth analysis of the source event. For the case where the signal-to-noise ratio of active source events in the monitoring records is low, and the case where no effective active source event records are detected, the effect of the above-mentioned invention method is not ideal, and it cannot even be used. The observation method of shallow well seismic detectors has the advantages of low equipment cost, no dependence on adjacent wells, and suitability for long-term monitoring, and has broad application prospects. However, in the process of shallow well monitoring, there are often long-term continuous production activities on the surface. The stable background signals generated by them reduce the quality of monitoring data and the signal-to-noise ratio of active source events, making it difficult for existing technologies to meet technical requirements.

[0003] Therefore, the existing shallow well seismic detector horizontal component orientation method has the problem of inaccurate correction angle of the detector and low detection quality due to low monitoring data quality and large active source signal noise ratio. Summary of the invention

[0004] In order to solve the above-mentioned problems in the prior art, that is, the prior art shallow well seismic detector horizontal component orientation method has the problem that the calibration angle of the detector is inaccurate and the detection quality is not high due to the low quality of monitoring data and the small signal-to-noise ratio of the active source signal, the present invention provides a shallow well seismic detector horizontal component orientation method based on a stable source signal, the method comprising:

[0005] Calculate the relative position of each interference source signal with respect to each detector;

[0006] Obtain representative records of each interference source signal detected by each detector; calculate the waveform matching coefficient of each representative record, screen each representative record based on the waveform matching coefficient, and obtain each stability record detected by each detector after screening; the interference source signal corresponding to each stability record is recorded as a stable source signal;

[0007] All stability records acquired by each detector are superimposed to obtain a superimposed channel record with enhanced signal-to-noise ratio of each detector;

[0008] Based on the superposition channel records of each detector, the covariance matrix of each detector is constructed and covariance analysis is performed to obtain the covariance matrix analysis results; based on the covariance matrix analysis results of each detector and the relative azimuth of each stable source signal relative to each detector, the calibration angle of each detector relative to each stable source signal is calculated; based on the calibration angle of each detector relative to each stable source signal, the correction angle of the horizontal component of each detector is calculated.

[0009] The method for calculating the relative position of each interference source signal relative to each detector is:

[0010] Taking the north direction as the reference direction, clockwise is positive. m The interference source signal is relative to the n The orientation of the detectors is:

[0011] ;

[0012] in, m is the interference source number, n is the detector number, sx m is the horizontal position coordinate of the interference source signal, sy m is the longitudinal position coordinate of the interference source signal, rx m is the transverse position coordinate of the detector, ry m is the longitudinal position coordinate of the detector, Indicates m The interference source signal is relative to the n The position of the detector, Indicates direction.

[0013] In a preferred embodiment, the method for obtaining the stability record of each interference source signal detected by each detector includes:

[0014] First, calculate the waveform matching coefficient of each representative record; the waveform coefficient of the representative record of the mth interference source signal obtained by the nth detector is:

[0015] ;

[0016] Among them, c is the number of the three components, x, y, z are the three directions of the three components respectively, For the n The first detector obtains m The waveform coefficients of a representative record of an interferer signal, i Represents the sampling point number, j Indicates the sampling point number in the time window, Indicates that the sampling point number in the time window is j The representative record of the cth trisection corresponding to the sampling point in the time window, represents a representative record of the cth component obtained by the nth detector; Nwin represents the time window length, n Number the detector;

[0017] No. n The first detector obtains m Waveform matching coefficients of representative recorded waveform coefficients of interference source signals When it is greater than the threshold E, the current representative record is retained and used as the stability record obtained by the nth detector, and the interference source signal corresponding to each stability record is recorded as the stable source signal.

[0018] In a preferred embodiment, the method for obtaining the superimposed trace record with enhanced signal-to-noise ratio of each detector is:

[0019] For the n detectors, the number of stability records detected in a certain period of time is K; each stability record corresponds to a stable source signal;

[0020] No. n The kth stability record detected by the detector is denoted by ;

[0021] The signal-to-noise ratio enhanced stack trace record of the nth detector is:

[0022] ;

[0023] in, For the n The signal-to-noise ratio of the detector is enhanced by the stacked channel record. K is the number of stability records, j is the sampling point, k is the event number, c is the number of the three components, x , y , zThey are three directions of three parts respectively.

[0024] In a preferred embodiment, the covariance matrix analysis result includes the polarization angle and polarization degree recorded by the superposition trace of each detector;

[0025] The method for obtaining the covariance matrix analysis result is:

[0026] Based on n The superposition trace records of the geophones are used to construct the n The covariance matrix corresponding to the detectors is:

[0027] ;

[0028] in, and It starts from the first arrival time and the duration is N 2 The first arrival time is the first arrival time of the superposition channel record of the nth detector, and the time length N2 is the time length of the superposition channel record of the nth detector;

[0029] No. n The polarization angle recorded by the superposition channel of each detector is:

[0030] ;

[0031] No. n The polarization degree recorded by the superposition channel of each detector is:

[0032] ;

[0033] in, For the n The polarization angle recorded by the superposition channel of the detectors is For the n The polarization degree recorded by the superposition channel of each detector; is the first eigenvalue of the covariance matrix, is the second eigenvalue of the covariance matrix, is the eigenvalue The corresponding first eigenvector is is the eigenvalue The corresponding second eigenvector.

[0034] In a preferred embodiment, the method for calculating the calibration angle of each detector relative to each stable source signal is:

[0035] No. n The calibration angle of each detector relative to the reference direction is:

[0036] ;

[0037] in, For the n The detector is relative to the k The calibration angle of a stable source signal, For the k The stable source signal is relative to the n The position of the detector, For the n The polarization angle recorded by the superposition channel of the detectors is n is the detector number, k is the number of the stable source signal.

[0038] In a preferred embodiment, n The correction angle of the horizontal component of each detector is:

[0039] ;

[0040] ;

[0041] in, For the n The correction angle of the horizontal component of a detector is, K is the number of stable source signals, For the k The detector is relative to the m The calibration angle of a stable source signal, For the n The detector is relative to the k The weighting coefficient of a stable source signal, k is the number of the stable signal source, r is a positive number greater than 1, n is the detector number, k is the number of the stable source signal.

[0042] In a second aspect of the present invention, a shallow well seismic detector horizontal component orientation system based on a stable source signal is proposed, the system comprising:

[0043] A bearing calculation module calculates the relative bearing of each interference source signal relative to each detector;

[0044] The stable source screening module is used to obtain representative records of each interference source signal detected by each detector; calculate the waveform matching coefficient of each representative record, screen each representative record based on the waveform matching coefficient, and obtain each stability record detected by each detector after screening; the interference source signal corresponding to each stability record is recorded as a stable source signal;

[0045] A stacking channel stacking module stacks all stability records acquired by each detector to obtain a stacking channel record with enhanced signal-to-noise ratio of each detector;

[0046] The correction angle calculation module is used to construct the covariance matrix of each detector based on the superposition channel record of each detector and perform covariance analysis to obtain the covariance matrix analysis result; calculate the calibration angle of each detector relative to each stable source signal based on the covariance matrix analysis result of each detector and the relative azimuth of each stable source signal relative to each detector; and calculate the correction angle of the horizontal component of each detector based on the calibration angle of each detector relative to each stable source signal.

[0047] The third aspect of the present invention proposes an electronic device, comprising: at least one processor; and a memory communicatively connected to at least one of the processors; wherein the memory stores instructions executable by the processor, and the instructions are used to be executed by the processor to implement the above-mentioned shallow well seismic detector horizontal component orientation method based on stable source signals.

[0048] In a fourth aspect of the present invention, a computer-readable storage medium is proposed, wherein the computer-readable storage medium stores computer instructions, and the computer instructions are used to be executed by the computer to implement the above-mentioned shallow well seismic detector horizontal component orientation method based on stable source signals.

[0049] Beneficial effects of the present invention:

[0050] (1) The present invention improves the signal-to-noise ratio of stable source events and the quality of monitoring data by superimposing interference source signals;

[0051] (2) The present invention improves the quality of monitoring data, improves the accuracy of the correction angle of the detector, and thus improves the detection quality of the detector;

[0052] (3) The present invention can improve the accuracy of the horizontal component orientation of the detector and the azimuth analysis of the earthquake source event by improving the detection quality of the detector. BRIEF DESCRIPTION OF THE DRAWINGS

[0053] Other features, objects and advantages of the present application will become more apparent by reading the detailed description of non-limiting embodiments made with reference to the following drawings:

[0054] Figure 1 It is a schematic diagram of a method for orienting the horizontal component of a shallow well geophone based on a stable source signal according to an embodiment of the present invention;

[0055] Figure 2 is a schematic diagram of a shallow well three-component observation system in an embodiment of the present invention;

[0056] Figure 3 is a schematic diagram of 4 minutes of monitoring data at the 25th observation position in the shallow well three-component observation system in an embodiment of the present invention;

[0057] Figure 4 is a schematic diagram of a representative record at the 25th measurement position of the shallow well three-component observation system of an embodiment of the present invention;

[0058] Figure 5 is a schematic diagram of waveform matching coefficients of a representative record at the 25th measurement position in an embodiment of the present invention;

[0059] Figure 6 is the stability record of the first stable source signal at the 25th position in the embodiment of the present invention ( Figure 6 (a)) and superimposed trace records ( Figure 6 (b) Schematic diagram;

[0060] Figure 7 is the stability record of the second stable source signal at the 25th measurement position in the embodiment of the present invention ( Figure 7 (a)) and superimposed trace records ( Figure 7 (b) Schematic diagram;

[0061] Figure 8 is the polarization analysis result of the superposition trace record of the second stable source signal at the 25th measurement position in the embodiment of the present invention;

[0062] Fig. 9 It is a schematic diagram of the detector azimuth calibration angle and polarization degree obtained based on different stable source signals at different times in the 25th position measurement in the embodiment, and the final azimuth calibration angle obtained;

[0063] Fig.10 It is a structural diagram of a computer system of a server for implementing the method, system, and device embodiments of the present application. DETAILED DESCRIPTION

[0064] The present application will be further described in detail below in conjunction with the accompanying drawings and embodiments. It is to be understood that the specific embodiments described herein are only used to explain the relevant invention, rather than to limit the invention. It is also necessary to explain that, for ease of description, only the parts related to the relevant invention are shown in the accompanying drawings.

[0065] It should be noted that, in the absence of conflict, the embodiments and features in the embodiments of the present application can be combined with each other. The present application will be described in detail below with reference to the accompanying drawings and in combination with the embodiments.

[0066] The present invention provides a method for orienting the horizontal component of a shallow well seismic detector based on a stable source signal, the method comprising:

[0067] Calculate the relative position of each interference source signal with respect to each detector;

[0068] Obtain representative records of each interference source signal detected by each detector; calculate the waveform matching coefficient of each representative record, screen each representative record based on the waveform matching coefficient, and obtain each stability record detected by each detector after screening; the interference source signal corresponding to each stability record is recorded as a stable source signal;

[0069] All stability records acquired by each detector are superimposed to obtain a superimposed channel record with enhanced signal-to-noise ratio of each detector;

[0070] Based on the superposition channel records of each detector, the covariance matrix of each detector is constructed and covariance analysis is performed to obtain the covariance matrix analysis results; based on the covariance matrix analysis results of each detector and the relative azimuth of each stable source signal relative to each detector, the calibration angle of each detector relative to each stable source signal is calculated; based on the calibration angle of each detector relative to each stable source signal, the correction angle of the horizontal component of each detector is calculated.

[0071] In order to more clearly explain the shallow well seismic geophone horizontal component orientation method based on the stable source signal of the present invention, the following is combined with Figure 1 Each step in the embodiment of the present invention is described in detail.

[0072] The steps of the shallow well seismic detector horizontal component orientation method based on the stable source signal of the first embodiment of the present invention are described in detail as follows:

[0073] Calculate the relative position of each interference source signal with respect to each detector;

[0074] In this embodiment, the method for calculating the relative position of each interference source signal relative to each detector is:

[0075] Taking the north direction as the reference direction and clockwise as positive, the direction of the mth interference source signal relative to the nth detector is:

[0076] ;

[0077] in, m is the interference source number, n is the detector number, sx m is the horizontal position coordinate of the interference source signal, sy m is the longitudinal position coordinate of the interference source signal, rx m is the transverse position coordinate of the detector, ry mis the longitudinal position coordinate of the detector, Indicates m The interference source signal is relative to the n The position of the detector, Indicates direction.

[0078] Acquire representative records of each interference source signal detected by each detector; calculate the waveform matching coefficient of each representative record, screen the representative record based on the waveform matching coefficient, and obtain each stability record detected by each detector after screening; the interference source signal corresponding to each stability record is recorded as a stable source signal;

[0079] Methods for obtaining each stability record detected by each detector include:

[0080] First, calculate the waveform matching coefficient of each representative record; the waveform coefficient of the representative record of the mth interference source signal obtained by the nth detector is:

[0081] ;

[0082] Among them, c is the number of the three components, x, y, z are the three directions of the three components respectively, For the n The waveform coefficient of the representative record of the mth interference source signal obtained by the detector, i Represents the sampling point number, j Indicates the sampling point number in the time window, Indicates that the sampling point number in the time window is j The representative record of the cth trisection corresponding to the sampling point in the time window, represents a representative record of the cth component obtained by the nth detector; Nwin represents the time window length, n Number the detector;

[0083] No. n The waveform matching coefficient of the waveform coefficient of the representative record of the mth interference source signal obtained by the detector When it is greater than the threshold E, the current representative record is retained and used as the stability record obtained by the nth detector, and the interference source signal corresponding to each stability record is recorded as the stable source signal.

[0084] All stability records acquired by each detector are superimposed to obtain a superimposed channel record with enhanced signal-to-noise ratio of each detector;

[0085] In this embodiment, the method for obtaining the superimposed trace record with enhanced signal-to-noise ratio of each detector is:

[0086] For the nth detector, the number of stability records detected within a certain period of time is K; each stability record corresponds to a stable source signal;

[0087] The kth stability record detected by the nth detector is recorded as ;

[0088] The signal-to-noise ratio enhanced stack trace record of the nth detector is:

[0089] ;

[0090] in, For the n The signal-to-noise ratio of the detector is enhanced by the stacked channel record. K is the number of stability records, j is the sampling point, k is the event number, c is the three-component number, x , y , z They are three directions of three parts respectively.

[0091] Based on the superposition channel records of each detector, the covariance matrix of each detector is constructed and covariance analysis is performed to obtain the covariance matrix analysis results; based on the covariance matrix analysis results of each detector and the relative azimuth of each stable source signal relative to each detector, the calibration angle of each detector relative to each stable source signal is calculated; based on the calibration angle of each detector relative to each stable source signal, the correction angle of the horizontal component of each detector is calculated.

[0092] In this embodiment, the covariance matrix analysis result includes the polarization angle and polarization degree recorded by the superposition trace of each detector;

[0093] The method for obtaining the covariance matrix analysis result is:

[0094] Based on n The superposition trace records of the geophones are used to construct the n The covariance matrix corresponding to the detectors is:

[0095] ;

[0096] in, and It starts from the first arrival time and the duration is N 2 The first arrival time is the first arrival time of the superposition channel record of the nth detector, and the time length N2 is the time length of the superposition channel record of the nth detector;

[0097] No.n The polarization angle recorded by the superposition channel of each detector is:

[0098] ;

[0099] No. n The polarization degree recorded by the superposition channel of each detector is:

[0100] ;

[0101] in, For the n The polarization angle recorded by the superposition channel of the detectors is For the n The polarization degree recorded by the superposition channel of each detector; is the first eigenvalue of the covariance matrix, is the second eigenvalue of the covariance matrix, is the eigenvalue The corresponding first eigenvector is is the eigenvalue The corresponding second eigenvector.

[0102] In this embodiment, the method for calculating the calibration angle of each detector relative to each stable source signal is:

[0103] No. n The calibration angle of each detector relative to the reference direction is:

[0104] ;

[0105] in, For the n The detector is relative to the k The calibration angle of a stable source signal, For the k The stable source signal is relative to the n The position of the detector, For the n The polarization angle recorded by the superposition channel of the detectors is n is the detector number, k is the number of the stable source signal.

[0106] In this embodiment, n The correction angle of the horizontal component of each detector is:

[0107] ;

[0108]

[0109] in, For then The correction angle of the horizontal component of a detector is, K is the number of stable source signals, For the k The detector is relative to the m The calibration angle of a stable source signal, For the n The detector is relative to the k The weighting coefficient of a stable source signal, k is the number of the stable signal source, r is a positive number greater than 1, n is the detector number, k is the number of the stable source signal.

[0110] In order to better reflect the calibration of the geophone azimuth angle when there is no active source event in the current continuous monitoring record, this embodiment provides an example of a continuous seismic monitoring case in a water flooding front scenario of a water injection development reservoir. Figure 2 This is the observation system composed of 44 shallow well three-component detectors used in the example. In the figure, “○” indicates the target well location, “ ” is the location of the oil pump in the monitoring area, “▽” is the location of the detector, the depth is about 10m, and the sampling rate of the acquisition system is 250Hz. Due to the dense network of injection and production wells in the study area, the seismic monitoring data acquisition cannot avoid the oil wells being produced, resulting in a large number of long-duration background interference signals in the continuous monitoring records. These signals are mainly generated by the up and down reciprocating vibrations of the oil pumps in the study area. The active source event signal excited in this case is weak and is masked by the strong interference background noise, making it impossible for conventional methods to calibrate the azimuth of the horizontal components of each detector. According to the method of the present invention, the relative azimuth and distance can be calculated based on the stable source and the detector position coordinates, and the detector number with the smallest distance to each stable source can be obtained. For example, there is a stable source near the detector numbered 25. Figure 3 for Figure 2 An example of 4 minutes of monitoring data from the 25th observation location in the observation system. It can be clearly seen from the figure that there are periodic stable source signal records in the continuous monitoring record acquisition.

[0111] According to the method of the present invention, the waveform matching coefficient of each representative record is calculated, and the representative records are screened based on the waveform matching coefficient, and the stability records detected by each detector are obtained after screening; the interference source signal corresponding to each stability record is recorded as a stable source signal, and the representative record of the interference signal source is intercepted from the continuous record. Since the stable source signal interval is about 12.95s, Figure 4 As shown in Figure 1, the interception time length of the representative record is 2000 sampling points (8s). The waveform matching coefficient of the representative record is calculated with a threshold of 0.3, as shown in Figure 1. Figure 5The figure shows the waveform matching coefficients of the representative records at the 25th measurement position.

[0112] According to the method of the present invention, all stability records obtained by each detector are superimposed in the morning to obtain a superimposed channel record with enhanced signal-to-noise ratio of each detector. By identifying the signal moment, all stable signal source records can be intercepted from the continuous monitoring record. Figure 6 (a) is the stability record of the stable source signal received by the 25th detector; Figure 6 As shown in (b), by aligning the stable source signal records at the signal recognition time and superimposing the records to form a superimposed track, the signal-to-noise ratio of the signal can be enhanced, and the relative first arrival tp(25) of the stable source signal record superimposed track can be obtained. Similarly, the records and superimposed tracks of other stable signal sources at the measurement position 25 can be obtained, such as Figure 7 As shown in a, it is the stability record of another stable source in detector No. 25. The stable source signal cannot be clearly seen from the figure, such as Figure 7 As shown in b, the signal-to-noise ratio of the stable source signal recording can be significantly enhanced by superposition.

[0113] According to the method of the present invention, the covariance matrix corresponding to each detector is constructed based on the superposition channel record of each detector, the calibration angle of each detector relative to the stable source signal is calculated based on the covariance matrix corresponding to each detector and the relative orientation of each stable source signal relative to each detector, and the correction angle of the horizontal component of each detector is calculated based on the calibration angle of each detector relative to the stable source signal. Polarization analysis is performed on the superposition channel data set of the stable source signal record to obtain the polarization angle and polarization degree of the stable source signal recorded by the detector. Taking the record of the 25th position measurement as an example, Figure 8 for Figure 7 Polarization analysis results of the stable source signal superposition channel record. Considering the relative position relationship between the detector and each stable source in different time and space, combined with the polarization degree of each stable source signal, the calibration angle of the detector horizontal component can be calculated. Fig. 9 The detector azimuth calibration angle and polarization degree obtained for different stable source signals at different times in the 25th position measurement in the embodiment and the final azimuth calibration angle are -64.48°. According to the above method flow, the horizontal component orientation of all three-component detectors in the monitoring area can be achieved, while the conventional method relying on active source events cannot achieve the purpose of detector azimuth calibration.

[0114] Although the various steps in the above embodiment are described in the above-mentioned order, those skilled in the art can understand that in order to achieve the effect of this embodiment, different steps do not have to be executed in such an order. They can be executed simultaneously (in parallel) or in a reverse order. These simple changes are within the scope of protection of the present invention.

[0115] A shallow well geophone horizontal component orientation system based on a stable source signal according to a second embodiment of the present invention comprises:

[0116] A bearing calculation module calculates the relative bearing of each interference source signal relative to each detector;

[0117] A stable source screening module is used to obtain representative records of each interference source signal detected by each detector; calculate the waveform matching coefficient of each representative record, screen the representative record based on the waveform matching coefficient, and obtain each stability record detected by each detector after screening; the interference source signal corresponding to each stability record is recorded as a stable source signal;

[0118] A stacking channel stacking module stacks all stability records acquired by each detector to obtain a stacking channel record with enhanced signal-to-noise ratio of each detector;

[0119] The correction angle calculation module is used to construct the covariance matrix of each detector based on the superposition channel record of each detector and perform covariance analysis to obtain the covariance matrix analysis result; calculate the calibration angle of each detector relative to each stable source signal based on the covariance matrix analysis result of each detector and the relative azimuth of each stable source signal relative to each detector; and calculate the correction angle of the horizontal component of each detector based on the calibration angle of each detector relative to each stable source signal.

[0120] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working process and related instructions of the system described above can refer to the corresponding process in the aforementioned method embodiment, and will not be repeated here.

[0121] It should be noted that the shallow well seismic detector horizontal component directional system based on stable source signal provided in the above embodiment is only illustrated by the division of the above functional modules. In practical applications, the above functions can be assigned to different functional modules as needed, that is, the modules or steps in the embodiments of the present invention can be decomposed or combined. For example, the modules in the above embodiments can be combined into one module, or further divided into multiple sub-modules to complete all or part of the functions described above. The names of the modules and steps involved in the embodiments of the present invention are only for distinguishing the modules or steps, and are not regarded as improper limitations of the present invention.

[0122] An electronic device according to the third embodiment of the present invention comprises: at least one processor; and a memory communicatively connected to at least one of the processors; wherein the memory stores instructions executable by the processor, and the instructions are used to be executed by the processor to implement the above-mentioned shallow well seismic detector horizontal component orientation method based on stable source signals.

[0123] A fourth embodiment of the present invention is a computer-readable storage medium, wherein the computer-readable storage medium stores computer instructions, and the computer instructions are used to be executed by the computer to implement the above-mentioned shallow well seismic detector horizontal component orientation method based on stable source signals.

[0124] Those skilled in the art can clearly understand that, for the convenience and brevity of description, the specific working process and related instructions of the electronic device and storage medium described above can refer to the corresponding process in the aforementioned method embodiment and will not be repeated here.

[0125] Those skilled in the art should be able to appreciate that the modules and method steps of each example described in conjunction with the embodiments disclosed herein can be implemented with electronic hardware, computer software, or a combination of the two, and the programs corresponding to the software modules and method steps can be placed in random access memory (RAM), memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disks, removable disks, CD-ROMs, or any other form of storage medium known in the technical field. In order to clearly illustrate the interchangeability of electronic hardware and software, the composition and steps of each example have been generally described in the above description according to the function. Whether these functions are performed in electronic hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered to exceed the scope of the present invention.

[0126] Reference below Fig.10 , which shows a schematic diagram of the structure of a computer system of a server for implementing the method, system, and device embodiments of the present application. Fig.10 The server shown is merely an example and should not bring any limitation to the functions and scope of use of the embodiments of the present application.

[0127] like Fig.10 As shown, the computer system includes a central processing unit (CPU) 601, which can perform various appropriate actions and processes according to the program stored in the read-only memory (ROM) 602 or the program loaded from the storage part 608 to the random access memory (RAM) 603. Various programs and data required for system operation are also stored in the RAM 603. The CPU 601, ROM 602 and RAM 603 are connected to each other through a bus 604. An input / output (I / O) interface 605 is also connected to the bus 604.

[0128] The following components are connected to the I / O interface 605: an input section 606 including a keyboard, a mouse, etc.; an output section 607 including a cathode ray tube (CRT), a liquid crystal display (LCD), etc., and a speaker; a storage section 608 including a hard disk, etc.; and a communication section 609 including a network interface card such as a LAN (Local Area Network) card, a modem, etc. The communication section 609 performs communication processing via a network such as the Internet. A drive 610 is also connected to the I / O interface 605 as needed. A removable medium 611, such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, etc., is installed on the drive 610 as needed so that a computer program read therefrom is installed into the storage section 608 as needed.

[0129] In particular, according to an embodiment of the present disclosure, the process described above with reference to the flowchart can be implemented as a computer software program. For example, an embodiment of the present disclosure includes a computer program product, which includes a computer program carried on a computer-readable medium, and the computer program includes a program code for executing the method shown in the flowchart. In such an embodiment, the computer program can be downloaded and installed from the network through the communication part 609, and / or installed from the removable medium 611. When the computer program is executed by the central processing unit (CPU) 601, the above-mentioned functions defined in the method of the present application are executed. It should be noted that the above-mentioned computer-readable medium of the present application can be a computer-readable signal medium or a computer-readable storage medium or any combination of the above two. The computer-readable storage medium can be, for example, - but not limited to - an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or device, or any combination of the above. More specific examples of computer-readable storage media may include, but are not limited to, an electrical connection with one or more conductors, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above. In the present application, a computer-readable storage medium may be any tangible medium containing or storing a program that can be used by or in combination with an instruction execution system, an apparatus or a device. In the present application, a computer-readable signal medium may include a data signal propagated in a baseband or as part of a carrier wave, which carries a computer-readable program code. Such propagated data signals may take a variety of forms, including but not limited to electromagnetic signals, optical signals, or any suitable combination of the above. A computer-readable signal medium may also be any computer-readable medium other than a computer-readable storage medium, which may send, propagate, or transmit a program for use by or in combination with an instruction execution system, an apparatus or a device. The program code contained on the computer-readable medium may be transmitted using any appropriate medium, including but not limited to: wireless, wire, optical cable, RF, etc., or any suitable combination of the above.

[0130] Computer program code for performing the operations of the present application may be written in one or more programming languages ​​or a combination thereof, including object-oriented programming languages, such as Java, Smalltalk, C++, and conventional procedural programming languages, such as "C" or similar programming languages. The program code may be executed entirely on the user's computer, partially on the user's computer, as a separate 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 a remote computer, the remote computer may 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 may be connected to an external computer (e.g., via the Internet using an Internet service provider).

[0131] The flow chart and block diagram in the accompanying drawings illustrate the possible architecture, function and operation of the system, method and computer program product according to various embodiments of the present application. In this regard, each square box in the flow chart or block diagram can represent a module, a program segment or a part of a code, and the module, the program segment or a part of the code contains one or more executable instructions for realizing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the square box can also occur in a sequence different from that marked in the accompanying drawings. For example, two square boxes represented in succession can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each square box in the block diagram and / or flow chart, and the combination of the square boxes in the block diagram and / or flow chart can be implemented with a dedicated hardware-based system that performs a specified function or operation, or can be implemented with a combination of dedicated hardware and computer instructions.

[0132] The terms "first", "second", etc. are used to distinguish similar objects rather than to describe or indicate a particular order or sequence.

[0133] The term "comprise" or any other similar term is intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus / device that includes a list of elements includes not only those elements but also other elements not expressly listed, or also includes elements inherent to such process, method, article, or apparatus / device.

[0134] So far, the technical solutions of the present invention have been described in conjunction with the preferred embodiments shown in the accompanying drawings. However, it is easy for those skilled in the art to understand that the protection scope of the present invention is obviously not limited to these specific embodiments. Without departing from the principle of the present invention, those skilled in the art can make equivalent changes or substitutions to the relevant technical features, and the technical solutions after these changes or substitutions will fall within the protection scope of the present invention.

Claims

1. A method for horizontal component orientation of shallow well seismic detectors based on stable source signals, characterized in that: The method comprises: Calculate the relative position of each interference source signal with respect to each detector; Obtain representative records of each interference source signal detected by each detector; calculate the waveform matching coefficient of each representative record, screen each representative record based on the waveform matching coefficient, and obtain each stability record detected by each detector after screening; the interference source signal corresponding to each stability record is recorded as a stable source signal; Methods for obtaining each stability record detected by each detector include: First, calculate the waveform matching coefficient of each representative record; the waveform coefficient of the representative record of the mth interference source signal obtained by the nth detector is: Where c is the number of the three components, x, y, and z are the three directions of the three components, SI(i) is the waveform coefficient of the representative record of the mth interference source signal obtained by the nth detector, i represents the sampling point number, j represents the sampling point number in the time window, and p c (j) represents the representative record of the cth trisection corresponding to the sampling point numbered j in the time window, s c (n, j+i) represents the representative record of the cth component obtained by the nth detector; Nwin represents the time window length, and n is the detector number; When the waveform matching coefficient SI(i) of the waveform coefficient of the representative record of the mth interference source signal obtained by the nth detector is greater than the threshold E, the current representative record is retained and used as the stability record obtained by the nth detector, and the interference source signal corresponding to each stability record is recorded as a stable source signal; All stability records acquired by each detector are superimposed to obtain a superimposed channel record with enhanced signal-to-noise ratio of each detector; Based on the superposition channel records of each detector, the covariance matrix of each detector is constructed and covariance analysis is performed to obtain the covariance matrix analysis results; based on the covariance matrix analysis results of each detector and the relative azimuth of each stable source signal relative to each detector, the calibration angle of each detector relative to each stable source signal is calculated; based on the calibration angle of each detector relative to each stable source signal, the correction angle of the horizontal component of each detector is calculated.

2. The method for shallow well seismic geophone horizontal component orientation based on stable source signal according to claim 1, characterized in that: The method for calculating the relative position of each interference source signal relative to each detector is: Taking the north direction as the reference direction and clockwise as positive, the direction of the mth interference source signal relative to the nth detector is: Where m is the interference source number, n is the detector number, sx m is the horizontal position coordinate of the interference source signal, sy m is the longitudinal position coordinate of the interference source signal, rx n is the transverse position coordinate of the detector, ry n is the longitudinal position coordinate of the detector, α(m, n) represents the azimuth of the mth interference source signal relative to the nth detector, and α represents the azimuth.

3. The method for shallow well seismic detector horizontal component orientation based on stable source signal according to claim 2, characterized in that: The method to obtain the stacked trace record with enhanced signal-to-noise ratio for each detector is: For the nth detector, the number of stability records detected in a certain period of time is K; each stability record corresponds to a stable source signal; The kth stability record detected by the nth detector is denoted by ps c (k, j, n); The signal-to-noise ratio enhanced stack trace record of the nth detector is: Among them, Sps c (j, n) is the superposition channel record of the signal-to-noise ratio enhanced superposition channel record of the nth detector, K is the number of stability records, j is the sampling point, k is the event number, c is the three-component number, x, y, z are the three directions of the three components respectively.

4. The method for shallow well seismic geophone horizontal component orientation based on stable source signal according to claim 3, characterized in that: The covariance matrix analysis result includes the polarization angle and polarization degree recorded by the superposition channel of each detector; The method for obtaining the covariance matrix analysis result is: Construct the covariance matrix corresponding to the nth detector based on the superposition channel record of the nth detector: Among them, x i and i The two horizontal component records start from the first arrival time and have a time length of N2; the first arrival time is the first arrival time of the superposition channel record of the nth detector, and the time length N2 is the time length of the superposition channel record of the nth detector; The polarization angle recorded by the superposition channel of the nth detector is: The polarization degree recorded by the superposition channel of the nth detector is: Where θ(n) is the polarization angle of the superposition channel record of the nth detector, η(n) is the polarization degree of the superposition channel record of the nth detector; λ1 is the first eigenvalue of the covariance matrix, λ2 is the second eigenvalue of the covariance matrix, and u x is the first eigenvector corresponding to the first eigenvalue λ1, u y is the second eigenvector corresponding to the eigenvalue λ1.

5. The method for shallow well seismic geophone horizontal component orientation based on stable source signal according to claim 4, characterized in that: The method for calculating the calibration angle of each detector relative to each stable source signal is: The calibration angle of the nth detector relative to the kth stable source signal is: β(k,n)=α(k,n)-θ(n); Wherein, β(k, n) is the calibration angle of the nth detector relative to the kth stable source signal, α(k, n) is the azimuth of the kth stable source signal relative to the nth detector, θ(n) is the polarization angle recorded by the superposition channel of the nth detector, n is the detector number, and k is the number of the stable source signal.

6. The method for shallow well seismic detector horizontal component orientation based on stable source signal according to claim 5, characterized in that: The correction angle of the horizontal component of the nth detector is: ω(k,n)=η(n)r; Wherein, β′(n) is the correction angle of the horizontal component of the nth detector, K is the number of stable source signals, β(k, n) is the calibration angle of the kth detector relative to the mth stable source signal, ω(k, n) is the weighting coefficient of the nth detector relative to the kth stable source signal, k is the number of the stable signal source, r is a positive number greater than 1, n is the detector number, and k is the number of the stable source signal.

7. A shallow well seismic detector horizontal component directional system based on stable source signal, characterized in that: The system comprises: A bearing calculation module calculates the relative bearing of each interference source signal relative to each detector; The stable source screening module is used to obtain representative records of each interference source signal detected by each detector; calculate the waveform matching coefficient of each representative record, screen each representative record based on the waveform matching coefficient, and obtain each stability record detected by each detector after screening; the interference source signal corresponding to each stability record is recorded as a stable source signal; Methods for obtaining each stability record detected by each detector include: First, calculate the waveform matching coefficient of each representative record; the waveform coefficient of the representative record of the mth interference source signal obtained by the nth detector is: Where c is the number of the three components, x, y, and z are the three directions of the three components, SI(i) is the waveform coefficient of the representative record of the mth interference source signal obtained by the nth detector, i represents the sampling point number, j represents the sampling point number in the time window, and p c (j) represents the representative record of the cth trisection corresponding to the sampling point numbered j in the time window, s c (n, j+i) represents the representative record of the cth component obtained by the nth detector; Nwin represents the time window length, and n is the detector number; When the waveform matching coefficient SI(i) of the waveform coefficient of the representative record of the mth interference source signal obtained by the nth detector is greater than the threshold E, the current representative record is retained and used as the stability record obtained by the nth detector, and the interference source signal corresponding to each stability record is recorded as a stable source signal; A stacking channel stacking module stacks all stability records acquired by each detector to obtain a stacking channel record with enhanced signal-to-noise ratio of each detector; The correction angle calculation module is used to construct the covariance matrix of each detector based on the superposition channel record of each detector and perform covariance analysis to obtain the covariance matrix analysis result; calculate the calibration angle of each detector relative to each stable source signal based on the covariance matrix analysis result of each detector and the relative azimuth of each stable source signal relative to each detector; and calculate the correction angle of the horizontal component of each detector based on the calibration angle of each detector relative to each stable source signal.

8. An electronic device, characterized in that: include: at least one processor; as well as a memory communicatively connected to at least one of the processors; wherein, The memory stores instructions executable by the processor, and the instructions are used to be executed by the processor to implement the shallow well seismic detector horizontal component orientation method based on stable source signals as described in any one of claims 1-6.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores computer instructions, and the computer instructions are used to be executed by the computer to implement the shallow well seismic detector horizontal component orientation method based on stable source signals as described in any one of claims 1-6.

Images