A method for correlating active fault borehole stratigraphy with high resolution well log data

Abstract

The application discloses a kind of active fault borehole stratum comparison methods combined with high-resolution logging data, it is related to stratum comparison technical field, comprising: first respectively obtaining test fault region and test bottom layer region;Real-time data acquisition is carried out based on electrical imaging logging, and differential analysis model is established;The differential analysis model is optimized, and the ground fault comparison model is obtained;Active fault borehole stratum comparison method is obtained based on the data in ground fault comparison model;The application is used to solve the problem of low precision and efficiency in the existing active fault borehole stratum comparison method.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of stratum correlation, and particularly relates to a method for correlating a drilling stratum of an active fault by combining high-resolution logging data. BACKGROUND

[0002] An active fault is a trace left by seismic activity and records the history of seismic disasters; the activity degree of an active fault can be measured by the latest activity era and activity intensity; the later the latest activity era, the higher the activity intensity, and the higher the activity degree; a drilling stratum refers to a stratum sample and data obtained by drilling in geological exploration; drilling exploration is one of the basic exploration techniques for a buried active fault, and a core sample is obtained by drilling, and stratum division and correlation are performed to determine quantitative parameters such as the spatial position, dislocation mode, and dislocation period of the fault; the methods for drilling stratum division and correlation include core logging, magnetic susceptibility curve, and resistivity curve, which can help accurately divide and correlate strata.

[0003] The existing method for correlating a drilling stratum of an active fault is usually based on an existing logging curve sample set, encodes and decodes original data, builds a learning framework to perform data set experiments and error analysis, and thus realizes stratum horizon division and intelligent stratum correlation; although this improved method can accurately divide strata based on existing data, it only improves the analysis efficiency when multiple strata are simultaneously analyzed; if only two strata, such as an active fault and a drilling stratum, are correlated, the above method cannot quickly obtain the feature relationship between the strata according to the related parameters of the strata, and thus cannot effectively obtain a correlation method for correlating the strata, resulting in a problem of low stratum correlation efficiency; therefore, it is necessary to improve the existing method for correlating a drilling stratum of an active fault. SUMMARY

[0004] The present application aims to at least solve one of the technical problems in the prior art, and provides a method for correlating a drilling stratum of an active fault by combining high-resolution logging data, which can compensate for stratum loss and stratum order reversal caused by low drilling core sampling rate and mistakes of drilling workers, compensate for subjective errors in direct visual stratum layer logging, improve the resolution of drilling stratum layering, and improve the accuracy of active fault detection and identification, thereby solving the problem of low accuracy and efficiency of active fault detection and identification in the existing method for correlating a drilling stratum of an active fault.

[0005] To achieve the above-mentioned purpose, the present application provides a method for correlating a drilling stratum of an active fault by combining high-resolution logging data, which comprises the following steps:

[0006] respectively, and the test fault region and the test bottom region are respectively recorded as a test fault region and a test bottom region;

[0007] Real-time data acquisition is performed on the test fault region and the test bottom region based on the electrical imaging logging, and a differential analysis model is established based on the data obtained through the real-time data acquisition;

[0008] A plurality of test fault regions and test bottom regions at different depths are obtained, and the differential analysis model is optimized based on the data obtained through the electrical imaging logging, and the optimized differential analysis model is recorded as a ground fault comparison model;

[0009] An active fault borehole formation comparison method is obtained based on the data in the ground fault comparison model.

[0010] Further, the real-time data acquisition on the test fault region and the test bottom region based on the electrical imaging logging, and the establishment of the differential analysis model based on the data obtained through the real-time data acquisition, include:

[0011] The acquisition data of the real-time data acquisition on the test fault region and the test bottom region by the electrical imaging logging are respectively recorded as fault acquisition data and bottom acquisition data;

[0012] A resistivity analysis coordinate system, a density analysis coordinate system and a porosity analysis coordinate system are respectively established, wherein the unit of the X-axis of the resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system is m, the unit of the Y-axis of the resistivity analysis coordinate system is Ω·m, the unit of the Y-axis of the density analysis coordinate system is kg / m³, and the unit of the Y-axis of the porosity analysis coordinate system is %.

[0013] Further, the real-time data acquisition on the test fault region and the test bottom region based on the electrical imaging logging, and the establishment of the differential analysis model based on the data obtained through the real-time data acquisition, further include:

[0014] Based on the straight-line distance from the surface during the data acquisition by the electrical imaging logging in the fault acquisition data, and the resistivity, density and porosity in the fault acquisition data, corresponding curves in the resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system are plotted, and are respectively recorded as a fault resistivity curve, a fault density curve and a fault porosity curve;

[0015] Based on the analysis method of the fault acquisition data, corresponding curves of the bottom acquisition data in the resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system are obtained, and are respectively recorded as a bottom resistivity curve, a bottom density curve and a bottom porosity curve.

[0016] Further, the real-time data acquisition is performed on the test fault area and the test bottom layer area based on the electrical imaging logging, and the differential analysis model is established based on the data obtained through the real-time data acquisition, and the method further comprises:

[0017] The resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system are analyzed in sequence by using the coordinate system analysis method, and the multi-parameter contrast features are obtained based on the analysis results, wherein the multi-parameter contrast features comprise resistivity contrast features, density contrast features and porosity contrast features.

[0018] All curves corresponding to the fault acquisition data and the bottom layer acquisition data are analyzed by using the single-layer curve analysis method, and the single-layer contrast features corresponding to the active fault and the drilling bottom layer are obtained based on the analysis results.

[0019] Further, the coordinate system analysis method comprises:

[0020] For the coordinate system to be analyzed: the curve corresponding to the fault acquisition data of two curves in the coordinate system is recorded as curve UC, and the curve corresponding to the bottom layer acquisition data is recorded as curve IC;

[0021] The difference between the maximum value and the minimum value of the ordinate in the curve UC is recorded as Y1, and the difference between the maximum value and the minimum value of the ordinate in the curve IC is recorded as Y2; the point with the maximum absolute value of the slope in the curve UC and the curve IC is recorded as point K1 and point K2, respectively;

[0022] The value of α1 is set as Y1 / Y2; the value of α2 is set as the absolute value of the slope of point K1 in the curve UC divided by the absolute value of the slope of point K2 in the curve IC;

[0023] α3 is set as L1 / L2, wherein L1 is the distance between the point with the maximum ordinate in the curve UC and point K1, and L2 is the distance between the point with the maximum ordinate in the curve IC and point K2;

[0024] When the coordinate to be analyzed is the resistivity analysis coordinate system, [α1, α2, α3] is recorded as resistivity contrast features; when the coordinate to be analyzed is the density analysis coordinate system, [α1, α2, α3] is recorded as density contrast features; when the coordinate to be analyzed is the porosity analysis coordinate system, [α1, α2, α3] is recorded as porosity contrast features.

[0025] Further, the single-layer curve analysis method comprises:

[0026] For the fault acquisition data or the bottom layer acquisition data: the corresponding curves in the resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system are put into the same coordinate system and recorded as single-layer analysis curve DF1, single-layer analysis curve DF2 and single-layer analysis curve DF3 in sequence;

[0027] An interval occupied by the single-layer analysis curve DF1, the single-layer analysis curve DF2 and the single-layer analysis curve DF3 in the X-axis is denoted as [P1, P2]; Q points are uniformly obtained in X=P1 to X=P2, and all the Q points are denoted as single-layer analysis points;

[0028] For the abscissa r of any single-layer analysis point, distances between the point with the abscissa r in the single-layer analysis curve DF1, the single-layer analysis curve DF2 and the single-layer analysis curve DF3 and the X-axis are respectively denoted as a first analysis value, a second analysis value and a third analysis value of the single-layer analysis point;

[0029] The single-layer characteristic algorithm is used to obtain the single-layer comparison characteristic, and the single-layer characteristic algorithm is as follows: wherein, F is the single-layer comparison characteristic, q ji is the i-th analysis value of the j-th single-layer analysis point in all single-layer analysis points, q jsq is an average value of all analysis values of the j-th single-layer analysis point in all single-layer analysis points.

[0030] Further, real-time data acquisition is performed on the test fault region and the test bottom layer region based on the electrical imaging logging, and the difference analysis model is established based on the data obtained by the real-time data acquisition, and the difference analysis model further comprises:

[0031] The difference analysis model is established, wherein the difference analysis model is used to store a method of drawing curves in the resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system by the acquisition data of the electrical imaging logging, and the coordinate system analysis method and the single-layer curve analysis method are recorded in the difference analysis model;

[0032] The difference analysis model comprises a characteristic comparison library and a single-layer analysis library, wherein the characteristic comparison library is used to store the multi-parameter comparison characteristic, and the single-layer analysis library is used to store the single-layer comparison characteristic.

[0033] Further, the difference analysis model is optimized based on the data obtained by the electrical imaging logging, and the optimized difference analysis model is denoted as a fault comparison model, and the fault comparison model comprises:

[0034] A plurality of groups of test fault regions and test bottom layer regions at different depths are obtained, and based on the acquisition data obtained by the electrical imaging logging, the difference analysis model is used to obtain the multi-parameter comparison characteristic and the single-layer comparison characteristic corresponding to each group of test fault regions and test bottom layer regions;

[0035] When any new multi-parameter comparison characteristic β is obtained, a multi-parameter comparison characteristic with the same name as the multi-parameter comparison characteristic β in the characteristic comparison library is denoted as an original characteristic, and the original characteristic is updated as a union set of the original characteristic and the multi-parameter comparison characteristic β;

[0036] When any one new single layer comparison feature δ is obtained, the single layer comparison features in the single layer analysis library corresponding to the data of the single layer comparison feature δ are recorded as base features, and all the base features and the closed interval formed by the maximum value and the minimum value in the single layer comparison feature δ are recorded as the single layer feature interval of the data corresponding to the single layer comparison feature δ, wherein the single layer feature interval includes the single layer feature interval corresponding to the fault acquisition data and the single layer feature interval corresponding to the bottom layer acquisition data.

[0037] Further, the differentiated analysis model is optimized based on the data obtained by the electrical imaging logging, and the optimized differentiated analysis model is recorded as the fault-bottom layer comparison model.

[0038] When a set of new test fault regions and test bottom layer regions corresponding multi-parameter comparison features and single layer comparison features are obtained by using the differentiated analysis model, and the obtained original feature interval and single layer feature interval do not change, the differentiated analysis model is recorded as the fault-bottom layer comparison model.

[0039] Further, the active fault borehole stratum comparison method is obtained based on the data in the fault-bottom layer comparison model.

[0040] The data in the fault-bottom layer comparison model is analyzed, and the active fault borehole stratum comparison method is obtained based on the analysis result.

[0041] The active fault borehole stratum comparison method is: inputting the existing acquisition data based on the electrical imaging logging into the fault-bottom layer comparison model, and obtaining the multi-parameter comparison features and single layer comparison features respectively, the proportion in the original feature interval and single layer feature interval corresponding to the fault acquisition data and bottom layer acquisition data, and the fault acquisition data or bottom layer acquisition data corresponding to the maximum proportion is taken as the representation data of the acquisition data of the electrical imaging logging.

[0042] The active fault or borehole stratum corresponding to the representation data is recorded as the geological state of the area detected by the electrical imaging logging.

[0043] The application first respectively acquires the active fault and the bottom layer of the drill hole with the same depth from the ground surface, and records them as a test fault area and a test bottom layer area respectively; real-time data acquisition is performed on the test fault area and the test bottom layer area based on the electrical imaging logging, and a differential analysis model is established based on the data obtained through the real-time data acquisition, which has the advantages that the data acquired through the real-time data acquisition of the test fault area and the test bottom layer area is more accurate, which helps to improve the accuracy of the comparison of the active fault and the bottom layer of the drill hole in subsequent analysis; through the establishment of the differential analysis model, the method of data analysis on the acquired data can be stored, and the related features of the active fault and the bottom layer of the drill hole can be stored and further filled, which helps to improve the efficiency of the comparison result based on the acquired data when the active fault and the bottom layer of the drill hole are compared in the future, so as to realize that the feature relationship between the strata can be quickly acquired according to the corresponding related parameters of the strata, and then the comparison method of comparing the strata can be effectively acquired.

[0044] The application also acquires a plurality of test fault areas and test bottom layer areas with different depths, and optimizes the differential analysis model based on the data obtained through the electrical imaging logging, and records the optimized differential analysis model as a ground fault comparison model; finally, an active fault and drill hole stratum comparison method is acquired based on the data in the ground fault comparison model, which has the advantages that through the acquisition of the ground fault comparison model and the active fault and drill hole stratum comparison method, the active fault or drill hole stratum corresponding to the acquired data can be efficiently acquired based on the data analysis method in the ground fault comparison model and the stored multi-parameter comparison features and single-layer comparison features after the acquisition of the acquired data, and the stratum comparison efficiency is improved. BRIEF DESCRIPTION OF DRAWINGS

[0045] Figure 1 The method of the application is a step flowchart;

[0046] Figure 2 The acquisition schematic diagram of the porosity comparison features of the application is shown in the figure;

[0047] Figure 3 The structure schematic diagram of the electronic device of the application is shown in the figure. DETAILED DESCRIPTION

[0048] The technical solutions in the embodiments of the application will be clearly and completely described below with reference to the drawings in the embodiments of the application. Obviously, the described embodiments are only a part of the embodiments of the application, rather than all the embodiments of the application. Based on the embodiments in the application, all other embodiments obtained by those skilled in the art without creative labor fall within the scope of protection of the application.

[0049] Embodiment 1, please refer toFigure 1 As shown, the application provides an active fault drilling formation correlation method combined with high-resolution logging data, characterized in that it comprises the following steps:

