A method, apparatus, equipment, and medium for determining the area of ​​crack current anomalies.

Abstract

This invention relates to the field of oil and gas exploration and development, and discloses a method, apparatus, equipment, and medium for determining the area of ​​fracture current anomalies. It can generate conductivity imaging data and trajectory curve expressions for multiple fractures based on raw electrical imaging logging data from fractured oil and gas reservoirs. The conductivity imaging data and each trajectory curve expression are mapped to a predetermined three-dimensional coordinate system to generate a fracture trajectory curve corresponding to each trajectory curve expression and multiple discrete data points corresponding to the conductivity imaging data in the three-dimensional coordinate system. The three-dimensional coordinate system includes a depth axis, a wellbore expansion direction axis, and a conductivity axis. The background conductivity of the three-dimensional coordinate system is determined based on each discrete data point. Based on the background conductivity and each fracture trajectory curve, the current anomaly area of ​​each fracture is determined. This invention can effectively determine the area of ​​fracture current anomalies, enriching the methods for determining fracture current anomalies.

Description

Technical Field

[0001] This invention relates to the field of oil and gas exploration and development, and in particular to a method, apparatus, equipment and medium for determining the area of ​​fracture current anomalies. Background Technology

[0002] With the deepening of exploration and development, fractured oil and gas reservoirs, represented by buried hills and tight sandstone, have become an important area for increasing oil and gas reserves and production. The width and distribution characteristics of fractures have a decisive impact on the estimation of oil and gas reserves, production forecasting, and development strategies. Fracture width not only directly affects the physical properties of the formation but also plays a crucial role in the permeability and reservoir performance of oil and gas. Therefore, accurately evaluating fracture width is a core task for achieving efficient exploration and development of fractured oil and gas reservoirs.

[0003] Electrical imaging logging, with its high resolution and high precision, can intuitively reflect the distribution of underground fractures and has become the most commonly used method for evaluating fracture width. Among these methods, the fracture current anomaly area is a key parameter for evaluating fracture width based on electrical imaging logging.

[0004] However, domestic research in this field remains relatively weak, resulting in reliance on foreign logging interpretation software for accurate calculation of fracture current anomaly area and fracture width, which restricts the precise logging evaluation of fractured oil and gas reservoirs. Therefore, breaking through key technological bottlenecks and independently developing quantitative calculation methods for fracture current anomaly area is crucial for the efficient exploration and development of fractured oil and gas reservoirs. Summary of the Invention

[0005] This invention provides a method, apparatus, device, and medium for determining the area of ​​crack current anomalies, thereby addressing the deficiency of limited methods for determining the area of ​​crack current anomalies in related technologies and enriching the methods for determining the area of ​​crack current anomalies.

[0006] In a first aspect, the present invention provides a method for determining the area of ​​crack current anomalies, comprising: Based on the raw electrical imaging logging data of fractured oil and gas reservoirs, electrical conductivity imaging data and trajectory curve expressions for multiple fractures are generated. The conductivity imaging data and each of the trajectory curve expressions are mapped to a set three-dimensional coordinate system to generate a fracture trajectory curve corresponding to each trajectory curve expression and multiple discrete data points corresponding to the conductivity imaging data in the three-dimensional coordinate system; wherein, the three-dimensional coordinate system includes a depth axis, a wellbore development direction axis and a conductivity axis; The background conductivity of the three-dimensional coordinate system is determined based on each of the discrete data points. Based on the background conductivity and the trajectory curve of each crack, the current anomaly area of ​​each crack is determined.

[0007] Optionally, the step of generating electrical conductivity imaging data and trajectory curve expressions for multiple fractures based on the raw electrical imaging logging data of fractured oil and gas reservoirs includes: The raw data from electrical imaging logging of fractured oil and gas reservoirs are preprocessed to generate dynamic and static image data. The static image data is calibrated for conductivity to obtain the conductivity imaging data. An imaging operation is performed on the dynamic image data, the static image data, or the conductivity imaging data to obtain the corresponding target image; Crack identification and feature extraction are performed on the target image to obtain feature data for each crack, including crack depth, orientation, and dip angle. Based on the characteristic data of each crack, a trajectory curve expression for each crack is generated.

[0008] Optionally, determining the background conductivity of the three-dimensional coordinate system based on each of the discrete data points includes: The coordinate interval on the wellbore unfolding direction axis is divided into multiple sub-intervals; For any of the sub-intervals, each target data point is found among the plurality of discrete data points, and the depth coordinate value of each target data point is located within the sub-interval. A dynamic threshold is generated based on the conductivity coordinate value of each target data point. Each background data point is determined among each target data point, and the conductivity coordinate value of each background data point is less than the dynamic threshold. The average conductivity coordinate value of each background data point is calculated to obtain the calculation result, which is used as the background conductivity corresponding to the sub-interval. The background conductivity of each sub-interval is taken as the background conductivity of the three-dimensional coordinate system.

