CN111239822A - Fault lower speed modeling method for eliminating well seismic closure difference and processing terminal - Google Patents
Fault lower speed modeling method for eliminating well seismic closure difference and processing terminal Download PDFInfo
- Publication number
- CN111239822A CN111239822A CN202010126380.9A CN202010126380A CN111239822A CN 111239822 A CN111239822 A CN 111239822A CN 202010126380 A CN202010126380 A CN 202010126380A CN 111239822 A CN111239822 A CN 111239822A
- Authority
- CN
- China
- Prior art keywords
- model
- depth
- velocity
- velocity model
- well
- 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
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. analysis, for interpretation, for correction
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- 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. analysis, for interpretation, for correction
- G01V1/282—Application of seismic models, synthetic seismograms
-
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
Abstract
The invention relates to a fault lower speed modeling method for eliminating well seismic closure difference and a processing terminal, wherein the method comprises the following steps: step 1: establishing a depth domain initial layer velocity model according to the seismic data; step 2: establishing a construction model; and step 3: picking up the residual time difference of the common imaging point gather; and 4, step 4: updating and iterating the initial layer velocity model of the depth domain; and 5: performing prestack depth migration calculation to obtain a prestack depth migration profile and a common imaging point gather of prestack depth migration; step 6: the well mark and the horizon are correlated to establish a structural model containing well information and well layers; and 7: maintaining travel time tomography; and 8: performing prestack depth migration calculation on the updated and iterated layer velocity model again; and step 9: and repeating the optimization iteration of the depth domain layer velocity model to obtain a final imaging result profile. The invention improves the imaging precision and the precision of the speed model and effectively reduces the exploration risk.
Description
Technical Field
The invention relates to the technical field of seismic exploration, in particular to a method for modeling the speed under a fault by eliminating well seismic closure difference and a processing terminal.
Background
In seismic exploration, particularly in oil and gas reservoir exploration, complex fault block imaging and structure implementation become key factors, and a series of geological risk factors corresponding to structure trapping implementation, reservoir carving, accurate drilling target entering and the like are all urgent to be considered in actual exploration, so that the geological exploration is usually faced with two problems: imaging accuracy and structural distortion.
The prestack depth migration is a currently accepted method for effectively solving seismic imaging, but the imaging precision of the prestack depth migration depends on the rationality of a depth-velocity model to a great extent, so that the depth-velocity model is usually optimized by means of a tomography method (namely, a tomography inversion method) and data driving, the tomography method needs to explain a horizon by modeling along the horizon, and the influence of a fracture system on the model is not considered; for relying on data-driven, it makes full use of the acquired seismic data information to perform velocity modeling. Due to the limitation of the principle of the tomography method, uncertainty is easy to occur in tomography inversion speed modeling during inversion solving, so that the inverted speed model is difficult to be matched with a real geological structure, particularly, under the condition of complex fault block development, the speed modeling by using the tomography inversion method is difficult to meet the exploration evaluation requirement of a target area, and one important reason is that well seismic closure difference cannot be eliminated, so that the seismic depth and the well depth have errors.
At present, no fault lower speed modeling method which is good and effective for eliminating well seismic closure difference exists. The relevant references are as follows:
1. the application of the complex structure modeling technology in seismic inversion [ J ]. geophysical exploration of petroleum, 2018,53 (supplement & supplement): 270-class 275.
2. Scheffer, Dumin, Wenpeng, etc. the application of the grid chromatography speed inversion method in quasi-three-dimensional Xisha hydrate [ J ]. geophysical prospecting and chemical prospecting, 2017,41(5): 846-.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a fault lower speed modeling method for eliminating well seismic closure difference, which can solve the problem of eliminating the well seismic closure difference;
the second object of the present invention is to provide a processing terminal, which can solve the problem of eliminating the poor closing of the well earthquake.
The technical scheme for realizing one purpose of the invention is as follows: a method for modeling the speed under a fault by eliminating the well seismic closure difference comprises the following steps:
step 1: acquiring seismic data, obtaining a common midpoint gather according to the seismic data, and establishing a depth domain initial layer velocity model according to the common midpoint gather;
step 2: defining fault and stratum properties, loading fault and stratum information to establish a fault grid and fault blocks, then establishing a layer grid, and constructing a geological grid according to the fault grid, the fault blocks and the layer grid so as to establish a structural model;
and step 3: picking up the residual time difference of the common imaging point gather;
and 4, step 4: performing at least two times of update iteration processing on the depth domain initial layer velocity model by adopting a tomography inversion method of structural model constraint to obtain a final updated depth domain layer velocity model, wherein the updated depth domain layer velocity model comprises velocity three-component information, the velocity three-component information comprises inclination angle, continuity and azimuth angle information,
the residual time difference is used as soft constraint of geological grid information input, and the soft constraint controls the deviation of a current depth domain layer velocity model and a last depth domain layer velocity model in each updating iteration process;
and 5: performing prestack depth migration calculation according to the depth domain layer velocity model obtained in the step 4 to obtain a prestack depth migration profile and a common imaging point gather of prestack depth migration;
step 6: associating the well marks with the horizons, reestablishing a structural model containing well information and well layers according to the method in the step 2, and calculating well seismic closure differences corresponding to each well layer;
and 7: obtaining a depth domain layer velocity model after well seismic closure difference updating iteration by keeping travel time tomography, thereby obtaining a depth domain layer velocity model after well seismic closure difference updating iteration under fault control, velocity three-component information of the depth domain layer velocity model and an updated anisotropic parameter model;
and 8: simultaneously updating the obtained three-component information of the velocity into the construction model in the step 6, recalculating a new round of well seismic closure difference for the new construction model, and performing prestack depth migration calculation on the updated and iterated layer velocity model obtained in the step 8 again according to the newly calculated well seismic closure difference and the updated anisotropic parameter model in the step 7 to obtain a corresponding prestack depth migration gather and a prestack depth migration profile;
and step 9: and (5) repeating the steps 6 to 8 to perform optimization iteration on the depth domain layer velocity model until the depth domain layer velocity model after the optimization iteration meets the imaging geological requirements, and performing prestack depth migration volume migration according to the final depth domain layer velocity model to obtain a final imaging result profile.
Further, the establishing of the depth domain initial layer velocity model according to the common midpoint gather includes the following steps:
and (3) picking out the root mean square velocity of the time domain from the common midpoint gather, establishing a time domain initial layer velocity model which accords with the geological velocity change rule by the root mean square velocity of the time domain through a CVI constraint velocity inversion method, and obtaining a depth domain initial layer velocity model by time domain-depth domain conversion.
Further, the common midpoint gather is processed, wherein the processing comprises one or more of preprocessing, noise suppression, horizontal superposition and superposition offset, and the root mean square velocity of the time domain is picked out from the processed common midpoint gather.
Further, in step 5, the tomography inversion method is any one of grid tomography, structural model-based tomography, and travel-time-keeping tomography.
The second technical scheme for realizing the aim of the invention is as follows: a processing terminal, comprising,
a memory for storing program instructions;
a processor for executing the program instructions to perform the steps in the method for modeling downhole velocities for eliminating well seismic closure differentials.
The invention has the beneficial effects that: according to the invention, through fault control, the velocity model and the imaging precision are ensured, the problem of errors of seismic depth and well depth is effectively solved, reliable seismic data are provided for target evaluation, the exploration risk is effectively reduced, and the precision of the velocity model is improved.
Drawings
FIG. 1 is a schematic flow diagram of the present invention;
FIG. 2 is a schematic diagram of a depth domain initial layer velocity model;
FIG. 3 is a schematic diagram of a tectonic model established under fault control;
FIG. 4 is a schematic illustration of a residual time difference;
FIG. 5 is an updated depth domain layer velocity model obtained by a tomographic inversion method with formation model constraints;
FIG. 6(a) is a velocity comparison diagram of the depth domain initial layer velocity model and the depth domain layer velocity model obtained through tomography at the point A of the fault footwall;
FIG. 6(B) is a velocity comparison diagram of the depth domain initial layer velocity model and the tomographic depth domain layer velocity model at the B point on the slice;
FIG. 7 is a schematic diagram showing the superposition of a depth domain interval velocity model and a seismic prestack section;
FIG. 8 is a schematic diagram of the dip angle of the velocity three-component information of the depth domain layer velocity model obtained by the processing of this step;
FIG. 9 is an azimuth schematic diagram of the velocity three-component information of the depth domain layer velocity model obtained by the processing of this step;
FIG. 10 is a schematic view of the joint display of four wells (numbered W01-W04) with fault and horizon model information;
FIG. 11 is a three-dimensional visualization display of well seismic closure differences;
FIG. 12 is a schematic diagram showing a well layer velocity model and a well layer display before travel time tomography inversion (the well depth and the horizon depth after velocity inversion have a certain depth error);
FIG. 13 is a schematic view showing a depth domain interval velocity model and a well stratification after travel time tomography inversion (the well depth and the depth of the interval after velocity inversion have a certain depth error of zero); fig. 14 is a schematic structural diagram of a processing terminal according to the present invention.
Detailed description of the preferred embodiments
The invention will be further described with reference to the accompanying drawings and specific embodiments:
as shown in fig. 1-13, a method for modeling the velocity of a fault by eliminating the closing difference of well-seisms comprises the following steps:
step 1: acquiring seismic data including shot line data, performing conventional denoising processing on the seismic data in a sequence time domain to obtain a common midpoint gather, and establishing a depth domain initial layer velocity model according to the common midpoint gather. The root mean square velocity of the time domain can be used for establishing a time domain initial layer velocity model which accords with the geological velocity change rule through a CVI constraint velocity inversion method, and the time-depth (namely the time domain-depth domain) conversion is utilized to obtain a depth domain initial layer velocity model. The rms velocity in the time domain may be picked from the common midpoint gather after conventional processing, which generally includes preprocessing, various noise suppression (linear or non-linear noise, multiples, etc.), horizontal stacking, and stack shifting. And then, carrying out travel time calculation based on ray tracking on the depth domain initial layer velocity model, and carrying out prestack depth migration calculation on the depth domain initial layer velocity model according to the calculated travel time so as to obtain a prestack depth migration profile and a common imaging point gather (CIG) of prestack depth migration. FIG. 2 is a schematic diagram of a depth-domain initial layer velocity model, wherein the velocity values range from 3651m/s to 6500 m/s.
Step 2: and establishing a construction model, wherein the construction model comprises a fault model and a horizon model. Firstly, defining fault and stratum properties, loading fault and stratum information to establish a fault grid and fault blocks, then establishing a layer grid, and establishing a geological grid according to the fault grid, the fault blocks and the layer grid, thereby establishing a structural model. Fig. 3 is a schematic diagram of a structural model established under fault control, which includes fault grids and main horizon information.
And step 3: and picking up the residual time difference of the common imaging point gather, wherein the residual time difference can be used as the estimation of the non-levelness degree of the common imaging point gather of the pre-stack depth deviation, namely leveling the common imaging point gather. The method for calculating the residual time difference is the prior art and is not described herein again. Fig. 4 is a schematic diagram of the residual time difference.
And 4, step 4: and (3) taking the depth domain initial layer velocity model obtained in the step (1), the construction model obtained in the step (2) and the residual time difference obtained in the step (3) as a first comprehensive database, wherein the first comprehensive database is used for providing basic data for subsequent chromatographic inversion.
And 5: and (3) performing updating iteration processing on the depth domain initial layer velocity model in the step (1) for at least two times by adopting a tomography inversion method based on the structural model constraint in the step (2) to obtain a final updated depth domain layer velocity model. In the updating iteration processing process, the residual time difference is used as soft constraint of geological grid information input, the constraint item controls the deviation of the current layer velocity model and the previous layer velocity model in each iteration process, so that the updated and iterated depth domain layer velocity model and velocity three-component information of the depth domain layer velocity model are obtained, and the velocity three-component information comprises inclination angle, continuity and azimuth angle information. The depth domain initial layer velocity model which is not updated does not have velocity three-component information, and only the depth domain layer velocity model after at least one tomography iteration based on the structural model constraint. FIG. 5 is an updated depth domain layer velocity model obtained by a tomography inversion method constrained by a structural model, and the updated depth domain layer velocity model can accurately depict underground structural features and can truly represent complex underground geological features. FIG. 6(a) is a velocity comparison diagram of the depth domain initial layer velocity model and the depth domain layer velocity model obtained through tomography at the point A of the fault footwall; fig. 6(B) is a velocity comparison diagram of the depth domain initial layer velocity model and the tomography depth domain layer velocity model at the disk B point on the fault, wherein the left graph is the depth domain initial layer velocity model, and the right graph is the tomography depth domain layer velocity model. In fig. 6 (including (a) and (b)), taking the inline 720 as an example, it is found by comparison that: the speed of the initial layer speed model in the depth domain at the point A of the lower wall of the fault is 4308.24m/s, the speed of the point B of the upper wall of the fault is 4306.7m/s, and the speed difference at the upper wall and the lower wall of the fault (point BA) is 1.54 m/s; the speed of the depth domain layer velocity model obtained through tomography at the point A of the lower disc of the fault is 4346.79m/s, the speed of the point B of the upper disc of the fault is 4504.5m/s, and the difference of the speeds at the upper disc and the lower disc of the fault (point BA) is 157.71 m/s. Therefore, the depth domain layer velocity model obtained after tomography based on the structural model can accurately depict the velocity difference between the upper and lower plates of the broken layer.
FIG. 7 is a schematic diagram of superposition display of a depth domain layer velocity model and a seismic prestack section, and from FIG. 7, it can be seen that the depth domain layer velocity model obtained through structural model constraint and a corresponding structural model keep high geological consistency, a velocity field is high in matching degree with real underground geology, fault imaging is effectively depicted, and the accurate position of fault trapping is realized.
In this step, the tomography inversion method generally includes methods such as grid tomography, tomography based on a structural model, and travel-time-keeping tomography, and the tomography inversion method is substantially an optimization problem, and aims to find a layer velocity model to minimize an objective function (or called a cost function), where the objective function mathematically quantifies the unflattened degree of a gather under the conditions of satisfying specific soft constraints (geological information input), hard constraints (well seismic closure difference: the maximum matching between a migration result and well data), and geological factor constraints, and the tomography inversion method belongs to the prior art and is not described in detail.
Step 6: and (5) taking the speed three-component information of the depth domain layer speed model obtained through the first round of updating iterative processing in the step (5) and the geological grid information as input to participate in the building of the structural model in the step (2), so as to obtain a new structural model. Meanwhile, pre-stack depth migration calculation is carried out on the depth domain layer velocity model obtained through the first round of updating iteration processing, and a new pre-stack depth migration profile and a new pre-stack depth migration common imaging point gather are obtained. FIG. 8 is a schematic diagram of the dip angles of the velocity three-component information of the depth domain layer velocity model obtained by the processing of this step, wherein the dip angle range of the velocity three-component is approximately about 0-40 °; fig. 9 is an azimuth schematic diagram of the three-component velocity information of the depth domain layer velocity model obtained by the processing in this step, wherein the azimuth range of the three-component velocity is approximately 10-320 ° in azimuth.
And then, repeating the steps 2-6 to perform optimization iteration of the depth domain layer velocity model until the objective function C (m) reaches the minimum value, and otherwise, continuing the optimization iteration. When the objective function C (m) reaches the minimum value, the optimized and iterated depth domain layer velocity model accords with imaging geological requirements, a final depth domain layer velocity model is correspondingly obtained, pre-stack depth migration calculation is carried out on the final depth domain layer velocity model, a final imaging result profile is obtained, namely, a high-precision layer velocity model under fault control is obtained, and velocity modeling is completed.
Through the processing of the steps 1 to 6, a high-precision layer velocity model under fault control can be obtained, the geological velocity modeling problem is solved, a set of targeted solution is provided for the problems (structure trapping and closing, reservoir carving and accurate drilling target entering) faced by the current exploration and development, the structural complexity of an underground medium is fully considered, the change characteristics of the velocities on two sides of a fault block are accurately carved, and the velocity model with the constraint of the geological structure model and the corresponding structure model keep high geological consistency.
And 7: and (5) directly entering the step, and performing prestack depth migration calculation according to the depth domain layer velocity model obtained in the step (5) to obtain a prestack depth migration profile and a common imaging point gather of prestack depth migration.
And 8: and (3) associating the well mark (marker) with the horizon, reestablishing a structural model containing well information and well layers according to the method in the step (2), and calculating the well seismic closure difference corresponding to each well layer. FIG. 10 is a schematic diagram of the joint display of four wells (numbered W01-W04) and fault and horizon model information. FIG. 11 is a three-dimensional visualization of well seismic closure differences, which are shown to have a depth range of-200-0.
And step 9: and (5) establishing a second comprehensive database, wherein the second comprehensive database comprises the construction model in the step (8), the well seismic closure error and the depth domain layer velocity model subjected to the chromatographic optimization of the construction model in the step (5).
Step 10: and keeping the travel time tomography to obtain well seismic closure error and updating the iterated depth domain layer velocity model. The method of preserving travel-time tomography is essentially applied by a process of building and solving matrices in which the first thomson parameter δ and the second thomson parameter e are equal in value. In this embodiment, data of the second integrated database are fused, so that a depth domain layer velocity model after fault control downhole seismic closure error update iteration, velocity three-component information of the depth domain layer velocity model, and an updated anisotropic parameter model are obtained. Wherein, the travel-time tomography process is actually a process Ax ═ b of a linear equation set solved by tomographic velocity inversion, wherein A is a matrix with dimension m multiplied by n, x is an array with dimension n multiplied by 1, and the physical meaning represented by the process is a slowness value (slowness change value) in each grid; b is an m x 1 dimensional array representing the physical meaning of the travel time (travel time difference) of the ray; dimension m represents the total number of rays and n is the number of meshes of the velocity model. The updated anisotropic parameter model will be used as input information for the second integrated database in the subsequent travel-time tomography iteration process. Since the tomographic imaging method while keeping traveling belongs to the prior art, it will not be described herein. Fig. 12 is a schematic diagram showing a depth domain layer velocity model and a well layer display before travel time tomography inversion (a certain depth error exists between the well depth and the depth of a layer after velocity inversion), and it can be seen from fig. 12 that a certain depth error exists between the well depth of four wells numbered W01-W04 and the depth of a layer after velocity inversion, wherein the position of a vertical third well marker of each well is not coincident with the depth of a nearby horizontal layer (a solid line curve), that is, a certain depth error exists between the well depth of four wells numbered W01-W04 and the depth of a layer after velocity inversion; and fig. 13 is a schematic diagram showing a depth domain layer velocity model and a well layer after travel time tomography inversion (the error between the well depth and the depth of a position after velocity inversion is zero), and the position of a vertical third well marker of each well in the diagram is overlapped with the depth (a solid curve) of a position nearby and transverse, so that the problem of the error between the seismic depth and the well depth is effectively solved, and a direct data basis can be provided for operations such as accurate drilling into a target.
Step 11: and simultaneously updating the three-component information of the velocity obtained in the step 10 into the structural model in the step 8, recalculating the new well-seismic closure difference for the new structural model, updating the newly calculated well-seismic closure difference and the updated anisotropic parameter model in the step 10 into a second comprehensive database to serve as a new second comprehensive database, and performing prestack depth migration calculation on the updated and iterated layer velocity model obtained in the step 10 to obtain a corresponding prestack depth migration gather and a prestack depth migration profile.
Step 12: and after the step 11 is completed, repeating the step 8 to the step 11 to perform optimization iteration on the depth domain layer velocity model until the depth domain layer velocity model after the optimization iteration meets the imaging geological requirements, and performing prestack depth migration volume migration according to the final depth domain layer velocity model to obtain a final imaging result profile.
Through the processing of the steps 1-5 and 7-12, a set of targeted solution can be provided for a series of problems of structural confinement implementation, reservoir carving, accurate target entering of drilling and the like faced by exploration and development at present, the problems of errors of seismic depth and well depth are effectively solved by ensuring a velocity model and imaging precision through fault control, reliable seismic data are provided for target evaluation, exploration risks are effectively reduced, and the precision of the velocity model is improved.
As shown in fig. 14, the present invention also provides an entity implementation processing terminal 100 of a high-precision velocity modeling method under fault control, which includes,
a memory 101 for storing program instructions;
a processor 102 for executing the program instructions to perform the steps of the high precision velocity modeling method under fault control.
The embodiments disclosed in this description are only an exemplification of the single-sided characteristics of the invention, and the scope of protection of the invention is not limited to these embodiments, and any other functionally equivalent embodiments fall within the scope of protection of the invention. Various other changes and modifications to the above-described embodiments and concepts will become apparent to those skilled in the art from the above description, and all such changes and modifications are intended to be included within the scope of the present invention as defined in the appended claims.
Claims (5)
1. A method for modeling the speed under a fault by eliminating well seismic closure difference is characterized by comprising the following steps:
step 1: acquiring seismic data, obtaining a common midpoint gather according to the seismic data, and establishing a depth domain initial layer velocity model according to the common midpoint gather;
step 2: defining fault and stratum properties, loading fault and stratum information to establish a fault grid and fault blocks, then establishing a layer grid, and constructing a geological grid according to the fault grid, the fault blocks and the layer grid so as to establish a structural model;
and step 3: picking up the residual time difference of the common imaging point gather;
and 4, step 4: performing at least two times of update iteration processing on the depth domain initial layer velocity model by adopting a tomography inversion method of structural model constraint to obtain a final updated depth domain layer velocity model, wherein the updated depth domain layer velocity model comprises velocity three-component information, the velocity three-component information comprises inclination angle, continuity and azimuth angle information,
the residual time difference is used as soft constraint of geological grid information input, and the soft constraint controls the deviation of a current depth domain layer velocity model and a last depth domain layer velocity model in each updating iteration process;
and 5: performing prestack depth migration calculation according to the depth domain layer velocity model obtained in the step 4 to obtain a prestack depth migration profile and a common imaging point gather of prestack depth migration;
step 6: associating the well marks with the horizons, reestablishing a structural model containing well information and well layers according to the method in the step 2, and calculating well seismic closure differences corresponding to each well layer;
and 7: obtaining a depth domain layer velocity model after well seismic closure difference updating iteration by keeping travel time tomography, thereby obtaining a depth domain layer velocity model after well seismic closure difference updating iteration under fault control, velocity three-component information of the depth domain layer velocity model and an updated anisotropic parameter model;
and 8: simultaneously updating the obtained three-component information of the velocity into the construction model in the step 6, recalculating a new round of well seismic closure difference for the new construction model, and performing prestack depth migration calculation on the updated and iterated layer velocity model obtained in the step 8 again according to the newly calculated well seismic closure difference and the updated anisotropic parameter model in the step 7 to obtain a corresponding prestack depth migration gather and a prestack depth migration profile;
and step 9: and (5) repeating the steps 6 to 8 to perform optimization iteration on the depth domain layer velocity model until the depth domain layer velocity model after the optimization iteration meets the imaging geological requirements, and performing prestack depth migration volume migration according to the final depth domain layer velocity model to obtain a final imaging result profile.
2. The method for modeling the downhole velocity of a fault with elimination of the well-to-seismic closure difference as claimed in claim 1, wherein the establishing of the initial interval velocity model of the depth domain according to the common midpoint gather comprises the following steps:
and (3) picking out the root mean square velocity of the time domain from the common midpoint gather, establishing a time domain initial layer velocity model which accords with the geological velocity change rule by the root mean square velocity of the time domain through a CVI constraint velocity inversion method, and obtaining a depth domain initial layer velocity model by time domain-depth domain conversion.
3. The method of claim 2, wherein the common midpoint gather is processed, wherein the processing comprises one or more of preprocessing, noise suppression, horizontal stacking, and stack offset, and wherein the root mean square velocity in the time domain is picked from the processed midpoint gather.
4. The method for modeling downhole velocities for eliminating well-to-seismic closure differences according to claim 1, wherein in the step 5, the tomography inversion method is any one of grid tomography, structural model-based tomography, and travel-time-keeping tomography.
5. A processing terminal, characterized in that it comprises,
a memory for storing program instructions;
a processor for executing the program instructions to perform the steps in the method of subsurface velocity modeling with elimination of well seismic closure differences as claimed in any of claims 1-4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010126380.9A CN111239822B (en) | 2020-02-27 | 2020-02-27 | Fault lower speed modeling method for eliminating well seismic closure difference and processing terminal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010126380.9A CN111239822B (en) | 2020-02-27 | 2020-02-27 | Fault lower speed modeling method for eliminating well seismic closure difference and processing terminal |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111239822A true CN111239822A (en) | 2020-06-05 |
| CN111239822B CN111239822B (en) | 2022-08-23 |
Family
ID=70875062
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010126380.9A Active CN111239822B (en) | 2020-02-27 | 2020-02-27 | Fault lower speed modeling method for eliminating well seismic closure difference and processing terminal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111239822B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115877449A (en) * | 2021-09-28 | 2023-03-31 | 中国石油化工股份有限公司 | Computer-implemented method for obtaining images of subsurface stacks within a survey area |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050038605A1 (en) * | 2003-08-14 | 2005-02-17 | Marin Irina A. | Method and apparatus for contemporaneous utilization of a higher order probe in pre-stack and post-stack seismic domains |
| CN103852789A (en) * | 2014-03-12 | 2014-06-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Nonlinear chromatography method and device for seismic data |
| CN107355205A (en) * | 2017-08-08 | 2017-11-17 | 广州海洋地质调查局 | A kind of hydraulic slotted liner technique method of semisubmersible drilling platform or drill ship |
| CN109738945A (en) * | 2018-11-08 | 2019-05-10 | 成都捷科思石油天然气技术发展有限公司 | A method of structural map is directly generated using pre-stack depth migration achievement |
-
2020
- 2020-02-27 CN CN202010126380.9A patent/CN111239822B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050038605A1 (en) * | 2003-08-14 | 2005-02-17 | Marin Irina A. | Method and apparatus for contemporaneous utilization of a higher order probe in pre-stack and post-stack seismic domains |
| CN103852789A (en) * | 2014-03-12 | 2014-06-11 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Nonlinear chromatography method and device for seismic data |
| CN107355205A (en) * | 2017-08-08 | 2017-11-17 | 广州海洋地质调查局 | A kind of hydraulic slotted liner technique method of semisubmersible drilling platform or drill ship |
| CN109738945A (en) * | 2018-11-08 | 2019-05-10 | 成都捷科思石油天然气技术发展有限公司 | A method of structural map is directly generated using pre-stack depth migration achievement |
Non-Patent Citations (2)
| Title |
|---|
| 薛花，等: "层析成像反演速度建模方法研究", 《CT理论与应用研究》 * |
| 韩令贺，等: "井震联合网格层析各向异性速度建模研究及应用", 《岩性油气藏》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115877449A (en) * | 2021-09-28 | 2023-03-31 | 中国石油化工股份有限公司 | Computer-implemented method for obtaining images of subsurface stacks within a survey area |
| CN115877449B (en) * | 2021-09-28 | 2024-02-13 | 中国石油化工股份有限公司 | Computer-implemented method for obtaining subsurface superimposed images within a survey area |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111239822B (en) | 2022-08-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104597490B (en) | Multi-wave AVO reservoir elastic parameter inversion method based on accurate Zoeppritz equations | |
| Radwan | Three-dimensional gas property geological modeling and simulation | |
| CN110361781B (en) | Sedimentary-phase ancient landform restoration method based on seismic reflection amplitude analysis | |
| CN107783187B (en) | Method for establishing three-dimensional velocity field by combining logging velocity and seismic velocity | |
| CN104459782B (en) | Horizon chromatography inversion earthquake speed modeling method with thin layers used and modeling unit | |
| US10732310B2 (en) | Seismic attributes derived from the relative geological age property of a volume-based model | |
| CN105319589B (en) | A kind of fully automatic stereo chromatography conversion method using local lineups slope | |
| CN109541685B (en) | River channel sand body identification method | |
| CN110244021B (en) | Stratum layering method based on anisotropic interpolation | |
| CN113031068B (en) | Reflection coefficient accurate base tracking prestack seismic inversion method | |
| CN109799540B (en) | Volcanic rock type uranium deposit magnetic susceptibility inversion method based on geological information constraint | |
| CN110531410A (en) | A kind of least square reverse-time migration gradient Preconditioning method based on through wave field | |
| Barthwal et al. | Passive seismic tomography using recorded microseismicity: Application to mining-induced seismicity | |
| CN115877449A (en) | Computer-implemented method for obtaining images of subsurface stacks within a survey area | |
| CN108957554B (en) | Seismic inversion method in geophysical exploration | |
| CN111239822B (en) | Fault lower speed modeling method for eliminating well seismic closure difference and processing terminal | |
| Al-Khazraji et al. | Uncertainty analysis to assess depth conversion accuracy: A case study of Subba Oilfield, Southern Iraq | |
| CN111257969B (en) | High-precision speed modeling method under fault control and processing terminal | |
| CN111198400B (en) | Speed modeling method capable of eliminating well seismic closure difference and processing terminal | |
| CN107179553A (en) | Desert area surface static correction method based on double Compaction Laws | |
| CN111208558B (en) | Method and device for establishing ultra-deep low-amplitude three-dimensional geological structure | |
| CN111999768B (en) | Construction elevation determination method based on construction model in drilling layered design | |
| CN105093324A (en) | Residual static correction method | |
| CN112649893B (en) | Multi-data and multi-parameter fusion modeling method and system for thin reservoir | |
| CN114861515A (en) | Method, device, equipment and medium for calculating layer speed data volume |
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 | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20201014 Address after: 510000, compound 477, Ring East Road, Guangzhou, Guangdong Applicant after: GUANGZHOU MARINE GEOLOGICAL SURVEY Applicant after: Guangdong Provincial Laboratory of marine science and engineering of South China (Guangzhou) Address before: 510000, compound 477, Ring East Road, Guangzhou, Guangdong Applicant before: GUANGZHOU MARINE GEOLOGICAL SURVEY |
|
| TA01 | Transfer of patent application right | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |