Natural potential differential measurement integral inversion method, device, equipment and medium

Abstract

The invention discloses a natural potential difference measurement integral inversion method, device and equipment and a medium, and relates to the technical field of oilfield development logging, and the method comprises the steps: obtaining a fixed electrode distance potential difference measurement curve of a to-be-logged well changing with the depth of a borehole; integrating the fixed electrode distance potential difference measurement curve to obtain a natural potential measurement curve, and folding each natural potential value in the natural potential measurement curve according to a preset period to obtain a periodic analysis chart; performing data remainder taking, data deletion and correction curve drawing on the periodic analysis chart to obtain a correction curve, and correcting the two-dimensional array storing each natural potential value in the periodic analysis chart by using the correction curve; and generating a continuous natural potential curve based on the corrected two-dimensional array, and taking the continuous natural potential curve as a natural potential difference measurement integral inversion result of the to-be-logged well, thereby realizing natural potential difference measurement integral inversion of the to-be-logged well even if the to-be-logged well does not contain a stable stratum, and improving the accuracy of the natural potential difference measurement integral inversion.

Description

Technical Field

[0001] This invention relates to the field of well logging technology for oilfield development, and in particular to a method, apparatus, equipment and medium for integral inversion of natural potential difference measurement. Background Technology

[0002] With the technological upgrade in oilfield development, development wells mostly adopt horizontal well designs. However, due to well conditions, cable delivery cannot be completed in horizontal wells. Currently, horizontal wells and complex wells mostly use direct-push logging. However, the ideal state for the spontaneous potential reference electrode in direct-push logging is to use the drill string to achieve a reliable, low-resistance connection between the reference electrode and the ground. In reality, due to the age and severe corrosion of the drill string, a large amount of insulating thread oil is applied to the joints to avoid damage to the threads, resulting in unstable contact resistance of the drill string. If the contact resistance is high, it is equivalent to the reference electrode being only connected to the instrument housing or only to one or two drill strings. This will lead to a decrease in the quality of the spontaneous potential measurement curve, mainly manifested in low curve amplitude and unstable mudstone baseline.

[0003] As can be seen from the above, how to achieve integral inversion of spontaneous potential difference measurement in non-stationary strata and improve the accuracy of spontaneous potential difference measurement integral inversion is a problem to be solved in this field. Summary of the Invention

[0004] In view of this, the purpose of this invention is to provide a method, apparatus, device, and medium for integrating and inverting the spontaneous potential difference measurement, which can realize the integrating and inverting of the spontaneous potential difference measurement in non-stationary strata and improve the accuracy of the integrating and inverting of the spontaneous potential difference measurement. The specific solution is as follows:

[0005] In a first aspect, this application discloses a method for integral inversion of natural potential difference measurement, comprising:

[0006] Obtain the fixed electrode distance potential difference measurement curve of the well to be tested as the wellbore depth changes;

[0007] The potential difference measurement curve of the fixed electrode is integrated to obtain the natural potential measurement curve. The natural potential values ​​in the natural potential measurement curve are folded according to a preset period to obtain a periodic analysis graph.

[0008] The periodic analysis graph is subjected to data remainder, data deletion, and correction curve plotting to obtain a correction curve. The correction curve is then used to correct the two-dimensional array in the periodic analysis graph used to store each of the natural potential values ​​to obtain the corrected two-dimensional array.

[0009] A continuous spontaneous potential curve is generated based on the corrected two-dimensional array, and the continuous spontaneous potential curve is used as the result of the spontaneous potential differential measurement integration inversion of the well to be measured.

[0010] Optionally, acquiring the fixed electrode distance potential difference measurement curve of the well to be tested as a function of wellbore depth includes:

[0011] A differential measurement logging instrument is used to measure the well to be logged, so as to directly obtain the fixed electrode distance potential difference measurement curve of the well as the wellbore depth from the differential measurement logging instrument.

[0012] Optionally, integrating the potential difference measurement curve of the fixed electrode to obtain the natural potential measurement curve includes:

[0013] The natural potential measurement curve is obtained by integrating all potential points in the fixed electrode-to-potential difference measurement curve using the natural potential value calculation formula; wherein, the natural potential value calculation formula is:

[0014] ;

[0015] in, Let L be the natural potential value at the m-th potential point, and L be the electrode distance. Let E(m) be the natural potential value of a point that is a distance of the electrode distance from the m-th potential point, and let E(m) be the potential difference of the m-th potential point.

[0016] Optionally, folding the natural potential values ​​in the natural potential measurement curve according to a preset period includes:

[0017] The natural potential values ​​in the natural potential measurement curve are folded according to a preset period using a two-dimensional array calculation formula; the two-dimensional array calculation formula is as follows:

[0018] ;

[0019] in, The x-coordinate of the two-dimensional array Let T be the i-th period, and T be the period. Let be the natural potential value corresponding to the horizontal coordinate of the two-dimensional array in the i-th period.

[0020] Optionally, the step of performing data modulo, data deletion, and correction curve plotting on the periodic analysis graph to obtain the correction curve includes:

[0021] Determine the total number of data points in the two-dimensional array of the periodic analysis graph, and take the remainder between the total number of data points and the preset period to obtain the remainder;

[0022] The periodic analysis graph is then deleting data based on the remainder, and a correction curve is plotted on the periodic analysis graph after data deletion to obtain the correction curve.

[0023] Optionally, the step of deleting data from the periodicity analysis graph according to the remainder, and plotting a correction curve on the periodicity analysis graph after data deletion to obtain the correction curve, includes:

[0024] Determine whether the remainder is 0. If the remainder is not 0, delete data from the two-dimensional array according to the remainder to obtain the periodic analysis chart after data deletion.

[0025] The natural potential values ​​of all phase points in different periods of the periodic analysis graph after data deletion are averaged to obtain the longitudinal average values, and the correction curves are plotted using the longitudinal average values.

[0026] Optionally, the step of correcting the two-dimensional array in the periodic analysis graph using the correction curve to obtain the corrected two-dimensional array; and generating a continuous natural potential curve based on the corrected two-dimensional array, includes:

[0027] The two-dimensional array and the correction curve in the periodic analysis graph are subtracted to obtain the corrected two-dimensional array;

[0028] The corrected two-dimensional array is expanded and spliced ​​to generate a continuous natural potential curve.

[0029] Secondly, this application discloses a natural potential difference measurement integration and inversion device, comprising:

[0030] The curve acquisition module is used to acquire the fixed electrode distance potential difference measurement curve of the well to be tested, which varies with the wellbore depth.

[0031] The curve integration and potential value folding module is used to integrate the potential difference measurement curve of the fixed electrode to obtain the natural potential measurement curve, and fold each natural potential value in the natural potential measurement curve according to a preset period to obtain a periodic analysis graph.

[0032] The correction module is used to perform data remainder, data deletion and correction curve plotting on the periodic analysis graph to obtain the correction curve, and to use the correction curve to correct the two-dimensional array in the periodic analysis graph used to store each of the natural potential values ​​to obtain the corrected two-dimensional array.

[0033] The spontaneous potential curve generation module is used to generate a continuous spontaneous potential curve based on the corrected two-dimensional array, and to use the continuous spontaneous potential curve as the result of the spontaneous potential differential measurement integration inversion of the well to be measured.

[0034] Thirdly, this application discloses an electronic device, including:

[0035] Memory, used to store computer programs;

[0036] A processor is used to execute the computer program to implement the aforementioned natural potential difference measurement integral inversion method.

[0037] Fourthly, this application discloses a computer storage medium for storing a computer program; wherein, when the computer program is executed by a processor, it implements the steps of the aforementioned disclosed natural potential difference measurement integration inversion method.

[0038] As can be seen, this application provides a spontaneous potential difference measurement integration inversion method, including acquiring a fixed electrode distance potential difference measurement curve of the well to be measured as it varies with wellbore depth; integrating the fixed electrode distance potential difference measurement curve to obtain a spontaneous potential measurement curve; folding each spontaneous potential value in the spontaneous potential measurement curve according to a preset period to obtain a periodic analysis diagram; performing data remainder, data deletion, and correction curve plotting on the periodic analysis diagram to obtain a correction curve; using the correction curve to correct the two-dimensional array used to store each spontaneous potential value in the periodic analysis diagram to obtain the corrected two-dimensional array; generating a continuous spontaneous potential curve based on the corrected two-dimensional array; and using the continuous spontaneous potential curve as the spontaneous potential difference measurement integration inversion result of the well to be measured. After obtaining the fixed electrode distance potential difference measurement curve of the well to be logged, this application performs curve integration and spontaneous potential value folding to obtain a periodic analysis diagram. Combining physical principles and mathematical methods, the periodic analysis diagram is subjected to data remainder, data deletion, and correction curve plotting. The correction curve is used to correct the two-dimensional array in the periodic analysis diagram, thereby generating a continuous spontaneous potential curve, improving the quality and accuracy of the spontaneous potential measurement curve. It can calculate the trend of change even without a stable reference formation, without the need for a stable reference formation, realizing the integral inversion of spontaneous potential difference measurement in non-stable formations, and improving the accuracy of spontaneous potential difference measurement integral inversion. Attached Figure Description

[0039] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on the provided drawings without creative effort.

[0040] Figure 1 This is a flowchart of a natural potential difference measurement integration inversion method disclosed in this application;

[0041] Figure 2 This is a graph showing the measurement curve of the potential difference between a fixed electrode and the electrode disclosed in this application.

[0042] Figure 3 This application discloses a measurement curve of integrated natural potential values.

[0043] Figure 4 This is a schematic diagram of the principle of a fixed electrode distance potential difference measurement curve unfolded periodically according to the present application.

[0044] Figure 5 This is a periodic analysis diagram disclosed in this application;

[0045] Figure 6 This is a modified curve diagram disclosed in this application;

[0046] Figure 7 This application discloses a periodic analysis graph including a corrected two-dimensional array;

[0047] Figure 8 This application discloses a continuous natural potential curve.

[0048] Figure 9 This application discloses a comparative full-well curve diagram;

[0049] Figure 10 This is a magnified view of a portion of the overall well curve disclosed in this application;

[0050] Figure 11 This is a schematic diagram of the structure of a natural potential difference measurement integration and inversion device disclosed in this application;

[0051] Figure 12 This application provides a structural diagram of an electronic device. Detailed Implementation

[0052] The technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0053] With the technological upgrades in oilfield development, horizontal wells are increasingly being designed for development. However, the well conditions of horizontal wells often prevent cable delivery from being completed. Currently, horizontal and complex wells mostly use direct-push logging. However, ideally, the spontaneous potential (SP) reference electrode in direct-push logging should be reliably and with low resistance connected to the ground using the drill string. In reality, due to the age and severe corrosion of the drill string, and the large amount of insulating thread oil applied to the joints to prevent thread damage, the contact resistance of the drill string becomes unstable. If the contact resistance is high, it's as if the reference electrode is only connected to the instrument housing or only to one or two drill strings. This leads to a decrease in the quality of the SP measurement curve, mainly manifested in low curve amplitude and unstable mudstone baseline. Therefore, how to achieve integral inversion of SP differential measurements even without stable formations, and improve the accuracy of SP differential measurement integral inversion, is a problem that needs to be solved in this field.

[0054] See Figure 1 As shown, this embodiment of the invention discloses a method for integral inversion of natural potential difference measurement, which specifically includes:

[0055] Step S11: Obtain the fixed electrode distance potential difference measurement curve of the well to be tested as the wellbore depth changes.

[0056] In this embodiment, a differential measurement logging instrument is used to measure the well to be logged, so as to directly obtain the fixed electrode distance potential difference measurement curve of the well as the wellbore depth from the differential measurement logging instrument. Specifically, the fixed electrode distance potential difference measurement curve of the potential difference E as the wellbore depth h is directly obtained by the differential measurement logging instrument as shown below. Figure 2 As shown.

[0057] Step S12: Integrate the potential difference measurement curve of the fixed electrode to obtain the natural potential measurement curve, and fold each natural potential value in the natural potential measurement curve according to a preset period to obtain a periodic analysis graph.

[0058] In this embodiment, the natural potential value is calculated by integrating all potential points in the fixed electrode-to-potential difference measurement curve using the natural potential value calculation formula to obtain the natural potential measurement curve; wherein, the natural potential value calculation formula is:

