The present disclosure relates to the field of optical fibre sensing technology, in particular to a method and device for calibrating the depth of an optical fibre in a well, a processor and a computer-readable storage medium.
With the continuous development of science and technology, more and more new technologies are constantly applied in various fields of life and production. Since the data acquired by optical fibres has significant advantages such as high density sampling, high temperature and high pressure resistance and whole well section acquisition, the optical fibres are widely used in the field of borehole geophysics.
Reliable formation information can be obtained by directly coupling the optical fibres with the formation. However, in the process of processing the optical fibres into an armoured optical cable, there is a certain redundant length in order to adapt to the stretching of the cable, that is, the length of the optical cable is different from the length of the fibre core therein, and it is necessary to know the actual depth position of each receiving point of optical fibre sensing data when optical fibre sensing is applied in wells. In the prior art, firstly, the wellhead position is taped for depth zero setting in DAS output channels, and then the optical fibre detection tool OTDR (optical time-domain reflectometer) is used to measure the length of whole optical cable, and the length is compared with the length read by the cable counter to obtain the length correction coefficient, so as to calculate the accurate depth position of each optical fibre when the data is output.
However, in the actual application process, the existing depth calibration method is greatly influenced by the accuracy of the optical fibre detection tool OTDR and the counting accuracy of the above-mentioned cable counter, so there is a technical problem of insufficient calculation accuracy, and the actual needs of current technicians cannot be satisfied.
In order to overcome the above technical problems existing in the prior art, embodiments of the present disclosure provide a method for calibrating the depth of an optical fibre in a well, which performs de-noising processing on the actual seismic wave field data acquired by the optical fibre, and compares the existing logging data and the optical fibre seismic wave field data for analysis, so that the depth calibration information of the optical fibre is accurately determined, and the depth calibration accuracy is improved.
In order to achieve the above object, an embodiment of the present disclosure provides a method for calibrating the depth of an optical fibre in a well, comprising: acquiring borehole seismic wave field data based on an optical fibre acoustic wave sensor; determining first arrival time information based on the borehole seismic wave field data; determining a downgoing wave first arrival amplitude based on the first arrival time information and the borehole seismic wave field data; determining an optical fibre amplitude feature point based on the downgoing wave first arrival amplitude; determining a wellhead initial position and a receiving point spacing of the optical fibre in the well based on the optical fibre amplitude feature point and logging curve feature points; and determining depth calibration information of the optical fibre in the well based on the wellhead initial position and the receiving point spacing.
Preferably, the method further comprises: before determining the downgoing wave first arrival amplitude, optimising the borehole seismic wave field data based on a preset filtering algorithm to obtain an optimised downgoing wave field; and determining the downgoing wave first arrival amplitude based on the optimised downgoing wave field and the first arrival time information.
Preferably, determining the downgoing wave first arrival amplitude based on the first arrival time information and the borehole seismic wave field data comprises: determining a preset short time window based on the first arrival time information; determining specific seismic wave field data within the preset short time window; and determining a downgoing wave first arrival amplitude based on the specific seismic wave field data.
Preferably, determining the amplitude feature point based on the downgoing wave first arrival amplitude comprises: performing a spherical spreading amplitude compensation operation on the downgoing wave first arrival amplitude to obtain a compensated amplitude; generating a relative impedance curve based on the compensated amplitude; acquiring preset logging data; and determining the corresponding optical fibre amplitude feature point based on the preset logging data and the relative impedance curve.
Preferably, the method further comprises: determining a borehole sampling parameter of the optical fibre in the well based on the depth calibration information; performing a depth resampling operation based on the borehole sampling parameter to obtain resampled data; and generating updated depth calibration information based on the resampled data.
Correspondingly, the present disclosure also provides a device for calibrating the depth of an optical fibre in a well, comprising: a seismic wave field data acquisition unit, used for acquiring borehole seismic wave field data based on an optical fibre acoustic wave sensor; a first arrival time determination unit, used for determining first arrival time information based on the borehole seismic wave field data; a downgoing wave first arrival amplitude determination unit, used for determining a downgoing wave first arrival amplitude based on the first arrival time information and the borehole seismic wave field data; a feature point determination unit, used for determining an optical fibre amplitude feature point based on the downgoing wave first arrival amplitude; an optical fibre parameter determination unit, used for determining a wellhead initial position and a receiving point spacing of the optical fibre in the well based on the optical fibre amplitude feature point; and a depth calibration unit, used for determining depth calibration information of the optical fibre in the well based on the wellhead initial position and the receiving point spacing.
Preferably, the device further comprises an optimisation unit, used for: optimising the borehole seismic wave field data based on a preset filtering algorithm to obtain an optimised downgoing wave field; the downgoing wave first arrival amplitude determination unit is further used for determining the downgoing wave first arrival amplitude based on the optimised downgoing wave field and the first arrival time information.
Preferably, the downgoing wave first arrival amplitude determination unit comprises: a time window determination module, used for determining a preset short time window based on the first arrival time information; a data determination module, used for determining specific seismic wave field data within the preset short time window; and a first arrival amplitude determination module, used for determining a downgoing wave first arrival amplitude based on the specific seismic wave field data.
Preferably, the feature point determination unit comprises: a compensation module, used for performing a spherical spreading amplitude compensation operation on the downgoing wave first arrival amplitude to obtain a compensated amplitude; an impedance information generation module, used for generating a relative impedance curve based on the compensated amplitude; a logging data acquisition module, used for acquiring preset logging data; and a feature point generation module, used for determining the corresponding optical fibre amplitude feature point based on the preset logging data and the relative impedance curve.
Preferably, the device further comprises a resampling unit, used for: determining a borehole sampling parameter of the optical fibre in the well based on the depth calibration information; performing a depth resampling operation based on the borehole sampling parameter to obtain resampled data; and generating updated depth calibration information based on the resampled data.
In another aspect, the present disclosure also provides a processor, configured to perform the method for calibrating the depth of an optical fibre in a well provided by the present disclosure.
In another aspect, the present disclosure also provides a computer-readable storage medium, storing a computer program therein, wherein the computer program, when executed by a processor, implements the method for calibrating the depth of an optical fibre in a well provided by the present disclosure.
Through the technical solutions provided by the present disclosure, the present disclosure has at least the following technical effects:
by improving the existing optical fibre depth correction method, de-noising processing is performed on the seismic wave field data acquired by the field optical fibre, and then the wellhead initial position and receiving point spacing of the optical fibre in the well are accurately determined in combination with the existing logging data, so that the depth calibration position of the downhole optical fibre can be accurately determined, and the depth calibration accuracy is improved.
Other features and advantages of embodiments of the present disclosure will be described in detail in the Detailed Description section that follows.
The accompanying drawings are intended to provide a further understanding of the embodiments of the present disclosure and constitute a part of this specification, and together with the detailed description below, serve to explain, but not to limit, the embodiments of the present disclosure. In the drawings:
The specific implementations of the embodiments of the present disclosure will be described in detail with reference to the accompanying drawings. It should be understood that the specific implementations described here are only used to illustrate and explain the embodiments of the present disclosure, and are not used to limit the embodiments of the present disclosure.
The terms “system” and “network” in the embodiments of the present disclosure can be used interchangeably. “plurality of” means two or more, and in view of this, “plurality of” can also be understood as “at least two” in the embodiments of the present disclosure. “and/or”, which describes the association relationship of associated objects, indicates that three relationships can exist, for example, A and/or B, which can mean that A exists alone, A and B exist together, and B exists alone. In addition, the character “/”, unless otherwise specified, generally indicates that the associated objects before and after the “/” are in an “OR” relationship. In addition, it should be understood that in the description of the embodiments of the present disclosure, words such as “first” and “second” are only used for the purpose of distinguishing descriptions, and cannot be understood as indicating or implying relative importance, nor can they be understood as indicating or implying the order.
The background of the present disclosure will be first described below.
When the formation is observed, seismic wave field and other data of the formation are acquired by disposing an optical fibre in a well, connecting the optical fibre to a DAS (distributed acoustic sensor), and collecting information of the optical fibre by the DAS. However, in the actual application process, the prior art cannot accurately detect the exact position of the optical fibre in the well, thereby causing great trouble for the subsequent analysis of the optical fibre sensing data.
In order to solve the above technical problem, referring to
In one possible implementation, borehole seismic wave field data is first acquired by excitation on the ground or in the well with an explosive source or an artificial vibroseis, and then the borehole seismic wave field data is collected by the optical fibre. At this time, first arrival time information is determined according to the borehole seismic wave field data, the first arrival time information can be obtained by a preset first arrival picking method, which includes, but is not limited to, a maximum energy method, a short-term average/long-term average method, a cross-correlation method, a phase tracking method, and the like. At this time, the downgoing wave first arrival amplitude can be determined based on the first arrival time information and the borehole seismic wave field data.
However, in the actual application process, there is a lot of noise in the borehole seismic wave field data that is initially acquired, and therefore, in order to improve the accuracy of subsequent processing, it is necessary to optimise the borehole seismic wave field data.
In the embodiment of the present disclosure, the method further comprises: before determining the downgoing wave first arrival amplitude, optimising the borehole seismic wave field data based on a preset filtering algorithm to obtain an optimised downgoing wave field; and determining the downgoing wave first arrival amplitude based on the optimised downgoing wave field and the first arrival time information.
After the above borehole seismic wave field data is acquired and before the downgoing wave first arrival amplitude is determined, the borehole seismic wave field data can be optimised by methods such as a median filtering method or a mean filtering method, so as to remove upgoing waves in the borehole seismic wave field data, eliminate the interference effects of upgoing waves in the original borehole seismic wave field data, and improve the accuracy of subsequent data processing, and at this time, the downgoing wave first arrival amplitude can be determined based on the optimised information and the acquired first arrival time information.
Referring to
In one possible implementation, after the first arrival time information is acquired, a preset short time window is determined according to the first arrival time information, the preset short time window can be a dominant period in the first arrival time information or a time window having a length of 20-50 ms, and corresponding specific seismic wave field data is extracted within the preset short time window, a downgoing wave first arrival amplitude is determined according to the specific seismic wave field data, a root mean square amplitude or a maximum amplitude of the specific seismic wave field data can be determined as the downgoing wave first arrival amplitude information Ai, wherein i is a corresponding channel number for different channels in the well, and each channel has a corresponding first arrival amplitude. At this time, the optical fibre amplitude feature point is further determined according to the amplitude information.
Referring to
In one possible implementation, in order to eliminate the influence caused by the seismic waves on the seismic amplitude information during propagation, a spherical spreading amplitude compensation operation is first performed on the downgoing wave first arrival amplitude of the optical fibre seismic wave field to eliminate the influence of amplitude variation caused by spherical spreading, and obtain the compensated information. The relative impedance curve is then obtained by reciprocal calculation based on the compensated information, the relative impedance curve is the relative impedance Ir, referring to
At this time, the wellhead initial position and receiving point spacing of the optical fibre in the well are determined according to the optical fibre amplitude feature points. For example, in the embodiment of the present disclosure, the initial position i0 and the receiving point spacing ΔZ of the optical fibre in the well can be determined by the following calculation rules respectively:
ΔZ=(D2−D1)/(i2−i1);
i0=(i1*D2−i2*D1)/(D2−D1);
at this time, the depth calibration information of the optical fibre in the well can be determined according to the wellhead initial position and the receiving point spacing, for example, the accurate downhole depth calibration information of the optical fibre in the well can be determined by performing depth resampling on the seismic channel below the depth zero point with ΔZ according to the initial position i0 and the receiving point spacing ΔZ obtained by calculation.
In the embodiment of the present disclosure, the existing logging data and the borehole seismic wave field data detected by the optical fibre are analysed, and the borehole seismic wave field data is optimised, so that accurate seismic wave field data is obtained, and then accurate optical fibre depth calibration information is further obtained, which provides accurate data support for subsequent analysis of data collected through the optical fibre, and satisfies the actual demands of technicians. In order to further improve the accuracy of the data, in the embodiment of the present disclosure, the method further comprises: determining a borehole sampling parameter of the optical fibre in the well based on the depth calibration information; and performing a depth resampling operation based on the borehole sampling parameter to obtain resampled data.
For example, after the depth calibration information is determined for the first time, the borehole sampling parameter of the optical fibre in the well can be determined again according to the depth calibration information, a depth resampling operation is performed according to the well sampling parameter to obtain resampled data, the data processing process is performed according to the resampled data to obtain more accurate optical fibre depth calibration information, and the depth calibration information acquired for the first time is updated based on the more accurate optical fibre depth calibration information.
In the embodiment of the present disclosure, by repeating sampling for a plurality of times and comparing the feature points, the accuracy of data sampling and data analysis can be further improved, and the technical problem of large deviation existing in single data sampling and analysis can be effectively avoided, and thus the accuracy of data analysis can be further improved.
A device for calibrating the depth of an optical fibre in a well according to an embodiment of the present disclosure will now be described with reference to the accompanying drawings.
Referring to
In an embodiment of the present disclosure, the device further comprises an optimisation unit, used for optimising the borehole seismic wave field data based on a preset filtering algorithm to obtain an optimised downgoing wave field; the downgoing wave first arrival amplitude determination unit is further used for determining the downgoing wave first arrival amplitude based on the optimised downgoing wave field and the first arrival time information.
In an embodiment of the present disclosure, the downgoing wave first arrival amplitude determination unit comprises a time window determination module, used for determining a preset short time window based on the first arrival time information; a data determination module, used for determining specific seismic wave field data within the preset short time window; and a first arrival amplitude determination module, used for determining a downgoing wave first arrival amplitude based on the specific seismic wave field data.
In an embodiment of the present disclosure, the feature point determination unit comprises a compensation module, used for performing a spherical spreading amplitude compensation operation on the downgoing wave first arrival amplitude to obtain compensated information; an impedance information generation module, used for generating a relative impedance curve based on the compensated amplitude; a logging data acquisition module, used for acquiring preset logging data; and a feature point generation module, used for determining corresponding optical fibre amplitude feature points based on the preset logging data and the relative impedance curve.
In an embodiment of the present disclosure, the device further comprises a resampling unit, used for determining a borehole sampling parameter of the optical fibre in the well based on the depth calibration information; performing a depth resampling operation based on the borehole sampling parameter to obtain resampled data; and generating updated depth calibration information based on the resampled data.
Further, embodiments of the present disclosure also provide a processor configured to perform the method for calibrating the depth of an optical fibre in a well according to the present disclosure.
Further, embodiments of the present disclosure also provide a computer-readable storage medium, storing a computer program therein, wherein the computer program, when executed by a processor, implements the method for calibrating the depth of an optical fibre in a well according to the present disclosure.
Optional implementations of embodiments of the present disclosure are described in detail above in conjunction with the accompanying drawings, however, the embodiments of the present disclosure are not limited to the specific details in the above implementations, and many simple variations can be made to the technical solutions of the embodiments of the present disclosure within the scope of the technical idea of the embodiments of the present disclosure, and these simple variations all belong to the protection scope of the embodiments of the present disclosure.
It needs to be further noted that the various specific features described in the above detailed description can be combined in any suitable manner without contradiction. In order to avoid unnecessary repetition, embodiments of the present disclosure will not explain any possible combinations.
Those skilled in the art will appreciate that all or part of the steps in the method for realising the embodiments described above may be performed by instructing related hardware through a program, the program is stored in a storage medium and includes instructions for causing a single-chip microcomputer, chip or processor to execute all or part of the steps of the method according to the various embodiments of the present application. The aforementioned storage medium includes a USB flash disk, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk or an optical disk and other media that can store program codes.
In addition, various different implementations of the embodiments of the present disclosure can also be arbitrarily combined, as long as they do not violate the idea of the embodiments of the present disclosure, which should also be considered as the content disclosed in the embodiments of the present disclosure.
Number | Date | Country | Kind |
---|---|---|---|
202111648600.5 | Dec 2021 | CN | national |
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CN2022/102880 | 6/30/2022 | WO |