CN113356834A - Borehole acoustic time difference correction method, correction device and electronic equipment - Google Patents
Borehole acoustic time difference correction method, correction device and electronic equipment Download PDFInfo
- Publication number
- CN113356834A CN113356834A CN202110797553.4A CN202110797553A CN113356834A CN 113356834 A CN113356834 A CN 113356834A CN 202110797553 A CN202110797553 A CN 202110797553A CN 113356834 A CN113356834 A CN 113356834A
- Authority
- CN
- China
- Prior art keywords
- borehole
- acoustic
- time difference
- invasion
- correction value
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 181
- 238000000034 method Methods 0.000 title claims abstract description 38
- 230000009545 invasion Effects 0.000 claims abstract description 132
- 239000012535 impurity Substances 0.000 claims abstract description 55
- 230000015572 biosynthetic process Effects 0.000 claims description 42
- 230000001902 propagating effect Effects 0.000 claims description 3
- 239000000356 contaminant Substances 0.000 claims description 2
- 238000005755 formation reaction Methods 0.000 description 39
- 230000008901 benefit Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The embodiment of the application discloses a method and a device for correcting borehole acoustic time difference and electronic equipment. The method for correcting the borehole acoustic time difference comprises the following steps: acquiring the invasion radius of impurities in the borehole and the theoretical value of the borehole sound wave time difference; establishing an invasion correction model based on the invasion radius of impurities in the borehole and the theoretical value of the borehole sound wave time difference; acquiring a slope correction value of the borehole acoustic time difference; and inputting the inclination correction value of the borehole acoustic time difference into an invasion correction model to obtain the corrected borehole acoustic time difference. According to the method and the device, on the basis of correcting the influence of the inclination angle of the stratum interface on the well sound wave, the influence of invasion of impurities in the well on the well sound wave time difference is further corrected, the accuracy of the corrected well sound wave time difference is ensured, the geological conditions around the well can be truly reflected by the well sound wave time difference, and the correction efficiency of the well sound wave time difference can be improved.
Description
Technical Field
The embodiment of the application relates to the technical field of borehole sound wave time difference correction, in particular to a borehole sound wave time difference correction method, a borehole sound wave time difference correction device and electronic equipment.
Background
Shale gas has gradually changed from a potential resource to a clean energy resource which is being developed and utilized on a large scale, and the shale gas becomes an effective substitute for the conventional oil and gas resources. The longitudinal distribution, the anisotropic change, the structure and the like of the shale in different areas are different, and how to more accurately acquire parameters such as the anisotropy of the shale in different areas becomes more important.
In the related technology, the acoustic time difference is obtained in an acoustic logging mode and is used as basic data of the following processes of well seismic calibration, shale physical modeling, fracturing optimization design and the like. However, the acoustic moveout is affected by the surrounding environment of the borehole, such as invasion of mud in the borehole, the inclination angle of the stratum, and other factors, and the acoustic moveout is deviated, thereby affecting the accuracy of the geological condition of shale around the borehole obtained by the acoustic moveout.
Disclosure of Invention
In order to solve at least one of the above technical problems, embodiments of the present application provide a method for correcting a borehole acoustic time difference, a device for correcting a borehole acoustic time difference, and an electronic apparatus.
In a first aspect, an embodiment of the present application provides a method for correcting a borehole acoustic time difference, including: acquiring the invasion radius of impurities in the borehole and the theoretical value of the borehole sound wave time difference; establishing an invasion correction model based on the invasion radius of impurities in the borehole and the theoretical value of the borehole sound wave time difference; acquiring a slope correction value of the borehole acoustic time difference; and inputting the inclination correction value of the borehole acoustic time difference into an invasion correction model to obtain the corrected borehole acoustic time difference.
In one possible embodiment, the step of establishing an invasion correction model based on the invasion radius of the impurity in the borehole and the theoretical value of the borehole acoustic time difference comprises: acquiring an invasion acoustic wave time difference based on the invasion radius of impurities in the borehole; and establishing an invasion correction model based on the theoretical values of the invasion acoustic wave time difference and the borehole acoustic wave time difference.
In one possible embodiment, the step of obtaining the difference in the invasion acoustic wave time based on the invasion radius of the impurity in the borehole includes: according to the formula
Acquiring an invasion sound wave time difference; wherein y is the difference in time of invasion sound waves, x is the invasion radius of impurities in the borehole, y0Is constant, A is constant, w is constant, and xc is constant.
In one possible embodiment, the step of establishing an invasion correction model based on the theoretical values of the invasion acoustic time difference and the borehole acoustic time difference comprises: calculating the ratio of the theoretical values of the time difference of the invasion sound waves and the time difference of the borehole sound waves; and establishing an invasion correction model according to the ratio and the invasion radius of the impurities in the borehole.
In one possible embodiment, the step of obtaining a slope correction value for the borehole acoustic moveout comprises: and acquiring a slope correction value of the acoustic time difference of the borehole according to at least one of the radius of the borehole, the position of the acoustic source, the position of the acoustic receiver, the propagation speed of the acoustic wave in the impurity, the propagation speed of the acoustic wave in the stratum and the inclination angle of the stratum.
In one possible embodiment, the step of obtaining a slope correction value for the borehole acoustic moveout based on at least one of a radius of the borehole, a location of the acoustic source, a location of the acoustic receiver, a speed of sound propagation through the contaminant, a speed of sound propagation through the formation, and a dip angle of the formation comprises: according to the formula
Acquiring a slope correction value of the borehole acoustic time difference; it is composed ofWhere t is the slope correction for the borehole acoustic time difference, r is the radius of the borehole, L is the linear distance between the acoustic source and the acoustic receiver, b is the linear distance between the acoustic receiver and the intersection of the borehole and the formation, vfV is the propagation velocity of the acoustic wave in the impurity, C is a variable related to the angle of inclination of the formation, D is a variable related to the angle of inclination of the formation, and E is a variable related to the angle of inclination of the formation.
In one possible embodiment, the sonic receivers include a first sonic receiver and a second sonic receiver, and the step of obtaining a slope correction value for the borehole sonic moveout comprises: acquiring a first inclination correction value through the position of the first sound wave receiver; acquiring a second slope correction value according to the position of the second sound wave receiver; a slope correction value for the borehole acoustic time difference is obtained based on the first slope correction value and the second slope correction value.
In one possible embodiment, the step of obtaining a slope correction value for the borehole acoustic moveout based on the first slope correction value and the second slope correction value comprises: according to the formula
Acquiring a slope correction value of the borehole acoustic time difference; where T is the slope correction for the borehole acoustic time difference, L1Is the linear distance, L, of the first acoustic receiver from the acoustic source2Is the linear distance, t, of the second acoustic receiver from the acoustic source1Is a first slope correction value, t2Is a second slope correction value.
In a second aspect, embodiments of the present application provide a device for correcting borehole acoustic time differences, comprising: the acquisition device is used for acquiring the invasion radius of impurities in the borehole and the theoretical value of the borehole acoustic time difference; the model establishing device establishes an invasion correction model based on the invasion radius of impurities in the borehole and the theoretical value of the borehole sound wave time difference; the inclination correction device is used for acquiring inclination correction values of the borehole acoustic time difference; and the input device is used for inputting the inclination correction value of the borehole acoustic time difference into the invasion correction model to obtain the corrected borehole acoustic time difference.
In a third aspect, embodiments of the present application provide an electronic device, which includes a processor, a memory, and a program or instructions stored in the memory and executable on the processor, where the program or instructions, when executed by the processor, implement the steps of the method for correcting borehole acoustic time difference according to the first aspect.
The embodiment of the application has the following beneficial effects:
the inclination correction value of the borehole sound wave time difference is corrected through the invasion correction model, namely on the basis of correcting the influence of the inclination angle of the stratum on the borehole sound wave time difference, the influence of invasion of impurities in the borehole on the borehole sound wave time difference is further corrected, the accuracy of the corrected borehole sound wave time difference is ensured, the borehole sound wave time difference can truly reflect the geological conditions around the borehole, and reliable data are provided for subsequent geological modeling, fracturing optimization and the like.
In addition, the inclination correction value of the borehole acoustic wave time difference is corrected by establishing an invasion correction model, so that the correction efficiency of the borehole acoustic wave time difference is further improved.
Additional aspects and advantages of embodiments of the present application will be set forth in part in the description which follows and, in part, will be obvious from the description, or may be learned by practice of the present application.
Drawings
Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the application. Also, like reference numerals are used to refer to like parts throughout the drawings. In the drawings:
FIG. 1 is a flow chart illustrating one embodiment of a method for correcting borehole acoustic moveout according to the present disclosure;
FIG. 2 is a second flowchart illustrating the steps of a method for correcting for borehole acoustic moveout according to an embodiment of the present disclosure;
FIG. 3 is a schematic diagram illustrating the time difference of the arrival sound wave according to an embodiment of the present disclosure;
FIG. 4 is a third flowchart illustrating the steps of a method for correcting borehole acoustic moveout according to an embodiment of the present disclosure;
FIG. 5 is a graph illustrating an intrusion correction model according to an embodiment of the present disclosure;
FIG. 6 is a schematic illustration of a borehole acoustic time difference measurement according to an embodiment provided herein;
FIG. 7 is a flow chart illustrating the steps of a method for correcting borehole acoustic moveout according to an embodiment of the present disclosure;
FIG. 8 is a schematic block diagram illustrating an exemplary embodiment of a borehole acoustic time difference correction apparatus;
FIG. 9 is a waveform illustrating a corrected borehole acoustic time difference according to an embodiment of the present disclosure.
Wherein, the correspondence between the reference numbers and the part names in fig. 6 and 8 is:
t: sound source, R1: first acoustic receiver, R2: second acoustic receiver, 100: device for correcting borehole acoustic time difference, 110: acquisition means, 120: model building apparatus, 130: slope correction device, 140: an input device.
Detailed Description
In order that the above objects, features and advantages of the present application can be more clearly understood, the present application will be described in further detail with reference to the accompanying drawings and detailed description. It should be noted that the embodiments and features of the embodiments of the present application may be combined with each other without conflict.
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present application, however, the present application may be practiced in other ways than those described herein, and therefore the scope of the present application is not limited by the specific embodiments disclosed below.
In a first aspect, as shown in fig. 1, an embodiment of the present application provides a method for correcting a borehole acoustic time difference, including:
s101, acquiring the invasion radius of impurities in a borehole and the theoretical value of borehole sound wave time difference;
s102, establishing an invasion correction model based on the invasion radius of impurities in the borehole and the theoretical value of the borehole sound wave time difference;
step S103, obtaining a slope correction value of the borehole sound wave time difference;
and step S104, inputting the inclination correction value of the borehole acoustic wave time difference into an invasion correction model to obtain the corrected borehole acoustic wave time difference.
It will be appreciated that the acoustic moveout of the borehole is the difference in time between the acoustic source location and the traveling longitudinal wave propagating along the borehole wall at the acoustic receiver. And acquiring the invasion radius of the impurities in the well hole and the theoretical value of the sound wave time difference of the well hole, wherein the impurities in the well hole can be mud understandably. The theoretical value of the sound wave time difference of the borehole is the sound wave time difference of the borehole under the condition of not being influenced by environmental factors.
In some examples, the invasion radius of the impurities in the borehole and the theoretical value of the borehole acoustic time difference can be obtained by consulting data, performing field investigation, calculating or simulating, and the like.
And establishing an invasion correction model based on the invasion radius of impurities in the borehole and the theoretical value of the borehole sound wave time difference. It will be appreciated that the invasion correction model is used to correct for the effects of impurities in the borehole on borehole acoustic moveout.
In some examples, different invasion correction models can be established according to different impurity species invading into the well bore, and the accuracy and reliability of the invasion correction models for correcting the well bore acoustic wave time difference are improved.
A slope correction of the borehole acoustic time difference is obtained, which is understood to be the borehole acoustic time difference after correcting for the effects of the angle of inclination of the formation. In some examples, slope corrections for borehole acoustic moveouts may be obtained by reviewing data, calculations, field surveys, or simulation.
In some examples, the slope correction value of the borehole acoustic moveout varies with changes in the angle of inclination of the formation.
The inclination correction value of the borehole sound wave time difference is input into the invasion correction model, so that the inclination correction value of the borehole sound wave time difference is corrected through the invasion correction model, namely, on the basis of correcting the influence of the formation inclination angle on the borehole sound wave time difference, the influence of invasion of impurities in the borehole on the borehole sound wave time difference is further corrected, the accuracy of the corrected borehole sound wave time difference is ensured, the borehole sound wave time difference can truly reflect the geological conditions around the borehole, and reliable data are provided for subsequent geological modeling, fracturing optimization and the like. In addition, the inclination correction value of the borehole acoustic wave time difference is corrected by establishing an invasion correction model, so that the correction efficiency of the borehole acoustic wave time difference is further improved.
In some examples, as shown in FIG. 2, the step of establishing an invasion correction model based on the invasion radius of the impurity in the borehole and the theoretical value of the borehole acoustic time difference comprises:
step S201, acquiring invasion sound wave time difference based on invasion radius of impurities in a borehole;
and S202, establishing an invasion correction model based on the theoretical values of the invasion acoustic wave time difference and the borehole acoustic wave time difference.
It is understood that the invasion acoustic moveout is the acoustic moveout of the borehole under the influence of different invasion radii. In some examples, the invasion acoustic moveout of the borehole at different invasion radii can be obtained by consulting data, performing field exploration, calculating or simulating, etc.
The method comprises the steps of establishing an invasion correction model based on the theoretical values of the invasion acoustic wave time difference and the borehole acoustic wave time difference, further improving the accuracy of the invasion correction model, improving the correction effect of the invasion correction model on the inclination correction value of the borehole acoustic wave time difference, ensuring the accuracy of the corrected borehole acoustic wave time difference, enabling the corrected borehole acoustic wave time difference to accurately reflect the geological parameters of the borehole, and improving the reliability of the correction method of the borehole acoustic wave time difference.
In some examples, the step of obtaining the difference in the invasion acoustic wave time based on the invasion radius of the impurity in the borehole comprises:
according to the formula
Acquiring an invasion sound wave time difference;
wherein y is the difference in time of invasion sound waves, x is the invasion radius of impurities in the borehole, y0Is constant, A is constant, w is constant, and xc is constant.
It can be understood that, as shown in fig. 3, there is a corresponding variation in the difference in the invasion acoustic wave time of the borehole with the invasion radius of the foreign matter in the borehole. And fitting the time difference of the invasion sound waves of the borehole under different invasion radiuses, so as to obtain a variation curve of the time difference of the invasion sound waves along with the invasion radiuses, and further obtain a calculation formula of the time difference of the invasion sound waves.
In particular, the formula
Wherein y is the difference in time of invasion sound waves, x is the invasion radius of impurities in the borehole, and y0A, w and xc are constants.
The time difference y of the invasion sound wave can be calculated through a fitting formula, and the accuracy of the time difference of the invasion sound wave is ensured, so that the accuracy of the established invasion correction model is improved, and the reliability of the correction method of the well sound wave time difference is improved.
It will be appreciated that y may be set according to different circumstances of the wellbore0A, w and xc. In some examples, y0Is 260.827, xc is 0.01, w is 0.94, and a is 1.815.
In some examples, as shown in FIG. 4, the step of building an invasion correction model based on the theoretical values of the invasion acoustic time difference and the borehole acoustic time difference comprises:
step S301, calculating the ratio of the theoretical values of the time difference of the invasion sound waves and the time difference of the borehole sound waves;
and step S302, establishing an invasion correction model according to the ratio and the invasion radius of impurities in the borehole.
It can be understood that different geological parameters correspond to different theoretical values of the borehole acoustic time difference, so that the ratio of the time difference of the invasion acoustic wave to the theoretical value is different under different geological parameters. As shown in fig. 5, different curves represent the variation relationship between the ratio and the invasion radius under different geological parameters, and the applicability of the established invasion model is improved.
And calculating the ratio of the theoretical values of the time difference of the invasion sound waves and the time difference of the well sound waves, and establishing an invasion correction model according to the ratio of the time difference of the invasion sound waves and the time difference of the well sound waves and the invasion radius, so that the accuracy of the invasion model is improved, the accuracy of the corrected sound wave time difference is improved, the sound wave time difference can reflect the geological conditions more accurately, and the reliability of the correction method of the well sound wave time difference is improved. Meanwhile, different invasion models can be established according to different geological parameters, the applicability of the invasion models is improved, the method is simple and easy to implement, and the efficiency of the correction method of the borehole sound waves is improved.
In some examples, the step of obtaining a slope correction value for the borehole acoustic moveout comprises: and acquiring a slope correction value of the acoustic time difference of the borehole according to at least one of the radius of the borehole, the position of the acoustic source, the position of the acoustic receiver, the propagation speed of the acoustic wave in the impurity, the propagation speed of the acoustic wave in the stratum and the inclination angle of the stratum.
It will be appreciated that when measuring the acoustic moveout of the borehole, the radius of the borehole, the location of the acoustic source, the location of the acoustic receiver, the speed of propagation of the acoustic wave and the angle of inclination of the formation all contribute to the acoustic moveout of the borehole. Therefore, the borehole sound wave time difference is corrected according to at least one of the radius of the borehole, the position of the sound source, the position of the sound wave receiver, the propagation speed of the sound wave in the impurities, the propagation speed of the sound wave in the stratum and the inclination angle of the stratum, the inclination correction value of the borehole sound wave time difference is obtained, the influence of the inclination angle of the stratum on the borehole sound wave time difference is reduced, the accuracy of the obtained inclination correction value of the borehole sound wave time difference is improved, and the accuracy of the correction method of the borehole sound wave time difference is improved.
In some examples, there are a variety of different geologies around the borehole, and the propagation velocities of sound waves in different formations are not the same. The borehole acoustic time difference is corrected through the propagation speed of the acoustic wave in different stratums, the accuracy of the inclination correction value of the borehole acoustic time difference can be further improved, and the influence of the inclination angle of the stratum on the borehole acoustic time difference is reduced.
In some examples, the step of obtaining a slope correction value for the borehole acoustic moveout based on at least one of a radius of the borehole, a location of the acoustic source, a location of the acoustic receiver, a speed of propagation of the acoustic wave in the impurity, a speed of propagation of the acoustic wave in the formation, and a dip angle of the formation comprises:
according to the formula
Acquiring a slope correction value of the borehole acoustic time difference;
wherein t is a slope correction value of the borehole acoustic time difference, r is a radius of the borehole, L is a linear distance between the acoustic source and the acoustic receiver, b is a linear distance between the acoustic receiver and an intersection of the borehole and the formation, vfV is the propagation velocity of the acoustic wave in the impurity, C is a variable related to the angle of inclination of the formation, D is a variable related to the angle of inclination of the formation, and E is a variable related to the angle of inclination of the formation.
When calculating the slope correction value of the borehole acoustic time difference, r is the radius of the borehole, L is the linear distance between the acoustic source and the acoustic receiver, and b is the linear distance between the acoustic receiver and the intersection of the borehole and the formation. It will be appreciated that the intersection of the borehole and the formation may be the intersection between the central axis of the borehole and the formation. v. offThe propagation velocity of the acoustic wave in the impurity, which may be mud in some examples. v is the propagation velocity of the acoustic wave in the formation, which in some examples may be a number of different values due to the presence of a number of different geological environments around the borehole. C. D and E are variables related to the inclination angle of the formation.
Specifically, as shown in fig. 6, the included angle θ is the inclined angle of the formation interface, and after the sound source T emits the sound wave, the sound wave refracts and propagates between different formations, and is finally received by the first sound wave receiver R1. During the refraction and propagation of the sound wave, an included angle theta is formed1、θ2、θ3And theta4. The inclination correction value of the borehole acoustic time difference is calculated by the formula
Where t is the slope correction for the borehole acoustic time difference, r is the radius of the borehole, vfIs the propagation velocity of sound waves in impurities, L1Is the linear distance between the sound source T and the first sonic receiver R1, b1Is the linear distance between the first acoustic receiver R1 and the intersection of the borehole and the formation. v. of1And v2Is the propagation velocity of sound waves in different formations.
By the formula
And calculating and acquiring a slope correction value of the borehole sound wave time difference, and further improving the accuracy of the slope correction value of the borehole sound wave time difference, so that the accuracy and the reliability of the correction method of the borehole sound wave time difference are improved.
In some examples, the number of the acoustic receivers can be multiple, and the multiple acoustic receivers compensate each other for correction, so that the accuracy of the slope correction value of the acoustic time difference of the borehole is further improved.
In some examples, as shown in FIG. 7, the sonic receivers include a first sonic receiver and a second sonic receiver, and the step of obtaining a slope correction value for the borehole sonic moveout includes:
step S401, acquiring a first inclination correction value according to the position of the first sound wave receiver;
step S402, acquiring a second inclination correction value according to the position of the second sound wave receiver;
in step S403, a slope correction value of the borehole acoustic time difference is obtained based on the first slope correction value and the second slope correction value.
It is to be understood that, as shown in fig. 6, the acoustic wave receiver includes a first acoustic wave receiver R1 and a second acoustic wave receiver R2. After the sound source T emits sound waves, the sound waves refractively propagate between different strata and are finally received by the first sound receiver R1 and the second sound receiver R2. During the refraction and propagation of the sound wave, an included angle theta is formed1、θ2、θ3、θ4、θ1’、θ2’、θ3' and theta4’。
In particular, it can be based on a formula
Calculating a first slope correction value, where t1Is a first slope correction value, r is the radius of the borehole, vfIs the propagation velocity of sound waves in impurities, L1Is the linear distance between the sound source T and the first sonic receiver R1, b1Is the linear distance between the acoustic receiver and the intersection of the borehole and the formation. v. of1And v2Is the propagation velocity of sound waves in different formations.
According to the formula
Calculating a second slope correction value, where t2For the second slope correction value, r is the radius of the borehole, vfIs the propagation velocity of sound waves in impurities, L2Is the linear distance between the sound source T and the second sound receiver R2, b2Is the linear distance between the second acoustic receiver R2 and the intersection of the borehole and the formation. v. of1And v2Is the propagation velocity of sound waves in different formations.
And a first inclination correction value is acquired through the position of the first sound wave receiver R1, a second inclination correction value is acquired through the position of the second sound wave receiver R2, correction compensation can be carried out between the second inclination correction value and the first inclination correction value, the accuracy of the acquired inclination correction value of the borehole sound wave time difference is improved, and the reliability of the correction method of the borehole sound wave time difference is further improved.
In some examples, the step of obtaining a slope correction value for the borehole acoustic moveout based on the first slope correction value and the second slope correction value comprises:
according to the formula
Acquiring a slope correction value of the borehole acoustic time difference;
where T is the slope correction for the borehole acoustic time difference, L1Is the linear distance, L, of the first acoustic receiver from the acoustic source2Is the linear distance, t, of the second acoustic receiver from the acoustic source1Is a first slope correction value, t2Is a second slope correction value.
By the formula
And calculating a slope correction value of the borehole sound wave time difference, so that correction compensation can be performed between the first slope correction value and the second slope correction value, the accuracy of the obtained slope correction value of the borehole sound wave time difference is improved, the influence of the formation inclination angle on the sound wave time difference is reduced, and the reliability of the correction method of the borehole sound wave time difference is improved.
In a second aspect, as shown in fig. 8, based on the same inventive concept of the above method, the embodiment of the present application provides a device 100 for correcting borehole acoustic wave time difference, which includes an obtaining device 110, a model building device 120, a slope correcting device 130 and an input device 140. The acquisition device 110 is used for acquiring the invasion radius of impurities in the borehole and the theoretical value of the sound wave time difference of the borehole. The model building means 120 builds an invasion correction model based on the invasion radius of the impurities in the borehole and the theoretical value of the borehole acoustic time difference. The slope correction device 130 is used to obtain slope correction values for the borehole acoustic moveout. The input device 140 is configured to input the slope correction value of the borehole acoustic time difference into the intrusion correction model to obtain a corrected borehole acoustic time difference. The device for correcting borehole acoustic wave time difference in the embodiment of the present application can implement the steps of the method for correcting borehole acoustic wave time difference in the first aspect, and therefore has all the beneficial effects of the first aspect, and is not described herein again.
In a third aspect, based on the same inventive concept of the above method, an embodiment of the present application provides an electronic device, which includes a processor, a memory, and a program or instructions stored in the memory and executable on the processor, where the program or instructions, when executed by the processor, implement the steps of the method for correcting borehole acoustic time difference according to the first aspect, and therefore, all the advantages of the first aspect are achieved, and details are not repeated herein.
In one embodiment of the present application, the borehole acoustic time difference is corrected by first correcting for the effect of the formation dip angle on the borehole acoustic time difference, as shown in FIG. 6.
The sound wave emitted from the sound source T is refracted by different formation media, and then received by the first sound receiver R1 and the second sound receiver R2. The sound wave forms an included angle theta in the refraction process1、θ2、θ3、θ4、θ1’、θ2’、θ3' and theta4’。
From the set position of the first acoustic wave receiver R1, a first inclination correction value can be calculated
Wherein t is1Is a first slope correction value, r is the radius of the borehole, vfIs the propagation velocity of sound waves in impurities, L1Is the linear distance between the sound source T and the first sonic receiver R1, b1Is the linear distance between the acoustic receiver and the intersection of the borehole and the formation. v. of1And v2Is the propagation velocity of sound waves in different formations.
From the set position of the second sound wave receiver R2, a second slope correction value can be calculated
Wherein t is2For the second slope correction value, r is the radius of the borehole, vfIs the propagation velocity of sound waves in impurities, L2Is the linear distance between the sound source T and the second sound receiver R2, b2Is the linear distance between the second acoustic receiver R2 and the intersection of the borehole and the formation. v. of1And v2Is the propagation velocity of sound waves in different formations.
According to the formula
A slope correction of the borehole acoustic moveout may be calculated, i.e., the inclination angle of the formation corrected.
Where T is the slope correction for the borehole acoustic time difference, L1Is the linear distance, L, of the first acoustic receiver from the acoustic source2Is the linear distance, t, of the second acoustic receiver from the acoustic source1Is a first slope correction value, t2Is a second slope correction value.
According to the formula
Calculating the difference of the invasion sound wave time, wherein y is the difference of the invasion sound wave time, x is the invasion radius of the impurities in the borehole, y0A, w and xc are constants. In some examples, y0Is 260.827, xc is 0.01, w is 0.94, and a is 1.815.
And calculating the ratio of the theoretical values of the time difference of the invasion sound waves and the time difference of the borehole sound waves, and establishing an invasion correction model according to the ratio and the invasion radius. As shown in fig. 5, the theoretical value of the borehole acoustic wave time difference is different under different geological parameters, so the ratio of the difference of the invasion acoustic wave time difference to the theoretical value of the borehole acoustic wave time difference is different.
And inputting the inclination correction value of the borehole acoustic wave time difference into an invasion correction model to further correct the influence of invasion radius on the borehole acoustic wave time difference. Specifically, a slope correction value of the borehole sound wave time difference can be taken as the invasion sound wave time difference and brought into an invasion correction model, a theoretical value of the borehole sound wave time difference under the corresponding invasion radius is searched, the borehole sound wave time difference is corrected, the corrected borehole sound wave time difference is obtained, and the influence of impurity invasion and formation inclination angle on the borehole sound wave time difference is reduced. As shown in fig. 9, by correcting, deviation of the borehole acoustic time difference can be reduced, accuracy of the borehole acoustic time difference is improved, and reliable data is provided for subsequent geological modeling, fracturing optimization and the like.
In the present invention, the terms "first", "second", and "third" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance; the term "plurality" means two or more unless expressly limited otherwise. The terms "mounted," "connected," "fixed," and the like are to be construed broadly, and for example, "connected" may be a fixed connection, a removable connection, or an integral connection; "coupled" may be direct or indirect through an intermediary. The specific meanings of the above terms in the present invention can be understood by those skilled in the art according to specific situations.
In the description of the present invention, it is to be understood that the terms "upper", "lower", "left", "right", "front", "rear", and the like indicate orientations or positional relationships based on those shown in the drawings, and are only for convenience of description and simplification of description, but do not indicate or imply that the referred device or unit must have a specific direction, be constructed in a specific orientation, and be operated, and thus, should not be construed as limiting the present invention.
In the description of the present application, the description of the terms "one embodiment," "some embodiments," "specific embodiments," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this application, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
The above is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and various modifications and changes will occur to those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Claims (10)
1. A method of correcting for borehole acoustic moveout, comprising:
acquiring the invasion radius of impurities in the borehole and the theoretical value of the borehole sound wave time difference;
establishing an invasion correction model based on the invasion radius of impurities in the borehole and the theoretical value of the borehole acoustic time difference;
acquiring a slope correction value of the borehole acoustic time difference;
and inputting the inclination correction value of the borehole acoustic time difference into the invasion correction model to obtain the corrected borehole acoustic time difference.
2. The method of correcting for borehole acoustic moveout of claim 1, wherein the step of building the invasion correction model based on the radius of invasion of the impurity in the borehole and the theoretical value of the borehole acoustic moveout of the borehole comprises:
acquiring an invasion acoustic wave time difference based on the invasion radius of the impurities in the borehole;
and establishing the invasion correction model based on the invasion sound wave time difference and the theoretical value of the borehole sound wave time difference.
3. The method of correcting for borehole acoustic moveout of claim 2, wherein the step of deriving the invasion acoustic moveout based on the invasion radius of the impurity in the borehole comprises:
according to the formula
Acquiring the invasion acoustic wave time difference;
wherein y is the invasion acoustic time difference, x is the invasion radius of the impurity in the borehole, y0Is constant, A is constant, w is constant, and xc is constant.
4. The method of correcting for borehole acoustic time differences according to claim 2, wherein said step of building said intrusion correction model based on said theoretical values of said intrusion acoustic time differences and said borehole acoustic time differences comprises:
calculating the ratio of the invasion acoustic time difference to the theoretical value of the borehole acoustic time difference;
and establishing the invasion correction model according to the ratio and the invasion radius of the impurities in the borehole.
5. The method of correcting for borehole acoustic moveout of claim 1, wherein the step of obtaining a slope correction for the borehole acoustic moveout comprises:
and acquiring a slope correction value of the borehole acoustic wave time difference according to at least one of the radius of the borehole, the position of the acoustic source, the position of the acoustic receiver, the propagation speed of the acoustic wave in the impurity, the propagation speed of the acoustic wave in the stratum and the inclination angle of the stratum.
6. The method of correcting for borehole acoustic moveout of claim 5, wherein said step of obtaining a slope correction for the borehole acoustic moveout based on at least one of a radius of the borehole, a location of the acoustic source, a location of the acoustic receiver, a velocity of the acoustic wave propagating in the contaminant, a velocity of the acoustic wave propagating in the formation, and a slope angle of the formation comprises:
according to the formula
Obtaining a slope correction value of the borehole acoustic time difference;
wherein t is a slope correction value of the borehole acoustic time difference, r is a radius of the borehole, L is a linear distance between the acoustic source and the acoustic receiver, b is a linear distance between the acoustic receiver and an intersection of the borehole and the formation, vfThe propagation speed of the sound wave in the impurities is represented as v, the propagation speed of the sound wave in the stratum is represented as C, the variable related to the inclination angle of the stratum is represented as D, and the variable related to the inclination angle of the stratum is represented as E.
7. The method of correcting for borehole acoustic moveout of claim 6, wherein the acoustic receivers comprise a first acoustic receiver and a second acoustic receiver, and wherein the step of obtaining slope correction values for the borehole acoustic moveout comprises:
acquiring a first inclination correction value according to the position of the first sound wave receiver;
acquiring a second slope correction value according to the position of the second sound wave receiver;
obtaining a slope correction value for the borehole acoustic time difference based on the first slope correction value and the second slope correction value.
8. The method of correcting for borehole acoustic moveout of claim 7, wherein the step of obtaining a slope correction value for the borehole acoustic moveout based on the first slope correction value and the second slope correction value comprises:
according to the formula
Obtaining a slope correction value of the borehole acoustic time difference;
wherein T is a slope correction value of the borehole acoustic time difference, L1Is the linear distance, L, of the first acoustic receiver from the acoustic source2Is the linear distance, t, of the second acoustic receiver from the acoustic source1For said first slope correction value, t2Is the second slope correction value.
9. A borehole acoustic time difference correction apparatus, comprising:
the acquisition device is used for acquiring the invasion radius of impurities in the borehole and the theoretical value of the borehole sound wave time difference;
the model establishing device is used for establishing the invasion correction model based on the invasion radius of impurities in the borehole and the theoretical value of the borehole acoustic time difference;
a slope correction device for obtaining a slope correction value for the borehole acoustic moveout;
and the input device is used for inputting the inclination correction value of the borehole acoustic time difference into the invasion correction model to obtain the corrected borehole acoustic time difference.
10. An electronic device comprising a processor, a memory and a program or instructions stored on the memory and executable on the processor, the program or instructions when executed by the processor implementing the steps of the method of correcting borehole acoustic time differences according to any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110797553.4A CN113356834B (en) | 2021-07-14 | 2021-07-14 | Borehole acoustic time difference correction method, correction device and electronic equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110797553.4A CN113356834B (en) | 2021-07-14 | 2021-07-14 | Borehole acoustic time difference correction method, correction device and electronic equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113356834A true CN113356834A (en) | 2021-09-07 |
| CN113356834B CN113356834B (en) | 2023-04-18 |
Family
ID=77539506
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110797553.4A Active CN113356834B (en) | 2021-07-14 | 2021-07-14 | Borehole acoustic time difference correction method, correction device and electronic equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113356834B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101787884A (en) * | 2010-01-28 | 2010-07-28 | 中国石油集团川庆钻探工程有限公司 | Reservoir fluid type discrimination method based on difference value of acoustic porosity and neutron porosity |
| US20150267527A1 (en) * | 2014-03-21 | 2015-09-24 | Conocophillips Company | Method for analysing pore pressure in shale formations |
| CN105089663A (en) * | 2015-07-31 | 2015-11-25 | 中国石油天然气集团公司 | High-angle well/horizontal well formation resistivity anisotropy correcting method |
| US20160109605A1 (en) * | 2014-10-20 | 2016-04-21 | Schlumberger Technology Corporation | Fast Model Based Inversion of Acoustic Impedance of Annulus Behind Casing |
| CN106951660A (en) * | 2017-04-05 | 2017-07-14 | 中国石油天然气股份有限公司 | Sea facies clastic rock horizontal well reservoir logging interpretation method and device |
| CN107678064A (en) * | 2016-08-02 | 2018-02-09 | 中石化石油工程技术服务有限公司 | A kind of interval transit time real time extracting method |
| CN110320562A (en) * | 2018-03-30 | 2019-10-11 | 中国石油化工股份有限公司 | The bearing calibration of the shale gas reservoir-level sound wave in borehole time difference |
| CN110826249A (en) * | 2019-11-18 | 2020-02-21 | 核工业北京地质研究院 | Uranium ore comprehensive logging curve borehole environment correction method |
-
2021
- 2021-07-14 CN CN202110797553.4A patent/CN113356834B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101787884A (en) * | 2010-01-28 | 2010-07-28 | 中国石油集团川庆钻探工程有限公司 | Reservoir fluid type discrimination method based on difference value of acoustic porosity and neutron porosity |
| US20150267527A1 (en) * | 2014-03-21 | 2015-09-24 | Conocophillips Company | Method for analysing pore pressure in shale formations |
| US20160109605A1 (en) * | 2014-10-20 | 2016-04-21 | Schlumberger Technology Corporation | Fast Model Based Inversion of Acoustic Impedance of Annulus Behind Casing |
| CN105089663A (en) * | 2015-07-31 | 2015-11-25 | 中国石油天然气集团公司 | High-angle well/horizontal well formation resistivity anisotropy correcting method |
| CN107678064A (en) * | 2016-08-02 | 2018-02-09 | 中石化石油工程技术服务有限公司 | A kind of interval transit time real time extracting method |
| CN106951660A (en) * | 2017-04-05 | 2017-07-14 | 中国石油天然气股份有限公司 | Sea facies clastic rock horizontal well reservoir logging interpretation method and device |
| CN110320562A (en) * | 2018-03-30 | 2019-10-11 | 中国石油化工股份有限公司 | The bearing calibration of the shale gas reservoir-level sound wave in borehole time difference |
| CN110826249A (en) * | 2019-11-18 | 2020-02-21 | 核工业北京地质研究院 | Uranium ore comprehensive logging curve borehole environment correction method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113356834B (en) | 2023-04-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3433643B1 (en) | Method and device for estimating sonic slowness in a subterranean formation | |
| EP2093591B1 (en) | Method for Three Dimensional Seismic Travel Time Tomography in Transversely Isotropic Media | |
| US8670288B2 (en) | Velocity model for well time-depth conversion | |
| US20070162249A1 (en) | Traveltime calculation in three dimensional transversely isotropic (3D TTI) media by the fast marching method | |
| US20150049585A1 (en) | Acoustic logging systems and methods employing multi-mode inversion for anisotropy and shear slowness | |
| RU2611096C2 (en) | Calculation of anisotropy parameter | |
| WO2017210231A1 (en) | Estimation of horizontal stresses and nonlinear constants in anisotropic formations such as interbedded carbonate layers in organic-shale reservoirs | |
| CN103827696A (en) | Multi-well anisotropy inversion | |
| CN111257946B (en) | Geophysical drilling guiding method and method for updating stratum seismic velocity | |
| US11237288B2 (en) | Verifying measurements of elastic anisotropy parameters in an anisotropic wellbore environment | |
| CN103018778A (en) | Method and equipment for reservoir prediction by correcting acoustic logging curve | |
| CN112034517B (en) | Ground stress data prediction method and device, storage medium and electronic equipment | |
| CN101285381A (en) | Process for inversing soft ground horizontal wave velocity by leaky mode waves | |
| Tsang et al. | Modeling shear rigidity of stratified bedrock in site response analysis | |
| CA2731780C (en) | Estimating formation stresses using radial profiles of three shear moduli | |
| CN113356834B (en) | Borehole acoustic time difference correction method, correction device and electronic equipment | |
| CN111257932B (en) | Method for updating stratum seismic velocity | |
| US10989824B2 (en) | Depth-dependent mud density determination and processing for horizontal shear slowness in vertical transverse isotropy environment using full-waveform sonic data | |
| CN114740532B (en) | Cement acoustic impedance and slurry acoustic impedance calculation method | |
| US20140350857A1 (en) | Method Of Mapping A Subterranean Formation Based Upon Wellbore Position And Seismic Data And Related System | |
| US20240077630A1 (en) | A method of and apparatus for processing seismic signals | |
| CN115201915A (en) | Ocean wide azimuth data deep dessert prediction method, storage medium and equipment | |
| CN109061737A (en) | A kind of method for predicting reservoir and device based on synthetic seismogram | |
| US20240044239A1 (en) | A method of and apparatus for optimizing a drilling strategy | |
| CN113123780B (en) | Method and device for determining thickness of hydrocarbon source rock at each level of single well and storage medium |
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 |