CN113933901A - Depth domain seismic data splicing method and device - Google Patents
Depth domain seismic data splicing method and device Download PDFInfo
- Publication number
- CN113933901A CN113933901A CN202010602759.2A CN202010602759A CN113933901A CN 113933901 A CN113933901 A CN 113933901A CN 202010602759 A CN202010602759 A CN 202010602759A CN 113933901 A CN113933901 A CN 113933901A
- Authority
- CN
- China
- Prior art keywords
- depth
- block
- data
- time
- speed
- 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
- 238000000034 method Methods 0.000 title claims abstract description 47
- 238000006243 chemical reaction Methods 0.000 claims abstract description 69
- 238000013508 migration Methods 0.000 claims abstract description 39
- 230000005012 migration Effects 0.000 claims abstract description 39
- 238000012937 correction Methods 0.000 claims description 44
- 238000012545 processing Methods 0.000 claims description 32
- 230000003068 static effect Effects 0.000 claims description 25
- 230000006870 function Effects 0.000 claims description 22
- 238000004590 computer program Methods 0.000 claims description 13
- 238000009499 grossing Methods 0.000 claims description 6
- 230000007704 transition Effects 0.000 claims description 5
- 238000001914 filtration Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 16
- 238000010276 construction Methods 0.000 description 8
- 238000003384 imaging method Methods 0.000 description 8
- 238000004891 communication Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 6
- 239000000872 buffer Substances 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000008021 deposition Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000002349 favourable effect Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012958 reprocessing Methods 0.000 description 2
- 238000000638 solvent extraction Methods 0.000 description 2
- 230000001413 cellular effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid 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/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
- G01V2210/512—Pre-stack
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 discloses a depth domain seismic data splicing method and device, wherein the method comprises the following steps: obtaining depth domain seismic data and corresponding well logging data of a plurality of blocks, wherein the depth domain seismic data comprises: a depth migration velocity, the well log data comprising: logging speed; converting the depth domain seismic data of the block to which the seismic data belong into time domain data based on the depth migration velocity of each block, and respectively correcting the time domain data of each block; splicing the time domain data of each corrected block based on a preset splicing rule; and obtaining a time-depth conversion speed body according to the depth deviation speed and the logging speed of each block, and performing time-depth conversion operation on the spliced time domain data according to the time-depth conversion speed body to obtain the spliced depth domain data. The depth domain seismic data splicing scheme is low in cost, short in period and easy to implement.
Description
Technical Field
The invention relates to the field of data processing, in particular to a depth domain seismic data splicing method and device.
Background
With the continuous improvement of oil and gas exploration degree, the three-dimensional seismic data of a plurality of oil fields basically cover a main power oil production zone, but due to the limitation of the cost of the oil fields, the acquisition of the seismic data is carried out in a partitioning and partitioning mode, so that the data are obviously inconsistent in imaging, energy, phase and frequency due to the difference of construction and processing technologies. The early seismic data splicing is realized by adopting a time shifting method, and the method has large workload and low efficiency when three-dimensional seismic data are connected and cannot be realized on depth domain data. In 2009, plum relay light and the like adopt a cross-correlation time difference quantitative identification method in a victory oil field to greatly improve the efficiency of a time shift method, but still can be realized only in a time domain. The current mainstream prestack depth migration slice connection method is large-area and multi-block prestack slice connection reprocessing, the method can solve the structure, frequency and phase difference among different blocks, but with the increasing of the three-dimensional earthquake area, the prestack reprocessing cost is higher and higher, and the period is longer and longer, so that a method which is low in cost and short in period and can realize depth domain earthquake data splicing is needed to meet the requirement of regional structure interpretation.
The mature interpretation technology is based on time domain data, so the research on the problem of time domain data splicing is more, and the processing technology of energy, frequency and phase is mature. The prestack depth migration data has obvious advantages on accurate imaging of complex structures and transverse heterogeneous stratums and plays an important role in seismic exploration, but due to the fact that the influence of speed ambiguity is larger than time domain data due to the difference of different blocks, splicing difficulty of depth domain data is undoubtedly larger than that of time domain data, and at present, no mature method is available for solving the problem of splicing of depth domain data.
Disclosure of Invention
In view of the above, the present invention provides a depth domain seismic data stitching method and apparatus, so as to solve at least one of the above-mentioned problems.
According to a first aspect of the invention, there is provided a method of depth domain seismic data stitching, the method comprising: obtaining depth domain seismic data and corresponding well log data for a plurality of blocks, the depth domain seismic data comprising: a depth migration velocity, the well log data comprising: logging speed; converting the depth domain seismic data of the block to which the seismic data belong into time domain data based on the depth migration velocity of each block, and respectively correcting the time domain data of each block; splicing the time domain data of each corrected block based on a preset splicing rule; and obtaining a time-depth conversion speed body according to the depth deviation speed and the logging speed of each block, and performing time-depth conversion operation on the spliced time domain data according to the time-depth conversion speed body to obtain the spliced depth domain data.
According to a second aspect of the invention, there is provided a depth domain seismic data stitching apparatus, the apparatus comprising: a data acquisition unit for acquiring depth domain seismic data and corresponding well log data for a plurality of blocks, the depth domain seismic data comprising: a depth migration velocity, the well log data comprising: logging speed; the time domain conversion unit is used for converting the depth domain seismic data of the block to which the seismic data belong into time domain data based on the depth migration velocity of each block; a correction unit for performing correction processing on the time domain data of each block; the splicing unit is used for splicing the time domain data of each corrected block based on a preset splicing rule; the time-depth conversion speed body determining unit is used for obtaining a time-depth conversion speed body according to the depth deviation speed and the logging speed of each block; and the time-depth conversion unit is used for performing time-depth conversion operation on the spliced time domain data according to the time-depth conversion speed body to obtain the spliced depth domain data.
According to a third aspect of the present invention, there is provided an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the steps of the method when executing the program.
According to a fourth aspect of the invention, the invention provides a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, carries out the steps of the above-mentioned method.
According to the technical scheme, the depth domain seismic data of the block to which the seismic data belong are converted into time domain data based on the depth migration velocity of each block, the time domain data of each block are corrected respectively, the time domain data of each block after correction are spliced, a time-depth conversion velocity body is obtained according to the depth migration velocity and the well logging velocity of each block, and time-depth conversion operation is performed on the spliced time domain data according to the time-depth conversion velocity body, so that the spliced depth domain data can be obtained and reprocessed relative to a tandem stack.
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 introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a flow diagram of a method of depth domain seismic data stitching according to an embodiment of the invention;
FIG. 2 is a detailed flow chart of a depth domain seismic data stitching method according to an embodiment of the invention
FIG. 3 is a schematic cross-sectional view of two exemplary tiles being directly tiled;
FIG. 4 is a cross-sectional schematic diagram of two example tile splices according to embodiments of the invention;
FIG. 5 is a diagram of new seismic velocities for time-depth transitions in accordance with an embodiment of the present invention;
FIG. 6 is a block diagram of a depth domain seismic data stitching device according to an embodiment of the invention;
fig. 7 is a schematic block diagram of a system configuration of an electronic apparatus 600 according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, 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 some, but not all, embodiments of the present invention. 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 the prior art, the processing period of the area pre-stack continuous piece is long, the cost is high, and many data details are covered due to the need of balancing the differences (signal-to-noise ratio, azimuth angle and covering times) of different block data, so that compared with the processing effect of a single piece of data, the accuracy of the continuous pre-stack processed data is reduced, and part of geological information is lost. Based on the scheme, the scheme is low in cost and short in period, and the spliced data basically realizes the unification of the structure depth, the phase and the frequency of each block, so that the method has strong practicability.
In actual operation, the depth domain data of different blocks not only have differences of energy, frequency and phase, but also often have differences of structural forms, and the differences of the structural depths at different depths are also large, so that it is not practical to directly splice the depth domain data.
Because the image velocity ambiguity of prestack depth migration is stronger than prestack time migration, and the image depth is more influenced by the seismic interval velocity, the embodiment of the invention provides that the seismic data in the depth domain is firstly converted into the time domain for splicing, so that the influence of velocity ambiguity is eliminated, and then the spliced time domain is converted back into the depth domain, so that the splicing connection of the seismic data in the depth domain is realized.
Embodiments of the present invention are similarly described below with reference to the accompanying drawings.
FIG. 1 is a flow chart of a method of depth domain seismic data stitching according to an embodiment of the invention, as shown in FIG. 1, the method comprising:
And 102, converting the depth domain seismic data of the block to which the seismic data belong into time domain data based on the depth migration velocity of each block, and respectively correcting the time domain data of each block.
In one embodiment, the performing the correction process on the time domain data of each block comprises: respectively carrying out unified operation of frequency and phase on the time domain data of each block; and performing static correction operation on the repeated areas among the blocks.
The performing static correction operation on the repeated areas among the blocks specifically comprises: firstly, selecting one block from all blocks as a basic block, and determining the time difference between the other blocks and the basic block by picking up the extreme value of cross correlation; and then, performing smoothing operation on the time difference determined by each block, and performing correction operation on the block in a static correction mode.
In practical operation, when depth domain seismic data is converted into time domain seismic data, differences of different data are concentrated on differences of shallow processing and differences of frequency and phase. The processing technology of frequency and phase is mature, and the frequency and phase of different blocks can be unified by adopting the existing processing technology. In shallow processing, the most common method is static correction, and the data difference caused by the static correction is not changed in the time direction, so that the difference values caused by static correction of different blocks can be obtained by adopting a cross-correlation method, and the difference caused by the static correction can be basically solved by taking the difference values as the residual static correction value.
And 103, splicing the time domain data of each corrected block based on a preset splicing rule.
Specifically, the repeated areas between the blocks adopt a weighted overlap mode, and the boundaries of the blocks are subjected to predetermined filtering processing.
And 104, obtaining a time-depth conversion speed body according to the depth deviation speed and the logging speed of each block, and performing time-depth conversion operation on the spliced time domain data according to the time-depth conversion speed body to obtain the spliced depth domain data.
Specifically, control layer pickup operation may be performed on the spliced time domain data, and the initial speed may be filled in the control layer according to the depth offset speed of each block to obtain an initial speed body; and then, taking the logging speed as a target function to carry out inversion operation on the initial speed body to obtain a time-depth conversion speed body corresponding to the logging speed.
In practical operation, after time domain data splicing is completed, the biggest problem is that a uniform seismic velocity is needed to be corrected to a depth domain, the prestack depth migration velocity of each block is often different greatly and cannot be used directly, so that the low frequency quantity of the prestack depth migration velocity of each block can be used as an initial velocity, then a new connected seismic velocity is obtained through well velocity inversion, and the spliced time domain data is converted back to the depth domain by using the new velocity, so that the spliced connected seismic data of the depth domain is realized.
Because the time-depth conversion speed body is a new speed obtained by inverting the logging speed as a target function, the corresponding condition of the spliced data and the well data is better, and the imaging depth is more accurate.
According to the depth domain seismic data splicing method and device, the depth domain seismic data of the block are converted into the time domain data based on the depth migration velocity of the block, the time domain data of the block are corrected respectively, then the time domain data of the block after correction are spliced, a time-depth conversion velocity body is obtained according to the depth migration velocity and the logging velocity of the block, and time-depth conversion operation is conducted on the spliced time domain data according to the time-depth conversion velocity body, so that the spliced depth domain data can be obtained, and the splicing scheme is reprocessed relative to a prestack piece.
FIG. 2 is a detailed flow chart of a depth domain seismic data stitching method according to an embodiment of the invention, as shown in FIG. 2, the flow chart includes:
and step 211, performing time-depth conversion on the spliced time domain large area data according to the time-depth conversion speed body to obtain the final spliced depth domain large area data.
For a better understanding of the invention, an example is given below.
In this example, taking two blocks as an example, which are located in the western oil field of china, block 1 and block 2 have a repetition area of about one kilometer in the east-west direction. Fig. 3 is a cross-sectional view of directly splicing the block 1 and the block 2, and as shown in fig. 3, there is a large difference between the structure depth, the frequency, and the phase of the two repeated regions, including the energy intensity. The early application effect proves that the imaging of the work area 1 is more real, but the range is too small, so that the depth data of the work area 2 is spliced to the work area 1 by taking the work area 1 as a reference.
According to the embodiment of the invention, fig. 4 is a schematic cross-sectional view of two block data after splicing, splicing traces are basically not seen on the spliced cross-sectional view, and the time-depth conversion speed used in the embodiment of the invention is a new speed obtained by taking the logging speed as an objective function for inversion, specifically referring to the schematic new speed view shown in fig. 5, so that the spliced data and the well data have better corresponding conditions, and the imaging depth is more accurate. The specific implementation steps are as follows:
1. expanding the grid of the work area 1 to include all data of the work areas 1 and 2, loading depth domain achievements and speeds of the work areas 1 and 2, and loading well data in the two work areas;
2. for the work area 1, converting the prestack depth migration data volume into a time domain by adopting the prestack depth migration speed, and performing the same operation on the work area 2;
3. correcting the working area 2 as a reference surface to make the working area consistent with the working area 1;
4. carrying out frequency and phase processing on the work area 2 to make the work area 2 consistent with the work area 1;
5. performing cross correlation on the repeated areas of the work area 1 and the work area 2, and picking up a cross correlation extreme value;
6. performing spatial smoothing on the acquired cross-correlation extreme value, and applying the cross-correlation extreme value as a static correction value to the data of the work area 2 to enable the construction depth of the data to be consistent with that of the work area 1;
7. splicing the data of the work area 1 and the processed data of the work area 2, wherein the repeated area adopts weighted superposition, and the boundary gives a reasonable transition zone;
8. picking up a control layer on the spliced time domain data, wherein the density of the control layer can reflect the construction trend of a large area;
9. filling the control layers with initial speeds, wherein the initial speed of each layer is the low frequency of the depth migration speed of the two work areas, and processing the repeated areas according to the splicing mode of the step 7;
10. taking the well speed as a target function to invert the initial velocity body to obtain a time-depth conversion velocity body which is better matched with the well speed, and particularly referring to fig. 5;
11. and performing time-depth conversion on the spliced time domain large-area data by using the velocity volume obtained by inversion to obtain the final spliced depth domain large-area data, which is shown in fig. 4 specifically.
The time-depth conversion speed body is a new speed obtained by inverting the logging speed as a target function, so that the corresponding condition of spliced data and well data is better, and the imaging depth is more accurate.
According to the embodiment of the invention, the spliced data realizes the consistency of the construction depth, the frequency and the phase, avoids the one-sided understanding of scattered block data, provides indispensable favorable conditions for the overall promotion and summary of oil and gas geological understanding of regional construction, deposition, reservoir, formation and the like from the region, greatly reduces the multi-solution of geological understanding, and provides a powerful basis for oil and gas exploration and development decisions.
Based on similar inventive concepts, the embodiment of the invention also provides a depth domain seismic data splicing device, which is preferably used for realizing the flow in the method embodiment.
FIG. 6 is a block diagram of the structure of the depth domain seismic data stitching apparatus, as shown in FIG. 6, which includes: a data acquisition unit 61, a time domain conversion unit 62, a correction unit 63, a stitching unit 64, a time-depth conversion rate body determination unit 65, and a time-depth conversion unit 66, wherein:
a data obtaining unit 61, configured to obtain depth domain seismic data of a plurality of blocks and corresponding well logging data, where the depth domain seismic data includes: a depth migration velocity, the well log data comprising: logging speed;
a time domain conversion unit 62 for converting the depth domain seismic data of the belonging block into time domain data based on the depth migration velocity of each block;
a correction unit 63, configured to perform correction processing on the time domain data of each block;
and a splicing unit 64, configured to perform splicing processing on the time domain data of each corrected block based on a predetermined splicing rule.
Specifically, the stitching unit 64 performs stitching processing on the time domain data of each corrected block, wherein the repeated areas between the blocks adopt a weighted overlap mode, and each block boundary performs predetermined filtering processing.
A time-depth conversion velocity body determining unit 65, configured to obtain a time-depth conversion velocity body according to the depth deviation velocity and the logging velocity of each block;
and a time-depth conversion unit 66, configured to perform time-depth conversion operation on the spliced time domain data according to the time-depth conversion speed body, so as to obtain spliced depth domain data.
According to the depth domain seismic data splicing method and device, the time domain conversion unit 62 converts the depth domain seismic data of the block to which the time domain seismic data belong into the time domain data based on the depth migration velocity of each block, the correction unit 63 corrects the time domain data of each block respectively, the splicing unit 64 splices the time domain data of each block after correction, the time-depth conversion velocity body determining unit 65 obtains the time-depth conversion velocity body according to the depth migration velocity and the logging velocity of each block, and the time-depth conversion unit 66 conducts time-depth conversion operation on the spliced time domain data according to the time-depth conversion velocity body, so that the spliced depth domain data can be obtained and reprocessed relative to a prestack connected piece.
The time-depth conversion rate body determination unit 65 includes: an initial velocity body determination module 651 and a time-depth transition velocity body determination module 652, wherein:
an initial velocity body determining module 651, configured to perform a control layer pickup operation on the spliced time domain data, and fill an initial velocity in the control layer according to the depth offset velocity of each block to obtain an initial velocity body;
and the time-depth conversion velocity body determining module 652 is configured to perform an inversion operation on the initial velocity body by using the logging velocity as an objective function, so as to obtain a time-depth conversion velocity body corresponding to the logging velocity.
Specifically, the correction unit 63 includes: a frequency phase processing module 631 and a static correction operation module 632, wherein:
the frequency and phase processing module 631 is configured to perform unified operations of frequency and phase on the time domain data of each block respectively;
and a static correction operation module 632, configured to perform static correction operation on the repeated area between the blocks.
The static correction operation module 632 includes: a time difference determination sub-module 6321 and a static correction operation sub-module 6322, where:
a time difference determination submodule 6321, configured to select one block from the blocks as a basic block, and determine the time difference between the other blocks and the basic block by picking up an extremum of the cross correlation;
the static correction operation sub-module 6322 is configured to perform a smoothing operation on the time difference determined by each block, and perform a correction operation on the block in a static correction manner.
For specific execution processes of the units, the modules, and the sub-modules, reference may be made to the description in the foregoing method embodiments, and details are not described here again.
In practical operation, the units, the modules and the sub-modules may be combined or may be arranged singly, and the present invention is not limited thereto.
The present embodiment also provides an electronic device, which may be a desktop computer, a tablet computer, a mobile terminal, and the like, but is not limited thereto. In this embodiment, the electronic device may be implemented by referring to the above method embodiment and the depth domain seismic data splicing apparatus embodiment, and the contents thereof are incorporated herein, and repeated details are not repeated.
Fig. 7 is a schematic block diagram of a system configuration of an electronic apparatus 600 according to an embodiment of the present invention. As shown in fig. 7, the electronic device 600 may include a central processor 100 and a memory 140; the memory 140 is coupled to the central processor 100. Notably, this diagram is exemplary; other types of structures may also be used in addition to or in place of the structure to implement telecommunications or other functions.
In one embodiment, the depth domain seismic data stitching functionality may be integrated into the central processor 100. The central processor 100 may be configured to control as follows:
obtaining depth domain seismic data and corresponding well log data for a plurality of blocks, the depth domain seismic data comprising: a depth migration velocity, the well log data comprising: logging speed;
converting the depth domain seismic data of the block to which the seismic data belong into time domain data based on the depth migration velocity of each block, and respectively correcting the time domain data of each block;
splicing the time domain data of each corrected block based on a preset splicing rule;
and obtaining a time-depth conversion speed body according to the depth deviation speed and the logging speed of each block, and performing time-depth conversion operation on the spliced time domain data according to the time-depth conversion speed body to obtain the spliced depth domain data.
As can be seen from the above description, according to the electronic device provided in the embodiment of the present application, the depth domain seismic data of the block to which the electronic device belongs is converted into the time domain data based on the depth migration velocity of each block, the time domain data of each block is corrected, the time domain data of each block after the correction is spliced, a time-depth conversion velocity body is obtained according to the depth migration velocity and the logging velocity of each block, and the time-depth conversion operation is performed on the spliced time domain data according to the time-depth conversion velocity body, so that the spliced depth domain data can be obtained, and the splicing scheme of the depth domain seismic data is low in cost, short in period, and easy to implement, compared with the pre-stack connection piece.
In another embodiment, the depth domain seismic data stitching device may be configured separately from the central processor 100, for example, the depth domain seismic data stitching device may be configured as a chip connected to the central processor 100, and the depth domain seismic data stitching function is realized by the control of the central processor.
As shown in fig. 7, the electronic device 600 may further include: communication module 110, input unit 120, audio processing unit 130, display 160, power supply 170. It is noted that the electronic device 600 does not necessarily include all of the components shown in fig. 7; furthermore, the electronic device 600 may also comprise components not shown in fig. 7, which may be referred to in the prior art.
As shown in fig. 7, the central processor 100, sometimes referred to as a controller or operational control, may include a microprocessor or other processor device and/or logic device, the central processor 100 receiving input and controlling the operation of the various components of the electronic device 600.
The memory 140 may be, for example, one or more of a buffer, a flash memory, a hard drive, a removable media, a volatile memory, a non-volatile memory, or other suitable device. The information relating to the failure may be stored, and a program for executing the information may be stored. And the central processing unit 100 may execute the program stored in the memory 140 to realize information storage or processing, etc.
The input unit 120 provides input to the cpu 100. The input unit 120 is, for example, a key or a touch input device. The power supply 170 is used to provide power to the electronic device 600. The display 160 is used to display an object to be displayed, such as an image or a character. The display may be, for example, an LCD display, but is not limited thereto.
The memory 140 may be a solid state memory such as Read Only Memory (ROM), Random Access Memory (RAM), a SIM card, or the like. There may also be a memory that holds information even when power is off, can be selectively erased, and is provided with more data, an example of which is sometimes called an EPROM or the like. The memory 140 may also be some other type of device. Memory 140 includes buffer memory 141 (sometimes referred to as a buffer). The memory 140 may include an application/function storage section 142, and the application/function storage section 142 is used to store application programs and function programs or a flow for executing the operation of the electronic device 600 by the central processing unit 100.
The memory 140 may also include a data store 143, the data store 143 for storing data, such as contacts, digital data, pictures, sounds, and/or any other data used by the electronic device. The driver storage portion 144 of the memory 140 may include various drivers of the electronic device for communication functions and/or for performing other functions of the electronic device (e.g., messaging application, address book application, etc.).
The communication module 110 is a transmitter/receiver 110 that transmits and receives signals via an antenna 111. The communication module (transmitter/receiver) 110 is coupled to the central processor 100 to provide an input signal and receive an output signal, which may be the same as in the case of a conventional mobile communication terminal.
Based on different communication technologies, a plurality of communication modules 110, such as a cellular network module, a bluetooth module, and/or a wireless local area network module, may be provided in the same electronic device. The communication module (transmitter/receiver) 110 is also coupled to a speaker 131 and a microphone 132 via an audio processor 130 to provide audio output via the speaker 131 and receive audio input from the microphone 132 to implement general telecommunications functions. Audio processor 130 may include any suitable buffers, decoders, amplifiers and so forth. In addition, an audio processor 130 is also coupled to the central processor 100, so that recording on the local can be enabled through a microphone 132, and so that sound stored on the local can be played through a speaker 131.
Embodiments of the present invention also provide a computer-readable storage medium, on which a computer program is stored, where the computer program is executed by a processor to implement the steps of the depth domain seismic data stitching method.
In summary, the present invention provides a scheme for implementing splicing of consecutive blocks after stacking by using depth domain achievement data and well data, which can implement fast splicing of data of adjacent blocks collected and processed in different ages, and has low cost, short period and easy implementation compared with the preprocessing of consecutive blocks. The time-depth conversion speed used in the embodiment of the invention is a new speed obtained by taking the logging speed as the target function for inversion, so that the corresponding situation of spliced data and well data is better, and the imaging depth is more accurate. Moreover, the spliced data realizes the consistency of the construction depth, the frequency and the phase, avoids the one-sided recognition of scattered block data, provides indispensable favorable conditions for the integral promotion and summary of oil and gas geological recognition such as regional structure, deposition, reservoir formation and the like from the region, greatly reduces the multi-solution of the geological recognition and provides powerful basis for the decision of oil and gas exploration and development.
The preferred embodiments of the present invention have been described above with reference to the accompanying drawings. The many features and advantages of the embodiments are apparent from the detailed specification, and thus, it is intended by the appended claims to cover all such features and advantages of the embodiments which fall within the true spirit and scope thereof. Further, since numerous modifications and changes will readily occur to those skilled in the art, it is not desired to limit the embodiments of the invention to the exact construction and operation illustrated and described, and accordingly, all suitable modifications and equivalents may be resorted to, falling within the scope thereof.
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 principle and the implementation mode of the invention are explained by applying specific embodiments in the invention, and the description of the embodiments is only used for helping to understand the method and the core idea of the invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present invention.
Claims (12)
1. A method of depth domain seismic data stitching, the method comprising:
obtaining depth domain seismic data and corresponding well log data for a plurality of blocks, the depth domain seismic data comprising: a depth migration velocity, the well log data comprising: logging speed;
converting the depth domain seismic data of the block to which the seismic data belong into time domain data based on the depth migration velocity of each block, and respectively correcting the time domain data of each block;
splicing the time domain data of each corrected block based on a preset splicing rule;
and obtaining a time-depth conversion speed body according to the depth deviation speed and the logging speed of each block, and performing time-depth conversion operation on the spliced time domain data according to the time-depth conversion speed body to obtain the spliced depth domain data.
2. The method of claim 1, wherein obtaining a time-depth conversion velocity volume according to the depth migration velocity and the logging velocity of each block comprises:
performing control layer pickup operation on the spliced time domain data, and filling an initial speed into a control layer according to the depth migration speed of each block to obtain an initial speed body;
and carrying out inversion operation on the initial velocity body by taking the logging speed as a target function to obtain a time-depth conversion velocity body corresponding to the logging speed.
3. The method of claim 1, wherein the performing the calibration process on the time domain data of each block respectively comprises:
respectively carrying out unified operation of frequency and phase on the time domain data of each block;
and performing static correction operation on the repeated areas among the blocks.
4. The method of claim 3, wherein performing static correction operations on the repeated regions between the blocks comprises:
selecting one block from the blocks as a basic block, and determining the time difference between the other blocks and the basic block by picking up the extreme value of the cross correlation;
and performing smoothing operation on the time difference determined by each block, and performing correction operation on the block in a static correction mode.
5. The method of claim 1, wherein the splicing the time domain data of each corrected block based on a predetermined splicing rule comprises:
and splicing the time domain data of each corrected block, wherein the repeated area between the blocks adopts a weighted superposition mode, and the boundary of each block is subjected to preset filtering processing.
6. A depth domain seismic data stitching apparatus, comprising:
a data acquisition unit for acquiring depth domain seismic data and corresponding well log data for a plurality of blocks, the depth domain seismic data comprising: a depth migration velocity, the well log data comprising: logging speed;
the time domain conversion unit is used for converting the depth domain seismic data of the block to which the seismic data belong into time domain data based on the depth migration velocity of each block;
a correction unit for performing correction processing on the time domain data of each block;
the splicing unit is used for splicing the time domain data of each corrected block based on a preset splicing rule;
the time-depth conversion speed body determining unit is used for obtaining a time-depth conversion speed body according to the depth deviation speed and the logging speed of each block;
and the time-depth conversion unit is used for performing time-depth conversion operation on the spliced time domain data according to the time-depth conversion speed body to obtain the spliced depth domain data.
7. The apparatus according to claim 6, wherein the time-depth transition speed body determining unit includes:
the initial velocity body determining module is used for carrying out control layer pickup operation on the spliced time domain data and filling the initial velocity into the control layer according to the depth migration velocity of each block so as to obtain an initial velocity body;
and the time-depth conversion velocity body determining module is used for carrying out inversion operation on the initial velocity body by taking the logging speed as a target function to obtain a time-depth conversion velocity body corresponding to the logging speed.
8. The apparatus of claim 6, wherein the correction unit comprises:
the frequency phase processing module is used for respectively carrying out unified operation of frequency and phase on the time domain data of each block;
and the static correction operation module is used for performing static correction operation on the repeated areas among the blocks.
9. The apparatus of claim 8, wherein the static correction operation module comprises:
the time difference determining submodule is used for selecting one block from the blocks as a basic block, and the other blocks determine the time difference between the other blocks and the basic block by picking up the extreme value of cross correlation;
and the static correction operation sub-module is used for performing smoothing operation on the time difference determined by each block and performing correction operation on the block in a static correction mode.
10. The apparatus according to claim 6, wherein the splicing unit is specifically configured to:
and splicing the time domain data of each corrected block, wherein the repeated area between the blocks adopts a weighted superposition mode, and the boundary of each block is subjected to preset filtering processing.
11. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the steps of the method according to any of claims 1 to 5 are implemented when the processor executes the program.
12. A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010602759.2A CN113933901B (en) | 2020-06-29 | 2020-06-29 | Depth domain seismic data splicing method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010602759.2A CN113933901B (en) | 2020-06-29 | 2020-06-29 | Depth domain seismic data splicing method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113933901A true CN113933901A (en) | 2022-01-14 |
| CN113933901B CN113933901B (en) | 2024-01-30 |
Family
ID=79272592
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010602759.2A Active CN113933901B (en) | 2020-06-29 | 2020-06-29 | Depth domain seismic data splicing method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113933901B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116699684A (en) * | 2023-06-09 | 2023-09-05 | 中海石油(中国)有限公司深圳分公司 | Reservoir prediction-oriented seismic data splicing method, device, equipment and medium |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103226210A (en) * | 2013-04-01 | 2013-07-31 | 中国石油天然气股份有限公司 | Method for realizing closure |
| CN105259581A (en) * | 2015-10-22 | 2016-01-20 | 中国石油化工股份有限公司 | Seismic data time-depth conversion method |
| CN109188522A (en) * | 2018-10-09 | 2019-01-11 | 中国石油天然气股份有限公司 | Velocity field construction method and device |
-
2020
- 2020-06-29 CN CN202010602759.2A patent/CN113933901B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103226210A (en) * | 2013-04-01 | 2013-07-31 | 中国石油天然气股份有限公司 | Method for realizing closure |
| CN105259581A (en) * | 2015-10-22 | 2016-01-20 | 中国石油化工股份有限公司 | Seismic data time-depth conversion method |
| CN109188522A (en) * | 2018-10-09 | 2019-01-11 | 中国石油天然气股份有限公司 | Velocity field construction method and device |
Non-Patent Citations (3)
| Title |
|---|
| 吕铁良 等: "稀疏脉冲反演技术在井间地震反演中的应用", 石油物探, vol. 46, no. 1, pages 59 * |
| 孟宪军 等: "井间地震资料时间域波阻抗反演研究", 石油学报, vol. 26, no. 1, pages 47 - 49 * |
| 张志明 等: "时深转换中的井震联合速度建模方法研究与应用现状", 地球物理学进展, vol. 31, no. 5, pages 2276 - 2282 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116699684A (en) * | 2023-06-09 | 2023-09-05 | 中海石油(中国)有限公司深圳分公司 | Reservoir prediction-oriented seismic data splicing method, device, equipment and medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113933901B (en) | 2024-01-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104181599B (en) | A kind of refractive wave static correction processing method based near surface formation and system | |
| CN112415591A (en) | Diffracted wave imaging method and device, electronic equipment and storage medium | |
| CN111596347A (en) | Method and device for rapidly calculating surface layer longitudinal and transverse wave velocity ratio | |
| CN109696705A (en) | Deep-sea OBS uplink and downlink Seismic imaging data correcting method and device | |
| CN113933901A (en) | Depth domain seismic data splicing method and device | |
| CN114384584B (en) | Crack modeling method and device | |
| CN112711064B (en) | Initial layer velocity modeling method and device | |
| CN115480296A (en) | Method and device for predicting the bottom boundary of gas hydrate stability zone | |
| CN112255679B (en) | Method and device for processing diffraction depth deviation of seismic data | |
| CN114578428B (en) | High-efficiency construction method and device for speed field | |
| CN112462421B (en) | Reservoir information prediction method and device, electronic equipment and storage medium | |
| CN114114417B (en) | Earthquake imaging method and device | |
| CN112817042B (en) | Method and device for identifying oil and gas in sedimentary reservoir | |
| CN112505748B (en) | Gun domain reflected wave pickup method and device | |
| CN112462426B (en) | Transverse wave vector static correction method and device | |
| CN112363225B (en) | Transverse wave micro-logging interpretation method and device | |
| CN112462423A (en) | Method and device for predicting thickness of biological reef | |
| CN113075747B (en) | Reservoir fracture development area prediction method and device | |
| CN114578445A (en) | Method and device for determining fracture position based on gravity data | |
| CN113960655B (en) | Method and system for updating seismic data samples | |
| CN112285782B (en) | Near-surface seismic wave absorption attenuation investigation method and device | |
| CN118330739B (en) | Depth domain data volume acquisition method, system and storage medium based on well-seismic joint constraint | |
| CN114578418A (en) | Partition denoising method and device based on near-surface type | |
| CN113970783B (en) | Carbonate reservoir prediction method and device | |
| CN112764107B (en) | Seismic data static correction method and device for high-speed coverage area of ground surface |
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 |