[0050] Step S1, respectively acquiring an active fault and a drilling bottom layer with the same distance from the surface depth, and respectively recording them as a test fault area and a test bottom layer area;

[0051] Based on the electrical imaging logging, real-time data acquisition is performed on the test fault area and the test bottom layer area, and a differential analysis model is established based on the data obtained by real-time data acquisition;

[0052] Step S1 comprises: step S101, respectively recording the acquisition data of the test fault area and the test bottom layer area acquired by the electrical imaging logging for real-time data acquisition as fault acquisition data and bottom layer acquisition data;

[0053] Step S102, respectively establishing a resistivity analysis coordinate system, a density analysis coordinate system and a porosity analysis coordinate system, wherein the unit of the X axis of the resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system is m, the unit of the Y axis of the resistivity analysis coordinate system is Ω·m, the unit of the Y axis of the density analysis coordinate system is kg / m³, and the unit of the Y axis of the porosity analysis coordinate system is %;

[0054] Step S103, based on the straight-line distance from the surface when the electrical imaging logging in the fault acquisition data is acquired, and the resistivity, density and porosity in the fault acquisition data, corresponding curves are drawn in the resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system, and are respectively recorded as a fault resistivity curve, a fault density curve and a fault porosity curve;

[0055] In the specific implementation process, the resistivity, density and porosity can be added, deleted or modified according to the data parameters that can be acquired by the electrical imaging logging and that the active fault and the drilling bottom layer exist at the same time, such as increasing the magnetic susceptibility to improve the diversification of data analysis, which helps to improve the accuracy of subsequent formation correlation. In this embodiment, only resistivity, density and porosity are analyzed;

[0056] Step S104, based on the analysis method of the fault acquisition data, the corresponding curves of the bottom layer acquisition data in the resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system are obtained, and are respectively recorded as a bottom layer resistivity curve, a bottom layer density curve and a bottom layer porosity curve;

[0057] Step S105: The resistivity analysis coordinate system, density analysis coordinate system and porosity analysis coordinate system are analyzed sequentially using the coordinate system analysis method, and multi-parameter comparison features are obtained based on the analysis results. The multi-parameter comparison features include resistivity comparison features, density comparison features and porosity comparison features.

[0058] The coordinate system analysis method includes: step S1051, for the coordinate system to be analyzed: the curve corresponding to the fault acquisition data of the two curves in the coordinate system is denoted as curve UC, and the curve corresponding to the bottom acquisition data is denoted as curve IC.

[0059] Step S1052: Record the difference between the maximum and minimum values ​​of the ordinate in curve UC as Y1, and the difference between the maximum and minimum values ​​of the ordinate in curve IC as Y2; record the points with the largest absolute values ​​of the slopes in curves UC and IC as points K1 and K2, respectively.

[0060] In specific implementation processes, for example, during a data analysis, the coordinate system obtained for analysis is the porosity analysis coordinate system, and such as Figure 2 As shown, through analysis, the value of Y1 is 20%, and Y2 is 10%, so the value of α1 in the subsequent analysis is 2; in addition, the absolute values ​​of the slopes of points K1 and K2 are 4 and 2 respectively, so the value of α2 is 2; since α3 is the proportional relationship between L1 and L2, the actual values ​​of L1 and L2 do not affect the acquisition of α3. Figure 2 If the values ​​of L1 and L2 are 5cm and 6cm respectively, then α3 is 0.83, and the porosity comparison characteristics are [2,2,0.83].

[0061] Step S1053: Set the value of α1 to Y1 / Y2; set the value of α2 to the absolute value of the slope of point K1 in curve UC divided by the absolute value of the slope of point K2 in curve IC.

[0062] Step S1054: Set α3 to L1 / L2, where L1 is the distance between the point with the largest ordinate in curve UC and point K1, and L2 is the distance between the point with the largest ordinate in curve IC and point K2.

[0063] Step S1055: When the coordinates to be analyzed are resistivity analysis coordinates, [α1, α2, α3] are recorded as resistivity comparison features; when the coordinates to be analyzed are density analysis coordinates, [α1, α2, α3] are recorded as density comparison features; when the coordinates to be analyzed are porosity analysis coordinates, [α1, α2, α3] are recorded as porosity comparison features.

[0064] Step S106, using single layer curve analysis method to analyze all curves corresponding to fault acquisition data and bottom layer acquisition data, and obtaining single layer comparison features corresponding to active fault and bottom layer of the borehole based on the analysis result;

[0065] The single layer curve analysis method comprises the following steps: in step S1061, for the fault acquisition data or the bottom layer acquisition data: placing the corresponding curves in the resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system into the same coordinate system, and sequentially recording them as single layer analysis curve DF1, single layer analysis curve DF2 and single layer analysis curve DF3;

[0066] In the specific implementation process, since the analysis object in the subsequent analysis is the distance between the points in the single layer analysis curve DF1, the single layer analysis curve DF2 and the single layer analysis curve DF3 and the X axis, the unit of the corresponding ordinate of the single layer analysis curve DF1, the single layer analysis curve DF2 and the single layer analysis curve DF3 in the coordinate system does not need to be considered;

[0067] In step S1062, the interval occupied by the single layer analysis curve DF1, the single layer analysis curve DF2 and the single layer analysis curve DF3 in the X axis is recorded as [P1, P2]; Q points are uniformly obtained in X=P1 to X=P2, and are all recorded as single layer analysis points;

[0068] In the specific implementation process, the value of Q can be set according to the data analysis capability. When the actual data analysis capability is strong, the value of Q can be increased to obtain more accurate single layer comparison features; when the actual data analysis capability is weak, the value of Q can be reduced to reduce the data analysis pressure of the equipment;

[0069] In step S1063, for the horizontal coordinate r of any single layer analysis point, the distance between the point with the horizontal coordinate r in the single layer analysis curve DF1, the single layer analysis curve DF2 and the single layer analysis curve DF3 and the X axis is recorded as the first analysis value, the second analysis value and the third analysis value of the single layer analysis point, respectively;

[0070] In step S1064, the single layer comparison features are obtained by using a single layer feature algorithm, and the single layer feature algorithm is as follows: wherein F is the single layer comparison feature, q ji is the i-th analysis value of the j-th single layer analysis point in all single layer analysis points, q jsq is the average value of all analysis values of the j-th single layer analysis point in all single layer analysis points;21.33344.6667 3.5378

[0071] In the implementation process, for example, in a data analysis, the value of Q is 5, and the first analysis values corresponding to each single-layer analysis point are 2, 4, 2, 3 and 4 respectively, the second analysis values are 6, 5, 7, 6 and 8 respectively, and the third analysis values are 4, 5, 5, 5 and 4 respectively. Through calculation, the single-layer comparison feature is 0.404. In the implementation process, by obtaining the multi-parameter comparison feature and the single-layer comparison feature, the feature relationship between the strata can be quickly obtained according to the related parameters corresponding to the strata when the active fault and the drill hole bottom layer are compared in the subsequent comparison, that is, the relationship between the active fault and the drill hole bottom layer in the multi-parameter comparison feature and the single-layer comparison feature, and then the comparison method for comparing the strata can be effectively obtained.

[0072] Step S1 further includes step S107 of establishing a differential analysis model, wherein the differential analysis model is used to store a method of drawing curves in a resistivity analysis coordinate system, a density analysis coordinate system and a porosity analysis coordinate system by the acquisition data of the electrical imaging logging, and the coordinate system analysis method and the single-layer curve analysis method are recorded in the differential analysis model.

[0073] Step S108, the differential analysis model includes a feature comparison library and a single-layer analysis library, wherein the feature comparison library is used to store the multi-parameter comparison feature, and the single-layer analysis library is used to store the single-layer comparison feature.

[0074] Step S2, a plurality of test fault regions and test bottom layer regions of different depths are obtained, and the differential analysis model is optimized based on the data obtained by the electrical imaging logging, and the optimized differential analysis model is denoted as a ground fault comparison model;

[0075] Step S2 includes step S201 of obtaining a plurality of test fault regions and test bottom layer regions of different depths, and obtaining the multi-parameter comparison feature and the single-layer comparison feature corresponding to each test fault region and test bottom layer region using the differential analysis model based on the acquisition data obtained by the electrical imaging logging;

[0076] Step S202, when any one new multi-parameter comparison feature β is obtained, the multi-parameter comparison feature with the same name as the multi-parameter comparison feature β in the feature comparison library is denoted as an original feature, and the original feature is updated to the union set of the original feature and the multi-parameter comparison feature β;

[0077] In the implementation process, for example, a new porosity comparison feature is [3, 3, 1], and the original feature is [2, 2, 0.83]. The original feature can be updated to [(2, 3), (2, 3), (0.83, 1)];

[0078] Step S203, when any one new single-layer matching feature delta is obtained, the single-layer matching features with the same data as the single-layer matching feature delta corresponding to the single-layer matching feature delta in the single-layer analysis library are recorded as base features, and a closed interval formed by the maximum value and the minimum value in the single-layer matching feature delta is recorded as a single-layer feature interval of the data corresponding to the single-layer matching feature delta, wherein the single-layer feature interval includes a single-layer feature interval corresponding to the fault acquisition data and a single-layer feature interval corresponding to the bottom acquisition data;

[0079] In the specific implementation process, for example, in one data analysis, the single-layer matching feature corresponding to one new fault acquisition data obtained is 0.505, and the base feature is 0.404, so the single-layer feature interval can be set as [0.404, 0.505];

[0080] Step S204, when a plurality of parameter matching features and single-layer matching features corresponding to a group of new test fault regions and test bottom regions are obtained by using the differential analysis model, and the original features and the single-layer feature interval obtained do not change, the differential analysis model is recorded as a fault contrast model.

[0081] Step S3, an active fault borehole formation contrast method is obtained based on the data in the fault contrast model.

[0082] Step S3 includes: step S301, using the data in the fault contrast model for analysis, and obtaining an active fault borehole formation contrast method based on the analysis result;

[0083] Step S302, the active fault borehole formation contrast method is: inputting the existing acquisition data based on the electrical imaging logging into the fault contrast model, and obtaining the plurality of parameter matching features and the single-layer matching features respectively, the proportion in the original features and the single-layer feature interval corresponding to the fault acquisition data and the bottom acquisition data, and taking the fault acquisition data or the bottom acquisition data corresponding to the maximum proportion as the representation data of the acquisition data of the electrical imaging logging;

[0084] In the specific implementation process, for example, the single-layer matching feature based on the acquisition data is 0.504, the single-layer feature intervals corresponding to the fault acquisition data and the bottom acquisition data are [0.404, 0.505] and [0.505, 0.606] respectively, so the proportion of the single-layer matching feature in the single-layer feature interval of the fault acquisition data is 1, and the proportion of the single-layer matching feature in the single-layer feature interval of the bottom acquisition data is 0; if the plurality of parameter matching features corresponding to the acquisition data have the same proportion in the plurality of parameter matching features corresponding to the fault acquisition data and the bottom acquisition data, the fault acquisition data can be taken as the representation data of the acquisition data of the electrical imaging logging;

[0085] Step S303, the active fault or the borehole stratum corresponding to the characterization data is recorded as the geological state of the area detected by the electrical imaging logging.

[0086] Embodiment 2, please refer to Figure 3 as shown, Figure 3 An example of a structural schematic diagram of an electronic device can include: a processor, a communication interface, a memory and a communication bus, wherein the processor, the communication interface, the memory complete the communication among each other through the communication bus. The memory stores computer readable instructions, and the processor can call the instructions in the memory, when the computer readable instructions are executed by the processor, run as steps in a kind of active fault borehole stratum comparison method combined with high-resolution logging data, to realize the following functions: first, respectively obtain the active fault and the borehole bottom layer of the same distance from the ground depth, and respectively record as test fault area and test bottom layer area;Real-time data acquisition is carried out on test fault area and test bottom layer area based on electrical imaging logging, and differential analysis model is established based on the data obtained by real-time data acquisition;Then obtain a plurality of groups of test fault area and test bottom layer area of different depths, and optimize differential analysis model based on the data obtained by electrical imaging logging, and the optimized differential analysis model is recorded as ground fault comparison model;Finally, based on the data in ground fault comparison model, active fault borehole stratum comparison method is obtained.

[0087] In addition, the logical instructions in the memory described above can be realized in the form of software functional units and sold or used as independent products, which can be stored in a computer readable storage medium. Based on such understanding, the technical solutions of the present application essentially or the parts that contribute to the prior art or parts of the technical solutions can be embodied in the form of software products, and the computer software product is stored in a storage medium, including a plurality of instructions to make a computer device (which can be a personal computer, server, or network device, etc.) execute all or part of the steps of the method described in various embodiments of the present application. The foregoing storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), magnetic disk or optical disk and various program code storage media.

[0088] In embodiment 3, the application further provides a computer program product, which comprises a computer program stored on a computer readable storage medium, and the computer program comprises program instructions, when the program instructions are executed by a computer, the computer can execute the method of the above-mentioned method for comparing the active fault and the drilled formation, which comprises the following steps: first, respectively acquiring the active fault and the drilled bottom layer with the same depth from the ground surface, and respectively recording as the test fault area and the test bottom layer area; based on the electrical imaging logging, respectively collecting real-time data of the test fault area and the test bottom layer area, and establishing a differential analysis model based on the data obtained by the real-time data collection; then, acquiring a plurality of groups of test fault areas and test bottom layer areas with different depths, and optimizing the differential analysis model based on the data obtained by the electrical imaging logging, and recording the optimized differential analysis model as a ground fault comparison model; and finally, obtaining the method for comparing the active fault and the drilled formation based on the data in the ground fault comparison model.

[0089] In embodiment 4, the application further provides a computer readable storage medium, and the application provides a storage medium, which stores a computer program, when the computer program is executed by a processor, the steps of the method for comparing the active fault and the drilled formation based on the high-resolution logging data are executed to realize the following functions: first, respectively acquiring the active fault and the drilled bottom layer with the same depth from the ground surface, and respectively recording as the test fault area and the test bottom layer area; based on the electrical imaging logging, respectively collecting real-time data of the test fault area and the test bottom layer area, and establishing a differential analysis model based on the data obtained by the real-time data collection; then, acquiring a plurality of groups of test fault areas and test bottom layer areas with different depths, and optimizing the differential analysis model based on the data obtained by the electrical imaging logging, and recording the optimized differential analysis model as a ground fault comparison model; and finally, obtaining the method for comparing the active fault and the drilled formation based on the data in the ground fault comparison model.

[0090] Through the above description of the embodiments, the embodiments of the application can be provided as a method, a system or a computer program product. Based on such understanding, the above-mentioned technical solutions essentially or say the part of the prior art which makes a contribution can be embodied in the form of a software product, which can be stored in a computer readable storage medium, such as ROM / RAM, a magnetic disk, an optical disk, etc., and comprises a plurality of instructions to make a computer device (which can be a personal computer, a server, or a network device, etc.) execute the method described in each embodiment or some part of the embodiment.

[0091] In the embodiments of the present application, it should be understood that the disclosed system or method can be implemented in other ways. The embodiments described above are only illustrative, for example, the division of modules or units is only a logical function division, and other division manners can be used in actual implementation, for example, a plurality of modules or units can be combined or integrated into another system, or some features can be ignored or not executed. In addition, the coupling or direct coupling or communication connection between the displayed or discussed each other can be through some communication interface, the indirect coupling or communication connection between the system, the module and the unit can be electrical, mechanical or other forms.

[0092] Finally, it should be noted that: the above embodiments are only used to illustrate the technical solutions of the present application, and not to limit them; although the present application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand: it can still modify the technical solutions recorded in the foregoing embodiments, or make equivalent replacement for some technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for correlating active faulted borehole stratigraphy with high resolution logging data, characterized by, The method comprises the following steps: Respectively acquire the active fault and the bottom layer of the drill hole with the same depth from the ground surface, and mark them as the test fault area and the test bottom layer area respectively; Based on the electrical imaging logging, real-time data acquisition is performed on the test fault area and the test bottom layer area, and a differential analysis model is established based on the data obtained by the real-time data acquisition, which specifically comprises: Based on the linear distance from the ground surface when the electrical imaging logging is performed on the fault acquisition data, and the resistivity, density and porosity in the fault acquisition data, corresponding curves are plotted in the resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system, and are marked as the fault resistivity curve, the fault density curve and the fault porosity curve respectively; Based on the analysis method of the fault acquisition data, corresponding curves of the bottom layer acquisition data in the resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system are obtained, and are marked as the bottom layer resistivity curve, the bottom layer density curve and the bottom layer porosity curve respectively; The coordinate system analysis method is used to analyze the resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system in sequence, and multi-parameter comparison features are obtained based on the analysis results, wherein the multi-parameter comparison features include resistivity comparison features, density comparison features and porosity comparison features; The single-layer curve analysis method is used to analyze all the curves corresponding to the fault acquisition data and the bottom layer acquisition data, and single-layer comparison features corresponding to the active fault and the bottom layer of the drill hole are obtained based on the analysis results; The single-layer curve analysis method comprises: For the fault acquisition data or the bottom layer acquisition data, the corresponding curves in the resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system are placed in the same coordinate system, and are sequentially marked as single-layer analysis curve DF1, single-layer analysis curve DF2 and single-layer analysis curve DF3; The interval occupied by DF1, DF2 and DF3 in the X-axis is marked as [P1, P2]; Q points are uniformly obtained in X=P1 to X=P2, and are all marked as single-layer analysis points; For the horizontal coordinate r of any single-layer analysis point, the distances between the points with the horizontal coordinate r in DF1, DF2 and DF3 and the X-axis are respectively marked as the first analysis value, the second analysis value and the third analysis value of the single-layer analysis point; The single-layer comparison features are acquired by using a single-layer feature algorithm, and the single-layer feature algorithm is: wherein F is a single-layer comparison feature, q ji is the i-th analysis value of the j-th single-layer analysis point in all single-layer analysis points, q jsq is an average value of all analysis values of the j-th single-layer analysis point in all single-layer analysis points. A plurality of test fault areas and test bottom layer areas with different depths are acquired, and the differential analysis model is optimized based on the data obtained by the electrical imaging logging, and the optimized differential analysis model is marked as the fault comparison model; The existing acquisition data based on the electrical imaging logging is input into the fault comparison model, and the multi-parameter comparison features and the single-layer comparison features obtained are respectively acquired, and the proportion in the original feature and single-layer feature interval corresponding to the fault acquisition data or the bottom layer acquisition data is calculated, and the fault acquisition data or the bottom layer acquisition data corresponding to the maximum proportion is taken as the representation data of the acquisition data of the electrical imaging logging; The active fault or the bottom layer of the drill hole corresponding to the representation data is marked as the geological state of the area detected by the electrical imaging logging.

2. The method of claim 1, wherein the method is characterized by, The real-time data collection is performed on the test fault area and the test bottom layer area based on the electrical imaging logging, and a differential analysis model is established based on the data obtained by the real-time data collection, and the differential analysis model includes: The collected data of the real-time data collection performed on the test fault area and the test bottom layer area by the electrical imaging logging is recorded as fault collected data and bottom layer collected data, respectively; A resistivity analysis coordinate system, a density analysis coordinate system and a porosity analysis coordinate system are established, respectively, wherein the unit of the X axis of the resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system is m, the unit of the Y axis of the resistivity analysis coordinate system is Ω·m, the unit of the Y axis of the density analysis coordinate system is kg / m³, and the unit of the Y axis of the porosity analysis coordinate system is %.

3. The method of claim 2, wherein the method further comprises: The coordinate system analysis method includes: A curve corresponding to the fault collected data of two curves in the coordinate system is recorded as UC, and a curve corresponding to the bottom layer collected data is recorded as IC; The difference between the maximum value and the minimum value of the ordinate in UC is recorded as Y1, and the difference between the maximum value and the minimum value of the ordinate in IC is recorded as Y2; the point with the maximum absolute value of the slope in UC is recorded as point K1, and the point with the maximum absolute value of the slope in IC is recorded as point K2; The value of α1 is set as Y1 / Y2, and the value of α2 is set as the absolute value of the slope of point K1 in UC divided by the absolute value of the slope of point K2 in IC; α3 is set as L1 / L2, wherein L1 is the distance between the point with the maximum ordinate in UC and point K1, and L2 is the distance between the point with the maximum ordinate in IC and point K2; According to the differences of the coordinate systems, [α1, α2, α3] is recorded as resistivity contrast features or density contrast features or porosity contrast features.

4. The method of claim 3, wherein the method further comprises: The real-time data collection is performed on the test fault area and the test bottom layer area based on the electrical imaging logging, and a differential analysis model is established based on the data obtained by the real-time data collection, and the differential analysis model further includes: The differential analysis model is established, wherein the differential analysis model is used to store the method of drawing curves in the resistivity analysis coordinate system, the density analysis coordinate system and the porosity analysis coordinate system by the collected data of the electrical imaging logging, and the coordinate system analysis method and the single layer curve analysis method are recorded in the differential analysis model; The differential analysis model includes a feature contrast library and a single layer analysis library, wherein the feature contrast library is used to store multi-parameter contrast features, and the single layer analysis library is used to store single layer contrast features.

5. The method of claim 4, wherein the method further comprises: The differential analysis model is optimized based on the data obtained by the electrical imaging logging, and the optimized differential analysis model is recorded as a fault contrast model, and the fault contrast model includes: A plurality of groups of test fault areas and test bottom layer areas at different depths are obtained, and based on the collected data obtained by the electrical imaging logging, the differential analysis model is used to obtain multi-parameter comparison features and single layer comparison features corresponding to each group of test fault areas and test bottom layer areas; When any new multi-parameter comparison feature β is obtained, a multi-parameter comparison feature with the same name as the multi-parameter comparison feature β in the feature comparison library is recorded as an original feature, and the original feature is updated as the union of the original feature and the multi-parameter comparison feature β; When any one new single-layer comparison feature δ is obtained, the single-layer comparison features with the same data as the single-layer comparison feature δ in the single-layer analysis library are recorded as base features, and all the base features and the closed interval formed by the maximum value and the minimum value in the single-layer comparison feature δ are recorded as a single-layer feature interval of the data corresponding to the single-layer comparison feature δ, wherein the single-layer feature interval includes a single-layer feature interval corresponding to the fault acquisition data and a single-layer feature interval corresponding to the bottom acquisition data.

6. The method of claim 5, wherein the method further comprises: The differential analysis model is optimized based on the data obtained by the electrical imaging logging, and the optimized differential analysis model is recorded as a fault contrast model. When a group of new test fault regions and test bottom regions are obtained by using the differential analysis model, and the multi-parameter comparison features and the single-layer comparison features corresponding to the test fault regions and the test bottom regions are obtained, and the original features and the single-layer feature intervals obtained do not change, the differential analysis model is recorded as the fault contrast model.

Images