[0059] ;

[0060] in, Let L be the natural potential value at the m-th potential point, and L be the electrode distance. Let E(m) be the natural potential value of a point that is a distance of the electrode distance from the m-th potential point, and let E(m) be the potential difference of the m-th potential point.

[0061] Then, the natural potential values ​​in the natural potential measurement curve are folded according to a preset period using a two-dimensional array calculation formula; the two-dimensional array calculation formula is:

[0062] ;

[0063] in, The x-coordinate of the two-dimensional array Let T be the i-th period, and T be the period. Let be the natural potential value corresponding to the horizontal coordinate of the two-dimensional array in the i-th period.

[0064] Specifically, a preliminary integration is performed on the potential difference measurement curve at a fixed electrode distance. The integration method is as follows: assuming that the potential value at all phase points within the first electrode distance is 0, the integration is performed on the potential points corresponding to the same position within the electrode distance to obtain the natural potential value at the m-th potential point. The formula for calculating the natural potential value is shown above. The measurement curve obtained after integration, including each natural potential value, is shown below. Figure 3 As shown.

[0065] Determine the electrode distance L, with a set distance length as one period T. , To achieve resolution, there are several phase points (i.e., AK) within each period. The number of periods, n = h / T, is calculated. The principle of expanding the fixed electrode-potential difference measurement curve according to the period is as follows: Figure 4 As shown, the natural potential value is then folded according to the period T to obtain the following result: Figure 5 The periodic analysis diagram shown.

[0066] Step S13: Perform data remainder, data deletion, and correction curve plotting on the periodic analysis graph to obtain a correction curve. Use the correction curve to correct the two-dimensional array in the periodic analysis graph used to store each of the natural potential values ​​to obtain the corrected two-dimensional array.

[0067] In this embodiment, the total number of data points in the two-dimensional array of the periodic analysis graph is determined, and the remainder is obtained by taking the remainder between the total number of data points and a preset period. Data is deleted from the periodic analysis graph according to the remainder, and a correction curve is drawn on the periodic analysis graph after data deletion to obtain a correction curve. The two-dimensional array in the periodic analysis graph and the correction curve are subtracted to obtain the corrected two-dimensional array. The corrected two-dimensional array is expanded and spliced ​​to generate a continuous natural potential curve.

[0068] The data deletion and correction curve drawing process is as follows: determine whether the remainder is 0; if the remainder is not 0, delete the data in the two-dimensional array according to the remainder to obtain the periodic analysis diagram after data deletion; calculate the average value of the natural potential values ​​of all phase points in different periods of the periodic analysis diagram after data deletion to obtain each longitudinal average value; and draw the correction curve using each longitudinal average value.

[0069] Specifically, to determine whether the total number of data points in the two-dimensional array of the periodic analysis graph can form a complete period (by taking the remainder between the total number of data points and the preset period), the formula is as follows: q = N%T; where N is the total number of data points, T is the period, and q is the remainder. If a complete period cannot be formed (the remainder q is not 0), the last q data points are deleted, and then the average value of the natural potential values ​​of the corresponding phases of different periods is taken (or histogram analysis is performed first to find the data with the highest probability and then the average value is taken), for example:

[0070] ;

[0071] ;

[0072] …

[0073] ;

[0074] in, This is the longitudinal average. This represents the natural potential value of phase A in the first period. This represents the natural potential value of phase A in the second period. This represents the natural potential value of phase B in the 3rd period, and so on; n is the number of periods, and then the average value is... arrive Connect the curves to obtain the correction curve, as shown in the figure. Figure 6 As shown, the x-axis of the curve represents the phase point (AK), and the y-axis represents... The amplitude is calculated by subtracting the corrected curve from the two-dimensional array in the periodicity analysis graph to obtain the corrected two-dimensional array. The periodicity analysis graph including the corrected two-dimensional array is shown below. Figure 7 As shown:

[0075] Phase A:

[0076] ;

[0077] ;

[0078] …

[0079] ;

[0080] …

[0081] Phase K:

[0082] ;

[0083] ;

[0084] …

[0085] .

[0086] Step S14: Generate a continuous spontaneous potential curve based on the corrected two-dimensional array, and use the continuous spontaneous potential curve as the result of the spontaneous potential differential measurement integration inversion of the well to be measured.

[0087] The corrected two-dimensional array is unfolded and stitched together to form a continuous natural potential curve, such as Figure 8 As shown,

[0088] ;

[0089] in, This is the continuous natural potential curve after splicing, where T is the period, % is the modulo operation, and / is the integer division, taking the smaller nearest integer.

[0090] Comparing this patent with existing technologies, the prior art points out that during measurement, it is necessary to determine the spontaneous potential curve of a specified reference layer. If there is no reference layer with stable variation characteristics, measurement is impossible. Furthermore, if there are formations in the reference formation with spacing equal to or close to the first and second electrodes, "non-stable variation" will be misjudged as "stable variation," and the smaller the electrode spacing, the higher the probability of misjudgment, leading to reduced accuracy. Under the same testing environment, the original real curve, the spontaneous potential curves of this patent, and the prior art are compared and plotted. The comparison of the whole-well curves is shown below. Figure 9 As shown in the enlarged view of the full well curve, see the comparison. Figure 10 As shown. (Through) Figure 9 and Figure 10 It can be seen that when the natural potential curve does not have a stable flat segment (or the stable flat segment is not stable and flat enough), the comparative patent cannot truly reproduce the natural potential curve, while this application can truly reproduce the natural potential curve.

[0091] The innovation of this application is as follows: (1) Integral inversion calculation method: After differential measurement of the potential change curve, integration is required to restore the potential curve. During restoration, the potential change between adjacent points cannot be measured, so it is necessary to determine the potential relationship between adjacent points to determine the initial value of integration; (2) Calculation method of initial value of integration: The calculation method of the initial value of integration in this patent is calculated by using a probability method. First, it is assumed that the initial value of integration is zero. At the points with an interval of electrode distance throughout the well, the value obtained by integration will produce the same deviation. The deviation is added together and the average value is taken. The deviation from the zero value is the offset. The error of the initial value is eliminated by subtracting the average value from the value of all points.

[0092] This invention utilizes physical principles and mathematical methods to summarize the characteristics of spontaneous potential curves for differential measurement and integral reconstruction. It is applicable not only to spontaneous potential curves with "partially flat" characteristics, but also to curves that are completely "irregular" (noise curves), and can perfectly achieve integral reconstruction. Furthermore, it does not require a stable reference formation, and can use data analysis from the entire well to calculate the changing trend of non-stable formations, and use the changing trend to compensate for or offset the influence of unstable changing formations.

[0093] In this embodiment, a fixed electrode distance potential difference measurement curve of the well to be tested, varying with wellbore depth, is obtained; the fixed electrode distance potential difference measurement curve is integrated to obtain a spontaneous potential measurement curve; the spontaneous potential values ​​in the spontaneous potential measurement curve are folded according to a preset period to obtain a periodic analysis chart; the periodic analysis chart is subjected to data remainder, data deletion, and correction curve plotting to obtain a correction curve; the correction curve is used to correct the two-dimensional array used to store each spontaneous potential value in the periodic analysis chart to obtain the corrected two-dimensional array; a continuous spontaneous potential curve is generated based on the corrected two-dimensional array, and the continuous spontaneous potential curve is used as the spontaneous potential differential measurement integration inversion result of the well to be tested. After obtaining the fixed electrode distance potential difference measurement curve of the well to be logged, this application performs curve integration and spontaneous potential value folding to obtain a periodic analysis diagram. Combining physical principles and mathematical methods, the periodic analysis diagram is subjected to data remainder, data deletion, and correction curve plotting. The correction curve is used to correct the two-dimensional array in the periodic analysis diagram, thereby generating a continuous spontaneous potential curve, improving the quality and accuracy of the spontaneous potential measurement curve. It can calculate the changing trend of non-stationary formations without a stable reference formation, realizing the integral inversion of spontaneous potential difference measurement even without stable formations, and improving the accuracy of spontaneous potential difference measurement integral inversion.

[0094] See Figure 11 As shown, this embodiment of the invention discloses a natural potential difference measurement integration inversion device, which specifically includes:

[0095] The curve acquisition module 11 is used to acquire the fixed electrode distance potential difference measurement curve of the well to be measured as the wellbore depth changes.

[0096] The curve integration and potential value folding module 12 is used to integrate the potential difference measurement curve of the fixed electrode to obtain the natural potential measurement curve, and fold each natural potential value in the natural potential measurement curve according to a preset period to obtain a periodic analysis diagram.

[0097] The correction module 13 is used to perform data remainder, data deletion and correction curve drawing on the periodic analysis graph to obtain the correction curve, and to use the correction curve to correct the two-dimensional array in the periodic analysis graph used to store each of the natural potential values ​​to obtain the corrected two-dimensional array.

[0098] The spontaneous potential curve generation module 14 is used to generate a continuous spontaneous potential curve based on the corrected two-dimensional array, and to use the continuous spontaneous potential curve as the result of the spontaneous potential differential measurement integration inversion of the well to be measured.

[0099] In this embodiment, a fixed electrode distance potential difference measurement curve of the well to be tested, varying with wellbore depth, is obtained; the fixed electrode distance potential difference measurement curve is integrated to obtain a spontaneous potential measurement curve; the spontaneous potential values ​​in the spontaneous potential measurement curve are folded according to a preset period to obtain a periodic analysis chart; the periodic analysis chart is subjected to data remainder, data deletion, and correction curve plotting to obtain a correction curve; the correction curve is used to correct the two-dimensional array used to store each spontaneous potential value in the periodic analysis chart to obtain the corrected two-dimensional array; a continuous spontaneous potential curve is generated based on the corrected two-dimensional array, and the continuous spontaneous potential curve is used as the spontaneous potential differential measurement integration inversion result of the well to be tested. After obtaining the fixed electrode distance potential difference measurement curve of the well to be logged, this application performs curve integration and spontaneous potential value folding to obtain a periodic analysis diagram. Combining physical principles and mathematical methods, the periodic analysis diagram is subjected to data remainder, data deletion, and correction curve plotting. The correction curve is used to correct the two-dimensional array in the periodic analysis diagram, thereby generating a continuous spontaneous potential curve, improving the quality and accuracy of the spontaneous potential measurement curve. It can calculate the changing trend of non-stationary formations without a stable reference formation, realizing the integral inversion of spontaneous potential difference measurement even without stable formations, and improving the accuracy of spontaneous potential difference measurement integral inversion.

[0100] In some specific embodiments, the curve acquisition module 11 may specifically include:

[0101] The measurement module is used to measure the well to be logged using a differential measurement logging instrument, so as to directly obtain the fixed electrode distance potential difference measurement curve of the well to be logged as a function of wellbore depth from the differential measurement logging instrument.

[0102] In some specific embodiments, the curve integral and potential value folding module 12 may specifically include:

[0103] The integration module is used to integrate all potential points in the fixed electrode-to-potential difference measurement curve using the natural potential value calculation formula to obtain the natural potential measurement curve; wherein, the natural potential value calculation formula is:

[0104] ;

[0105] in, Let L be the natural potential value at the m-th potential point, and L be the electrode distance. Let E(m) be the natural potential value of a point that is a distance of the electrode distance from the m-th potential point, and let E(m) be the potential difference of the m-th potential point.

[0106] In some specific embodiments, the curve integral and potential value folding module 12 may specifically include:

[0107] The folding module is used to fold the natural potential values ​​in the natural potential measurement curve according to a preset period using a two-dimensional array calculation formula; the two-dimensional array calculation formula is:

[0108] ;

[0109] in, The x-coordinate of the two-dimensional array Let T be the i-th period, and T be the period. Let be the natural potential value corresponding to the horizontal coordinate of the two-dimensional array in the i-th period.

[0110] In some specific embodiments, the correction module 13 may specifically include:

[0111] The remainder module is used to determine the total number of data in the two-dimensional array in the periodic analysis graph, and to take the remainder by taking the remainder between the total number of data and the preset period.

[0112] The data deletion and curve plotting module is used to delete data from the periodic analysis graph according to the remainder, and to plot a correction curve on the periodic analysis graph after data deletion to obtain the correction curve.

[0113] In some specific embodiments, the correction module 13 may specifically include:

[0114] The data deletion module is used to determine whether the remainder is 0. If the remainder is not 0, the data in the two-dimensional array is deleted according to the remainder to obtain the periodic analysis chart after data deletion.

[0115] The average value calculation module is used to calculate the average value of the natural potential values ​​of all phase points of different periods in the periodic analysis graph after data deletion, so as to obtain the longitudinal average value, and to draw the correction curve using the longitudinal average value.

[0116] In some specific embodiments, the natural potential curve generation module 14 may specifically include:

[0117] The subtraction correction module is used to subtract the two-dimensional array and the correction curve in the periodic analysis graph to obtain the corrected two-dimensional array;

[0118] The unfolding and splicing module is used to unfold and splice the corrected two-dimensional array to generate a continuous natural potential curve.

[0119] Figure 12 This is a schematic diagram of an electronic device provided in an embodiment of this application. The electronic device 20 may specifically include: at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 stores a computer program, which is loaded and executed by the processor 21 to implement the relevant steps in the natural potential difference measurement integration inversion method disclosed in any of the foregoing embodiments.

[0120] In this embodiment, the power supply 23 is used to provide operating voltage for each hardware device on the electronic device 20; the communication interface 24 can create a data transmission channel between the electronic device 20 and external devices, and the communication protocol it follows can be any communication protocol applicable to the technical solution of this application, and is not specifically limited here; the input / output interface 25 is used to acquire external input data or output data to the outside world, and its specific interface type can be selected according to specific application needs, and is not specifically limited here.

[0121] In addition, the memory 22, as a carrier for resource storage, can be a read-only memory, random access memory, disk or optical disk, etc. The resources stored on it include operating system 221, computer program 222 and data 223, etc., and the storage method can be temporary storage or permanent storage.

[0122] The operating system 221 manages and controls the various hardware devices on the electronic device 20 and the computer program 222 to enable the processor 21 to perform calculations and processing on the data 223 in the memory 22. It can be Windows, Unix, Linux, etc. The computer program 222, in addition to including a computer program capable of performing the natural potential difference measurement integration inversion method executed by the electronic device 20 as disclosed in any of the foregoing embodiments, may further include computer programs capable of performing other specific tasks. The data 223 may include data received by the natural potential difference measurement integration inversion device from external devices, and may also include data collected by its own input / output interface 25.

[0123] The steps of the methods or algorithms described in conjunction with the embodiments disclosed herein can be implemented directly by hardware, a software module executed by a processor, or a combination of both. The software module can be located in random access memory (RAM), main memory, read-only memory (ROM), electrically programmable ROM, electrically erasable programmable ROM, registers, hard disk, removable disk, CD-ROM, or any other form of storage medium known in the art.

[0124] Furthermore, embodiments of this application also disclose a computer-readable storage medium storing a computer program. When the computer program is loaded and executed by a processor, it implements the steps of the natural potential difference measurement integration inversion method disclosed in any of the foregoing embodiments.

[0125] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.

[0126] The present invention provides a detailed description of a natural potential difference measurement integration inversion method, apparatus, device, and storage medium. Specific examples have been used to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of the present invention. Therefore, the content of this specification should not be construed as a limitation of the present invention.

Claims

1. A method for integral inversion of natural potential difference measurement, characterized in that, include: Obtain the fixed electrode distance potential difference measurement curve of the well to be tested as the wellbore depth changes; The potential difference measurement curve of the fixed electrode is integrated to obtain the natural potential measurement curve. The natural potential values ​​in the natural potential measurement curve are folded according to a preset period to obtain a periodic analysis graph. The periodic analysis graph is subjected to data remainder, data deletion, and correction curve plotting to obtain a correction curve. The correction curve is then used to correct the two-dimensional array in the periodic analysis graph used to store each of the natural potential values ​​to obtain the corrected two-dimensional array. A continuous spontaneous potential curve is generated based on the corrected two-dimensional array, and the continuous spontaneous potential curve is used as the result of the spontaneous potential differential measurement integration inversion of the well to be measured.

2. The method for integral inversion of natural potential difference measurement according to claim 1, characterized in that, The acquisition of the fixed electrode distance potential difference measurement curve of the well to be tested as a function of wellbore depth includes: A differential measurement logging instrument is used to measure the well to be logged, so as to directly obtain the fixed electrode distance potential difference measurement curve of the well as the wellbore depth from the differential measurement logging instrument.

3. The method for integral inversion of natural potential difference measurement according to claim 1, characterized in that, The step of integrating the potential difference measurement curve of the fixed electrode to obtain the natural potential measurement curve includes: The natural potential measurement curve is obtained by integrating all potential points in the fixed electrode-to-potential difference measurement curve using the natural potential value calculation formula; wherein, the natural potential value calculation formula is: ； in, Let L be the natural potential value at the m-th potential point, and L be the electrode distance. Let E(m) be the natural potential value of a point that is a distance of the electrode distance from the m-th potential point, and let E(m) be the potential difference of the m-th potential point.

4. The method for integral inversion of natural potential difference measurement according to claim 1, characterized in that, The step of folding each natural potential value in the natural potential measurement curve according to a preset period includes: The natural potential values ​​in the natural potential measurement curve are folded according to a preset period using a two-dimensional array calculation formula; the two-dimensional array calculation formula is as follows: ； in, The x-coordinate of the two-dimensional array Let T be the i-th period, and T be the period. Let be the natural potential value corresponding to the horizontal coordinate of the two-dimensional array in the i-th period.

5. The method for integral inversion of natural potential difference measurement according to claim 1, characterized in that, The process of performing data remainder, data deletion, and correction curve plotting on the periodic analysis graph to obtain the correction curve includes: Determine the total number of data points in the two-dimensional array of the periodic analysis graph, and take the remainder between the total number of data points and the preset period to obtain the remainder; The periodic analysis graph is then deleting data based on the remainder, and a correction curve is plotted on the periodic analysis graph after data deletion to obtain the correction curve.

6. The method for integral inversion of natural potential difference measurement according to claim 5, characterized in that, The step of deleting data from the periodicity analysis graph according to the remainder, and then plotting a correction curve on the periodicity analysis graph after data deletion to obtain the correction curve, includes: Determine whether the remainder is 0. If the remainder is not 0, delete data from the two-dimensional array according to the remainder to obtain the periodic analysis chart after data deletion. The natural potential values ​​of all phase points in different periods of the periodic analysis graph after data deletion are averaged to obtain the longitudinal average values, and the correction curves are plotted using the longitudinal average values.

7. The method for integral inversion of natural potential difference measurement according to any one of claims 1 to 6, characterized in that, The two-dimensional array in the periodic analysis graph is corrected using the correction curve to obtain the corrected two-dimensional array; Generate a continuous natural potential curve based on the corrected two-dimensional array, including: The two-dimensional array and the correction curve in the periodic analysis graph are subtracted to obtain the corrected two-dimensional array; The corrected two-dimensional array is expanded and spliced ​​to generate a continuous natural potential curve.

8. A natural potential difference measurement integration inversion device, characterized in that, include: The curve acquisition module is used to acquire the fixed electrode distance potential difference measurement curve of the well to be tested, which varies with the wellbore depth. The curve integration and potential value folding module is used to integrate the potential difference measurement curve of the fixed electrode to obtain the natural potential measurement curve, and fold each natural potential value in the natural potential measurement curve according to a preset period to obtain a periodic analysis graph. The correction module is used to perform data remainder, data deletion and correction curve plotting on the periodic analysis graph to obtain the correction curve, and to use the correction curve to correct the two-dimensional array in the periodic analysis graph used to store each of the natural potential values ​​to obtain the corrected two-dimensional array. The spontaneous potential curve generation module is used to generate a continuous spontaneous potential curve based on the corrected two-dimensional array, and to use the continuous spontaneous potential curve as the result of the spontaneous potential differential measurement integration inversion of the well to be measured.

9. An electronic device, characterized in that, include: Memory, used to store computer programs; A processor for executing the computer program to implement the natural potential difference measurement integral inversion method as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, Used to store a computer program; wherein, when the computer program is executed by a processor, it implements the natural potential difference measurement integral inversion method as described in any one of claims 1 to 7.

Images