[0009] Optionally, generating a dynamic threshold based on the conductivity coordinate value of each target data point includes: The median of the conductivity coordinates for each target data point is obtained by statistically analyzing the median coordinate values. Calculate the absolute deviation between the conductivity coordinate value and the median of the coordinate values ​​for each of the target data points; For each absolute deviation value, the median absolute deviation is obtained by performing median statistics. The median absolute deviation is multiplied by a set factor to obtain a multiplication result, and the multiplication result is added to the median coordinate value to obtain the dynamic threshold.

[0010] Optionally, determining the current anomaly area of ​​each crack based on the background conductivity and the trajectory curve of each crack includes: For any of the fracture trajectory curves, data points are sampled along the wellbore unfolding direction at a set sampling interval to determine multiple sampling points; For any sampling point on the crack trajectory curve, a target data range is determined based on the depth coordinate value of the sampling point. The absolute value of the difference between the depth coordinate value of any point in the target data range and the depth coordinate value of the sampling point is not greater than a set threshold. Among each discrete data point included in the target data range, the discrete data point with the largest conductivity coordinate value is determined and used as a candidate data point. Each valid data point is determined among each candidate data point. The current anomaly area of ​​the crack is determined based on the conductivity coordinate value of each valid data point.

[0011] Optionally, determining each valid data point among each candidate data point includes: For any candidate data point, with the candidate data point as the center, draw a circle with a set radius on the plane enclosed by the depth axis and the wellbore unfolding direction axis of the three-dimensional coordinate system, and search for whether there is a blank strip with a conductivity coordinate value of 0 in the drawn circle. If the blank stripe exists in the drawn circular surface, then the candidate data point is determined to be a non-valid data point; If no blank stripe exists in the drawn circular surface, then the candidate data point is determined to be a valid data point.

[0012] Optionally, determining the current anomaly area of ​​the crack based on the conductivity coordinate value of each valid data point includes: For any valid data point on the crack trajectory curve, the circular surface drawn with the valid data point as the center is determined as the target integration region. The target variable is used as the integration object, which is the difference between the conductivity coordinate value of the coordinate point and the background conductivity corresponding to the coordinate point. The integration object in the target integration region is integrated to obtain the integration result. The integration result is multiplied by a set instrument coefficient to obtain the multiplication result, which is used as the current anomaly area corresponding to the valid data point. The current anomaly area corresponding to each valid data point on the crack trajectory curve is taken as the current anomaly area of ​​the crack corresponding to the crack trajectory curve.

[0013] In a second aspect, the present invention provides a crack current anomaly area determination device, applied to the crack current anomaly area determination method of the first aspect or any corresponding embodiment thereof, the device comprising: The first generation unit is used to generate electrical conductivity imaging data and trajectory curve expressions of multiple fractures based on the raw electrical imaging logging data of fractured oil and gas reservoirs. The second generation unit is used to map the conductivity imaging data and each of the trajectory curve expressions to a set three-dimensional coordinate system, so as to generate a fracture trajectory curve corresponding to each trajectory curve expression and multiple discrete data points corresponding to the conductivity imaging data in the three-dimensional coordinate system; wherein, the three-dimensional coordinate system includes a depth axis, a wellbore development direction axis and a conductivity axis; The first determining unit is used to determine the background conductivity of the three-dimensional coordinate system based on each of the discrete data points; The second determining unit is used to determine the current anomaly area of ​​each crack based on the background conductivity and the trajectory curve of each crack.

[0014] Thirdly, the present invention provides a computer device, comprising: a memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing computer instructions, and the processor executing the computer instructions to perform the crack current anomaly area determination method of the first aspect or any corresponding embodiment described above.

[0015] Fourthly, the present invention provides a computer-readable storage medium storing computer instructions for causing a computer to execute the crack current anomaly area determination method of the first aspect or any corresponding embodiment described above.

[0016] The present invention provides a method, apparatus, equipment, and medium for determining fracture current anomaly area. Based on raw electrical imaging logging data from fractured oil and gas reservoirs, it generates conductivity imaging data and trajectory curve expressions for multiple fractures. The conductivity imaging data and each trajectory curve expression are mapped to a defined three-dimensional coordinate system. This generates a fracture trajectory curve corresponding to each trajectory curve expression and multiple discrete data points corresponding to the conductivity imaging data within the three-dimensional coordinate system. The three-dimensional coordinate system includes a depth axis, a wellbore expansion direction axis, and a conductivity axis. The background conductivity of the three-dimensional coordinate system is determined based on each discrete data point. Based on the background conductivity and each fracture trajectory curve, the current anomaly area of ​​each fracture is determined. This invention can effectively achieve accurate quantitative evaluation of fracture current anomaly area, supporting efficient exploration and development of fractured oil and gas reservoirs. Attached Figure Description

[0017] To more clearly illustrate the technical solutions in this invention or related technologies, the accompanying drawings used in the description of the embodiments or related technologies will be briefly introduced below. Obviously, the accompanying drawings described below are some embodiments of this invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0018] Figure 1A flowchart illustrating a method for determining the fracture current anomaly area in a fractured oil and gas reservoir, provided as an embodiment of the present invention; Figure 2 This is a schematic diagram of picking up crack feature data provided in an embodiment of the present invention; Figure 3 A three-dimensional schematic diagram of electrical conductivity data at a crack provided in an embodiment of the present invention; Figure 4 A schematic diagram of discrete data points and crack trajectory curves provided in an embodiment of the present invention; Figure 5 This is a schematic diagram illustrating the effect of a dynamic threshold method for selecting a background region, provided in an embodiment of the present invention. Figure 6 A three-dimensional schematic diagram of a sampling circle and the data within the circle is provided as an embodiment of the present invention; Figure 7 A comparison chart of current anomaly area calculation results provided in an embodiment of the present invention; Figure 8 This is a schematic diagram of a device for determining the fracture current anomaly area of ​​a fractured oil and gas reservoir, provided in an embodiment of the present invention. Figure 9 This is a schematic diagram of the structure of a computer device provided in an embodiment of the present invention. Detailed Implementation

[0019] To make the objectives, technical solutions, and advantages of this invention clearer, the technical solutions of this invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of this invention. All other embodiments obtained by those skilled in the art based on the embodiments of this invention without creative effort are within the scope of protection of this invention.

[0020] The following is combined Figures 1-7 The present invention describes a method for determining the area of ​​crack current anomalies.

[0021] like Figure 1 As shown in this embodiment, a method for determining the area of ​​crack current anomalies is proposed. This method may include the following steps: S101. Based on the original electrical imaging logging data of fractured oil and gas reservoirs, generate electrical conductivity imaging data and trajectory curve expressions for multiple fractures.

[0022] Among them, the raw data of electrical imaging logging can be the data obtained by using downhole instruments to perform downhole electrical imaging logging on a well in a fractured oil and gas reservoir.

[0023] Specifically, the fracture refers to a fracture in a well within a fractured oil and gas reservoir.

[0024] Specifically, this embodiment can generate electrical conductivity imaging data and trajectory curve expressions of multiple fractures in fractured oil and gas reservoirs based on the original electrical imaging logging data.

[0025] Optionally, step S101 includes: The raw data from electrical imaging logging of fractured oil and gas reservoirs are preprocessed to generate dynamic and static image data. Conductivity calibration is performed on static image data to obtain conductivity imaging data; Imaging operations are performed on dynamic image data, static image data, or conductivity imaging data to obtain the corresponding target image; Crack identification and feature extraction are performed on the target image to obtain feature data for each crack, including crack depth, orientation, and dip angle. Based on the characteristic data of each crack, generate the trajectory curve expression for each crack.

[0026] Specifically, this embodiment can perform preprocessing on the raw data of electro-imaging logging of fractured oil and gas reservoirs, including acceleration correction, depth alignment and rearrangement, data normalization, emission voltage correction, and bad electrode correction, to obtain preprocessed data. The preprocessed data is then imaged and represented using images of different color levels to obtain dynamic and static images of electro-imaging logging. The image data is then resistivity calibrated, that is, the preprocessed data is calibrated using shallow resistivity data from conventional logging data to obtain data with the physical meaning of electrical conductivity. The longitudinal sampling interval of this data is 0.1 in, and there are 360 ​​data points around the well laterally, including some blank strip points.

[0027] Specifically, this embodiment can pick up fractures from dynamic and static images or calibrated images of electrical imaging logging, and extract characteristic data of the fractures such as well diameter, depth, azimuth, and dip angle, such as... Figure 2 As shown, after picking up the crack in the electro-imaging image, a tadpole diagram can be obtained. The vertical axis of the tadpole is the depth of the crack center point, the horizontal axis is the dip angle of the crack, and the direction of the tadpole indicates the azimuth angle of the crack.

[0028] Specifically, in this embodiment, the calibrated electrical imaging logging data of the well can be used as the bottom layer data for calculating the fracture width and coordinates can be set for the data points. The vertical coordinate is the depth of the data point, and the horizontal coordinate interval is calculated according to the well diameter data using formula (1).

[0029] -------------Formula (1); in, The horizontal sampling interval is expressed in meters. This represents the well diameter at the corresponding depth, in meters. The formula means dividing the well circumference equally among 360 data points.

[0030] No. n The x-coordinates of the data points are: ----------Formula (2); in, These are the coordinates for the wellbore deployment direction. Specifically, in this embodiment, the characteristic data of the crack can be integrated and formulated to obtain the corresponding trajectory curve expression, as shown in the following formula.

[0031] -------------Formula (3); in, This is the trajectory curve expression, in meters; The coordinates for the wellbore deployment direction are in meters. The depth at the center of the crack is in meters. The azimuth angle of the crack is in degrees. The angle of inclination of the crack is expressed in degrees. This refers to the well diameter, expressed in meters.

[0032] S102. Map the conductivity imaging data and each trajectory curve expression to a set three-dimensional coordinate system to generate the fracture trajectory curve corresponding to each trajectory curve expression and multiple discrete data points corresponding to the conductivity imaging data in the three-dimensional coordinate system; wherein, the three-dimensional coordinate system includes the depth axis, the wellbore development direction axis and the conductivity axis.

[0033] This embodiment can project conductivity imaging data and the expression for each trajectory curve into the same three-dimensional coordinate system, generating crack trajectory curves and multiple discrete data points, such as... Figure 3 The image shown is a three-dimensional image of conductivity imaging data near a certain fracture. It should be noted that, on the plane of the depth axis and the wellbore expansion direction axis (i.e., the horizontal axis), the fracture trajectory curve and the lattice formed by multiple discrete data points can be represented as follows: Figure 4 As shown, the dot matrix (sparsed for ease of display) represents conductivity imaging data, and the sine and cosine curves represent crack trajectory curves.

[0034] S103. Determine the background conductivity of the three-dimensional coordinate system based on each discrete data point.

[0035] It should be noted that the conductivity of the crack portion in the conductivity imaging data is significantly higher, while the background conductivity is significantly lower. Figure 3As shown, this is a three-dimensional image of conductivity imaging data near a crack. It can be seen that the conductivity at the crack is significantly higher than the background conductivity. In this embodiment, a threshold can be set to define the background area and the non-background area. Then, the conductivity data of the background area can be extracted and averaged to obtain the background conductivity.

[0036] Optionally, step S103 includes: The coordinate interval on the wellbore unfolding direction axis is divided into multiple sub-intervals; For any sub-interval, each target data point is found among multiple discrete data points. The depth coordinate value of each target data point is located within the sub-interval. A dynamic threshold is generated based on the conductivity coordinate value of each target data point. Each background data point is determined in each target data point. The conductivity coordinate value of each background data point is less than the dynamic threshold. The average conductivity coordinate value of each background data point is calculated to obtain the calculation result, which is used as the background conductivity corresponding to the sub-interval. The background conductivity of each sub-interval is taken as the background conductivity of the three-dimensional coordinate system.

[0037] Optionally, the above-mentioned generation of a dynamic threshold based on the conductivity coordinates of each target data point includes: The median of the conductivity coordinates for each target data point is obtained by statistically analyzing the median coordinate values. Calculate the absolute deviation between the conductivity coordinate value and the median of the coordinate values ​​for each target data point; For each absolute deviation value, the median is calculated to obtain the median absolute deviation. The median absolute deviation is multiplied by a set factor to obtain the multiplication result, and then the multiplication result is added to the median of the coordinate values ​​to obtain the dynamic threshold.

[0038] Specifically, in this embodiment, the coordinate interval on the wellbore unfolding direction axis can be divided into 36 sub-intervals. 10 columns of data are taken as a group, and 360 columns of data can be divided into 36 groups. Each group uses the calculated value in formula (4) as the dynamic threshold. Data less than the dynamic threshold is taken as background data, and data greater than the dynamic threshold is taken as non-background data. The average conductivity of the background data is determined as the background conductivity.

[0039] ----------Formula (4); in, The median of the data. As a multiple factor, This represents the absolute median difference. For a one-dimensional dataset... ,That Defined as the median of the dataset relative to each data point in the dataset. The median of the absolute deviation, i.e.: ; This is a multiplier factor, which controls the relative magnitude of the dynamic threshold. Choosing an appropriate multiplier factor is crucial. Afterwards, a more accurate background conductivity can be obtained. Through multiple sets of data tests, when... When the value is -0.6, a relatively reasonable background area division result can be obtained for most cracks, such as... Figure 5 As shown, The image shows the effect of distinguishing the background area of ​​a fracture in Well B at -0.6. The left image is the fracture electro-imaging image, and the right image shows the dark area as the non-background area and the light area as the background area. It can be seen that the light area is the area that is not affected by the fracture, and the distinction effect is better.

[0040] S104. Based on the background conductivity and the trajectory curve of each crack, determine the current anomaly area of ​​each crack.

[0041] Specifically, in this embodiment, the current anomaly area of ​​each crack can be calculated based on the background conductivity of each sub-interval and the trajectory curve of each crack.

[0042] Optionally, step S104 includes: For any fracture trajectory curve, data points are sampled along the wellbore development direction at a set sampling interval to determine multiple sampling points; For any sampling point on the crack trajectory curve, the target data range is determined based on the depth coordinate value of the sampling point. The absolute value of the difference between the depth coordinate value of any point in the target data range and the depth coordinate value of the sampling point is not greater than a set threshold. Among each discrete data point included in the target data range, the discrete data point with the largest conductivity coordinate value is determined and used as a candidate data point. Each valid data point is determined among each candidate data point, and the current anomaly area of ​​the crack is determined based on the conductivity coordinate value of each valid data point.

[0043] Optionally, the above process of determining each valid data point among each candidate data point includes: For any candidate data point, with the candidate data point as the center, draw a circle with a set radius on the plane enclosed by the depth axis and the wellbore unfolding direction axis of the three-dimensional coordinate system, and search for whether there is a blank strip with a conductivity coordinate value of 0 in the drawn circle. If blank stripes exist in the drawn circle, the candidate data point is determined to be a invalid data point; If there are no blank stripes in the drawn circle, then the candidate data point is determined to be a valid data point.

[0044] It should be noted that due to the limited wellbore coverage of electrical imaging logging instruments, some wellbore areas cannot be measured during instrument measurement. Therefore, data for these areas is missing during imaging, and they will appear as several blank stripes. In the 3D schematic diagram of fracture conductivity, this is shown as a long rectangle with a conductivity of 0, truncating the fracture.

[0045] Optionally, the determination of the current anomaly area of ​​the crack based on the conductivity coordinate value of each valid data point includes: For any valid data point on the crack trajectory curve, the circle drawn with the valid data point as the center is defined as the target integration region. The target variable is used as the integration object, which is the difference between the conductivity coordinate value of the coordinate point and the background conductivity corresponding to the coordinate point. The integration object in the target integration region is integrated to obtain the integration result. The integration result is multiplied by the set instrument coefficient to obtain the multiplication result, which is used as the current anomaly area corresponding to the valid data point. The current anomaly area corresponding to each valid data point on the crack trajectory curve is taken as the current anomaly area of ​​the crack corresponding to the crack trajectory curve.

[0046] Specifically, in this embodiment, for a given crack, sampling points are selected at 0.01m intervals along its crack trajectory curve. Since the electrical imaging data is discrete, the selected sampling points do not fall on specific data points. Therefore, it is necessary to search for the nearest data point near each sampling point. Because the actual crack trajectory does not completely coincide with the calculated crack trajectory curve, the point with the maximum conductivity within 5cm above and below the nearest data point is searched and used as a candidate data point. A circle is drawn on the plane of depth and horizontal coordinate with the candidate data point as the center and a radius of 0.02m. It is determined whether the data within the circle contains blank bands. If it does, the point is discarded, and the remaining data points are valid data points. A two-dimensional numerical integration is performed on the data within the circle corresponding to the valid data point by subtracting the background conductivity, thus obtaining the integral result corresponding to the valid data point. For example... Figure 6 As shown, in this embodiment, a circle can be drawn at the valid data points, and numerical integration can be performed on the data within the circle.

[0047] Specifically, this embodiment can create and use a formula to implement the above-mentioned two-dimensional numerical integration, in order to calculate the integration result corresponding to the valid data points. The formula is: -------Formula (5); -------Formula (6); in, For valid data points Abnormal conductivity, and These are the coordinates of the wellbore unfolding direction axis and depth axis in a three-dimensional coordinate system, respectively. To set the radius, For data points conductivity, This represents the background conductivity of the corresponding sub-interval. For instrument coefficients, The current anomaly area is the valid data point.

[0048] It should be noted that formulas (5) and (6) in this embodiment are improvements made by the inventor to the original Luthi-Souhaite two-dimensional integral formula. Specifically, the original Luthi-Souhaite two-dimensional integral formula is as follows: ; ; in, The width of the crack. For the area of ​​current anomaly, and These are constants related to the instrument. The resistivity of the drilling fluid. To flush the resistivity zone. This represents the plate potential value. The current value is obtained from electrical imaging logging. This represents the background current value.

[0049] The inventors of this invention analyzed actual electrical imaging logging raw data and discovered that the electrode potential value The current anomaly area varies at different depths; therefore, when calculating the current anomaly area of ​​an inclined crack, the depths at different integration points are different. Different; when integrating at the same integration point, The variation also exists over the integration interval. Therefore, the plate potential value needs to be considered when calculating the current anomaly area. The change in plate potential. The influence of depth variation is used to transform and derive the Luthi-Souhaite formula. First, the depth variation is transformed and derived. Move it inside the integral sign, then apply the conductance formula. G=I / U Substituting the values, we get: ; Further combining conductivity With conductivity G The relationship, that is From this, we can obtain the above formula (5).

[0050] Specifically, in this embodiment, the current anomaly area of ​​each valid data point on the crack trajectory curve can be determined, and the current anomaly area of ​​all valid data points on the crack trajectory curve can be taken as the current anomaly area of ​​the crack corresponding to the crack trajectory curve.

[0051] Taking the calculation results of 5 fractures in Well B as an example, this embodiment calculates the fracture current anomaly area at different effective data points for each fracture. Since the correctness of the calculated current anomaly area cannot be directly verified, the corresponding fracture width is further calculated, and the average fracture width of each fracture is obtained by taking the mean value. Simultaneously, the average width of each fracture is calculated using relevant software, and the relative error of the calculation results in this embodiment is calculated based on the software's calculation results to verify the correctness of the calculation results in this embodiment. Figure 7 As shown, the average crack width calculated in this embodiment is very close to the result calculated by the software, with a relative error of less than 15%, indicating that the crack current anomaly area calculated in this embodiment is highly accurate.

[0052] This embodiment uses calibrated electrical imaging logging data and an innovative fracture trace characterization method to obtain the geometric properties of fractures. Based on this, effective calculation points are selected, and a two-dimensional integral method is used to accurately calculate the area of ​​fracture current anomalies.

[0053] The method for determining the fracture current anomaly area proposed in this embodiment can generate conductivity imaging data and trajectory curve expressions for multiple fractures based on the raw data from electrical imaging logging of fractured oil and gas reservoirs. The conductivity imaging data and each trajectory curve expression are mapped to a defined three-dimensional coordinate system to generate a fracture trajectory curve corresponding to each trajectory curve expression and multiple discrete data points corresponding to the conductivity imaging data in the three-dimensional coordinate system. The three-dimensional coordinate system includes a depth axis, a wellbore expansion direction axis, and a conductivity axis. The background conductivity of the three-dimensional coordinate system is determined based on each discrete data point. Based on the background conductivity and each fracture trajectory curve, the current anomaly area of ​​each fracture is determined. This embodiment can effectively determine the fracture current anomaly area, ensure the accuracy of the fracture current anomaly area calculation, and enrich the methods for determining the fracture current anomaly area.

[0054] like Figure 8 As shown, this embodiment proposes a device for determining the area of ​​crack current anomalies, applicable to any of the aforementioned methods for determining the area of ​​crack current anomalies. This device may include: The first generation unit 101 is used to generate electrical conductivity imaging data and trajectory curve expressions of multiple fractures based on the original electrical imaging logging data of fractured oil and gas reservoirs. The second generation unit 102 is used to map the conductivity imaging data and each trajectory curve expression to a set three-dimensional coordinate system, so as to generate the fracture trajectory curve corresponding to each trajectory curve expression and multiple discrete data points corresponding to the conductivity imaging data in the three-dimensional coordinate system; wherein, the three-dimensional coordinate system includes a depth axis, a wellbore development direction axis and a conductivity axis; The first determining unit 103 is used to determine the background conductivity of the three-dimensional coordinate system based on each discrete data point; The second determining unit 104 is used to determine the current anomaly area of ​​each crack based on the background conductivity and the trajectory curve of each crack.

[0055] It should be noted that the processing procedures of the first generation unit 101, the second generation unit 102, the first determination unit 103, and the second determination unit 104, and their beneficial effects, can be referred to respectively. Figure 1 Steps S101 to S104 in the process will not be described again.

[0056] Optionally, the first generating unit 101 is also used for: The raw data from electrical imaging logging of fractured oil and gas reservoirs are preprocessed to generate dynamic and static image data. Conductivity calibration is performed on static image data to obtain conductivity imaging data; Imaging operations are performed on dynamic image data, static image data, or conductivity imaging data to obtain the corresponding target image; Crack identification and feature extraction are performed on the target image to obtain feature data for each crack, including crack depth, orientation, and dip angle. Based on the characteristic data of each crack, generate the trajectory curve expression for each crack.

[0057] Optionally, the first determining unit 103 is also used for: The coordinate interval on the wellbore unfolding direction axis is divided into multiple sub-intervals; For any sub-interval, each target data point is found among multiple discrete data points. The depth coordinate value of each target data point is located within the sub-interval. A dynamic threshold is generated based on the conductivity coordinate value of each target data point. Each background data point is determined in each target data point. The conductivity coordinate value of each background data point is less than the dynamic threshold. The average conductivity coordinate value of each background data point is calculated to obtain the calculation result, which is used as the background conductivity corresponding to the sub-interval. The background conductivity of each sub-interval is taken as the background conductivity of the three-dimensional coordinate system.

[0058] Optionally, the first determining unit 103 is also used for: The median of the conductivity coordinates for each target data point is obtained by statistically analyzing the median coordinate values. Calculate the absolute deviation between the conductivity coordinate value and the median of the coordinate values ​​for each target data point; For each absolute deviation value, the median is calculated to obtain the median absolute deviation. The median absolute deviation is multiplied by a set factor to obtain the multiplication result, and then the multiplication result is added to the median of the coordinate values ​​to obtain the dynamic threshold.

[0059] Optionally, the second determining unit 104 is also used for: For any fracture trajectory curve, data points are sampled along the wellbore development direction at a set sampling interval to determine multiple sampling points; For any sampling point on the crack trajectory curve, the target data range is determined based on the depth coordinate value of the sampling point. The absolute value of the difference between the depth coordinate value of any point in the target data range and the depth coordinate value of the sampling point is not greater than a set threshold. Among each discrete data point included in the target data range, the discrete data point with the largest conductivity coordinate value is determined and used as a candidate data point. Each valid data point is determined among each candidate data point, and the current anomaly area of ​​the crack is determined based on the conductivity coordinate value of each valid data point.

[0060] Optionally, the second determining unit 104 is also used for: For any candidate data point, with the candidate data point as the center, draw a circle with a set radius on the plane enclosed by the depth axis and the wellbore unfolding direction axis of the three-dimensional coordinate system, and search for whether there is a blank strip with a conductivity coordinate value of 0 in the drawn circle. If blank stripes exist in the drawn circle, the candidate data point is determined to be a invalid data point; If there are no blank stripes in the drawn circle, then the candidate data point is determined to be a valid data point.

[0061] Optionally, the second determining unit 104 is also used for: For any valid data point on the crack trajectory curve, the circle drawn with the valid data point as the center is defined as the target integration region. The target variable is used as the integration object, which is the difference between the conductivity coordinate value of the coordinate point and the background conductivity corresponding to the coordinate point. The integration object in the target integration region is integrated to obtain the integration result. The integration result is multiplied by the set instrument coefficient to obtain the multiplication result, which is used as the current anomaly area corresponding to the valid data point. The current anomaly area corresponding to each valid data point on the crack trajectory curve is taken as the current anomaly area of ​​the crack corresponding to the crack trajectory curve.

[0062] The fracture current anomaly area determination device proposed in this embodiment can generate conductivity imaging data and trajectory curve expressions for multiple fractures based on the raw electrical imaging logging data of fractured oil and gas reservoirs. The conductivity imaging data and each trajectory curve expression are mapped to a predetermined three-dimensional coordinate system to generate a fracture trajectory curve corresponding to each trajectory curve expression and multiple discrete data points corresponding to the conductivity imaging data in the three-dimensional coordinate system. The three-dimensional coordinate system includes a depth axis, a wellbore expansion direction axis, and a conductivity axis. The background conductivity of the three-dimensional coordinate system is determined based on each discrete data point. Based on the background conductivity and each fracture trajectory curve, the current anomaly area of ​​each fracture is determined. This embodiment can effectively achieve accurate quantitative evaluation of the fracture current anomaly area, supporting the efficient exploration and development of fractured oil and gas reservoirs.

[0063] The crack current abnormal area determination device in this embodiment is presented in the form of a functional unit. Here, a unit refers to an ASIC (Application Specific Integrated Circuit) circuit, a processor and memory that execute one or more software or fixed programs, and / or other devices that can provide the above functions.

[0064] This invention also provides a computer device having the above-described features. Figure 8 The device shown is for determining the area of ​​crack current anomalies.

[0065] Please see Figure 9 The present invention provides a schematic diagram of the structure of a computer device according to an optional embodiment. The computer device includes one or more processors 10, a memory 20, and interfaces for connecting the various components, including high-speed interfaces and low-speed interfaces. The various components are interconnected via different buses and can be mounted on a common motherboard or otherwise installed as needed. The processors can process instructions executed within the computer device, including instructions stored in or on memory to display graphical information of a GUI on an external input / output device (such as a display device coupled to the interface). In some optional embodiments, multiple processors and / or multiple buses can be used with multiple memories, if desired. Similarly, multiple computer devices can be connected, each providing some of the necessary operations (e.g., as a server array, a group of blade servers, or a multiprocessor system). Figure 9 Take a processor 10 as an example.

[0066] Processor 10 may be a central processing unit, a network processor, or a combination thereof. Processor 10 may further include a hardware chip. The hardware chip may be an application-specific integrated circuit (ASIC), a programmable logic device (PLD), or a combination thereof. The programmable logic device may be a complex programmable logic device (CAMP), a field-programmable gate array (FPGA), a general-purpose array logic (GDA), or any combination thereof.

[0067] The memory 20 stores instructions executable by at least one processor 10 to cause at least one processor 10 to perform the method shown in the above embodiments.

[0068] The memory 20 may include a program storage area and a data storage area. The program storage area may store the operating system and applications required for at least one function. The data storage area may store data created based on the use of the computer device. Furthermore, the memory 20 may include high-speed random access memory and may also include non-transitory memory, such as at least one disk storage device, flash memory device, or other non-transitory solid-state storage device. In some alternative embodiments, the memory 20 may optionally include memory remotely located relative to the processor 10, which can be connected to the computer device via a network. Examples of such networks include, but are not limited to, the Internet, intranets, local area networks, mobile communication networks, and combinations thereof.

[0069] Memory 20 may include volatile memory, such as random access memory. Memory may also include non-volatile memory, such as flash memory, hard disk, or solid-state drive. Memory 20 may also include combinations of the above types of memory.

[0070] The computer device also includes a communication interface 30 for communicating with other devices or communication networks.

[0071] This invention also provides a computer-readable storage medium. The methods described above according to embodiments of the invention can be implemented in hardware or firmware, or implemented as computer code that can be recorded on a storage medium, or implemented as computer code downloaded via a network and originally stored on a remote storage medium or a non-transitory machine-readable storage medium and then stored on a local storage medium. Thus, the methods described herein can be processed by software stored on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. The storage medium can be a magnetic disk, optical disk, read-only memory, random access memory, flash memory, hard disk, or solid-state drive, etc.; further, the storage medium can also include combinations of the above types of memory. It is understood that computers, processors, microprocessor controllers, or programmable hardware include storage components capable of storing or receiving software or computer code, which, when accessed and executed by the computer, processor, or hardware, implements the methods shown in the above embodiments.

[0072] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, and not to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and these modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for determining the area of ​​crack current anomalies, characterized in that, include: Based on the raw electrical imaging logging data of fractured oil and gas reservoirs, electrical conductivity imaging data and trajectory curve expressions for multiple fractures are generated. The conductivity imaging data and each of the trajectory curve expressions are mapped to a set three-dimensional coordinate system to generate a fracture trajectory curve corresponding to each trajectory curve expression and multiple discrete data points corresponding to the conductivity imaging data in the three-dimensional coordinate system; wherein, the three-dimensional coordinate system includes a depth axis, a wellbore development direction axis and a conductivity axis; The background conductivity of the three-dimensional coordinate system is determined based on each of the discrete data points. Based on the background conductivity and the trajectory curve of each crack, the current anomaly area of ​​each crack is determined.

