CN111458751A - Comprehensive time depth calibration method and device - Google Patents
Comprehensive time depth calibration method and device Download PDFInfo
- Publication number
- CN111458751A CN111458751A CN201910057366.5A CN201910057366A CN111458751A CN 111458751 A CN111458751 A CN 111458751A CN 201910057366 A CN201910057366 A CN 201910057366A CN 111458751 A CN111458751 A CN 111458751A
- Authority
- CN
- China
- Prior art keywords
- breakpoint
- seismic
- well
- distance
- goodness
- 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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
Abstract
The invention provides a comprehensive time depth calibration method and a device, wherein the method comprises the following steps: comparing the stratums among wells according to the logging curves of all wells to obtain the breakpoint position and the breakpoint distance of each well breakpoint; inputting the well logging curve, the breakpoint position and the break distance into an interpretation system to obtain a synthetic seismic record; carrying out well seismic calibration according to the wave resistance characteristics and the synthetic seismic record to obtain a time depth calibration result; inputting the time-depth calibration result, the breakpoint position and the fault distance into the seismic section, determining the goodness of fit of the breakpoints of all wells and the seismic breakpoints, and determining the goodness of fit of the fault distance of the breakpoints of all wells and the fault distance of the seismic upper and lower plates; and adjusting the synthetic seismic record according to the wave resistance characteristics until the coincidence degree of each well breakpoint and the seismic breakpoint and the coincidence degree of the fracture distance of the seismic upper and lower discs and the fracture distance of each well breakpoint reach preset values. The result of the invention can meet the demand of a geophysical operator for constructing and explaining reservoir prediction, and the subsequent result can provide an oil reservoir operator for analyzing the oil reservoir result.
Description
Technical Field
The invention relates to the field of geophysical exploration, in particular to a comprehensive time depth calibration method and a comprehensive time depth calibration device.
Background
In oil exploration, seismic exploration is the most important means for obtaining oil and gas discovery, seismic acquisition signals are time domain data, a well drilling acquisition receipt is depth domain data, time-depth calibration is needed in order to find out the corresponding relation between an earthquake and a well drilling, the quality of the calibration relation has great influence on structure interpretation and reservoir prediction results, and correct recognition of reservoir relation and successful implementation of subsequent drilling are finally influenced. The common geophysical methods are of two types: first, VSP data bridging, also known as corridor stacking; and secondly, calibrating the synthetic seismic records. The two methods can be carried out independently or can be calibrated in a cross-correlation mode. Because of the high cost of VSP data, there is relatively little VSP data in wells in a general work area, and it can be done separately based on wave resistance characteristics. According to the convolution principle, if the reflection coefficient sequence and the wavelets are known, a seismic channel can be synthesized by using an acoustic impedance curve (a one-dimensional forward modeling method), then the synthesized seismic channel is compared with the seismic channel around the well point, if the correlation between the synthesized seismic channel and the seismic channel around the well point is good, the synthesized seismic channel is considered reasonable, the reasonability of the synthesized seismic channel can be determined through manual observation or through computer correlation analysis, if the correlation is high, the time depth is considered reasonable, otherwise, the correlation is considered unreasonable.
The existing comparison method has clear wave resistance relation or special lithology, forms a stronger reflection interface and has no larger error when logging data is reasonable. However, in two cases, errors are easy to occur or the calibration uniqueness cannot be determined due to multiple possibilities, the first case is a fracture complex area, and the wave resistance characteristics of seismic data are often unclear due to the influence of faults; the second situation is that a plurality of sets of sand shale sections are generated due to frequent sand shale interbedding, the seismic reflection characteristics of the sand shale sections are extremely similar, the well seismic contrast is easy to generate the situation of 'stretch-hat plum wear' and up-down contrast error, and even if VSP data are added for correction, the problem of calibration uniqueness is difficult to eliminate because the VSP data are influenced by seismic processing and the VSP data are used and then the synthetic seismic records are required to be stretched and compressed. However, whether the time-depth calibration is correct or not is important for construction interpretation and reservoir prediction. Therefore, how to determine the uniqueness of the time-depth contrast is very important.
Disclosure of Invention
In order to solve the problems of inaccuracy, high error probability and the like of the conventional calibration method, the embodiment of the invention provides a comprehensive time depth calibration method, which comprises the following steps:
comparing the stratums among wells according to the logging curves of all wells to obtain the breakpoint position and the breakpoint distance of each well breakpoint;
inputting the well logging curve, the breakpoint position and the break distance into an interpretation system to obtain a synthetic seismic record;
carrying out well seismic calibration according to the wave resistance characteristics and the synthetic seismic record to obtain a time depth calibration result;
inputting the time-depth calibration result, the breakpoint position and the fault distance into the seismic section, determining the goodness of fit of the breakpoints of all wells and the seismic breakpoints, and determining the goodness of fit of the fault distance of the breakpoints of all wells and the fault distance of the seismic upper and lower plates;
and adjusting the synthetic seismic record according to the wave resistance characteristics until the coincidence degree of each well breakpoint and the seismic breakpoint and the coincidence degree of the fracture distance of the seismic upper and lower discs and the fracture distance of each well breakpoint reach preset values.
Optionally, in an embodiment of the present invention, the comparing the interwell strata according to the log of each well to obtain the breakpoint position and the break distance of each well breakpoint includes: acquiring the logging curve according to the logging data of each well; comparing the horizons and faults among multiple wells according to the logging curve, and determining a full section of the region; and comprehensively analyzing the stratum according to the region full-section, and determining the breakpoint position and the breakpoint distance of each well breakpoint.
Optionally, in an embodiment of the present invention, the adjusting the synthetic seismic record according to the wave impedance characteristics includes: and adjusting the position of the synthetic seismic record up and down according to the wave resistance characteristics to obtain an adjusted time-depth calibration result.
Optionally, in an embodiment of the present invention, the method further includes: and determining the goodness of fit between each well breakpoint and the earthquake breakpoint and the goodness of fit between the fault distance of the earthquake upper and lower plates and the fault distance of each well breakpoint after adjustment according to the adjusted time-depth calibration result.
The embodiment of the invention also provides a comprehensive time depth calibration device, which comprises:
the fracture data unit is used for comparing the stratums among the wells according to the logging curves of all the wells to obtain the fracture positions and the fracture distances of the well breakpoints;
the synthetic seismic unit is used for inputting the well logging curve, the breakpoint position and the break distance into an interpretation system to obtain a synthetic seismic record;
the time-depth calibration unit is used for carrying out well seismic calibration according to the wave resistance characteristics and the synthetic seismic record to obtain a time-depth calibration result;
the goodness of fit determining unit is used for inputting the time-depth calibration result, the breakpoint position and the break distance into the seismic section, determining the goodness of fit of each well breakpoint and the seismic breakpoint, and determining the goodness of fit of the break distance of each well breakpoint and the break distance of the seismic upper and lower plates;
and the goodness of fit adjusting unit is used for adjusting the synthetic seismic record according to the wave resistance characteristics until the goodness of fit of the well breakpoints and the seismic breakpoints and the goodness of fit of the fracture distance of the seismic upper and lower discs and the fracture distance of the well breakpoints reach preset values.
Optionally, in an embodiment of the present invention, the breakpoint data unit includes: the logging curve module is used for acquiring the logging curve according to logging data of each well; the regional full-profile module is used for comparing horizons and faults among multiple wells according to the logging curve and determining a regional full profile; and the breakpoint and fault distance module is used for comprehensively analyzing the stratum according to the full section of the region and determining the breakpoint position and fault distance of each well breakpoint.
Optionally, in an embodiment of the present invention, the goodness-of-fit adjusting unit includes: and the time depth adjusting module is used for adjusting the position of the synthetic seismic record up and down according to the wave resistance characteristics to obtain an adjusted time depth calibration result.
Optionally, in an embodiment of the present invention, the goodness of fit adjusting unit further includes: and the goodness of fit adjusting module is used for determining the goodness of fit of each adjusted well breakpoint and earthquake breakpoint and the goodness of fit of the fault distance of the earthquake upper and lower plates and the fault distance of each well breakpoint according to the adjusted time-depth calibration result.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and capable of running on the processor, and when the processor executes the computer program, the following steps are implemented:
comparing the stratums among wells according to the logging curves of all wells to obtain the breakpoint position and the breakpoint distance of each well breakpoint;
inputting the well logging curve, the breakpoint position and the break distance into an interpretation system to obtain a synthetic seismic record;
carrying out well seismic calibration according to the wave resistance characteristics and the synthetic seismic record to obtain a time depth calibration result;
inputting the time-depth calibration result, the breakpoint position and the fault distance into the seismic section, determining the goodness of fit of the breakpoints of all wells and the seismic breakpoints, and determining the goodness of fit of the fault distance of the breakpoints of all wells and the fault distance of the seismic upper and lower plates;
and adjusting the synthetic seismic record according to the wave resistance characteristics until the coincidence degree of each well breakpoint and the seismic breakpoint and the coincidence degree of the fracture distance of the seismic upper and lower discs and the fracture distance of each well breakpoint reach preset values.
An embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, implements the following steps:
comparing the stratums among wells according to the logging curves of all wells to obtain the breakpoint position and the breakpoint distance of each well breakpoint;
inputting the well logging curve, the breakpoint position and the break distance into an interpretation system to obtain a synthetic seismic record;
carrying out well seismic calibration according to the wave resistance characteristics and the synthetic seismic record to obtain a time depth calibration result;
inputting the time-depth calibration result, the breakpoint position and the fault distance into the seismic section, determining the goodness of fit of the breakpoints of all wells and the seismic breakpoints, and determining the goodness of fit of the fault distance of the breakpoints of all wells and the fault distance of the seismic upper and lower plates;
and adjusting the synthetic seismic record according to the wave resistance characteristics until the coincidence degree of each well breakpoint and the seismic breakpoint and the coincidence degree of the fracture distance of the seismic upper and lower discs and the fracture distance of each well breakpoint reach preset values.
The invention is calibrated by combining geological geophysical methods, and the method principle is that the geology and the earthquake can be perfectly combined under the correct condition, and the geology and the earthquake can be mutually checked. And a complete quality control basis is provided, the result can meet the requirement of a geophysical operator for constructing and explaining reservoir prediction, and the subsequent result can provide an oil reservoir operator for analyzing the oil reservoir result. The calibration result is obtained and used for construction interpretation, reservoir prediction and oil reservoir analysis, so that the accuracy of the result can be improved, and a more accurate basis is provided for oil field deployment.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without inventive labor.
FIG. 1 is a flowchart of a comprehensive time depth calibration method according to an embodiment of the present invention;
FIG. 2A is a cross-sectional top stratigraphic representation of an embodiment of the present invention;
FIG. 2B is a cross-sectional view of a lower portion of a fault in an embodiment of the invention;
FIG. 3 is a time depth calibration diagram when wave drag characteristics are matched according to an embodiment of the present invention;
FIG. 4 is a diagram illustrating matching of well seismic breakpoint locations when breakpoints do not coincide in an embodiment of the present invention;
FIG. 5 is a time depth calibration graph after adjustment according to an embodiment of the present invention;
FIG. 6 is a well seismic breakpoint position matching graph when breakpoints are anastomosed in accordance with an embodiment of the present invention;
fig. 7 is a schematic structural diagram of an integrated time depth calibration apparatus according to an embodiment of the present invention.
Detailed Description
The embodiment of the invention provides a comprehensive time depth calibration method and device.
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Fig. 1 is a flowchart of an integrated time depth calibration method according to an embodiment of the present invention, where the method includes: step S1, cross-well stratum comparison is carried out according to the logging curve of each well, and the breakpoint position and the breakpoint distance of each well breakpoint are obtained;
step S2, inputting the logging curve, the breakpoint position and the break distance into an interpretation system to obtain a synthetic seismic record;
step S3, well seismic calibration is carried out according to the wave resistance characteristics and the synthetic seismic record, and a time depth calibration result is obtained;
step S4, inputting the time-depth calibration result, the breakpoint position and the fault distance into an earthquake section, determining the goodness of fit of each well breakpoint and the earthquake breakpoint, and determining the goodness of fit of the fault distance of each well breakpoint and the fault distance of the upper and lower seismic plates;
and step S5, adjusting the synthetic seismic record according to the wave resistance characteristics until the coincidence degree of the breakpoints of the wells and the seismic breakpoints and the coincidence degree of the fracture distances of the seismic upper and lower discs and the fracture distances of the breakpoints of the wells reach preset values.
As an embodiment of the present invention, the comparing the interwell strata according to the log of each well to obtain the breakpoint position and the break distance of each well breakpoint comprises: acquiring the logging curve according to the logging data of each well; comparing the horizons and faults among multiple wells according to the logging curve, and determining a full section of the region; and comprehensively analyzing the stratum according to the region full-section, and determining the breakpoint position and the breakpoint distance of each well breakpoint.
As an embodiment of the invention, said adjusting said synthetic seismic record according to said wave impedance characteristics comprises: and adjusting the position of the synthetic seismic record up and down according to the wave resistance characteristics to obtain an adjusted time-depth calibration result.
In this embodiment, according to the adjusted time-depth calibration result, the goodness of fit between each well breakpoint and the seismic breakpoint and the goodness of fit between the fracture distance of the seismic upper and lower trays and the fracture distance of each well breakpoint are determined.
The invention is based on well logging layering and seismic wave resistance, and is a comprehensive seismic comparison method, the depth and the distance of a break point obtained by comparison in a well are important rings of calibration, the calibration is considered to be proper only when the break point of the well is matched with the seismic break point and the wave resistance is consistent, otherwise, the reason needs to be found, and the unification of the break point and the seismic break point is achieved through adjustment.
Firstly, the synthetic seismic record is calibrated, but the calibration quality of the synthetic seismic record is not the only standard of the calibration result, and the most important link is that the calibration result needs to be put on a seismic profile to match the goodness of fit of a well breakpoint and a seismic breakpoint, so that the full unification of seismic-geological stratification and breakpoint data is considered, and the principle of combining well medium wave impedance with the seismic profile is also considered. The calibration accuracy is greatly improved, and the non-uniqueness of the calibration result is eliminated. Meanwhile, the used early data are considered to be reliable results.
The invention mainly adopts the following steps:
1) and (3) performing interwell stratum comparison, wherein multiple wells are required for comparison, a comparison sequence of each well is established, and the depth point and the fracture distance of each well breakpoint in the well are required to be compared in a key way.
2) Inputting curve data, layered data, break points and fault distance data into an interpretation system, making a synthetic seismic record, adjusting the position of the synthetic seismic record according to wave resistance characteristics, projecting the calibrated time depth to a corresponding seismic section according to the well point position, observing the goodness of fit of the well break points and the seismic section break points, and estimating the goodness of fit of the seismic upper and lower plate fault distances and the well fault distance.
3) The results are consistent, and the time-depth relation is considered to be reasonable, so that the time-depth relation can be used for construction explanation and reservoir prediction; moving the synthetic seismic record up and down to observe wave resistance characteristics when the synthetic seismic record is not matched, if the wave resistance characteristics are matched, putting the synthetic seismic record to a seismic section, and observing the matching degree of each breakpoint in a well and a seismic breakpoint, if the synthetic seismic record is matched, considering that the calibration is reasonable; if the time depth is not consistent, the time depth needs to be adjusted again.
4) And (5) repeating the operation of the step 3 until the wave resistance is reasonable and the breakpoints are matched, wherein the time depth is reasonable, and the result is used for constructing explanation and reservoir prediction.
In a specific embodiment, in the method, a geological stratification data table, a breakpoint position and a fault distance size stratification data table are confirmed according to logging data, the data tables are loaded into an interpretation system and then can be calibrated conventionally, and the matching relation between the well seismic breakpoint position and the fault distance size is checked after the time depth is determined. And finally, when the calibration is in a proper condition, the method is deeply used for construction explanation and reservoir prediction. The specific process is as follows:
1) interwell contrast determination of horizons and faults
The method is characterized by multi-well comparison, determination of the full section of the region, and comprehensive analysis of the missing stratum to determine the position of each well breakpoint and the well to be compared with the breakpoint. FIGS. 2A and 2B are graphs of the stratification and breakpoint positions and the fracture distance sizes obtained by comparing the Phoenix w-1 well with the surrounding wells.
By multi-well comparison, the characteristics of well logging curve convolution and rhythm are considered, and the similarity of combination curves among multiple wells is considered. And dividing horizon, break point and break distance by using the logging curve.
Taking Phoenix W-1 well as an example, the longitudinal breakpoint position and the distance of the fracture are determined by comparing the Phoenix W-1 well with different wells. If the well 518-602 meters is missing in the 461m interval compared to the Mimosa E-2 well, it is inferred that the well is 84 meters off-line compared to the Mimosa E-2 well, the breakpoint is 416 meters, and the fault is defined as F1. According to the same method, the well is compared with a PhoenixS-1 well for a fault distance of 462 m, the position of a breakpoint is 964 m, and a fault is defined as F2, so the comparison principles are the same, and the description is omitted. The well is divided accordingly and includes a hierarchical data table of breakpoint positions and corresponding horizons (as shown in table 1).
TABLE 1 hierarchical data table (including breakpoint position, break distance size)
2) Synthetic seismic calibration
After reasonable hierarchical data are determined, the hierarchical data (including breakpoint position and fault distance size data) and well logging curve data (including sound wave and density curves) are loaded into an interpretation system, the loaded data are shown in table 1, but the format needs to be adjusted according to different interpretation systems. Well seismic calibration is performed according to wave drag characteristics (fig. 3) and the time depth results are injected into the seismic profile (fig. 4), which is consistent with conventional time depth calibration. Wherein, A1 and B1 in FIG. 4 are the points of non-coincidence.
3) Seismic and logging breakpoint position verification
According to the calibrated time-depth relation, the horizon and the break point are all input into the seismic section, the possible positions of the upper and lower discs are determined according to the seismic wave resistance characteristics, and the seismic fault distance is calculated according to the layer velocity obtained by calibration and compared with the fault distance obtained by logging layering. And meanwhile, confirming whether the well breaking point position is matched with the earthquake breaking point position or not. Through the inspection of fig. 4, the well-seismic breakpoint positions are not matched, and the time-depth calibration result is considered to be wrong. The depth of readjustment (fig. 5) confirmed the matching of the break point, the break distance and the wave resistance (fig. 6), wherein a2 and B2 in fig. 6 are the break point matching points, and the result is considered to be reliable. In rare cases, breakpoints and wave resistance do not coincide, and the earlier-stage results need to be checked as follows:
a. the logging layering is reliable, the breakpoints are used as the standard, and the quality of the logging curve is required to be checked and corrected when the wave resistance characteristics are not matched.
b. If the wave drag characteristics are considered to be reliable, more well analyses need to be added, and the position of the logging breakpoint and the size of the fault distance are further confirmed.
c. And checking the seismic processing speed to ensure that the processing speed is normal.
4) Time depth result determination
Upon further confirmation, the time-depth relationship is determined for use in constructing interpretations and reservoir predictions. And the earthquake, well logging and time depth after the third step of inspection can be used for reservoir prediction work, and a better result can be obtained.
The calibration method provided by the invention is used for calibrating by combining geological geophysical methods, and the method principle is that the geology and the earthquake can be perfectly combined under the correct condition, and the geology and the earthquake can be mutually checked. On the basis of conventional calibration, a complete quality control basis is provided, the result can meet the requirement of a geophysical operator for constructing and explaining reservoir prediction, and the subsequent result can be provided for an oil reservoir operator to analyze the oil reservoir result. The calibration result is obtained and used for construction interpretation, reservoir prediction and oil reservoir analysis, so that the accuracy of the result can be improved, and a more accurate basis is provided for oil field deployment.
Fig. 7 is a schematic structural diagram of an integrated time depth calibration apparatus according to an embodiment of the present invention, where the apparatus includes: the breakpoint data unit 10 is configured to perform interwell stratum comparison according to the logging curve of each well to obtain a breakpoint position and a breakpoint distance of each well breakpoint;
the synthetic seismic unit 20 is used for inputting the logging curve, the breakpoint position and the fault distance into an interpretation system to obtain a synthetic seismic record;
the time-depth calibration unit 30 is used for carrying out well seismic calibration according to the wave resistance characteristics and the synthetic seismic record to obtain a time-depth calibration result;
the goodness of fit determining unit 40 is used for inputting the time-depth calibration result, the breakpoint positions and the breaking distances into the seismic section, determining the goodness of fit of each well breakpoint and the seismic breakpoint, and determining the goodness of fit of the breaking distances of each well breakpoint and the breaking distances of the seismic upper and lower plates;
and the goodness of fit adjusting unit 50 is used for adjusting the synthetic seismic record according to the wave resistance characteristics until the goodness of fit of the well breakpoints and the seismic breakpoints and the goodness of fit of the fracture distance of the seismic upper and lower discs and the fracture distance of the well breakpoints reach preset values.
As an embodiment of the present invention, the breakpoint data unit includes: the logging curve module is used for acquiring the logging curve according to logging data of each well; the regional full-profile module is used for comparing horizons and faults among multiple wells according to the logging curve and determining a regional full profile; and the breakpoint and fault distance module is used for comprehensively analyzing the stratum according to the full section of the region and determining the breakpoint position and fault distance of each well breakpoint.
As an embodiment of the present invention, the goodness-of-fit adjustment unit includes: and the time depth adjusting module is used for adjusting the position of the synthetic seismic record up and down according to the wave resistance characteristics to obtain an adjusted time depth calibration result.
In this embodiment, the goodness of fit adjustment unit further includes: and the goodness of fit adjusting module is used for determining the goodness of fit of each adjusted well breakpoint and earthquake breakpoint and the goodness of fit of the fault distance of the earthquake upper and lower plates and the fault distance of each well breakpoint according to the adjusted time-depth calibration result.
Based on the same application concept as the comprehensive time depth calibration method, the invention also provides the comprehensive time depth calibration device. Because the principle of solving the problems of the comprehensive time depth calibration device is similar to that of a comprehensive time depth calibration method, the implementation of the comprehensive time depth calibration device can refer to the implementation of the comprehensive time depth calibration method, and repeated parts are not repeated.
The calibration device provided by the invention is used for calibrating by combining geological geophysical methods, and the principle is that geology and earthquake can be perfectly combined under the correct condition, and the geology and earthquake can be mutually checked. On the basis of conventional calibration, well seismic breakpoint inspection and fault distance matching after time-depth calibration are added, a complete quality control basis is provided, the result can meet the requirement of a geophysical operator for constructing and explaining reservoir prediction, and a subsequent result can provide an oil reservoir operator for analyzing the oil reservoir result. The calibration result is obtained and used for construction interpretation, reservoir prediction and oil reservoir analysis, so that the accuracy of the result can be improved, and a more accurate basis is provided for oil field deployment.
The embodiment of the present invention further provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and capable of running on the processor, and when the processor executes the computer program, the following steps are implemented:
comparing the stratums among wells according to the logging curves of all wells to obtain the breakpoint position and the breakpoint distance of each well breakpoint;
inputting the well logging curve, the breakpoint position and the break distance into an interpretation system to obtain a synthetic seismic record;
carrying out well seismic calibration according to the wave resistance characteristics and the synthetic seismic record to obtain a time depth calibration result;
inputting the time-depth calibration result, the breakpoint position and the fault distance into the seismic section, determining the goodness of fit of the breakpoints of all wells and the seismic breakpoints, and determining the goodness of fit of the fault distance of the breakpoints of all wells and the fault distance of the seismic upper and lower plates;
and adjusting the synthetic seismic record according to the wave resistance characteristics until the coincidence degree of each well breakpoint and the seismic breakpoint and the coincidence degree of the fracture distance of the seismic upper and lower discs and the fracture distance of each well breakpoint reach preset values.
An embodiment of the present invention further provides a computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, implements the following steps:
comparing the stratums among wells according to the logging curves of all wells to obtain the breakpoint position and the breakpoint distance of each well breakpoint;
inputting the well logging curve, the breakpoint position and the break distance into an interpretation system to obtain a synthetic seismic record;
carrying out well seismic calibration according to the wave resistance characteristics and the synthetic seismic record to obtain a time depth calibration result;
inputting the time-depth calibration result, the breakpoint position and the fault distance into the seismic section, determining the goodness of fit of the breakpoints of all wells and the seismic breakpoints, and determining the goodness of fit of the fault distance of the breakpoints of all wells and the fault distance of the seismic upper and lower plates;
and adjusting the synthetic seismic record according to the wave resistance characteristics until the coincidence degree of each well breakpoint and the seismic breakpoint and the coincidence degree of the fracture distance of the seismic upper and lower discs and the fracture distance of each well breakpoint reach preset values.
Based on the same application concept as the comprehensive time depth calibration method, the invention also provides the computer equipment and the computer readable storage medium. Since the principle of solving the problem of the computer device and the computer-readable storage medium is similar to that of the comprehensive time-depth calibration method, the implementation of the computer device and the computer-readable storage medium can refer to the implementation of the comprehensive time-depth calibration method, and repeated parts are not repeated.
The computer device and the computer readable storage medium of the invention are used for calibration by combining geological geophysical methods, and the principle is that the geology and the earthquake can be perfectly combined under the correct condition, and the geology and the earthquake can be mutually checked. On the basis of conventional calibration, well seismic breakpoint inspection and fault distance matching after time-depth calibration are added, a complete quality control basis is provided, the result can meet the requirement of a geophysical operator for constructing and explaining reservoir prediction, and a subsequent result can provide an oil reservoir operator for analyzing the oil reservoir result. The calibration result is obtained and used for construction interpretation, reservoir prediction and oil reservoir analysis, so that the accuracy of the result can be improved, and a more accurate basis is provided for oil field deployment.
It will be understood by those skilled in the art that all or part of the steps in the method for implementing the above embodiments may be implemented by relevant hardware instructed by a program, and the program may be stored in a computer readable storage medium, such as ROM/RAM, magnetic disk, optical disk, etc.
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 comprehensive time depth calibration method is characterized by comprising the following steps:
comparing the stratums among wells according to the logging curves of all wells to obtain the breakpoint position and the breakpoint distance of each well breakpoint;
inputting the well logging curve, the breakpoint position and the break distance into an interpretation system to obtain a synthetic seismic record;
carrying out well seismic calibration according to the wave resistance characteristics and the synthetic seismic record to obtain a time depth calibration result;
inputting the time-depth calibration result, the breakpoint position and the fault distance into the seismic section, determining the goodness of fit of the breakpoints of all wells and the seismic breakpoints, and determining the goodness of fit of the fault distance of the breakpoints of all wells and the fault distance of the seismic upper and lower plates;
and adjusting the synthetic seismic record according to the wave resistance characteristics until the coincidence degree of each well breakpoint and the seismic breakpoint and the coincidence degree of the fracture distance of the seismic upper and lower discs and the fracture distance of each well breakpoint reach preset values.
2. The method of claim 1, wherein the cross-well stratigraphic comparison according to the log of each well to obtain the breakpoint position and the break distance of each well breakpoint comprises:
acquiring the logging curve according to the logging data of each well;
comparing the horizons and faults among multiple wells according to the logging curve, and determining a full section of the region;
and comprehensively analyzing the stratum according to the region full-section, and determining the breakpoint position and the breakpoint distance of each well breakpoint.
3. The method of claim 1, wherein said adjusting the synthetic seismic records according to the wave impedance characteristics comprises: and adjusting the position of the synthetic seismic record up and down according to the wave resistance characteristics to obtain an adjusted time-depth calibration result.
4. The method of claim 3, further comprising: and determining the goodness of fit between each well breakpoint and the earthquake breakpoint and the goodness of fit between the fault distance of the earthquake upper and lower plates and the fault distance of each well breakpoint after adjustment according to the adjusted time-depth calibration result.
5. An integrated time depth calibration device, comprising:
the fracture data unit is used for comparing the stratums among the wells according to the logging curves of all the wells to obtain the fracture positions and the fracture distances of the well breakpoints;
the synthetic seismic unit is used for inputting the well logging curve, the breakpoint position and the break distance into an interpretation system to obtain a synthetic seismic record;
the time-depth calibration unit is used for carrying out well seismic calibration according to the wave resistance characteristics and the synthetic seismic record to obtain a time-depth calibration result;
the goodness of fit determining unit is used for inputting the time-depth calibration result, the breakpoint position and the break distance into the seismic section, determining the goodness of fit of each well breakpoint and the seismic breakpoint, and determining the goodness of fit of the break distance of each well breakpoint and the break distance of the seismic upper and lower plates;
and the goodness of fit adjusting unit is used for adjusting the synthetic seismic record according to the wave resistance characteristics until the goodness of fit of the well breakpoints and the seismic breakpoints and the goodness of fit of the fracture distance of the seismic upper and lower discs and the fracture distance of the well breakpoints reach preset values.
6. The apparatus of claim 5, wherein the breakpoint data unit comprises:
the logging curve module is used for acquiring the logging curve according to logging data of each well;
the regional full-profile module is used for comparing horizons and faults among multiple wells according to the logging curve and determining a regional full profile;
and the breakpoint and fault distance module is used for comprehensively analyzing the stratum according to the full section of the region and determining the breakpoint position and fault distance of each well breakpoint.
7. The apparatus according to claim 5, wherein the goodness-of-fit adjustment unit includes: and the time depth adjusting module is used for adjusting the position of the synthetic seismic record up and down according to the wave resistance characteristics to obtain an adjusted time depth calibration result.
8. The apparatus according to claim 7, wherein the goodness-of-fit adjustment unit further comprises: and the goodness of fit adjusting module is used for determining the goodness of fit of each adjusted well breakpoint and earthquake breakpoint and the goodness of fit of the fault distance of the earthquake upper and lower plates and the fault distance of each well breakpoint according to the adjusted time-depth calibration result.
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 when executing the computer program implements the steps of:
comparing the stratums among wells according to the logging curves of all wells to obtain the breakpoint position and the breakpoint distance of each well breakpoint;
inputting the well logging curve, the breakpoint position and the break distance into an interpretation system to obtain a synthetic seismic record;
carrying out well seismic calibration according to the wave resistance characteristics and the synthetic seismic record to obtain a time depth calibration result;
inputting the time-depth calibration result, the breakpoint position and the fault distance into the seismic section, determining the goodness of fit of the breakpoints of all wells and the seismic breakpoints, and determining the goodness of fit of the fault distance of the breakpoints of all wells and the fault distance of the seismic upper and lower plates;
and adjusting the synthetic seismic record according to the wave resistance characteristics until the coincidence degree of each well breakpoint and the seismic breakpoint and the coincidence degree of the fracture distance of the seismic upper and lower discs and the fracture distance of each well breakpoint reach preset values.
10. A computer-readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of:
comparing the stratums among wells according to the logging curves of all wells to obtain the breakpoint position and the breakpoint distance of each well breakpoint;
inputting the well logging curve, the breakpoint position and the break distance into an interpretation system to obtain a synthetic seismic record;
carrying out well seismic calibration according to the wave resistance characteristics and the synthetic seismic record to obtain a time depth calibration result;
inputting the time-depth calibration result, the breakpoint position and the fault distance into the seismic section, determining the goodness of fit of the breakpoints of all wells and the seismic breakpoints, and determining the goodness of fit of the fault distance of the breakpoints of all wells and the fault distance of the seismic upper and lower plates;
and adjusting the synthetic seismic record according to the wave resistance characteristics until the coincidence degree of each well breakpoint and the seismic breakpoint and the coincidence degree of the fracture distance of the seismic upper and lower discs and the fracture distance of each well breakpoint reach preset values.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910057366.5A CN111458751B (en) | 2019-01-22 | 2019-01-22 | Comprehensive time depth calibration method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910057366.5A CN111458751B (en) | 2019-01-22 | 2019-01-22 | Comprehensive time depth calibration method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111458751A true CN111458751A (en) | 2020-07-28 |
| CN111458751B CN111458751B (en) | 2023-02-28 |
Family
ID=71683058
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910057366.5A Active CN111458751B (en) | 2019-01-22 | 2019-01-22 | Comprehensive time depth calibration method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111458751B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112147717A (en) * | 2020-09-27 | 2020-12-29 | 北京中恒利华石油技术研究所 | Pre-stack and post-stack well-seismic combination time depth calibration method |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2191414C1 (en) * | 2001-12-28 | 2002-10-20 | Копилевич Ефим Абрамович | Method of geophysical prospecting for low-amplitude tectonic diastrophism in oil- and gas-efficient rocks |
| CN102121995A (en) * | 2010-12-03 | 2011-07-13 | 中国石油天然气集团公司 | Seismic inversion reservoir stratum forecasting method for overthrust-containing complex structure |
| WO2015053659A1 (en) * | 2013-10-10 | 2015-04-16 | Светлана Викторовна ШОЛОХОВА | Method of producing an a priori hodograph for carrying out lithostratigraphic correlation |
| CN105626054A (en) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | Method for subtly depicting broken edges of thick-layer fault block oil reservoir |
| CN105626053A (en) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | Method for conducting fault zone research by drilling fault zone through cored well |
| CN105759312A (en) * | 2016-03-01 | 2016-07-13 | 中国海洋石油总公司 | Reservoir characteristic method well-seismic calibration method |
| US20170075007A1 (en) * | 2014-06-04 | 2017-03-16 | Halliburton Energy Services, Inc. | Analyzing geomechanical properties of subterranean rock based on seismic data |
-
2019
- 2019-01-22 CN CN201910057366.5A patent/CN111458751B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2191414C1 (en) * | 2001-12-28 | 2002-10-20 | Копилевич Ефим Абрамович | Method of geophysical prospecting for low-amplitude tectonic diastrophism in oil- and gas-efficient rocks |
| CN102121995A (en) * | 2010-12-03 | 2011-07-13 | 中国石油天然气集团公司 | Seismic inversion reservoir stratum forecasting method for overthrust-containing complex structure |
| WO2015053659A1 (en) * | 2013-10-10 | 2015-04-16 | Светлана Викторовна ШОЛОХОВА | Method of producing an a priori hodograph for carrying out lithostratigraphic correlation |
| US20170075007A1 (en) * | 2014-06-04 | 2017-03-16 | Halliburton Energy Services, Inc. | Analyzing geomechanical properties of subterranean rock based on seismic data |
| CN105626054A (en) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | Method for subtly depicting broken edges of thick-layer fault block oil reservoir |
| CN105626053A (en) * | 2014-11-07 | 2016-06-01 | 中国石油化工股份有限公司 | Method for conducting fault zone research by drilling fault zone through cored well |
| CN105759312A (en) * | 2016-03-01 | 2016-07-13 | 中国海洋石油总公司 | Reservoir characteristic method well-seismic calibration method |
Non-Patent Citations (4)
| Title |
|---|
| DAJUN SUN等: "Array geometry calibration for underwater compact arrays", 《APPLIED ACOUSTICS》 * |
| 姜宏章等: "复杂断块油藏的精细构造解释方法", 《大庆石油地质与开发》 * |
| 王东明: "井震联合技术在断层精细解释中的应用", 《石油天然气学报 （江汉石油学院学报）》 * |
| 甘利灯等: "测井-地震-油藏模拟一体化技术及其在老油田挖潜中的应用", 《石油物探》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112147717A (en) * | 2020-09-27 | 2020-12-29 | 北京中恒利华石油技术研究所 | Pre-stack and post-stack well-seismic combination time depth calibration method |
| CN112147717B (en) * | 2020-09-27 | 2021-11-16 | 北京中恒利华石油技术研究所 | Pre-stack and post-stack well-seismic combination time depth calibration method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111458751B (en) | 2023-02-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9097821B2 (en) | Integrated workflow or method for petrophysical rock typing in carbonates | |
| CN100557464C (en) | Seismic prospecting horizon calibration method based on the prestack wave-field simulation | |
| CN108897066B (en) | Carbonate rock crack density quantitative prediction method and device | |
| US20150088424A1 (en) | Identifying geological formation depth structure using well log data | |
| EP2856217A2 (en) | Methods for generating depofacies classifications for subsurface oil or gas reservoirs or fields | |
| Cooper et al. | Quantitative biostratigraphy of the Taranaki Basin, New Zealand: A deterministic and probabilistic approach | |
| US10024991B2 (en) | Method and device for processing seismic signals | |
| CN105697002A (en) | Method for recognizing coal measure strata lithology | |
| Hirsche et al. | Avoiding pitfalls in geostatistical reservoir characterization: A survival guide | |
| AU2014386285B2 (en) | Validation of depth-depth curves using time-seismic depth interval velocity | |
| CN111458751B (en) | Comprehensive time depth calibration method and device | |
| CN114609675A (en) | Quantitative recovery method for carbonate rock stratum sedimentary micro-landform based on high-frequency cycle | |
| CN108254785B (en) | Reservoir determination method and device | |
| Kocon et al. | Quality assessment of microseismic event locations and traveltime picks using a multiplet analysis | |
| CN109630089B (en) | Horizontal well geological structure recognition method and device | |
| US6084825A (en) | Apparatus and methods for seismic data processing | |
| CN112528106A (en) | Volcanic lithology identification method | |
| Hogarth et al. | Controlled-source velocity calibration for real-time downhole microseismic monitoring | |
| Saberi | Rock-physics-assisted well-tie analysis for structural interpretation and seismic inversion | |
| CN112859172A (en) | Longitudinal and transverse wave micro-logging data processing method and device | |
| Walker et al. | Stochastic inversion for facies: A case study on the Schiehallion field | |
| CN112180443B (en) | Shale gas two-dimensional seismic dessert area optimization method and device | |
| CN112394392B (en) | Method and device for evaluating distribution condition of hydrocarbon source rock | |
| CN108260359A (en) | A kind of seismic data processing methods, devices and systems | |
| WO2018075046A1 (en) | Correcting biases in microseismic-event data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |