CN111736219A - Method and device for processing multi-component seismic signals - Google Patents
Method and device for processing multi-component seismic signals Download PDFInfo
- Publication number
- CN111736219A CN111736219A CN202010395096.1A CN202010395096A CN111736219A CN 111736219 A CN111736219 A CN 111736219A CN 202010395096 A CN202010395096 A CN 202010395096A CN 111736219 A CN111736219 A CN 111736219A
- Authority
- CN
- China
- Prior art keywords
- component
- data
- component seismic
- nsam
- seismic
- 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.)
- Granted
Links
- 238000012545 processing Methods 0.000 title claims abstract description 51
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000005070 sampling Methods 0.000 claims abstract description 21
- 238000004590 computer program Methods 0.000 claims description 16
- 238000010276 construction Methods 0.000 claims description 6
- 230000003190 augmentative effect Effects 0.000 claims 2
- 238000013523 data management Methods 0.000 abstract description 3
- 238000007726 management method Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 12
- 238000001514 detection method Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 238000003672 processing method Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000003345 natural gas Substances 0.000 description 2
- 239000003209 petroleum derivative Substances 0.000 description 2
- 239000013598 vector Substances 0.000 description 2
- 102100029469 WD repeat and HMG-box DNA-binding protein 1 Human genes 0.000 description 1
- 101710097421 WD repeat and HMG-box DNA-binding protein 1 Proteins 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
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), voltage generated by the three-component detectors in X, Y, Z three directions is subjected to discrete time sampling 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 survey data has three (X, Y, Z three components) or four seismic signals (X, Y, Z, H four components) at a single demodulation point, and is generally implemented by processing each component of the multi-wave multi-component seismic survey data separately in order to allow for rapid analysis and processing of the multi-wave multi-component seismic survey data using conventional seismic data processing software. However, since the particle vibration induced at the wave detection point is projected on X, Y, Z three components of the three-component wave detector when the real seismic wave propagates in the three-dimensional underground medium, one component is processed alone, the wave field of the seismic wave in the underground three-dimensional medium cannot be truly reflected, and the advantages of the multi-wave multi-component seismic exploration data cannot be fully utilized.
In addition, the data transmission of the traditional longitudinal wave seismic exploration data processing software is also in a single-channel form, each processing unit can only receive and output seismic data of one channel, cannot synchronously transmit seismic data of X, Y, Z three components or X, Y, Z, H four 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 has low efficiency and poor effect, thereby seriously hindering the 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.
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;
and the processing module is used for processing the multi-component coding summation signal.
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 description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings 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 raw X, Y, Z three-component seismic signal with sample time in milliseconds (ms) on the abscissa and a sample interval of 1 ms; 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, Z three-component seismic signal 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, and 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 derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within 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 this 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 orthogonal random code sequence may be used to reconstruct the signal of the multi-component code-summed signal, so as to obtain reconstructed multi-component seismic survey data.
1) Case of three-component seismic data:
step 1041: extracting X, Y, Z three-component seismic signals at any one demodulation point from the three-component seismic data, and recording as: { X0j,j＝1,2,……,nsam}、{Y0j,j＝1,2,……,nsam}、{Z0jJ ═ 1,2, … …, nsam }, where X is0jRepresents the value of the X component at the jth sample point; y is0jRepresents the value of the Y component at the jth sample point; z0jRepresents 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:
{Cj,j＝1,2,……,nsam}；
wherein, Cj＝sign(mod(100×|X0j+Y0j+Z0j|/(|X0j|+|Y0j|+|Z0j|),1)-0.5)，CjIs 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 value in parentheses according to the positive and negative values, and takes a value of +1 when the value in parentheses is greater than zero, otherwise takes a value of-1.
Step 1043: constructing an orthogonal code sequence according to the following formula;
H1＝[1,1,1,1]，H2＝[1,-1,1,-1]，H3＝[1,1,-1,-1]，H4＝[1,-1,-1,1]；
wherein, the sequence H1、H2、H3、H4Is composed of +1 and-1 elements and is orthogonal to each other two by two.
H is a 4x4 matrix, H1, H2, H3 and H4 are 4 row vectors of H, and the row vectors are orthogonal pairwise.
Step 1044: constructing an orthogonal random code sequence as follows;
using orthogonal code sequences H1, H2、H3、H4For random code sequence { CjJ is 1,2, … …, nsam to obtain orthogonal random code sequence CH1、CH2、CH3、CH4Is prepared from CH1、CH2、CH3、CH4Is uniformly denoted as { CHkjK is 1,2,3, 4; j ═ 1,2, … …, nsam }, where CHkj＝Cj×Hk，CHkjIs an array containing four elements. For example, when CjWhen is ═ 1, CH1j＝-1×H1＝[-1,-1,-1,-1]，CH2j＝-1×H2＝[-1,1,-1,1]，CH3j＝-1×H3＝[-1,-1,1,1]，CH4j＝-1×H4＝[-1,1,1,-1]。
Will construct the completed orthogonal random code sequence CHkjK is 1,2,3, 4; j is 1,2, … …, nsam, and is stored for standby.
Step 106: from orthogonal random code sequences CH1、CH2、CH3、CH4Sequentially 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}、{Y0j,j＝1,2,……,nsam}、{Z0jJ is one of 1,2, … …, nsam, to obtain encoded three-component seismic encoding data { CXj,j＝1,2,……,nsam}、{CYj,j＝1,2,……,nsam}、{CZjAnd j is 1,2, … …, nsam, and each component seismic coded data after being coded is an array containing four elements.
In particular, from { CHkjK is 1,2,3, 4; selecting CH from 1,2, … …, nsam1、CH2、CH3Respectively with { X0j,j＝1,2,……,nsam}、{Y0j,j＝1,2,……,nsam}、{Z0jJ is multiplied by 1,2, … …, nsam as follows:
CXj＝X0j×CH1jis an array containing four elements;
CYj＝Y0j×CH2jis an array containing four elements;
CZj＝Z0j×CH3jis an array containing four elements.
For example, when X0j＝1234.5678，CH1j＝[-1，-1，-1，-1]When is, CXj＝[-1234.5678,-1234.5678,-1234.5678,-1234.5678](ii) a When Y is0j＝1.56，CH2j＝[-1，1，-1，1]When, CYj＝[-1.56,1.56,-1.56,1.56](ii) a When X is present0j＝-2.6，CH1j＝[-1，-1，1，1]When, CZj＝[2.6,2.6,-2.6,-2.6]。
The three-component encoded signal may also be implemented by:
CXj＝X0j×CH2j；CYj＝Y0j×CH3j；CZj＝Z0j×CH4j。
or:
CXj＝X0j×CH4j；CYj＝Y0j×CH1j；CZj＝Z0j×CH2j。
when all three-component seismic signals in the three-component seismic data are encoded, three orthogonal random code sequences are used, and the use sequence is always consistent.
Each X, Y, Z seismic signal in the three-component seismic data is encoded with an orthogonal random code sequence, resulting in three-component seismic encoded data.
Step 108: component seismic encoding data { CX ] corresponding to each component seismic prospecting data in three-component seismic prospecting dataj,j＝1,2,……,nsam}、{CYj,j＝1,2,……,nsam}、{CZjJ is 1,2, … …, nsam, and the sum operation is performed according to the sampling time sequence to obtain multi-component coding sum data { CXYZ }jJ-1, 2, … …, nsam, wherein CXYZj＝CXj+CYj+CZjIs an array containing four elements.
For example, when CX isj＝[-1234.5678,-1234.5678,-1234.5678,-1234.5678]、CYj＝[-1.56,1.56,-1.56,1.56]、CZj＝[2.6,2.6,-2.6,-2.6]When, there is CXYZj＝[-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 { CXYZjAnd j is 1,2, … …, nsam, and the orthogonal random code sequence { CH } is completed by the above constructionkjK is 1,2,3, 4; j is 1,2, … …, nsam, and the signal is reconstructed to obtain a reconstructed three-component seismic signal { DXj,j＝1,2,……,nsam}，{DYj,j＝1,2,……,nsam}，{DZjJ ═ 1,2, … …, nsam }, where:
CXYZjand CH1jIs an array comprising four elements, CXYZj(1)、CH1j(1) Representing the first value in the array, and the rest being analogized;
CXYZjand CH2jIs an array comprising four elements, CXYZj(1)、CH2j(1) Representing the first number in an arrayValue, the rest is analogized;
CXYZjand CH3jIs an array comprising four elements, CXYZj(1)、CH3j(1) Representing the first value in the array, and the rest being analogized;
and after the operation is completed on all three-component coding summation signals in the three-component seismic coding summation data, three-component seismic reconstruction data are obtained.
2) When the multi-wave multi-component seismic exploration is acquired by a submarine streamer, 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: x, Y, Z, H four-component seismic signals at any one geophone junction are extracted from the four-component seismic data and are recorded as: { X0j,j＝1,2,……,nsam}、{Y0j,j＝1,2,……,nsam}、{Z0j,j＝1,2,……,nsam}、{H0jJ ═ 1,2, … …, nsam }, where X is0jRepresents the value of the X component at the jth sample point; y is0jRepresents the value of the Y component at the jth sample point; z0jRepresents the value of the Z component at the jth sample point; h0jRepresents the value of the H component at the jth sample point; j denotes 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:
{Cj,j＝1,2,……,nsam}；
wherein, Cj＝sign(mod(100×|X0j+Y0j+Z0j+H0j|/(|X0j|+|Y0j|+|Z0j|++|H0j|),1)-0.5)，CjIs 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 value in parentheses according to the positive and negative values, and takes a value of +1 when the value in parentheses is greater than zero, otherwise takes a value of-1.
Step 1043: a method for constructing orthogonal code sequences for three-component seismic signals.
Step 1044: a method for constructing orthogonal random code sequences of three-component seismic signals.
Step 106: using orthogonal random code sequences CH1、CH2、CH3、CH4Multiplying by four-component seismic survey data { X) respectively in sequence0j,j＝1,2,……,nsam}、{Y0j,j＝1,2,……,nsam}、{Z0j,j＝1,2,……,nsam}、{H0jJ is one of 1,2, … …, nsam, to obtain encoded quartering seismic data { CX }j,j＝1,2,……,nsam}、{CYj,j＝1,2,……,nsam}、{CZj,j＝1,2,……,nsam}、{CHjAnd j is 1,2, … …, nsam, wherein each component seismic encoding data after encoding is an array containing four elements.
When four-component seismic signals are coded, four orthogonal random code sequences CH1、CH2、CH3、CH4In 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, H seismic signal in the four-component seismic data by using a corresponding orthogonal random code sequence according to a fixed 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 dataj,j＝1,2,……,nsam}、{CYj,j＝1,2,……,nsam}、{CZj,j＝1,2,……,nsam}、{CHjJ is 1,2, … …, nsam, and the sum operation is carried out according to the sampling time sequence to obtain multi-component coding sum data { CXYZH }jJ-1, 2, … …, nsam, wherein CXYZj＝CXj+CYj+CZj+CHjIs an array containing four elements. After the four-component encoding signals in the four-component seismic encoding data are subjected to the addition operation, the four-component seismic encoding is obtainedAnd (4) adding the data.
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 a code summation signal at the position of a required detection point from four-component seismic code summation data, and recording the code summation signal as { CXYZjAnd j is 1,2, … …, nsam, and the orthogonal random code sequence { CH } is completed by the above constructionkjK is 1,2,3, 4; j is 1,2, … …, nsam, and the four-component seismic signal { DX after reconstruction is obtainedj,j＝1,2,……,nsam}，{DYj,j＝1,2,……,nsam}，{DZjJ ═ 1,2, … …, nsam }, where:
CXYZjand CH1jIs an array comprising four elements, CXYZj(1)、CH1j(1) Representing the first value in the array, and the rest being analogized;
CXYZjand CH2jIs an array comprising four elements, CXYZj(1)、CH2j(1) Representing the first value in the array, and the rest being analogized;
CXYZjand CH3jIs an array comprising four elements, CXYZj(1)、CH3j(1) Representing the first value in the array, and the rest being analogized;
CXYZjand CH4jIs an array comprising four elements, CXYZj(1)、CH4j(1) Representing the first value in the array, and the rest being analogized;
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 raw X, Y, Z three-component seismic signal with sample time in milliseconds (ms) on the abscissa and a sample interval of 1 ms; 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, Z three-component seismic signal 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, and DZ three-component seismic signal. Comparing and analyzing fig. 3, fig. 6, and fig. 7, it can be seen that the error between the reconstructed signal and the original signal is very small, and the larger error is distributed at the beginning and the ending of the recording time, and the time period is far away from the interested middle time period, so 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.
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 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 (15)
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.
2. The method of processing multicomponent seismic signals according to claim 1, wherein constructing a set of pairwise orthogonal random codes from the multicomponent seismic survey data comprises:
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.
3. A method of processing multicomponent seismic signals according to claim 2, wherein when the multicomponent 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:
{Cj,j＝1,2,……,nsam}；
wherein, Cj＝sign(mod(100×|X0j+Y0j+Z0j|/(|X0j|+|Y0j|+|Z0j|),1)-0.5)，CjIs 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; { X0jJ-1, 2, … …, nsam is X-component seismic survey data, X, of the three-component seismic survey data0jRepresents the value of the X component at the jth sample point; { Y0jJ-1, 2, … …, nsam is Y-component seismic survey data, Y, of the three-component seismic survey data0jRepresents the value of the Y component at the jth sample point; { Z0j,j＝1,2,……Nsam is Z component seismic survey data, Z, of three component seismic survey data0jRepresents the value of the Z component at the jth sample point; and nsam is the number of sampling points.
4. A method of processing multicomponent seismic signals according to claim 2, 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:
{Cj,j＝1,2,……,nsam}；
wherein, Cj＝sign(mod(100×|X0j+Y0j+Z0j+H0j|/(|X0j|+|Y0j|+|Z0j|++|H0j|),1)-0.5)，CjIs 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; { X0jJ-1, 2, … …, nsam is X-component seismic survey data, X, of the three-component seismic survey data0jRepresents the value of the X component at the jth sample point; { Y0jJ-1, 2, … …, nsam is Y-component seismic survey data, Y, of the three-component seismic survey data0jRepresents the value of the Y component at the jth sample point; { Z0jJ-1, 2, … …, nsam is Z-component seismic survey data, Z, of the three-component seismic survey data0jRepresents the value of the Z component at the jth sample point; { H0j,j＝1,2,……,nsam}，H0jRepresents the value of the H component at the jth sample point; and nsam is the number of sampling points.
5. A method of processing multicomponent seismic signals according to claim 3 or 4, wherein orthogonal code sequences are constructed as follows;
H1＝[1,1,1,1]，H2＝[1,-1,1,-1]，H3＝[1,1,-1,-1]，H4＝[1,-1,-1,1]；
wherein, the sequence H1、H2、H3、H4Are orthogonal with each other pairwise.
6. The method of processing multicomponent seismic signals of claim 5, 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 H1、H2、H3、H4For random code sequence { CjJ is 1,2, … …, nsam to obtain orthogonal random code sequence CH1、CH2、CH3、CH4Is prepared from CH1、CH2、CH3、CH4Is uniformly denoted as { CHkj,k＝1,2,3,4；j＝1,2,……,nsam}。
7. The method of processing multicomponent seismic signals of claim 6, wherein the random code sequences are separately augmented with orthogonal code sequences as follows:
CHkj＝Cj×Hk；
wherein, CHkjIs an array containing four elements.
8. The method of processing multicomponent seismic signals of claim 6, 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 CH1、CH2、CH3、CH4Sequentially 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}、{Y0j,j＝1,2,……,nsam}、{Z0jJ is one of 1,2, … …, nsam, to obtain codedThree-component seismic encoding data { CXj,j＝1,2,……,nsam}、{CYj,j＝1,2,……,nsam}、{CZjAnd j is 1,2, … …, nsam, and each component seismic coded data after being coded is an array containing four elements.
9. The method of processing multicomponent seismic signals of claim 6, wherein when the multicomponent seismic survey data is quarter-component seismic survey data, each of the component seismic survey data 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:
using orthogonal random code sequences CH1、CH2、CH3、CH4Multiplying by four-component seismic survey data { X) respectively in sequence0j,j＝1,2,……,nsam}、{Y0j,j＝1,2,……,nsam}、{Z0j,j＝1,2,……,nsam}、{H0jJ is one of 1,2, … …, nsam, to obtain encoded quartering seismic data { CX }j,j＝1,2,……,nsam}、{CYj,j＝1,2,……,nsam}、{CZj,j＝1,2,……,nsam}、{CHjAnd j is 1,2, … …, nsam, wherein each component seismic encoding data after encoding is an array containing four elements.
10. The method of processing multicomponent seismic signals of claim 1, wherein processing the multicomponent encoded summed 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.
11. 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 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;
and the processing module is used for processing the multi-component coding summation signal.
12. The apparatus for processing multicomponent seismic signals of claim 11, wherein the orthogonal random code construction 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.
13. The apparatus for processing multicomponent seismic signals of claim 11, 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.
14. 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 10 when executing the computer program.
15. 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 10.
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 true CN111736219A (en) | 2020-10-02 |
| CN111736219B 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 attribute |
| 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 attribute |
| 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 |
|---|
| CHEN JUNFENG等: "Signal Process Method Research about Multi-Component", 《 2016 IEEE/OES CHINA OCEAN ACOUSTICS (COA)》 * |
| 李向阳，王九拴: "多波地震勘探及裂缝储层预测研究进展", 《石油科学通报》 * |
| 王有新等: "地震勘探技术概述", 《油气地球物理》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111736219B (en) | 2023-04-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US20070274155A1 (en) | Coding and Decoding: Seismic Data Modeling, Acquisition and Processing | |
| US11828894B2 (en) | Multi-scale unsupervised seismic velocity inversion method based on autoencoder for observation data | |
| US20100161235A1 (en) | Imaging of multishot seismic data | |
| US9857488B2 (en) | Efficient wavefield compression in seismic imaging | |
| Ikelle | Coding and decoding: Seismic data: The concept of multishooting | |
| US11221426B2 (en) | Method for improved processing of data with time overlapping recordings of energy sources | |
| CN112285768B (en) | High-frequency marine acoustic guided wave frequency dispersion analysis device and method | |
| CN114839673B (en) | Separation method, separation system and computer equipment for multi-seismic-source efficient acquisition wave field | |
| CN111025386B (en) | Vertical and horizontal wave separation method without separation false image | |
| de Ridder | Passive seismic surface-wave interferometry for reservoir-scale imaging | |
| CN109031412B (en) | Elastic wave passive source data primary wave estimation method | |
| Wang et al. | An iterative zero-offset VSP wavefield separating method based on the error analysis of SVD filtering | |
| Nakayama et al. | Machine-learning-based data recovery and its contribution to seismic acquisition: Simultaneous application of deblending, trace reconstruction, and low-frequency extrapolation | |
| Wang et al. | A self-supervised deep learning method for seismic data deblending using a blind-trace network | |
| CN111045077A (en) | Full waveform inversion method of land seismic data | |
| CN111736219B (en) | Method and device for processing multi-component seismic signals | |
| Zhong et al. | Multi-Scale Encoder-Decoder Network for DAS Data Simultaneous Denoising and Reconstruction | |
| CN112882101B (en) | Random noise attenuation method and device for pre-stack seismic data | |
| Saengduean et al. | Multi-source wavefield reconstruction combining interferometry and compressive sensing: application to a linear receiver array | |
| CN107450101B (en) | Seismic data dividing method based on multiple wavelet | |
| CN115453621B (en) | Longitudinal and transverse wave decoupling separation false image removing method based on first-order speed-stress equation | |
| CN114509805B (en) | Vector convex set projection multi-component three-dimensional seismic data reconstruction method and device | |
| Ma et al. | Accurate Reconstruction of Short Duration Passive Seismic Data with Transformer Integrating Multi-Scale Dense Network | |
| Wang et al. | $\mathbf {\emph {DeFault}} $: $\mathbf {\emph {De}} $ ep-learning Based $\mathbf {\emph {Fault}} $ Delineation Using the IBDP Passive Seismic Data at the Decatur CO $ _2 $ Storage Site | |
| Brankovic | Towards an Improved DAS Workflow for Geothermal Resource Development |
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 |