2. The method according to claim 1, characterized in that, The process of generating electrical conductivity imaging data and trajectory curve expressions for multiple fractures based on raw electrical imaging logging data of fractured oil and gas reservoirs includes: The raw data from electrical imaging logging of fractured oil and gas reservoirs are preprocessed to generate dynamic and static image data. The static image data is calibrated for conductivity to obtain the conductivity imaging data. An imaging operation is performed on the dynamic image data, the static image data, or the conductivity imaging data to obtain the corresponding target image; Crack identification and feature extraction are performed on the target image to obtain feature data for each crack, including crack depth, orientation, and dip angle. Based on the characteristic data of each crack, a trajectory curve expression for each crack is generated.

3. The method according to claim 1, characterized in that, Determining the background conductivity of the three-dimensional coordinate system based on each discrete data point includes: The coordinate interval on the wellbore unfolding direction axis is divided into multiple sub-intervals; For any of the sub-intervals, each target data point is found among the plurality of discrete data points, and the depth coordinate value of each target data point is located within the sub-interval. A dynamic threshold is generated based on the conductivity coordinate value of each target data point. Each background data point is determined among each target data point, and the conductivity coordinate value of each background data point is less than the dynamic threshold. The average conductivity coordinate value of each background data point is calculated to obtain the calculation result, which is used as the background conductivity corresponding to the sub-interval. The background conductivity of each sub-interval is taken as the background conductivity of the three-dimensional coordinate system.

4. The method according to claim 3, characterized in that, The step of generating a dynamic threshold based on the conductivity coordinate value of each target data point includes: The median of the conductivity coordinates for each target data point is obtained by statistically analyzing the median coordinate values. Calculate the absolute deviation between the conductivity coordinate value and the median of the coordinate values ​​for each of the target data points; For each absolute deviation value, the median absolute deviation is obtained by performing median statistics. The median absolute deviation is multiplied by a set factor to obtain a multiplication result, and the multiplication result is added to the median coordinate value to obtain the dynamic threshold.

5. The method according to claim 3, characterized in that, The determination of the current anomaly area of ​​each crack based on the background conductivity and the trajectory curve of each crack includes: For any of the fracture trajectory curves, data points are sampled along the wellbore unfolding direction at a set sampling interval to determine multiple sampling points; For any sampling point on the crack trajectory curve, a target data range is determined based on the depth coordinate value of the sampling point. The absolute value of the difference between the depth coordinate value of any point in the target data range and the depth coordinate value of the sampling point is not greater than a set threshold. Among each discrete data point included in the target data range, the discrete data point with the largest conductivity coordinate value is determined and used as a candidate data point. Each valid data point is determined among each candidate data point. The current anomaly area of ​​the crack is determined based on the conductivity coordinate value of each valid data point.

6. The method according to claim 5, characterized in that, The step of determining each valid data point from each of the candidate data points includes: For any candidate data point, with the candidate data point as the center, draw a circle with a set radius on the plane enclosed by the depth axis and the wellbore unfolding direction axis of the three-dimensional coordinate system, and search for whether there is a blank strip with a conductivity coordinate value of 0 in the drawn circle. If the blank stripe exists in the drawn circular surface, then the candidate data point is determined to be a non-valid data point; If no blank stripe exists in the drawn circular surface, then the candidate data point is determined to be a valid data point.

7. The method according to claim 6, characterized in that, Determining the current anomaly area of ​​the crack based on the conductivity coordinate value of each valid data point includes: For any valid data point on the crack trajectory curve, the circular surface drawn with the valid data point as the center is determined as the target integration region. The target variable is used as the integration object, which is the difference between the conductivity coordinate value of the coordinate point and the background conductivity corresponding to the coordinate point. The integration object in the target integration region is integrated to obtain the integration result. The integration result is multiplied by a set instrument coefficient to obtain the multiplication result, which is used as the current anomaly area corresponding to the valid data point. The current anomaly area corresponding to each valid data point on the crack trajectory curve is taken as the current anomaly area of ​​the crack corresponding to the crack trajectory curve.

8. A device for determining the area of ​​crack current anomalies, characterized in that, The apparatus used in the crack current anomaly area determination method according to any one of claims 1 to 7 comprises: The first generation unit is used to generate electrical conductivity imaging data and trajectory curve expressions of multiple fractures based on the raw electrical imaging logging data of fractured oil and gas reservoirs. The second generation unit is used to map the conductivity imaging data and each of the trajectory curve expressions to a set three-dimensional coordinate system, so as to generate a fracture trajectory curve corresponding to each trajectory curve expression and multiple discrete data points corresponding to the conductivity imaging data in the three-dimensional coordinate system; wherein, the three-dimensional coordinate system includes a depth axis, a wellbore development direction axis and a conductivity axis; The first determining unit is used to determine the background conductivity of the three-dimensional coordinate system based on each of the discrete data points; The second determining unit is used to determine the current anomaly area of ​​each crack based on the background conductivity and the trajectory curve of each crack.

9. A computer device, characterized in that, include: A memory and a processor are communicatively connected, the memory stores computer instructions, and the processor executes the computer instructions to perform the crack current anomaly area determination method according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores computer instructions for causing the computer to perform the crack current anomaly area determination method according to any one of claims 1 to 7.

Images