CN111736219B - Method and device for processing multi-component seismic signals - Google Patents
Method and device for processing multi-component seismic signals Download PDFInfo
- Publication number
- CN111736219B CN111736219B CN202010395096.1A CN202010395096A CN111736219B CN 111736219 B CN111736219 B CN 111736219B CN 202010395096 A CN202010395096 A CN 202010395096A CN 111736219 B CN111736219 B CN 111736219B
- Authority
- CN
- China
- Prior art keywords
- component
- data
- component seismic
- nsam
- orthogonal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method and a device for processing multi-component seismic signals, wherein the method comprises the following steps: acquiring multi-component seismic exploration data; constructing a group of pairwise orthogonal random codes according to the multi-component seismic exploration data; respectively encoding each component seismic exploration data of the multi-component seismic exploration data according to the orthogonal random codes to obtain component seismic encoding data corresponding to each component seismic exploration data; summing the component seismic encoding data corresponding to each component seismic exploration data according to a sampling time sequence to obtain multi-component encoding summed data; the multi-component coded summed signal is processed. Because the multi-component seismic data acquired by the scheme are coded and added to obtain new coded and added data, the multi-component seismic data can be regarded as only one component, the data management can be conveniently carried out by using conventional longitudinal wave seismic data processing software, and the management efficiency of the multi-component seismic exploration data is improved.
Description
Technical Field
The invention relates to the technical field of multi-component seismic data processing, in particular to a multi-component seismic signal processing method and device.
Background
Geophones used in conventional longitudinal wave seismic exploration typically record only a single component seismic signal, i.e., at one geophone site, only the seismic signal at that site is vibrating in a direction perpendicular to the earth's surface. The invention and the application of the three-component geophone enable people to have the ability of simultaneously recording three-component vibration signals of ground geophone points in three-dimensional space caused by seismic waves, and thus promote the rapid development of a multi-wave multi-component seismic exploration technology. The three-component detector has three detection devices which are two-by-two vertical, and generally consists of a detection device which is vertical to a horizontal plane (generally called as a Z component) and two detection devices which are in the horizontal plane but vertical to each other (generally called as an X component and a Y component). In the multi-wave multi-component seismic exploration, a plurality of three-component detectors are generally adopted to simultaneously record vibration conditions of seismic waves caused by a plurality of detection point positions in the field (land or sea), discrete time sampling is carried out on voltages generated by the three-component detectors in the X direction, the Y direction and the Z direction according to sampling step lengths at equal intervals, analog signals are converted into digital signals, the digital signals are recorded and stored on a magnetic medium, and then the three-component seismic signals collected in the field are processed and analyzed indoors. In the multi-wave multi-component seismic exploration, a four-component detector is generally adopted for collecting data on the seabed, and the four-component detector is formed by adding a hydrophone (generally called an H component) on the basis of a three-component detector. The multi-wave multi-component seismic exploration is excited by adopting a longitudinal wave source or a transverse wave source, is received by a three-component detector, and detects underground geological information by comprehensively utilizing a plurality of types of seismic wave field information such as longitudinal waves, transverse waves, converted waves and the like, so that the multi-solution of seismic exploration results can be effectively reduced, the exploration precision of underground geological conditions is improved, and the recognition range of underground geological phenomena is widened.
However, conventional seismic data processing software is designed to process conventional longitudinal wave seismic survey data and can only be used to process the case where there is only one trace of seismic signal at the monitoring point. The multi-wave multi-component seismic exploration data has three channels (three components of X, Y and Z) or four channels (four components of X, Y, Z and H) at a wave detection point, and in order to rapidly analyze and process the multi-wave multi-component seismic exploration data by using conventional seismic data processing software, a method of individually processing each component of the multi-component seismic exploration data is generally adopted. However, when real seismic waves propagate in a three-dimensional underground medium, particle vibration induced at a wave detecting point is projected on three components of X, Y and Z of a three-component wave detector respectively, so that one component is processed independently, a wave field of the seismic waves in the underground three-dimensional medium cannot be reflected truly, and the advantages of multi-wave multi-component seismic exploration data cannot be fully utilized.
In addition, data transmission of the traditional longitudinal wave seismic exploration data processing software is also in a single-channel mode, each processing unit can only receive and output seismic data of one channel, cannot synchronously transmit seismic data of X, Y and Z components or X, Y, Z and H components, can only decompose multi-component data into data of a plurality of single components, and independently process each single-component data, and the mode of decomposing the multi-component seismic data into the single components for processing is low in efficiency and poor in effect, and seriously hinders popularization and application of the multi-wave multi-component seismic exploration technology.
Disclosure of Invention
The embodiment of the invention provides a method and a device for processing multi-component seismic signals, which solve the technical problem of low processing efficiency caused by the fact that only each single-component data in multi-component seismic exploration data can be processed independently in the prior art.
The embodiment of the invention provides a method for processing multi-component seismic signals, which comprises the following steps:
acquiring multi-component seismic exploration data;
constructing a group of pairwise orthogonal random codes according to the multi-component seismic exploration data;
respectively encoding each component seismic exploration data of the multi-component seismic exploration data according to the orthogonal random codes to obtain component seismic encoding data corresponding to each component seismic exploration data;
summing the component seismic encoding data corresponding to each component seismic exploration data according to a sampling time sequence to obtain multi-component encoding summed data;
processing the multi-component coded summed signal;
constructing a set of pairwise orthogonal random codes according to the multi-component seismic survey data, comprising:
extracting a plurality of component seismic survey data at one of the waypoints from the multi-component seismic survey data;
constructing a random code sequence from the plurality of component seismic survey data;
constructing an orthogonal code sequence;
and carrying out expansion operation on the random code sequence by using the orthogonal code sequence to obtain the orthogonal random code sequence.
The embodiment of the invention also provides a processing device of the multi-component seismic signal, which comprises the following components:
the multi-component seismic exploration data acquisition module is used for acquiring multi-component seismic exploration data;
the orthogonal random code construction module is used for constructing a group of pairwise orthogonal random codes according to the multi-component seismic exploration data;
the encoding module is used for respectively encoding each component seismic exploration data of the multi-component seismic exploration data according to the orthogonal random codes to obtain component seismic encoding data corresponding to each component seismic exploration data;
the summation module is used for summing the component seismic coded data corresponding to each component seismic exploration data according to the sampling time sequence to obtain multi-component coded summation data;
a processing module for processing the multi-component coded summed signal;
the orthogonal random code constructing module is specifically configured to:
extracting a plurality of component seismic survey data at one of the waypoints from the multi-component seismic survey data;
constructing a random code sequence from the plurality of component seismic survey data;
constructing an orthogonal code sequence;
and carrying out expansion operation on the random code sequence by using the orthogonal code sequence to obtain the orthogonal random code sequence.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the method when executing the computer program.
The embodiment of the invention also provides a computer readable storage medium, and the computer readable storage medium stores a computer program for executing the method.
In the embodiment of the invention, each component of the multi-component seismic signals is respectively encoded by constructing a group of pairwise orthogonal random codes, and then the encoded multi-component seismic signals are added into a new encoded sum signal according to the sampling time sequence, so that only one component can be regarded as the single component, the data management can be conveniently carried out by utilizing conventional longitudinal wave seismic data processing software, the management efficiency of the multi-component seismic exploration data can be improved, the processing period is shortened, and the wide application of the multi-wave multi-component seismic exploration technology is facilitated.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the embodiments or the prior art descriptions will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
Fig. 1 is a flow chart (i) of a method for processing a multi-component seismic signal according to an embodiment of the present invention.
Fig. 2 is a flow chart of an orthogonal random code structure according to an embodiment of the present invention.
FIG. 3 is a graph of the original X, Y, Z three-component seismic signal, with sample time on the abscissa, in milliseconds (ms), and a sample interval of 1ms; the ordinate is the voltage recorded by the detector in microvolts (μ V).
FIG. 4 shows the encoded signals CX, CY, and CZ obtained by encoding the X, Y, and Z three-component seismic signals in FIG. 3, respectively.
Fig. 5 is a signal CXYZ obtained by summing the three-component seismic signals CX, CY, CZ encoded in fig. 4.
Fig. 6 is a block diagram showing that the sum signals in fig. 5 are decoded respectively to obtain reconstructed three-component seismic signals DX, DY, and DZ.
FIG. 7 is a difference between an original X, Y, Z three-component seismic signal and a reconstructed DX, DY, DZ three-component seismic signal.
Fig. 8 is a block diagram of a multi-component seismic signal processing apparatus according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
In an embodiment of the present invention, a method for processing a multi-component seismic signal is provided, as shown in fig. 1, the method including:
step 102: acquiring multi-component seismic exploration data;
step 104: constructing a group of pairwise orthogonal random codes according to the multi-component seismic exploration data;
step 106: respectively encoding each component seismic exploration data of the multi-component seismic exploration data according to the orthogonal random codes to obtain component seismic encoding data corresponding to each component seismic exploration data;
step 108: summing the component seismic encoding data corresponding to each component seismic exploration data according to a sampling time sequence to obtain multi-component encoding summed data;
step 110: processing the multi-component coded summed signal.
In the embodiment of the present invention, as shown in fig. 2, step 104 specifically includes:
step 1041: extracting a plurality of component seismic survey data at one of the waypoints from the multi-component seismic survey data;
step 1042: constructing a random code sequence from the plurality of component seismic survey data;
step 1043: constructing an orthogonal code sequence;
step 1044: and carrying out expansion operation on the random code sequence by using the orthogonal code sequence to obtain the orthogonal random code sequence.
In the embodiment of the present invention, the processing of the multi-component code-summed signal in step 110 may be management, storage, and transmission, or may be analysis and processing of a three-component seismic signal, and at this time, the multi-component code-summed signal may be reconstructed by using an orthogonal random code sequence, so as to obtain reconstructed multi-component seismic survey data.
1) Case of three-component seismic data:
step 1041: extracting the seismic signals of X, Y and Z components at any demodulation point from the three-component seismic data, and recording as follows: { X 0j ,j＝1,2,……,nsam}、{Y 0j ,j＝1,2,……,nsam}、{Z 0j J =1,2, \8230;, nsam }, where X 0j Represents the value of the X component at the jth sample point; y is 0j Represents the value of the Y component at the jth sample point; z 0j Represents the value of the Z component at the jth sample point; and nsam is the number of sampling points. The sampling points are time-sequential sampling points.
Step 1042: a random code sequence is constructed according to the following formula:
{C j ,j＝1,2,……,nsam}；
wherein, C j ＝sign(mod(100×|X 0j +Y 0j +Z 0j |/(|X 0j |+|Y 0j |+|Z 0j |),1)-0.5)，C j Is the value of the jth sampling instant; mod (a, b) is a remainder operation, i.e., a remainder obtained by dividing a by b; sign () takes the value of a value in parentheses according to the signAnd operating, when the value in the brackets is more than zero, taking the value as +1, otherwise, taking the value as-1.
Step 1043: constructing an orthogonal code sequence according to the following formula;
H 1 ＝[1,1,1,1]，H 2 ＝[1,-1,1,-1]，H 3 ＝[1,1,-1,-1]，H 4 ＝[1,-1,-1,1]；
wherein, the sequence H 1 、H 2 、H 3 、H 4 Is composed of +1 and-1 elements and is orthogonal to each other two by two.
H is a 4x4 matrix, H1, H2, H3, H4 are 4 row vectors of H, and the row vectors are orthogonal pairwise.
Step 1044: constructing an orthogonal random code sequence as follows;
respectively using orthogonal code sequences H 1 、H 2 、H 3 、H 4 For random code sequence { C j J =1,2, \8230;, nsam } performs an expansion operation to obtain an orthogonal random code sequence CH 1 、CH 2 、CH 3 、CH 4 Is prepared from CH 1 、CH 2 、CH 3 、CH 4 Is uniformly denoted as { CH kj K =1,2,3,4; j =1,2, \8230;, nsam }, wherein CH kj ＝C j ×H k ，CH kj Is an array containing four elements. For example, when C j When =1, CH 1j ＝-1×H 1 ＝[-1,-1,-1,-1]，CH 2j ＝-1×H 2 ＝[-1,1,-1,1]，CH 3j ＝-1×H 3 ＝[-1,-1,1,1]，CH 4j ＝-1×H 4 ＝[-1,1,1,-1]。
Will construct the completed orthogonal random code sequence CH kj K =1,2,3,4; j =1,2, \8230;, nsam } is kept for later use.
Step 106: from orthogonal random code sequences CH 1 、CH 2 、CH 3 、CH 4 Sequentially selecting three orthogonal random code sequences, and sequentially multiplying the three orthogonal random code sequences by three-component seismic exploration data { X } 0j ,j＝1,2,……,nsam}、{Y 0j ,j＝1,2,……,nsam}、{Z 0j J =1,2, \8230;, nsam } after being codedThree-component seismic encoding data { CX j ,j＝1,2,……,nsam}、{CY j ,j＝1,2,……,nsam}、{CZ j J =1,2, \8230;, nsam }, each component seismic encoded data after encoding is an array containing four elements.
In particular, from { CH kj K =1,2,3,4; j =1,2, \8230;, nsam } selecting CH 1 、CH 2 、CH 3 Respectively with { X 0j ,j＝1,2,……,nsam}、{Y 0j ,j＝1,2,……,nsam}、{Z 0j J =1,2, \8230;, nsam } multiplication as follows:
CX j ＝X 0j ×CH 1j is an array containing four elements;
CY j ＝Y 0j ×CH 2j is an array containing four elements;
CZ j ＝Z 0j ×CH 3j is an array containing four elements.
For example, when X 0j ＝1234.5678，CH 1j ＝[-1，-1，-1，-1]When is, CX j ＝[-1234.5678,-1234.5678,-1234.5678,-1234.5678](ii) a When Y is 0j ＝1.56，CH 2j ＝[-1，1，-1，1]When, CY j ＝[-1.56,1.56,-1.56,1.56](ii) a When X is present 0j ＝-2.6，CH 1j ＝[-1，-1，1，1]When, CZ j ＝[2.6,2.6,-2.6,-2.6]。
The three-component coded signal may also be implemented by the following operations:
CX j ＝X 0j ×CH 2j ；CY j ＝Y 0j ×CH 3j ；CZ j ＝Z 0j ×CH 4j 。
or:
CX j ＝X 0j ×CH 4j ；CY j ＝Y 0j ×CH 1j ；CZ j ＝Z 0j ×CH 2j 。
when all three-component seismic signals in the three-component seismic data are coded, three fixed orthogonal random code sequences are used, and the using sequence is always consistent.
And coding each X, Y and Z seismic signal in the three-component seismic data by an orthogonal random code sequence to finally obtain three-component seismic coding data.
Step 108: component seismic encoding data { CX } corresponding to each component seismic survey data in three-component seismic survey data j ,j＝1,2,……,nsam}、{CY j ,j＝1,2,……,nsam}、{CZ j J =1,2, \8230;, nsam } carries out summation operation according to the sequence of sampling time to obtain multi-component coding summation data { CXYZ } j J =1,2, \8230;, nsam }, wherein CXYZ j ＝CX j +CY j +CZ j Is an array containing four elements.
For example, when CX is j ＝[-1234.5678,-1234.5678,-1234.5678,-1234.5678]、CY j ＝[-1.56,1.56,-1.56,1.56]、CZ j ＝[2.6,2.6,-2.6,-2.6]When, there is CXYZ j ＝[-1233.5278，-1230.4078，-1238.7278，-1235.6078]。
And after the summation operation is completed on all three-component coding signals in the three-component seismic coding data, the three-component seismic coding summation data is obtained.
Step 110: the three-component seismic coding sum data can be managed, stored and transmitted by using traditional longitudinal wave seismic data processing software. When the three-component seismic signals need to be analyzed and processed, the three-component seismic signals can be reconstructed by adding three-component seismic codes and data.
The specific reconstruction method is as follows:
extracting a code summation signal at the position of a required detection point from three-component seismic code summation data, and recording the code summation signal as { CXYZ j J =1,2, \8230;, nsam }, an orthogonal random code sequence { CH } completed with the above construction kj K =1,2,3,4; j =1,2, \8230;, nsam } carries out signal reconstruction to obtain a reconstructed three-component seismic signal { DX j ,j＝1,2,……,nsam}，{DY j ,j＝1,2,……,nsam}，{DZ j J =1,2, \8230;, nsam }, wherein:
and after the three-component code summation signals in the three-component seismic code summation data are subjected to the operation, three-component seismic reconstruction data are obtained.
2) When the multi-wave multi-component seismic exploration is acquired by a submarine towrope, a four-component detector is generally adopted, and the four-component seismic signal is encoded and reconstructed, so that the method and the technology can be popularized and applied.
Step 1041: extracting the seismic signals of X, Y, Z and H components at any demodulation point from the four-component seismic data, and recording as: { X 0j ,j＝1,2,……,nsam}、{Y 0j ,j＝1,2,……,nsam}、{Z 0j ,j＝1,2,……,nsam}、{H 0j J =1,2, \8230;, nsam }, where X 0j Represents the value of the X component at the jth sample point; y is 0j Represents the value of the Y component at the jth sample point; z is a linear or branched member 0j Represents the value of the Z component at the jth sample point; h 0j Represents the value of the H component at the jth sample point; j represents the jth sample point and nsam is the number of sample points.
Step 1042: a random code sequence is constructed according to the following formula:
{C j ,j＝1,2,……,nsam}；
wherein, C j ＝sign(mod(100×|X 0j +Y 0j +Z 0j +H 0j |/(|X 0j |+|Y 0j |+|Z 0j |+|H 0j |),1)-0.5)，C j Is the value of the jth sampling instant; mod (a, b) is a remainder operation, i.e., a remainder obtained by dividing a by b; sign () is an operation of taking values of the numerical values in parentheses according to the positive and negative, and takes a value of +1 when the numerical values in parentheses are greater than zero, otherwise takes a value of-1.
Step 1043: a method for constructing an orthogonal code sequence of a same three-component seismic signal.
Step 1044: a method for constructing orthogonal random code sequences of three-component seismic signals.
Step 106: using orthogonal random code sequences CH 1 、CH 2 、CH 3 、CH 4 Multiplying by four-component seismic survey data { X) respectively in sequence 0j ,j＝1,2,……,nsam}、{Y 0j ,j＝1,2,……,nsam}、{Z 0j ,j＝1,2,……,nsam}、{H 0j J =1,2, \8230;, nsam } resulting in encoded four-component seismic encoded data { CX j ,j＝1,2,……,nsam}、{CY j ,j＝1,2,……,nsam}、{CZ j ,j＝1,2,……,nsam}、{CH j J =1,2, \8230;, nsam }, where each component seismic encoded data after encoding is an array containing four elements.
When four-component seismic signals are coded, four orthogonal random code sequences CH 1 、CH 2 、CH 3 、CH 4 In any order. When all four-component seismic signals in the four-component seismic data are encoded, the four orthogonal random code sequences used should be used in the same order. And carrying out coding operation on each X, Y, Z and H seismic signal in the four-component seismic data according to a fixed sequence by using a corresponding orthogonal random code sequence to finally obtain four-component seismic coding data.
Step 108: component seismic encoding data { CX ] corresponding to each component seismic prospecting data in four component seismic prospecting data j ,j＝1,2,……,nsam}、{CY j ,j＝1,2,……,nsam}、{CZ j ,j＝1,2,……,nsam}、{CH j J =1,2, \8230;, nsam } carries out summation operation according to the sequence of sampling time to obtain multi-component coding summation data { CXYZH } j J =1,2, \8230;, nsam }, wherein CXYZ j ＝CX j +CY j +CZ j +CH j Is an array containing four elements. And after the four-component coded signals in the four-component seismic coded data are subjected to the addition operation, the four-component seismic coded addition data are obtained.
Step 110: the four-component seismic coding sum data can be managed, stored and transmitted by using traditional longitudinal wave seismic data processing software. When the four-component seismic signals need to be analyzed and processed, the four-component seismic signals can be reconstructed by adding the four-component seismic codes and the data.
The specific reconstruction method is as follows:
extracting the code summation signal at the position of the required detection point from the four-component seismic code summation data, and recording the code summation signal as { CXYZ j J =1,2, \8230;, nsam }, an orthogonal random code sequence { CH } completed with the above construction kj K =1,2,3,4; j =1,2, \8230;, nsam } carries out signal reconstruction to obtain a reconstructed four-component seismic signal { DX } j ,j＝1,2,……,nsam}，{DY j ,j＝1,2,……,nsam}，{DZ j J =1,2, \8230;, nsam }, wherein:
and (4) after the four-component coded summation signals in the four-component seismic coding summation data are subjected to the operation, obtaining four-component seismic reconstruction data.
The following description will take three-component seismic signal encoding and reconstruction as an example.
FIG. 3 is a graph of the original X, Y, Z three-component seismic signal, with sample time on the abscissa, in milliseconds (ms), and a sample interval of 1ms; the ordinate is the voltage recorded by the detector in microvolts (μ V). Fig. 4 shows the encoded signals CX, CY, and CZ obtained by encoding the three-component seismic signals X, Y, and Z shown in fig. 3. Fig. 5 is a signal CXYZ obtained by summing the three-component seismic signals CX, CY, CZ encoded in fig. 4. FIG. 6 is a block diagram of the reconstructed three-component seismic signals DX, DY, and DZ obtained by decoding the summed signals of FIG. 5. FIG. 7 is a difference between an original X, Y, Z three-component seismic signal and a reconstructed DX, DY, DZ three-component seismic signal. As can be seen from the comparison analysis of the images in FIGS. 3, 6 and 7, the error between the reconstructed signal and the original signal is extremely small, and the larger error is distributed at the beginning end and the ending end of the recording time, and the time period is far away from the interested middle time period, so that the influence on the multi-wave multi-component seismic exploration result can be ignored.
Based on the same inventive concept, the embodiment of the present invention further provides a processing apparatus for multi-component seismic signals, as described in the following embodiments. Because the principle of the multi-component seismic signal processing device for solving the problem is similar to the multi-component seismic signal processing method, the implementation of the multi-component seismic signal processing device can refer to the implementation of the multi-component seismic signal processing method, and repeated details are not repeated. As used hereinafter, the term "unit" or "module" may be a combination of software and/or hardware that implements a predetermined function. Although the means described in the embodiments below are preferably implemented in software, an implementation in hardware or a combination of software and hardware is also possible and contemplated.
Fig. 8 is a block diagram of a multi-component seismic signal processing apparatus according to an embodiment of the present invention, as shown in fig. 8, including:
a multi-component seismic exploration data acquisition module 02 for acquiring multi-component seismic exploration data;
an orthogonal random code construction module 04, configured to construct a set of pairwise orthogonal random codes according to the multi-component seismic exploration data;
the encoding module 06 is configured to encode each component seismic exploration data of the multi-component seismic exploration data according to the orthogonal random code, and obtain component seismic encoding data corresponding to each component seismic exploration data;
the summing module 08 is used for summing the component seismic coded data corresponding to each component seismic exploration data according to a sampling time sequence to obtain multi-component coded summed data;
a processing module 10, configured to process the multi-component coded summed signal.
In the embodiment of the present invention, the orthogonal random code constructing module 04 is specifically configured to:
extracting a plurality of component seismic survey data at one of the waypoints from the multi-component seismic survey data;
constructing a random code sequence from the plurality of component seismic survey data;
constructing an orthogonal code sequence;
and carrying out expansion operation on the random code sequence by using the orthogonal code sequence to obtain the orthogonal random code sequence.
In the embodiment of the present invention, the processing module 10 is specifically configured to:
and performing signal reconstruction on the multi-component code summation signal by using the orthogonal random code sequence to obtain reconstructed multi-component seismic exploration data.
An embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, where the processor implements the method when executing the computer program.
The embodiment of the invention also provides a computer readable storage medium, and the computer readable storage medium stores a computer program for executing the method.
In conclusion, after the multi-component seismic data are coded, the conventional longitudinal wave seismic data processing software can be used for carrying out data management, storage and transmission on the multi-component seismic data, so that the indoor analysis and processing efficiency of the multi-component seismic exploration data can be improved, the petroleum and natural gas exploration and development efficiency can be improved, and the wide application of the multi-wave multi-component seismic exploration technology can be promoted.
When the multi-component seismic exploration coding sum data is analyzed and processed, the stored orthogonal random codes are used for decoding and reconstructing the multi-component seismic data to obtain reconstructed multi-component seismic signals, the multi-component seismic signals can be synchronously analyzed and processed, the real situation that seismic waves propagate in an underground three-dimensional medium is met, the multi-solution property of exploring underground geological conditions can be reduced, and the risk of petroleum and natural gas exploration is reduced.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and various modifications and changes may be made to the embodiment of the present invention by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (13)
1. A method of processing a multi-component seismic signal, comprising:
acquiring multi-component seismic exploration data;
constructing a group of pairwise orthogonal random codes according to the multi-component seismic exploration data;
respectively encoding each component seismic exploration data of the multi-component seismic exploration data according to the orthogonal random codes to obtain component seismic encoding data corresponding to each component seismic exploration data;
summing the component seismic encoding data corresponding to each component seismic exploration data according to a sampling time sequence to obtain multi-component encoding summed data;
processing the multi-component coded summed signal;
constructing a set of pairwise orthogonal random codes from the multi-component seismic survey data, comprising:
extracting a plurality of component seismic survey data at one of the waypoints from the multi-component seismic survey data;
constructing a random code sequence from the plurality of component seismic survey data;
constructing an orthogonal code sequence;
and carrying out expansion operation on the random code sequence by using the orthogonal code sequence to obtain the orthogonal random code sequence.
2. A method of processing multi-component seismic signals according to claim 1, wherein when the multi-component seismic survey data is three-component seismic survey data, a random code sequence is constructed from the plurality of component seismic survey data as follows:
{C j ,j＝1,2,……,nsam}；
wherein, C j ＝sign(mod(100×|X 0j +Y 0j +Z 0j |/(|X 0j |+|Y 0j |+|Z 0j |),1)-0.5)，C j Is the value of the jth sampling instant; mod (a, b) is a remainder operation, i.e., a remainder obtained by dividing a by b; sign () is an operation of taking a value of a numerical value in brackets according to the positive and negative, wherein the value is +1 when the numerical value in brackets is more than zero, and is-1 otherwise; { X 0j J =1,2, \8230;, nsam } is three-component seismic survey dataSeismic survey data of medium X component, X 0j Represents the value of the X component at the jth sample point; { Y 0j J =1,2, \8230;, nsam } is Y-component seismic survey data, Y, of three-component seismic survey data 0j Represents the value of the Y component at the jth sample point; { Z 0j J =1,2, \8230;, nsam } is Z-component seismic survey data, of three-component seismic survey data 0j Represents the value of the Z component at the jth sample point; and nsam is the number of sampling points.
3. A method of processing multicomponent seismic signals according to claim 1, wherein when the multicomponent seismic survey data is a quarter-component seismic survey data, a random code sequence is constructed from the plurality of component seismic survey data as follows:
{C j ,j＝1,2,……,nsam}；
wherein, C j ＝sign(mod(100×|X 0j +Y 0j +Z 0j +H 0j |/(|X 0j |+|Y 0j |+|Z 0j |+|H 0j |),1)-0.5)，C j Is the value of the jth sampling instant; mod (a, b) is a remainder operation, i.e., a remainder obtained by dividing a by b; sign () is an operation of taking values of the numerical values in brackets according to the positive and negative, when the numerical values in brackets are more than zero, the numerical values are taken as +1, otherwise, the numerical values are taken as-1; { X 0j J =1,2, \8230;, nsam } is X-component seismic survey data, X, of three-component seismic survey data 0j Represents the value of the X component at the jth sample point; { Y 0j J =1,2, \8230;, nsam } is Y-component seismic survey data, Y, of three-component seismic survey data 0j Represents the value of the Y component at the jth sample point; { Z 0j J =1,2, \8230;, nsam } is Z-component seismic survey data, of three-component seismic survey data 0j Represents the value of the Z component at the jth sample point; { H 0j ,j＝1,2,……,nsam}，H 0j Represents the value of the H component at the jth sample point; and nsam is the number of sampling points.
4. A method of processing multicomponent seismic signals as claimed in claim 2 or 3, characterized by constructing orthogonal code sequences as follows;
H 1 ＝[1,1,1,1]，H 2 ＝[1,-1,1,-1]，H 3 ＝[1,1,-1,-1]，H 4 ＝[1,-1,-1,1]；
wherein, the sequence H 1 、H 2 、H 3 、H 4 Are orthogonal with each other pairwise.
5. The method of processing multicomponent seismic signals of claim 4, wherein the random code sequence is augmented with an orthogonal code sequence in a manner to obtain an orthogonal random code sequence;
respectively using orthogonal code sequences H 1 、H 2 、H 3 、H 4 For random code sequence { C j J =1,2, \8230;, nsam } performs an expansion operation to obtain an orthogonal random code sequence CH 1 、CH 2 、CH 3 、CH 4 Is prepared from CH 1 、CH 2 、CH 3 、CH 4 Is uniformly denoted as { CH kj ,k＝1,2,3,4；j＝1,2,……,nsam}。
6. The method of processing multicomponent seismic signals of claim 5, wherein the random code sequences are separately augmented with orthogonal code sequences as follows:
CH kj ＝C j ×H k ；
wherein, CH kj Is an array containing four elements.
7. The method of processing multicomponent seismic signals of claim 5, wherein when the multicomponent seismic survey data is three-component seismic survey data, each of the multicomponent seismic survey data is separately encoded according to the orthogonal random code in a manner to obtain component seismic encoded data corresponding to each of the component seismic survey data as follows:
from orthogonal random code sequences CH 1 、CH 2 、CH 3 、CH 4 Sequentially selecting three orthogonal random code sequences, and sequentially multiplying the three orthogonal random code sequences by three-component seismic exploration data { X } 0j ,j＝1,2,……,nsam}、{Y 0j ,j＝1,2,……,nsam}、{Z 0j J =1,2, \8230;, nsam } to obtain encoded three-component seismic encoding data { CX j ,j＝1,2,……,nsam}、{CY j ,j＝1,2,……,nsam}、{CZ j J =1,2, \8230;, nsam }, each component seismic encoded data after encoding is an array containing four elements.
8. The method of processing multicomponent seismic signals of claim 5, wherein when the multicomponent seismic survey data is a quarter-component seismic survey data, each of the multicomponent seismic survey data is separately encoded according to the orthogonal random code to obtain component seismic encoded data corresponding to each of the component seismic survey data as follows:
using orthogonal random code sequences CH 1 、CH 2 、CH 3 、CH 4 Multiplying the four-component seismic exploration data by the four-component seismic exploration data respectively in sequence { X } 0j ,j＝1,2,……,nsam}、{Y 0j ,j＝1,2,……,nsam}、{Z 0j ,j＝1,2,……,nsam}、{H 0j J =1,2, \8230;, nsam } resulting in encoded four-component seismic encoded data { CX j ,j＝1,2,……,nsam}、{CY j ,j＝1,2,……,nsam}、{CZ j ,j＝1,2,……,nsam}、{CH j J =1,2, \8230;, nsam }, where each component of encoded seismic encoded data is an array containing four elements.
9. The method of processing multicomponent seismic signals of claim 1, wherein processing the multicomponent coded sum signal comprises:
and performing signal reconstruction on the multi-component code summation signal by using the orthogonal random code sequence to obtain reconstructed multi-component seismic exploration data.
10. Apparatus for processing multicomponent seismic signals, comprising:
the multi-component seismic exploration data acquisition module is used for acquiring multi-component seismic exploration data;
the orthogonal random code construction module is used for constructing a group of pairwise orthogonal random codes according to the multi-component seismic exploration data;
the encoding module is used for respectively encoding each component seismic exploration data of the multi-component seismic exploration data according to the orthogonal random codes to obtain component seismic encoding data corresponding to each component seismic exploration data;
the summing module is used for summing component seismic encoding data corresponding to each component seismic exploration data according to a sampling time sequence to obtain multi-component encoding summed data;
a processing module for processing the multi-component coded summed signal;
the orthogonal random code constructing module is specifically configured to:
extracting a plurality of component seismic survey data at one of the waypoints from the multi-component seismic survey data;
constructing a random code sequence from the plurality of component seismic survey data;
constructing an orthogonal code sequence;
and carrying out expansion operation on the random code sequence by using the orthogonal code sequence to obtain the orthogonal random code sequence.
11. The apparatus for processing multicomponent seismic signals of claim 10, wherein the processing module is specifically configured to:
and performing signal reconstruction on the multi-component code summation signal by using the orthogonal random code sequence to obtain reconstructed multi-component seismic exploration data.
12. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the method of any of claims 1 to 9 when executing the computer program.
13. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 9.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010395096.1A CN111736219B (en) | 2020-05-12 | 2020-05-12 | Method and device for processing multi-component seismic signals |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010395096.1A CN111736219B (en) | 2020-05-12 | 2020-05-12 | Method and device for processing multi-component seismic signals |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111736219A CN111736219A (en) | 2020-10-02 |
| CN111736219B true CN111736219B (en) | 2023-04-07 |
Family
ID=72647107
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010395096.1A Active CN111736219B (en) | 2020-05-12 | 2020-05-12 | Method and device for processing multi-component seismic signals |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111736219B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1575424A (en) * | 2001-10-26 | 2005-02-02 | 维斯特恩格科地震控股有限公司 | A method of and an apparatus for processing seismic data |
| CN102257403A (en) * | 2008-08-17 | 2011-11-23 | 格库技术有限公司 | Estimating and correcting perturbations on seismic particle motion sensors employing seismic source signals |
| CN102854526A (en) * | 2011-06-30 | 2013-01-02 | 中国石油化工股份有限公司 | Multi-component seismic data processing method |
| CN104520733A (en) * | 2012-06-25 | 2015-04-15 | 普拉德研究及开发股份有限公司 | Seismic Orthogonal Decomposition Properties |
| CN105403912A (en) * | 2014-08-25 | 2016-03-16 | 中国石油化工股份有限公司 | Device and system for collecting seismic data |
| WO2018087648A2 (en) * | 2016-11-10 | 2018-05-17 | Seismic Apparition Gmbh | Seismic acquisition and processing method |
-
2020
- 2020-05-12 CN CN202010395096.1A patent/CN111736219B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1575424A (en) * | 2001-10-26 | 2005-02-02 | 维斯特恩格科地震控股有限公司 | A method of and an apparatus for processing seismic data |
| CN102257403A (en) * | 2008-08-17 | 2011-11-23 | 格库技术有限公司 | Estimating and correcting perturbations on seismic particle motion sensors employing seismic source signals |
| CN102854526A (en) * | 2011-06-30 | 2013-01-02 | 中国石油化工股份有限公司 | Multi-component seismic data processing method |
| CN104520733A (en) * | 2012-06-25 | 2015-04-15 | 普拉德研究及开发股份有限公司 | Seismic Orthogonal Decomposition Properties |
| CN105403912A (en) * | 2014-08-25 | 2016-03-16 | 中国石油化工股份有限公司 | Device and system for collecting seismic data |
| WO2018087648A2 (en) * | 2016-11-10 | 2018-05-17 | Seismic Apparition Gmbh | Seismic acquisition and processing method |
Non-Patent Citations (3)
| Title |
|---|
| Signal Process Method Research about Multi-Component;CHEN Junfeng等;《 2016 IEEE/OES China Ocean Acoustics (COA)》;全文 * |
| 地震勘探技术概述;王有新等;《油气地球物理》(第01期);第1-9、15页 * |
| 多波地震勘探及裂缝储层预测研究进展;李向阳，王九拴;《石油科学通报》;第1卷(第1期);第45-60页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111736219A (en) | 2020-10-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109164483B (en) | Multi-component seismic data vector denoising method and multi-component seismic data vector denoising device | |
| US11835672B2 (en) | Non-uniform optimal survey design principles | |
| US9857488B2 (en) | Efficient wavefield compression in seismic imaging | |
| US20070274155A1 (en) | Coding and Decoding: Seismic Data Modeling, Acquisition and Processing | |
| US11221426B2 (en) | Method for improved processing of data with time overlapping recordings of energy sources | |
| CN101487899B (en) | Method for wavefield separation in 3D dual sensor towed streamer data | |
| US20100054082A1 (en) | Reverse-time depth migration with reduced memory requirements | |
| WO2008112036A1 (en) | Imaging of multishot seismic data | |
| CN112285768B (en) | High-frequency marine acoustic guided wave frequency dispersion analysis device and method | |
| CN113885077B (en) | A multi-source seismic data separation method based on deep learning | |
| RU2737846C2 (en) | System for installing ground-based seismic sensors with pairs of adjacent multicomponent seismic sensors at an average distance of at least twenty meters | |
| Wang et al. | A self-supervised deep learning method for seismic data deblending using a blind-trace network | |
| CA1210493A (en) | Seismic exploration system including analog-to- digital converter using delta modulation | |
| Liu et al. | A multitone model-based seismic data compression | |
| Abedi et al. | A multidirectional deep neural network for self-supervised reconstruction of seismic data | |
| Sopher | Converting scanned images of seismic reflection data into SEG-Y format | |
| CN111736219B (en) | Method and device for processing multi-component seismic signals | |
| Han et al. | Seismic event and phase detection using deep learning for the 2016 Gyeongju earthquake sequence | |
| EP2487511A1 (en) | Seismic survey method of the subsoil | |
| CN107450101B (en) | Seismic data dividing method based on multiple wavelet | |
| Saengduean et al. | Multi-source wavefield reconstruction combining interferometry and compressive sensing: application to a linear receiver array | |
| Wang et al. | A physics-augmented deep learning method for seismic data deblending | |
| Hernandez-Rojas et al. | Design of undersampled seismic acquisition geometries via end-to-end optimization | |
| EA003422B1 (en) | Method and system for processing three-component seismic data | |
| CN112882101B (en) | Random noise attenuation method and device for pre-stack seismic data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |