CN112711064B - Initial layer velocity modeling method and device - Google Patents

Initial layer velocity modeling method and device Download PDF

Info

Publication number
CN112711064B
CN112711064B CN201911016588.9A CN201911016588A CN112711064B CN 112711064 B CN112711064 B CN 112711064B CN 201911016588 A CN201911016588 A CN 201911016588A CN 112711064 B CN112711064 B CN 112711064B
Authority
CN
China
Prior art keywords
data
horizon
layer
migration
stratum
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201911016588.9A
Other languages
Chinese (zh)
Other versions
CN112711064A (en
Inventor
首皓
曾同生
王兴
崔栋
秦楠
张征
宋雅莹
康敬程
黄荣善
张小乐
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201911016588.9A priority Critical patent/CN112711064B/en
Publication of CN112711064A publication Critical patent/CN112711064A/en
Application granted granted Critical
Publication of CN112711064B publication Critical patent/CN112711064B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/282Application of seismic models, synthetic seismograms
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30Analysis
    • G01V1/303Analysis for determining velocity profiles or travel times
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/622Velocity, density or impedance
    • G01V2210/6222Velocity; travel time
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/66Subsurface modeling

Abstract

The invention provides an initial layer velocity modeling method and a device, wherein the method comprises the following steps: acquiring time migration data, and picking up horizon data of a stratum according to the time migration data; performing time migration before the overlapping on the horizon data through a zero offset distance to obtain a horizon reverse migration data body; performing prestack depth migration on the horizon reverse migration data volume to obtain migration horizon data; obtaining theoretical depth domain horizons of corresponding stratums through logging layering and/or tectonic graph data; updating the layer speed of the corresponding stratum according to the depth difference between the offset layer data and the theoretical depth domain layer; and constructing and obtaining an initial layer velocity model through the layer velocity after stratum updating.

Description

Initial layer velocity modeling method and device
Technical Field
The application relates to the technical field of seismic data processing, in particular to an initial interval velocity modeling method and device.
Background
The prestack depth migration needs to be subjected to layer velocity modeling, the means for updating the velocity model comprise grid chromatography, full waveform inversion and the like, the method for updating the velocity model needs to have an initial layer velocity model, and the accuracy of the initial layer velocity model determines the velocity and the accuracy of subsequent velocity updating.
In order to obtain an initial layer velocity model, a conventional method includes converting prestack time migration velocity into layer velocity, or interpolating logging information to obtain layer velocity, and the like.
Disclosure of Invention
The invention aims to provide an initial layer velocity modeling method and device so as to shorten the time for manually establishing an initial layer velocity model and improve the accuracy of initial layer velocity modeling.
To achieve the above object, the method for modeling initial layer velocity provided by the present invention specifically comprises: acquiring time migration data, and picking up horizon data of a stratum according to the time migration data; carrying out time migration before the reverse folding on the horizon data through zero offset distance to obtain a horizon reverse migration data body; performing prestack depth migration on the horizon reverse migration data volume to obtain migration horizon data; obtaining theoretical depth domain horizons of corresponding stratums through logging layering and/or tectonic graph data; updating the layer speed of the corresponding stratum according to the depth difference between the offset layer data and the theoretical depth domain layer; and constructing and obtaining an initial layer velocity model through the layer velocity after stratum updating.
In the above initial layer velocity modeling method, preferably, picking horizon data of the formation from the time offset data includes: the time migration data is pre-stack time migration data or post-stack time migration data; according to the time offset data, picking up horizon data of the stratum from shallow to deep; and taking the offset floating surface with the pre-stack depth offset or the post-stack depth offset as the initial horizon of the shallowest layer, wherein the picked horizon is below the initial horizon.
In the above initial horizon velocity modeling method, preferably, performing time migration before eclipsing on the horizon data by a zero offset distance to obtain a horizon reverse migration data volume includes: taking the horizon data as input data of reverse migration, performing reverse pre-stack time migration after the cheap distance of the horizon data is assigned to be 0, and obtaining the speed for the pre-stack time migration; obtaining the speed of the prestack time migration from the speed used for the prestack time migration
In the above initial horizon velocity modeling method, preferably, the performing prestack depth migration on the horizon reverse migration data volume to obtain migration horizon data includes: obtaining the layer speed of the corresponding layer position by carrying out DIX formula conversion on the layer position reverse offset data volume; and obtaining offset horizon data according to the horizon speed.
In the above method for modeling initial interval velocity, preferably, obtaining the theoretical depth domain horizon of the corresponding formation from the well logging layering and/or formation map data includes: and obtaining theoretical depth domain horizons through logging layered interpolation and/or formation map data.
In the above initial layer velocity modeling method, preferably, the updating the layer velocity of the corresponding formation according to the depth difference between the offset horizon data and the theoretical depth domain horizon includes: the formation velocity of the formation is obtained by calculating the following formula: v2= v1 x h2/h1; in the above formula, v1 is the layer velocity in the offset layer data, v2 is the layer velocity of the stratum obtained by calculation, h1 is the height difference between the layer after depth offset and the layer adjacent to the stratum, and h2 is the height difference between the layer in the theoretical depth domain and the layer adjacent to the stratum; and the adjacent upper horizon of the horizon after the depth deviation and the adjacent upper horizon of the theoretical depth domain horizon are the same horizon.
In the above method for modeling initial interval velocity, preferably, the obtaining of the initial interval velocity model by constructing the interval velocity after the formation update includes: and taking the current theoretical depth domain horizon as an initial stratum, picking up the stratum below the current theoretical depth domain horizon, updating the stratum speed of each stratum repeatedly until the stratum speeds corresponding to all the stratums are updated, and obtaining an initial stratum speed model according to the stratum speeds after the stratum updating.
The invention also provides an initial layer velocity modeling device, which comprises a data acquisition module, an analysis module and a construction module; the data acquisition module is used for acquiring time migration data and picking up horizon data of a stratum according to the time migration data; performing time migration before the overlapping on the horizon data through a zero offset distance to obtain a horizon reverse migration data body; performing prestack depth migration on the horizon reverse migration data volume to obtain migration horizon data; the analysis module is used for obtaining theoretical depth domain horizons of corresponding stratums through logging layering and/or tectonic map data; the construction module is used for updating the layer speed of the corresponding stratum according to the depth difference between the offset layer data and the theoretical depth domain layer; and constructing and obtaining an initial layer velocity model through the layer velocity after stratum updating.
The invention also provides a computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method when executing the computer program.
The present invention also provides a computer-readable storage medium storing a computer program for executing the above method.
The invention has the beneficial technical effects that: the time domain horizon is reversely offset to obtain the travel time information of the simulated seismic data of the horizon, the depth offset of the horizon is performed on the reversely offset horizon to obtain the depth information of the horizon after the depth offset, at the moment, the theoretical target depth is obtained based on well data or construction diagram and other data, and the horizon velocity of the stratum is obtained through the difference between the theoretical target depth and the depth after the depth offset.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the principle of the invention. In the drawings:
FIG. 1 is a schematic flow chart diagram illustrating a method for modeling an initial layer velocity according to an embodiment of the present invention;
FIG. 2 is a schematic diagram illustrating a flow chart of generating an anti-migration data volume according to an embodiment of the present invention;
FIG. 3 is a schematic structural diagram of an initial layer velocity modeling apparatus according to an embodiment of the present invention;
fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present invention.
FIG. 5 is a diagram illustrating picking horizon data according to an embodiment of the invention.
FIG. 6 is a schematic diagram of picking up first-set horizon data according to an embodiment of the invention.
FIG. 7 is a diagram illustrating first set of anti-horizon offset data according to an embodiment of the invention.
FIG. 8 is a schematic diagram of a horizon and a well layer after a depth offset of a first set of horizons according to an embodiment of the present invention.
FIG. 9 is a diagram illustrating a first set of horizon velocities according to an embodiment of the invention.
Fig. 10 is a schematic diagram of a velocity model after completion of all horizon updating according to an embodiment of the present invention.
Detailed Description
The following detailed description will be given with reference to the accompanying drawings and examples to explain how to apply the technical means to solve the technical problems and to achieve the technical effects. It should be noted that, unless otherwise specified, the embodiments and features of the embodiments of the present invention may be combined with each other, and the technical solutions formed are within the scope of the present invention.
Additionally, the steps illustrated in the flowcharts of the figures may be performed in a computer system such as a set of computer-executable instructions, and while a logical order is illustrated in the flowcharts, in some cases, the steps illustrated or described may be performed in an order different than here.
Referring to fig. 1, the initial layer velocity modeling method provided by the present invention specifically includes: s101, time migration data are obtained, and horizon data of a stratum are picked up according to the time migration data; s102, performing time migration before the overlap on the horizon data through a zero offset distance to obtain a horizon reverse migration data body; s103, performing prestack depth migration on the horizon reverse migration data volume to obtain migration horizon data; s104, obtaining theoretical depth domain horizons of corresponding stratums through logging layering and/or tectonic graph data; s105, updating the layer speed of the corresponding stratum according to the depth difference between the offset layer data and the theoretical depth domain layer; s106, constructing and obtaining an initial layer velocity model through the layer velocity after stratum updating; therefore, the time-lapse information of the simulated seismic data of the horizon can be obtained by reversely offsetting the horizon of the time domain through the embodiment, the depth information of the horizon after depth offset can be obtained by performing depth offset on the reversely offset horizon, at the moment, the theoretical target depth is obtained based on well data or structural diagram and other data, and the horizon speed of the stratum is obtained through the difference between the theoretical target depth and the depth after depth offset.
In the above embodiment, the overall process of the initial layer velocity modeling method may be divided into the following contents in practical application: picking formation horizons on the pre-stack time migration data; performing time migration before the reverse folding on the picked horizon data according to a zero offset distance to obtain a horizon reverse migration data body; performing prestack depth migration on the horizon reverse migration data volume; forming theoretical depth domain horizons of corresponding stratums based on logging layering or tectonic graph data; updating the layer speed of the corresponding stratum according to the depth difference between the layer position after the prestack depth deviation and the layer position of the theoretical depth domain; and (5) picking up the next set of stratum, and repeating the steps until all the stratums are updated to obtain an initial stratum velocity model. Specifically, in the above embodiment, the step S101 of picking up horizon data of the formation according to the time offset data includes: the time migration data may be pre-stack time migration data or post-stack time migration data; according to the time offset data, picking up horizon data of the stratum from shallow to deep; and taking the offset floating surface with the pre-stack depth offset or the post-stack depth offset as an initial layer of the shallowest layer, wherein the picked layer is below the initial layer.
Referring to fig. 2, in the above embodiment, fig. 5 is a schematic diagram of horizon data, where horizons of 5 earth strata including the earth surface and the underground are respectively represented by (1) - (5), taking the first earth strata as an example, fig. 6 is a first earth strata picked up separately, and in step S102, performing time offset before underlap on the horizon data by using zero offset distance to obtain a horizon anti-offset data volume (as shown in fig. 7) includes: s201, taking the horizon data as input data of reverse migration, assigning the cheap distance of the horizon data as 0, and then performing reverse pre-stack time migration to obtain the speed for the pre-stack time migration; s202 obtains the speed of the pre-stack time migration from the speed used for the pre-stack time migration.
In the above embodiment, performing prestack depth migration on the horizon reverse-migration data volume in step S103, and obtaining migration horizon data (as shown in fig. 8) includes: obtaining the layer speed of the corresponding layer position through DIX formula conversion on the layer position reverse offset data volume; and obtaining offset horizon data according to the horizon speed.
In the above embodiment, as shown in fig. 8, the obtaining the theoretical depth domain horizon of the corresponding formation through the well logging stratigraphic and/or structural map data in step S104 includes: and obtaining theoretical depth domain horizons through logging layered interpolation and/or tectonic map data.
In the above embodiment, the step S105 of updating the layer velocity of the corresponding formation according to the depth difference between the offset horizon data and the theoretical depth domain horizon includes: the stratum velocity of the stratum is obtained through the following formula: v2= v1 x h2/h1; in the above formula, v1 is the layer velocity in the offset layer position data, v2 is the layer velocity of the stratum obtained by calculation, h1 is the height difference between the layer position after depth offset and the layer position adjacent to the stratum, and h2 is the height difference between the layer position in the theoretical depth domain and the layer position adjacent to the stratum; and the adjacent upper horizon of the horizon after the depth deviation and the adjacent upper horizon of the theoretical depth domain horizon are the same horizon.
In the above embodiment, the obtaining of the initial layer velocity model (as shown in fig. 9) by the layer velocity construction after formation updating in step S106 includes: taking the current theoretical depth domain horizon as an initial stratum, picking up the stratum below the current theoretical depth domain horizon, and updating the stratum speed of each stratum repeatedly until the stratum speeds corresponding to all the stratums are updated, and then obtaining an initial stratum speed model according to the stratum speeds after the stratum updating (as shown in fig. 10).
Referring to fig. 3, the present invention further provides an initial layer velocity modeling apparatus, which includes a data acquisition module, an analysis module, and a construction module; the data acquisition module is used for acquiring time offset data and picking up horizon data of the stratum according to the time offset data; carrying out time migration before the reverse folding on the horizon data through zero offset distance to obtain a horizon reverse migration data body; performing prestack depth migration on the horizon reverse migration data volume to obtain migration horizon data; the analysis module is used for obtaining theoretical depth domain horizons of corresponding stratums through logging layering and/or tectonic graph data; the construction module is used for updating the layer speed of the corresponding stratum according to the depth difference between the offset horizon data and the theoretical depth domain layer; and constructing and obtaining an initial layer velocity model through the layer velocity after stratum updating.
The invention also provides a computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method when executing the computer program.
The present invention also provides a computer-readable storage medium storing a computer program for executing the above method.
The beneficial technical effects of the invention are as follows: the time domain horizon is reversely offset to obtain the travel time information of the simulated seismic data of the horizon, the depth offset of the horizon is performed on the reversely offset horizon to obtain the depth information of the horizon after the depth offset, at the moment, the theoretical target depth is obtained based on well data or construction diagram and other data, and the horizon velocity of the stratum is obtained through the difference between the theoretical target depth and the depth after the depth offset.
As shown in fig. 4, 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. 4; furthermore, the electronic device 600 may also comprise components not shown in fig. 4, which may be referred to in the prior art.
As shown in fig. 4, 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 devices. The information relating to the failure may be stored, and a program for executing the information may be stored. And the cpu 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 referred to as 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 a communication function and/or for performing other functions of the electronic device (e.g., a messaging application, a directory 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, etc., 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.
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-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (10)

1. A method for initial layer velocity modeling, the method comprising:
obtaining time offset data, and picking horizon data of a stratum according to the time offset data;
performing time migration before the overlapping on the horizon data through a zero offset distance to obtain a horizon reverse migration data body;
performing prestack depth migration on the horizon reverse migration data volume to obtain migration horizon data;
obtaining theoretical depth domain horizons of corresponding stratums through logging layering and/or tectonic graph data;
updating the layer speed of the corresponding stratum according to the depth difference between the offset layer data and the theoretical depth domain layer;
and constructing and obtaining an initial layer velocity model through the layer velocity after stratum updating.
2. The method of initial interval velocity modeling according to claim 1, wherein picking horizon data for a formation from the time-shifted data comprises:
the time migration data is pre-stack time migration data or post-stack time migration data;
according to the time migration data, picking up horizon data of the stratum from shallow to deep; and taking the offset floating surface with the pre-stack depth offset or the post-stack depth offset as the initial horizon of the shallowest layer, wherein the picked horizon is below the initial horizon.
3. The method of modeling initial layer velocity of claim 1, wherein time shifting the layer data before the eclipsing by a zero offset distance to obtain a layer reverse shift data volume comprises:
taking the horizon data as input data of reverse migration, assigning the migration distance of the horizon data as 0, and then performing reverse pre-stack time migration to obtain the speed for the pre-stack time migration;
the speed of the pre-stack time migration is obtained from the speed used for the pre-stack time migration.
4. The method of initial interval velocity modeling according to claim 1, wherein performing prestack depth migration on said volume of horizon inverse migration data, obtaining migration horizon data comprises: obtaining the layer speed of the corresponding layer position by carrying out DIX formula conversion on the layer position reverse offset data volume; offset horizon data is obtained from the horizon velocities.
5. The method of modeling initial interval velocity of claim 1, wherein obtaining theoretical depth domain horizons for the corresponding formation from well log layering and/or formation map data comprises: and obtaining theoretical depth domain horizons through logging layered interpolation and/or formation map data.
6. The method of initial interval velocity modeling according to claim 1, wherein updating the interval velocities of the corresponding strata according to the depth difference between the offset horizon data and the theoretical depth domain horizons comprises:
the formation velocity of the formation is obtained by calculating the following formula:
v2=v1*h2/h1;
in the above formula, v1 is the layer velocity in the offset layer position data, v2 is the layer velocity of the stratum obtained by calculation, h1 is the height difference between the layer position after depth offset and the layer position adjacent to the stratum, and h2 is the height difference between the layer position in the theoretical depth domain and the layer position adjacent to the stratum; and the adjacent upper horizon of the horizon after the depth deviation is the same as the adjacent upper horizon of the theoretical depth domain horizon.
7. The method of modeling initial layer velocity according to claim 1, wherein obtaining the initial layer velocity model through layer velocity construction after formation updating comprises: and taking the current theoretical depth domain horizon as an initial stratum, picking up the stratum below the current theoretical depth domain horizon, updating the stratum speed of each stratum repeatedly until the stratum speeds corresponding to all the stratums are updated, and obtaining an initial stratum speed model according to the stratum speeds after the stratum updating.
8. The initial layer velocity modeling device is characterized by comprising a data acquisition module, an analysis module and a construction module;
the data acquisition module is used for acquiring time offset data and picking up horizon data of the stratum according to the time offset data; carrying out time migration before the reverse folding on the horizon data through zero offset distance to obtain a horizon reverse migration data body; performing prestack depth migration on the horizon reverse migration data volume to obtain migration horizon data;
the analysis module is used for obtaining theoretical depth domain horizons of corresponding stratums through logging layering and/or tectonic graph data;
the construction module is used for updating the layer speed of the corresponding stratum according to the depth difference between the offset horizon data and the theoretical depth domain layer; and constructing and obtaining an initial layer velocity model through the layer velocity after stratum updating.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any of claims 1 to 7 when executing the computer program.
10. 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 7.
CN201911016588.9A 2019-10-24 2019-10-24 Initial layer velocity modeling method and device Active CN112711064B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911016588.9A CN112711064B (en) 2019-10-24 2019-10-24 Initial layer velocity modeling method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911016588.9A CN112711064B (en) 2019-10-24 2019-10-24 Initial layer velocity modeling method and device

Publications (2)

Publication Number Publication Date
CN112711064A CN112711064A (en) 2021-04-27
CN112711064B true CN112711064B (en) 2022-11-04

Family

ID=75541349

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911016588.9A Active CN112711064B (en) 2019-10-24 2019-10-24 Initial layer velocity modeling method and device

Country Status (1)

Country Link
CN (1) CN112711064B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102565853A (en) * 2011-12-20 2012-07-11 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Method for modeling velocity model of geologic structure
CN102879814A (en) * 2011-07-15 2013-01-16 中国石油天然气集团公司 Accurate depth domain layer speed updating method
CN104360385A (en) * 2014-11-26 2015-02-18 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Method for comprehensively establishing initial depth interval velocity model by combining seismogeology understanding
CN104865597A (en) * 2014-02-20 2015-08-26 中国石油化工股份有限公司 Modeling method of depth domain interval velocity initial model
CN105974479A (en) * 2016-06-30 2016-09-28 恒泰艾普(北京)能源科技研究院有限公司 Chromatographic 2D/3D anisotropic depth domain velocity modeling method for GPU (Graphics Processing Unit) spatial grids
CN106094032A (en) * 2016-08-30 2016-11-09 中国石油集团川庆钻探工程有限公司地球物理勘探公司 A kind of method building formation velocity model
CN106597533A (en) * 2016-11-17 2017-04-26 中国石油化工股份有限公司 Depth domain velocity modeling method for piedmont zone seismic data processing
CN109581501A (en) * 2019-01-18 2019-04-05 中国石油化工股份有限公司 Method for desert area Depth Domain velocity modeling

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8902709B2 (en) * 2008-07-18 2014-12-02 William Marsh Rice University Methods for concurrent generation of velocity models and depth images from seismic data

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102879814A (en) * 2011-07-15 2013-01-16 中国石油天然气集团公司 Accurate depth domain layer speed updating method
CN102565853A (en) * 2011-12-20 2012-07-11 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Method for modeling velocity model of geologic structure
CN104865597A (en) * 2014-02-20 2015-08-26 中国石油化工股份有限公司 Modeling method of depth domain interval velocity initial model
CN104360385A (en) * 2014-11-26 2015-02-18 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Method for comprehensively establishing initial depth interval velocity model by combining seismogeology understanding
CN105974479A (en) * 2016-06-30 2016-09-28 恒泰艾普(北京)能源科技研究院有限公司 Chromatographic 2D/3D anisotropic depth domain velocity modeling method for GPU (Graphics Processing Unit) spatial grids
CN106094032A (en) * 2016-08-30 2016-11-09 中国石油集团川庆钻探工程有限公司地球物理勘探公司 A kind of method building formation velocity model
CN106597533A (en) * 2016-11-17 2017-04-26 中国石油化工股份有限公司 Depth domain velocity modeling method for piedmont zone seismic data processing
CN109581501A (en) * 2019-01-18 2019-04-05 中国石油化工股份有限公司 Method for desert area Depth Domain velocity modeling

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
速度建模技术在吉林探区复杂构造成像处理中的应用;邱敏 等;《中外能源》;20130315;第18卷(第03期);63-69 *

Also Published As

Publication number Publication date
CN112711064A (en) 2021-04-27

Similar Documents

Publication Publication Date Title
CN112415591A (en) Diffracted wave imaging method and device, electronic equipment and storage medium
CN104181599A (en) Refracted wave static correction handling method and system based on near-surface layer
CN105116448B (en) A kind of converted wave azimuthal anisotropy bearing calibration and device
RU2014136576A (en) SYSTEM AND METHOD OF VERTICAL SEISMIC PROFILING WITH REPRESENTATION OF EXPLORATION DATA AS A COMBINATION OF PARAMETRIZED COMPRESSION, SHEAR AND DISPERSIVE WAVE FIELDS
CN104199103A (en) Static correction method and static correction device for converted wave
CN105093319A (en) Ground micro-seismic static correction method based on three-dimensional seismic data
CN110376645A (en) A kind of determination method, device and equipment of Thin Sandbody top interface location
CN109696705A (en) Deep-sea OBS uplink and downlink Seismic imaging data correcting method and device
CN112711064B (en) Initial layer velocity modeling method and device
CN113433587B (en) Rock physical inversion method and system
CN114384584A (en) Crack modeling method and device
CN113933901B (en) Depth domain seismic data splicing method and device
CN114114392A (en) Layer velocity model establishing method and device
CN112255679B (en) Method and device for processing diffraction depth deviation of seismic data
CN112363225A (en) Transverse wave micro-logging interpretation method and device
CN112505748A (en) Shot domain reflected wave pickup method and device
CN114114417B (en) Earthquake imaging method and device
CN113671572A (en) Indoor sand box-based seismic data imaging method and device
CN113960655B (en) Method and system for updating seismic data samples
CN111474573A (en) Seismic source excitation control system, method and device for seismic exploration
CN111754611B (en) Automatic drawing method and system for inter-well sand body communication diagram based on geologic model constraint
CN110632650B (en) Time-sharing window extraction method and device for full-waveform inversion data
CN111366973A (en) Frequency domain noise generation and addition method and device of forward model
CN112817042B (en) Method and device for identifying oil and gas in sedimentary reservoir
CN112462426A (en) Shear wave vector static correction method and device

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