Electromagnetic method inversion result evaluation method, system, device and medium

Abstract

The invention discloses an electromagnetic method inversion result evaluation method, system and device and a medium. The method comprises the steps that logging resistivity data and electromagnetic inversion resistivity data of a calibration well in a research area and electromagnetic inversion resistivity profile layering data and modeling layering data of a marker bed in the research area are obtained, and the marker bed in the research area and upper and lower stratums of the marker bed are in the transverse direction; the average electromagnetic inversion resistivity change rate and the average modeling inter-formation resistivity change rate of each measuring point are obtained; calculating a longitudinal trend goodness of fit according to the logging resistivity data and the electromagnetic inversion resistivity data; calculating a transverse consistency parameter according to the electromagnetic inversion resistivity profile hierarchical data and the modeling hierarchical data; calculating a resistivity change matching degree according to the average electromagnetic inversion resistivity change rate and the average modeling inter-formation resistivity change rate; and calculating a comprehensive evaluation result based on the longitudinal trend goodness of fit, the transverse consistency parameter and the resistivity change matching degree. According to the invention, the evaluation precision of the accuracy of the inversion result can be improved.

Description

Technical Field

[0001] The present invention relates to the technical field related to oil and gas exploration and development, and in particular to an electromagnetic inversion results evaluation method, system, device and medium. Background Art

[0002] Electromagnetic inversion is a key technology for inferring electrical parameters of subsurface media, but its results are susceptible to factors such as data noise and model simplification, resulting in insufficient accuracy. Existing evaluation methods for electromagnetic inversion results are often limited to a single dimension and lack systematic and comprehensive verification, making it difficult to fully reflect the reliability of inversion results. Summary of the Invention

[0003] The present invention aims to solve at least one of the technical problems existing in the prior art. To this end, the present invention proposes an electromagnetic inversion results evaluation method that can improve the accuracy of the inversion results.

[0004] The present invention also provides an electromagnetic inversion results evaluation system, a control device for executing the above-mentioned electromagnetic inversion results evaluation method, and a computer-readable storage medium.

[0005] According to a first aspect of the present invention, the electromagnetic inversion results evaluation method includes: Obtaining logging resistivity data and electromagnetic inversion resistivity data of calibration wells in the study area, electromagnetic inversion resistivity profile layering data and modeling layering data of the marker layer in the study area, and the average electromagnetic inversion resistivity change rate and average modeling inter-stratum resistivity change rate of the marker layer and the strata above and below the marker layer in the horizontal direction at each measuring point in the study area; Calculating a vertical trend fit based on the well logging resistivity data and the electromagnetic inversion resistivity data, wherein the vertical trend fit is used to characterize the consistency between the vertical trend of the well logging resistivity and the trend of the electromagnetic inversion resistivity; Calculating a lateral consistency parameter based on the electromagnetic inversion resistivity profile layered data and the modeling layered data, wherein the lateral consistency parameter is used to characterize the consistency of the electromagnetic inversion resistivity profile top and bottom interface layered elevation values and the modeling top and bottom interface layered elevation values of each measuring point of each measuring line in the marker layer in the lateral direction; Calculate the resistivity change matching degree based on the average electromagnetic inversion resistivity change rate of each measuring point and the average modeled inter-stratum resistivity change rate of the marker layer and the strata above and below the marker layer in the study area in the horizontal direction. The resistivity change matching degree is used to characterize the consistency between the average electromagnetic inversion resistivity change rate and the average modeled inter-stratum resistivity change rate of the marker layer and the strata above and below the marker layer in the horizontal direction; A comprehensive evaluation result is calculated based on the longitudinal trend consistency, the transverse consistency parameter and the resistivity change matching degree.

[0006] The electromagnetic inversion results evaluation method according to the embodiment of the present invention has at least the following beneficial effects: The vertical trend consistency reflects the vertical consistency between the electromagnetic inversion results and the actual logging data, which plays a key role in judging the stratigraphic stratification and lithologic changes. The horizontal consistency parameter is used to evaluate the matching degree between the stratigraphic interface of the electromagnetic inversion resistivity profile and the modeling stratification data, which is of great significance for understanding the lateral continuity and structural characteristics of the stratigraphic formation. In geological structural research, horizontal stratigraphic information helps to determine the location and morphology of structures such as faults and folds. The resistivity change matching degree reflects the degree of conformity between the stratigraphic resistivity change in the inversion result and the actual geological conditions. The present invention introduces three types of data: logging, seismic, and electromagnetic results. The evaluation parameters are diverse. The evaluation features include vertical trend consistency, horizontal consistency parameters, and resistivity change matching degree. It is not limited to a single linear feature evaluation, but includes horizontal and vertical multi-dimensional evaluation. By verifying and evaluating the electromagnetic inversion results from multiple dimensions, the evaluation accuracy of the inversion results is improved, providing a reliable basis for geological exploration and related decision-making.

[0007] According to some embodiments of the present invention, the calculating the vertical trend consistency based on the well logging resistivity data and the electromagnetic inversion resistivity data includes: Normalizing the logging resistivity data according to a vertical resolution of the electromagnetic inversion resistivity data to obtain normalized resistivity data, wherein the electromagnetic inversion resistivity data includes a plurality of electromagnetic inversion resistivities at different altitudes, and the normalized resistivity data includes a plurality of normalized resistivities corresponding one-to-one to the plurality of electromagnetic inversion resistivities; calculating a difference between two electromagnetic inversion resistivities at adjacent altitudes to obtain a plurality of first differences; calculating a difference between two normalized resistivities at adjacent altitudes to obtain a plurality of second differences; The longitudinal trend consistency is calculated based on a plurality of the first differences and a plurality of the second differences.

[0008] According to some embodiments of the present invention, calculating the longitudinal trend consistency based on a plurality of first differences and a plurality of second differences includes: Recording the number of the first differences and the second differences as a total number of differences; Compare the first difference value and the corresponding second difference value of each group. If the first difference value and the corresponding second difference value are both greater than zero, equal to zero, or less than zero, then add one to the count; otherwise, do not count, and obtain the total count of all groups. The longitudinal trend consistency is calculated based on the total number of counts and the total number of differences.

[0009] According to some embodiments of the present invention, the electromagnetic inversion resistivity profile layering data includes the marker layer in the horizontal direction, the electromagnetic inversion resistivity profile top interface layering elevation value and the electromagnetic inversion resistivity profile bottom interface layering elevation value of each measuring point on each survey line, and the modeling layering data includes the marker layer in the horizontal direction, the modeling top interface layering elevation value and the modeling bottom interface layering elevation value of each measuring point on each survey line; The calculating of the lateral consistency parameter according to the electromagnetic inversion resistivity profile layered data and the modeling layered data comprises: Determine the thickness value corresponding to each measuring point according to the modeling top interface layer elevation value and the modeling bottom interface layer elevation value of each measuring point on each measuring line; Calculating the difference between the top interface layer elevation value of the electromagnetic inversion resistivity profile and the top interface layer elevation value of the modeling at each measuring point on each measuring line to obtain a plurality of third differences; Calculating the difference between the bottom interface layer elevation value of the electromagnetic inversion resistivity profile and the modeling bottom interface layer elevation value at each measuring point on each measuring line to obtain a plurality of fourth differences; Calculating an average error value of the top interface of each measuring line according to the multiple thickness values and the multiple third difference values of each measuring line; Calculating an average bottom interface error value of each measuring line according to the plurality of thickness values and the plurality of fourth difference values of each measuring line; The lateral consistency parameter is calculated according to the average error value of the top interface and the average error value of the bottom interface of each survey line.

[0010] According to some embodiments of the present invention, calculating the lateral consistency parameter according to the average error value of the top interface and the average error value of the bottom interface of each survey line includes: Calculate the average of the top interface average error value and the bottom interface average error value of each survey line to obtain a comprehensive error value; The lateral consistency parameter of each survey line is obtained by subtracting the comprehensive error value from 1.

[0011] According to some embodiments of the present invention, the average electromagnetic inversion resistivity change rate includes a first average electromagnetic inversion resistivity change rate of the marker layer and the marker layer, and a second average electromagnetic inversion resistivity change rate of the marker layer and the marker layer below the marker layer; the average modeled inter-stratum resistivity change rate includes a first average modeled inter-stratum resistivity change rate of the marker layer and the marker layer, and a second average modeled inter-stratum resistivity change rate of the marker layer and the marker layer below the marker layer; The calculating of the resistivity change matching degree according to the average electromagnetic inversion resistivity change rate of each measuring point and the average modeling inter-stratum resistivity change rate in the horizontal direction of the marker layer and the strata above and below the marker layer in the study area includes: The resistivity change matching degree of the upper part of the marker layer is calculated based on the first average electromagnetic inversion resistivity change rate of each measuring point and the first average modeling inter-stratum resistivity change rate; The resistivity change matching degree of the lower part of the marker layer is calculated based on the second average electromagnetic inversion resistivity change rate of each measuring point and the second average modeling inter-stratum resistivity change rate; The resistivity change matching degree is calculated based on the resistivity change matching degree of the upper portion of the marker layer and the resistivity change matching degree of the lower portion of the marker layer.

[0012] According to some embodiments of the present invention, the calculating of a comprehensive evaluation result based on the longitudinal trend consistency, the transverse consistency parameter, and the resistivity change matching degree includes: The comprehensive evaluation result is calculated based on the longitudinal trend matching degree and a preset first weight, the transverse consistency parameter and a preset second weight, and the resistivity change matching degree and a preset third weight.

[0013] According to the electromagnetic inversion results evaluation system of the second embodiment of the present invention, the system includes: a data acquisition unit, configured to acquire logging resistivity data and electromagnetic inversion resistivity data of calibration wells in a study area, electromagnetic inversion resistivity profile layering data and modeling layering data of a marker layer in the study area, and average electromagnetic inversion resistivity change rates of each measuring point and average modeling inter-stratum resistivity change rates of the marker layer and the strata above and below the marker layer in the horizontal direction; a vertical trend consistency calculation unit, configured to calculate a vertical trend consistency based on the well logging resistivity data and the electromagnetic inversion resistivity data, wherein the vertical trend consistency is used to characterize the consistency between the vertical trend of the well logging resistivity and the trend of the electromagnetic inversion resistivity; a lateral consistency parameter calculation unit, configured to calculate a lateral consistency parameter based on the electromagnetic inversion resistivity profile layered data and the modeling layered data, wherein the lateral consistency parameter is used to characterize the consistency between the top and bottom interface layered elevation values of the electromagnetic inversion resistivity profile and the modeling top and bottom interface layered elevation values of each measuring point of each measuring line in the horizontal direction of the marker layer; a resistivity change matching degree calculation unit, configured to calculate a resistivity change matching degree based on the average electromagnetic inversion resistivity change rate of each measuring point and the average modeled inter-stratum resistivity change rate of the marker layer and the strata above and below the marker layer in the study area in the horizontal direction, wherein the resistivity change matching degree is used to characterize the consistency between the average electromagnetic inversion resistivity change rate and the average modeled inter-stratum resistivity change rate of the marker layer and the strata above and below the marker layer in the horizontal direction; A comprehensive evaluation result calculation unit is used to calculate a comprehensive evaluation result based on the longitudinal trend matching degree, the transverse consistency parameter and the resistivity change matching degree.

[0014] The electromagnetic inversion results evaluation system according to the embodiment of the present invention has at least the following beneficial effects: The vertical trend consistency reflects the vertical consistency between the electromagnetic inversion results and the actual logging data, which plays a key role in judging the stratigraphic stratification and lithologic changes. The horizontal consistency parameter is used to evaluate the matching degree between the stratigraphic interface of the electromagnetic inversion resistivity profile and the modeling stratification data, which is of great significance for understanding the lateral continuity and structural characteristics of the stratigraphic formation. In geological structural research, horizontal stratigraphic information helps to determine the location and morphology of structures such as faults and folds. The resistivity change matching degree reflects the degree of conformity between the stratigraphic resistivity change in the inversion result and the actual geological conditions. The present invention introduces three types of data: logging, seismic, and electromagnetic results. The evaluation parameters are diverse. The evaluation features include vertical trend consistency, horizontal consistency parameters, and resistivity change matching degree. It is not limited to a single linear feature evaluation, but includes horizontal and vertical multi-dimensional evaluation. By verifying and evaluating the electromagnetic inversion results from multiple dimensions, the evaluation accuracy of the inversion results is improved, providing a reliable basis for geological exploration and related decision-making.

[0015] A control device according to a third embodiment of the present invention includes a memory, a processor, and a computer program stored in the memory and executable on the processor. When the processor executes the computer program, it implements the electromagnetic inversion results evaluation method described in the first embodiment. Because the control device utilizes all of the technical solutions of the electromagnetic inversion results evaluation method described in the above embodiment, it at least has all of the beneficial effects provided by the technical solutions of the above embodiment.

[0016] A computer-readable storage medium according to a fourth embodiment of the present invention stores computer-executable instructions for executing the electromagnetic inversion results evaluation method described in the first embodiment. Because the computer-readable storage medium incorporates all of the technical solutions of the electromagnetic inversion results evaluation method described in the aforementioned embodiment, it at least has all of the beneficial effects provided by the technical solutions of the aforementioned embodiment.

[0017] Other features and advantages of the present invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] The above and / or additional aspects and advantages of the present invention will become apparent and readily understood from the description of the embodiments with reference to the following drawings, in which: Figure 1 This is a flow chart of a method for evaluating electromagnetic inversion results according to one embodiment of the present invention; Figure 21 is a schematic diagram of normalization processing of well logging data according to an embodiment of the present invention; Figure 3 FIG. 4 is a schematic diagram of the vertical resolution of the unified electromagnetic inversion resistivity according to an embodiment of the present invention. DETAILED DESCRIPTION

[0019] The following describes embodiments of the present invention in detail. Examples of the embodiments are shown in the accompanying drawings, wherein the same or similar reference numerals throughout represent the same or similar elements or elements having the same or similar functions. The embodiments described below with reference to the accompanying drawings are exemplary and are intended only to explain the present invention and are not to be construed as limiting the present invention.

[0020] In the description of the present invention, if there is a description of first, second, etc., it is only for the purpose of distinguishing technical features, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features or implicitly indicating the order of the indicated technical features.

[0021] In the description of the present invention, it should be understood that descriptions involving orientation, such as the orientation or positional relationship indicated by up, down, etc., are based on the orientation or positional relationship shown in the accompanying drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operated in a specific orientation, and therefore cannot be understood as a limitation on the present invention.

[0022] In the description of the present invention, it should be noted that, unless otherwise clearly defined, terms such as setting, installing, and connecting should be understood in a broad sense, and technicians in the relevant technical field can reasonably determine the specific meanings of the above terms in the present invention based on the specific content of the technical solution.

[0023] Electromagnetic inversion is a key technology for inferring electrical parameters of subsurface media, but its results are susceptible to factors such as data noise and model simplification, resulting in insufficient accuracy. Existing evaluation methods for electromagnetic inversion results are often limited to a single dimension and lack systematic and comprehensive verification, making it difficult to fully reflect the reliability of inversion results.

[0024] For example, in the prior art, there is a method that uses the ratio of the logging Mahalanobis distance to the inversion Mahalanobis distance as the Mahalanobis distance ratio, calculates the difference between the ratio and the benchmark value, and determines whether the inversion result is reliable by comparing the difference with a preset value, but the evaluation parameters are single.

[0025] Other existing techniques use well logging data to extract reservoir sensitivity curve features, generating a sample set xi representing the reservoir at the wellpoint locations. Based on the waveform structure features of the seismic inversion data volume, a sample set yi representing the reservoir at the predicted points is obtained. Furthermore, the waveform structure features of the seismic data volume are extracted based on xi and yi, respectively, to generate waveform sample sets z1i and z2i representing the reservoir at the wellpoint locations and predicted points. Correlation coefficients are calculated between the sample sets xi and yi, and between the sample sets z1i and z2i, to quantitatively assess the confidence of the seismic inversion results based on the magnitude of the correlation coefficients. When the two correlation coefficients fall within the same correlation interval, the confidence is assessed by calculating their average; when they fall outside the same interval, the confidence is assessed based on the correlation coefficient between the sample sets z1i and z2i. This method calculates correlations based on measured data, comprehensively considering the seismic waveform characteristics and attribute curve characteristics of the reservoir, and achieves a quantitative assessment of the confidence of the inversion results. However, this quantitative assessment is a single, one-dimensional, linear assessment.

[0026] The above two methods are limited to a single dimension and lack systematic comprehensive verification, making it difficult to fully reflect the reliability of the inversion results.

[0027] In order to solve the problem that the electromagnetic inversion results evaluation methods in the prior art are mostly limited to a single dimension, lack systematic comprehensive verification, and are difficult to fully reflect the reliability of the inversion results, the embodiment of the present invention proposes an electromagnetic inversion results evaluation method, which can verify and evaluate the electromagnetic inversion results from multiple dimensions, thereby improving the assessment accuracy of the inversion results.

[0028] The following will be combined Figures 1 to 3 A clear and complete description is given of the electromagnetic inversion results evaluation method according to an embodiment of the present invention. Obviously, the embodiment described below is only a part of the embodiments of the present invention, not all of the embodiments.

[0029] refer to Figures 1 to 3 , Figure 1 This is a flow chart of a method for evaluating electromagnetic inversion results according to one embodiment of the present invention; Figure 2 1 is a schematic diagram of normalization processing of well logging data according to an embodiment of the present invention; Figure 3 FIG. 4 is a schematic diagram of the vertical resolution of the unified electromagnetic inversion resistivity according to an embodiment of the present invention.

[0030] According to the first embodiment of the present invention, the electromagnetic inversion results evaluation method includes: Obtain logging resistivity data and electromagnetic inversion resistivity data of calibration wells in the study area, electromagnetic inversion resistivity profile layer data and modeling layer data of marker layers in the study area, and the average electromagnetic inversion resistivity change rate and average modeling inter-stratum resistivity change rate of the marker layer and the strata above and below the marker layer in the horizontal direction at each measuring point in the study area; The vertical trend agreement is calculated based on the well logging resistivity data and the electromagnetic inversion resistivity data. The vertical trend agreement is used to characterize the consistency between the vertical trend of the well logging resistivity and the trend of the electromagnetic inversion resistivity. The lateral consistency parameter is calculated based on the electromagnetic inversion resistivity profile layer data and the modeling layer data. The lateral consistency parameter is used to characterize the consistency of the top and bottom interface layer elevation values of the electromagnetic inversion resistivity profile and the modeling top and bottom interface layer elevation values of the marker layer in the horizontal direction at each measuring point of each survey line; The resistivity change matching degree is calculated based on the average electromagnetic inversion resistivity change rate of each measuring point and the average resistivity change rate between the modeled strata in the horizontal direction of the marker layer and the strata above and below the marker layer in the study area. The resistivity change matching degree is used to characterize the consistency of the average electromagnetic inversion resistivity change rate and the average resistivity change rate between the modeled strata in the horizontal direction of the marker layer and the strata above and below the marker layer. The comprehensive evaluation results are calculated based on the longitudinal trend consistency, lateral consistency parameters and resistivity change matching.

[0031] It can be understood that a calibration well is a well within the study area whose wellhead coordinates coincide with or are close to the electromagnetic acquisition measurement points. Generally, full-well logging data from such vertical wells or data from above the buildup point for deviated wells are selected for electromagnetic inversion resistivity comparison at a unified location. A marker layer is a formation near the target layer in the study area that has relatively stable lithology, electrical properties, and thickness, and is easily identifiable.

[0032] First, the electromagnetic field information, seismic data, well logging data, geological data, geodetic coordinates and altitude values of each measuring point in the study area are obtained to determine the calibration wells and marker layers.

[0033] The position of the measuring point can be represented by the geodetic coordinates ( , , ) is expressed in the form of coordinate points, where , The geodetic coordinate system of the location 、 value, To correspond to the altitude depth. Electromagnetic recording method can also be used ( , )express, Indicates the line number during specific wide-area electromagnetic measurement. ( , ) indicates the specific location of the measurement, Indicates the distance between each measuring point on the plane and the starting point, as shown in Table 1 below, 0 is the starting point, Represents the distance from the starting point along the survey line Meter distance, and so on, Indicates the different vertical altitude depths of a certain measuring point. Indicated by 、 The electromagnetic inversion resistivity measured at the location determined by the value.

[0034] Table 1

[0035] Well locations and seismic records are usually expressed in geodetic coordinates. Seismic and well logging data often use geodetic coordinates to record the location of data points. The geodetic coordinates correspond to the location of each measuring point and are further converted into electromagnetic recording methods to uniformly correspond the measuring point coordinate information with the calibration well location coordinates and seismic coordinate information.

[0036] Based on electromagnetic field information, seismic data, and geological data, a 3D seismic geological framework is extracted. Well logging layers are used to calibrate and construct a geological model. The initial geoelectrical model is constructed using well logging resistivity as an electrical model. The initial geoelectrical model includes sub-layer thickness and initial resistivity. Sub-layer thickness is determined by well logging layers, while initial resistivity is determined by well logging resistivity.

[0037] Based on the initial geoelectric model, electromagnetic inversion of the known geological framework was carried out to obtain the inversion results data of each survey line in the study area.

[0038] The inversion results data for each survey line include ( , ) and the corresponding electromagnetic inversion resistivity . By the location of each measuring point : The distance between each measuring point on the measuring line and the starting point, and the different altitude depths of a certain measuring point in the vertical direction , and by 、 The electromagnetic inversion resistivity measured at the location determined by the value .

[0039] The above is a method of expressing the electromagnetic inversion resistivity at different altitudes of each measuring point using electromagnetic parameters or geodetic coordinates. Seismic and well logging data often use geodetic coordinates to record the location of data points. These two recording methods can be used to unify the data comparison dimensions and convert them into electromagnetic parameter recording methods or geodetic coordinate recording methods for comparison.

[0040] In some embodiments of the present invention, reference Figure 2 , calculate the vertical trend consistency based on the well logging resistivity data and the electromagnetic inversion resistivity data, including: Normalizing the logging resistivity data according to the vertical resolution of the electromagnetic inversion resistivity data to obtain normalized resistivity data, wherein the electromagnetic inversion resistivity data includes a plurality of electromagnetic inversion resistivities at different altitudes, and the normalized resistivity data includes a plurality of normalized resistivities corresponding one-to-one to the plurality of electromagnetic inversion resistivities; calculating a difference between two electromagnetic inversion resistivities at adjacent altitudes to obtain a plurality of first differences; calculating a difference between two normalized resistivities at adjacent altitudes to obtain a plurality of second differences; The longitudinal trend fit is calculated based on the plurality of first differences and the plurality of second differences.

[0041] According to the geodetic coordinate position of the calibration well, the measuring point at the same position or close to the position is confirmed, and the logging resistivity of the calibration well is compared with the electromagnetic inversion resistivity of the measuring point.

[0042] It is understandable that the vertical resolution of the electromagnetic inversion resistivity needs to be analyzed to obtain the top and bottom elevation data of each measurement value. 、 , and the corresponding altitude difference The vertical resolution of electromagnetic inversion resistivity differs from that of well logging and seismic data. When comparing multidimensional data, the vertical resolution must be normalized to unify it. Since the vertical resolution of electromagnetic inversion resistivity is the lowest among several types of data, they are all unified to the vertical resolution of electromagnetic inversion resistivity, and then the well logging resistivity data of the calibration well and the electromagnetic inversion resistivity data are compared.

[0043] The logging resistivity data of the calibration well contains the logging depth and the corresponding resistivity value, but the vertical resolution of the logging resistivity and the electromagnetic inversion resistivity are inconsistent, such as Figure 2 As shown, ~ The electromagnetic inversion resistivity of the altitude section is The logging resistivity corresponding to this altitude section contains multiple values. By taking weighted average of multiple groups of logging resistivity between altitudes, normalization is performed to unify the vertical resolution and normalized resistivity. The constraint formula is: ;Formula (1) in, represents the normalized resistivity, express ~ Layer thickness at altitude, express ~ The small layers in the altitude section are thick.

[0044] According to the vertical resolution of each set of electromagnetic inversion resistivity, all the logging resistivity are normalized accordingly. Each corresponds to one .

[0045] The difference between two electromagnetic inversion resistivities at all upper and lower adjacent altitudes in the vertical direction is calculated to obtain a plurality of first differences, and the difference between two normalized resistivities at all upper and lower adjacent altitudes in the vertical direction is calculated to obtain a plurality of second differences.

[0046] In some embodiments of the present invention, calculating the longitudinal trend consistency based on the plurality of first differences and the plurality of second differences includes: Record the number of first differences and second differences as the total number of differences; Compare the first difference value of each group with the corresponding second difference value. If the first difference value and the corresponding second difference value are both greater than zero, equal to zero, or less than zero at the same time, then add one to the count. Otherwise, do not count. Get the total count of all groups. The longitudinal trend was calculated based on the total number of counts and the total number of differences.

[0047] It can be understood that the sizes of the first difference and the second difference are divided into three cases: greater than zero, equal to zero, or less than zero. If the first difference and the corresponding second difference are both greater than zero, equal to zero, or less than zero at the same time, it means that the change trends are consistent, and the count is increased by one. Otherwise, the count is not made, and the total accumulated count is the trend consistency number.

[0048] Longitudinal trend consistency The constraint formula is: ;Formula (2) There may be multiple calibration wells in a study area. When calculating the vertical trend consistency, the following method is used: ;Formula (3) In some embodiments of the present invention, the electromagnetic inversion resistivity profile layering data includes the marker layer in the horizontal direction, the electromagnetic inversion resistivity profile top interface layering elevation value and the electromagnetic inversion resistivity profile bottom interface layering elevation value of each measuring point of each survey line, and the modeling layering data includes the marker layer in the horizontal direction, the modeling top interface layering elevation value and the modeling bottom interface layering elevation value of each measuring point of each survey line; Calculate lateral consistency parameters based on electromagnetic inversion resistivity profile layer data and modeling layer data, including: Determine the thickness value corresponding to each measuring point according to the modeling top interface layer elevation value and the modeling bottom interface layer elevation value of each measuring point on each measuring line; Calculate the difference between the top interface layer elevation value of the electromagnetic inversion resistivity profile and the top interface layer elevation value of the modeling at each measuring point on each measuring line to obtain multiple third difference values; Calculating the difference between the bottom interface layer elevation value of the electromagnetic inversion resistivity profile and the bottom interface layer elevation value of the modeling at each measuring point on each measuring line to obtain a plurality of fourth difference values; Calculate the average error value of the top interface of each measuring line according to multiple thickness values and multiple third difference values of each measuring line; Calculating an average bottom interface error value of each survey line according to multiple thickness values and multiple fourth difference values of each survey line; The lateral consistency parameter is calculated based on the average error value of the top interface and the average error value of the bottom interface of each survey line.

[0049] The electromagnetic inversion resistivity profile layer data of the marker layer is obtained through the following steps: According to the electromagnetic inversion resistivity profile and inversion results data, the first-order derivative of the inversion results data is performed to obtain the first-order difference curve. The first-order difference curve is used to characterize the difference in the average resistivity of adjacent depth formations. Then, according to the sawtooth points of the first-order difference curve and combined with geological information, the stratification position of the formation is determined to obtain the electromagnetic inversion resistivity profile stratification data.

[0050] It should be noted that the specific stratification principle of the above-mentioned electromagnetic inversion resistivity profile stratification data is an existing technology known to those skilled in the art, and the specific process and principle will not be described in detail here.

[0051] The top and bottom elevation values of the marker layer seismic stratification line at each measuring point are recorded as ( , ), and calculate the thickness of the marker layer corresponding to each measuring point The top and bottom elevation values of the marker layer at each measuring point in the electromagnetic inversion resistivity profile are recorded as ( , ).

[0052] Calculate the difference between the top interface layer elevation value of the electromagnetic inversion resistivity profile and the top interface layer elevation value of the modeling at each measuring point on each measuring line to obtain multiple third difference values. ( , ); Calculate the difference between the bottom interface layer elevation value of the electromagnetic inversion resistivity profile and the bottom interface layer elevation value of the modeling at each measuring point on each measuring line, and obtain multiple fourth difference values. ( , ).

[0053] The third difference ( , ) is: ;Formula (4) in, Indicates the measurement line number. Indicates the location of the measuring point. Indicates the layered elevation value of the top interface of the modeling, Indicates the elevation value of the top interface of the electromagnetic inversion resistivity profile. express Survey line No., distance from the starting point along the survey line The difference between the elevation value of the top interface layer of the modeling and the elevation value of the top interface layer of the electromagnetic inversion resistivity profile at the measuring point at 1000 meters.

[0054] Fourth difference ( , ) is: ;Formula (5) in, Indicates the measurement line number. Indicates the location of the measuring point. Indicates the layered elevation value of the modeling bottom interface, Indicates the bottom interface layer elevation value of the electromagnetic inversion resistivity profile, express Survey line No., distance from the starting point along the survey line At the measuring point at 100 meters, the difference between the modeled bottom interface layer elevation value and the electromagnetic inversion resistivity profile bottom interface layer elevation value.

[0055] Calculate the average error value of the top interface of each measuring line based on multiple thickness values and multiple third differences of each measuring line , calculate the average error value of the bottom interface of each measuring line based on multiple thickness values and multiple fourth differences of each measuring line , and then according to the average error value of the top interface and the average error value of the bottom interface Calculate the comprehensive error value of the marker layer .

[0056] Calculate the error value of the top interface of the layered marker layer at each measuring point of each measuring line marker layer and the bottom interface error value : = ( , ) *100%; formula (6) = ( , ) *100%; formula (7) Average error value of top interface The constraint formula is: ;Formula (8) in, Indicates measuring point The corresponding error value of the top interface of the marker layer, Indicates measuring point The thickness value of the corresponding marker layer.

[0057] Average error value of bottom interface The constraint formula is: ;Formula (9) in, Indicates measuring point The corresponding error value of the bottom interface of the marker layer, Indicates measuring point The thickness value of the corresponding marker layer.

[0058] In some embodiments of the present invention, calculating the lateral consistency parameter based on the average error value of the top interface and the average error value of the bottom interface of each survey line includes: Calculate the average of the top interface average error value and the bottom interface average error value of each survey line to obtain the comprehensive error value; Subtract the comprehensive error value from 1 to obtain the lateral consistency parameter of each survey line.

[0059] Comprehensive error value The constraint formula is: ;Formula (10) Horizontal consistency parameters The constraint formula is: ;Formula (11) In some embodiments of the present invention, reference Figure 3 The average electromagnetic inversion resistivity change rate includes a first average electromagnetic inversion resistivity change rate of the marker layer above and the marker layer, and a second average electromagnetic inversion resistivity change rate of the marker layer and the layer below the marker layer. The average modeled inter-stratum resistivity change rate includes a first average modeled inter-stratum resistivity change rate of the marker layer above and the marker layer, and a second average modeled inter-stratum resistivity change rate of the marker layer and the layer below the marker layer. The resistivity change matching degree is calculated based on the average electromagnetic inversion resistivity change rate of each measuring point and the average resistivity change rate between modeled strata in the horizontal direction of the marker layer and the strata above and below the marker layer in the study area, including: The resistivity change matching degree of the upper part of the marker layer is calculated based on the first average electromagnetic inversion resistivity change rate of each measuring point and the first average modeling inter-stratum resistivity change rate; The resistivity change matching degree of the lower part of the marker layer is calculated based on the second average electromagnetic inversion resistivity change rate of each measuring point and the second average modeling inter-stratum resistivity change rate; The resistivity change matching degree is calculated based on the resistivity change matching degree of the upper portion of the marker layer and the resistivity change matching degree of the lower portion of the marker layer.

[0060] Obtain the average resistivity value of each measuring point within the top and bottom range of the marker layer at each measuring point in the seismic inversion model, which is recorded as , and the average resistivity values of the upper and lower layers of the marker layer at each measuring point in the seismic inversion model 、 .

[0061] The first average modeled resistivity variation rate between layers The constraint formula is: ;Formula (12) Second average modeled resistivity change rate between formations The constraint formula is: ;Formula (13) Get the value of electromagnetic inversion resistivity contained in the top and bottom interfaces of the electromagnetic inversion resistivity marker layer , and the corresponding vertical resolution value , and a set of electromagnetic inversion resistivity values above and below the top and bottom interfaces of the marker layer 、 .

[0062] refer to Figure 3 , It is the weighted average value of multiple sets of electromagnetic inversion resistivity contained in the marker layer. The constraint formula is: ;Formula (14) in, Indicates the first The electromagnetic inversion resistivity corresponding to the small layer, Indicates the first The thickness of the small layer, Indicates the thickness of the marker layer.

[0063] The first average electromagnetic inversion resistivity change rate The constraint formula is: ;Formula (15) Second average electromagnetic inversion resistivity change rate The constraint formula is: ;Formula (16) Matching degree of resistivity change on the upper part of the marker layer The constraint formula is: ; Formula (17) in, Indicates measuring point The corresponding first average modeled inter-stratum resistivity change rate, Indicates measuring point The corresponding first average electromagnetic inversion resistivity change rate.

[0064] Matching degree of resistivity change in the lower part of the marker layer The constraint formula is: ; Formula (18) in, Indicates measuring point The corresponding second average modeled inter-stratum resistivity change rate, Indicates measuring point The corresponding second average electromagnetic inversion resistivity change rate.

[0065] Resistivity change matching The constraint formula is: ;Formula (19) In some embodiments of the present invention, a comprehensive evaluation result is calculated based on the longitudinal trend consistency, the lateral consistency parameter, and the resistivity change matching, including: A comprehensive evaluation result is calculated based on the longitudinal trend matching degree and a preset first weight, the transverse consistency parameter and a preset second weight, and the resistivity change matching degree and a preset third weight.

[0066] From a theoretical analysis, the vertical trend consistency reflects the vertical consistency between the electromagnetic inversion results and the actual logging data, which plays a key role in judging stratigraphic stratification and lithologic changes. According to relevant geological theories and exploration experience, the vertical trend consistency contributes significantly to the accuracy of the inversion results.

[0067] The lateral consistency parameter is used to assess the degree of match between the stratigraphic interfaces of the electromagnetic inversion resistivity profile and the modeling layer data. It is of great significance for understanding the lateral continuity and structural characteristics of the strata. In geological structural research, lateral stratigraphic information helps to determine the location and morphology of structures such as faults and folds. Although the lateral consistency parameter is also important, its direct impact on the accuracy of the inversion results is slightly weaker than that of the longitudinal trend fit.

[0068] The resistivity change matching degree of the marker layer can reflect the degree of conformity between the formation resistivity change in the inversion result and the actual geological conditions. In terms of sedimentary environment analysis, the resistivity change matching degree can provide important information. However, it is greatly affected by factors such as the selection of the marker layer and has a relatively high uncertainty.

[0069] The actual weight distribution should be determined based on the actual situation and experience of the study area, and depends on the completeness of the data. Generally speaking, the first, second, and third weights are generally arranged in the order: first weight > second weight > third weight. In some embodiments, the first weight can be set to 40%, the second weight to 35%, and the third weight to 25%.

[0070] The first weight, the second weight and the third weight are respectively denoted as 、 and , the longitudinal trend consistency, the lateral consistency parameter and the resistivity change matching degree are recorded as 、 、 Comprehensive evaluation results The constraint formula is: ;Formula (20) In some embodiments, when If it is less than 80%, it is unqualified and the model needs to be adjusted and re-modeled for electromagnetic inversion.

[0071] It should be noted that the specific value of the unqualified comprehensive evaluation result mentioned above is only an example to illustrate an embodiment and cannot be regarded as a limitation of the present invention.

[0072] According to the electromagnetic inversion results evaluation method of the embodiment of the present invention, the vertical trend fit reflects the vertical consistency of the electromagnetic inversion results and the actual logging data, which plays a key role in determining the stratigraphic stratification and lithologic changes. The lateral consistency parameter is used to evaluate the degree of match between the stratigraphic interface of the electromagnetic inversion resistivity profile and the modeled stratification data, which is of great significance for understanding the lateral continuity and structural characteristics of the stratigraphic formation. In geological structural research, lateral stratigraphic information helps to determine the location and morphology of structures such as faults and folds. The resistivity change match reflects the degree of conformity between the stratigraphic resistivity change in the inversion result and the actual geological conditions. The present invention introduces three types of data: well logging, seismic, and electromagnetic results, and has a variety of evaluation parameters. The evaluation features include vertical trend fit, lateral consistency parameter, and resistivity change match. It is not limited to single linear feature evaluation and includes horizontal and vertical multi-dimensional evaluation. By verifying and evaluating the electromagnetic inversion results from multiple dimensions, the accuracy of the inversion results is improved, providing a reliable basis for geological exploration and related decision-making.

[0073] According to the electromagnetic inversion results evaluation system of the second embodiment of the present invention, the system includes a data acquisition unit, a longitudinal trend matching calculation unit, a lateral consistency parameter calculation unit, a resistivity change matching calculation unit and a comprehensive evaluation result calculation unit.

[0074] A data acquisition unit is used to obtain logging resistivity data and electromagnetic inversion resistivity data of calibration wells in the study area, electromagnetic inversion resistivity profile layer data and modeling layer data of the marker layer in the study area, and the average electromagnetic inversion resistivity change rate and average modeling inter-stratum resistivity change rate of the marker layer and the strata above and below the marker layer in the horizontal direction at each measuring point in the study area; A vertical trend consistency calculation unit is used to calculate the vertical trend consistency based on the well logging resistivity data and the electromagnetic inversion resistivity data. The vertical trend consistency is used to characterize the consistency of the vertical trend of the well logging resistivity and the trend of the electromagnetic inversion resistivity. A lateral consistency parameter calculation unit is used to calculate the lateral consistency parameter based on the electromagnetic inversion resistivity profile layer data and the modeling layer data. The lateral consistency parameter is used to characterize the consistency of the electromagnetic inversion resistivity profile top and bottom interface layer elevation values and the modeling top and bottom interface layer elevation values of the marker layer in the horizontal direction at each measuring point of each measuring line; The resistivity change matching calculation unit is used to calculate the resistivity change matching degree based on the average electromagnetic inversion resistivity change rate of each measuring point and the average resistivity change rate between the modeled strata in the horizontal direction of the marker layer and the strata above and below the marker layer in the study area. The resistivity change matching degree is used to characterize the consistency of the average electromagnetic inversion resistivity change rate and the average resistivity change rate between the modeled strata in the horizontal direction of the marker layer and the strata above and below the marker layer; The comprehensive evaluation result calculation unit is used to calculate the comprehensive evaluation result based on the longitudinal trend consistency, the lateral consistency parameter and the resistivity change matching degree.

[0075] It can be understood that a calibration well is a well within the study area whose wellhead coordinates coincide with or are close to the electromagnetic acquisition measurement points. Generally, full-well logging data from such vertical wells or data from above the buildup point for deviated wells are selected for electromagnetic inversion resistivity comparison at a unified location. A marker layer is a formation near the target layer in the study area that has relatively stable lithology, electrical properties, and thickness, and is easily identifiable.

[0076] First, the electromagnetic field information, seismic data, well logging data, geological data, geodetic coordinates and altitude values of each measuring point in the study area are obtained to determine the calibration wells and marker layers.

[0077] The position of the measuring point can be represented by the geodetic coordinates ( , , ) is expressed in the form of coordinate points, where , The geodetic coordinate system of the location 、 value, To correspond to the altitude depth. Electromagnetic recording method can also be used ( , )express, Indicates the line number during specific wide-area electromagnetic measurement. ( , ) indicates the specific location of the measurement, Indicates the distance between each measuring point on the plane and the starting point, as shown in Table 1, 0 is the starting point, Represents the distance from the starting point along the survey line Meter distance, and so on, Indicates the different vertical altitude depths of a certain measuring point. Indicated by 、 The electromagnetic inversion resistivity measured at the location determined by the value.

[0078] Well locations and seismic records are usually expressed in geodetic coordinates. Seismic and well logging data often use geodetic coordinates to record the location of data points. The geodetic coordinates correspond to the location of each measuring point and are further converted into electromagnetic recording methods to uniformly correspond the measuring point coordinate information with the calibration well location coordinates and seismic coordinate information.

[0079] Based on electromagnetic field information, seismic data, and geological data, a 3D seismic geological framework is extracted. Well logging layers are used to calibrate and construct a geological model. The initial geoelectrical model is constructed using well logging resistivity as an electrical model. The initial geoelectrical model includes sub-layer thickness and initial resistivity. Sub-layer thickness is determined by well logging layers, while initial resistivity is determined by well logging resistivity.

[0080] Based on the initial geoelectric model, electromagnetic inversion of the known geological framework was carried out to obtain the inversion results data of each survey line in the study area.

[0081] The inversion results data for each survey line include ( , ) and the corresponding electromagnetic inversion resistivity . By the location of each measuring point : The distance between each measuring point on the measuring line and the starting point, and the different altitude depths of a certain measuring point in the vertical direction , and by 、 The electromagnetic inversion resistivity measured at the location determined by the value .

[0082] The above is a method of expressing the electromagnetic inversion resistivity at different altitudes of each measuring point using electromagnetic parameters or geodetic coordinates. Seismic and well logging data often use geodetic coordinates to record the location of data points. These two recording methods can be used to unify the data comparison dimensions and convert them into electromagnetic parameter recording methods or geodetic coordinate recording methods for comparison.

[0083] In some embodiments of the present invention, calculating the vertical trend consistency based on the well logging resistivity data and the electromagnetic inversion resistivity data includes: Normalizing the logging resistivity data according to the vertical resolution of the electromagnetic inversion resistivity data to obtain normalized resistivity data, wherein the electromagnetic inversion resistivity data includes a plurality of electromagnetic inversion resistivities at different altitudes, and the normalized resistivity data includes a plurality of normalized resistivities corresponding one-to-one to the plurality of electromagnetic inversion resistivities; calculating a difference between two electromagnetic inversion resistivities at adjacent altitudes to obtain a plurality of first differences; calculating a difference between two normalized resistivities at adjacent altitudes to obtain a plurality of second differences; Record the number of first differences and second differences as the total number of differences; Compare the first difference value of each group with the corresponding second difference value. If the first difference value and the corresponding second difference value are both greater than zero, equal to zero, or less than zero at the same time, then add one to the count. Otherwise, do not count. Get the total count of all groups. The longitudinal trend was calculated based on the total number of counts and the total number of differences.

[0084] According to the geodetic coordinate position of the calibration well, the measuring point at the same position or close to the position is confirmed, and the logging resistivity of the calibration well is compared with the electromagnetic inversion resistivity of the measuring point.

[0085] In some embodiments of the present invention, the electromagnetic inversion resistivity profile layering data includes the marker layer in the horizontal direction, the electromagnetic inversion resistivity profile top interface layering elevation value and the electromagnetic inversion resistivity profile bottom interface layering elevation value of each measuring point of each survey line, and the modeling layering data includes the marker layer in the horizontal direction, the modeling top interface layering elevation value and the modeling bottom interface layering elevation value of each measuring point of each survey line; Calculate lateral consistency parameters based on electromagnetic inversion resistivity profile layer data and modeling layer data, including: Determine the thickness value corresponding to each measuring point according to the modeling top interface layer elevation value and the modeling bottom interface layer elevation value of each measuring point on each measuring line; Calculate the difference between the top interface layer elevation value of the electromagnetic inversion resistivity profile and the top interface layer elevation value of the modeling at each measuring point on each measuring line to obtain multiple third difference values; Calculating the difference between the bottom interface layer elevation value of the electromagnetic inversion resistivity profile and the bottom interface layer elevation value of the modeling at each measuring point on each measuring line to obtain a plurality of fourth difference values; Calculate the average error value of the top interface of each measuring line according to multiple thickness values and multiple third difference values of each measuring line; Calculating an average bottom interface error value of each survey line according to multiple thickness values and multiple fourth difference values of each survey line; Calculate the average of the top interface average error value and the bottom interface average error value of each survey line to obtain the comprehensive error value; Subtract the comprehensive error value from 1 to obtain the lateral consistency parameter of each survey line.

[0086] In some embodiments of the present invention, the average electromagnetic inversion resistivity change rate includes a first average electromagnetic inversion resistivity change rate of the marker layer and the marker layer, and a second average electromagnetic inversion resistivity change rate of the marker layer and the layer below the marker layer. The average modeled inter-stratum resistivity change rate includes a first average modeled inter-stratum resistivity change rate of the marker layer and the marker layer, and a second average modeled inter-stratum resistivity change rate of the marker layer and the layer below the marker layer. The resistivity change matching degree is calculated based on the average electromagnetic inversion resistivity change rate of each measuring point and the average resistivity change rate between modeled strata in the horizontal direction of the marker layer and the strata above and below the marker layer in the study area, including: The resistivity change matching degree of the upper part of the marker layer is calculated based on the first average electromagnetic inversion resistivity change rate of each measuring point and the first average modeling inter-stratum resistivity change rate; The resistivity change matching degree of the lower part of the marker layer is calculated based on the second average electromagnetic inversion resistivity change rate of each measuring point and the second average modeling inter-stratum resistivity change rate; The resistivity change matching degree is calculated based on the resistivity change matching degree of the upper portion of the marker layer and the resistivity change matching degree of the lower portion of the marker layer.

[0087] In some embodiments of the present invention, a comprehensive evaluation result is calculated based on the longitudinal trend consistency, the lateral consistency parameter, and the resistivity change matching, including: A comprehensive evaluation result is calculated based on the longitudinal trend matching degree and a preset first weight, the transverse consistency parameter and a preset second weight, and the resistivity change matching degree and a preset third weight.

[0088] From a theoretical analysis, the vertical trend consistency reflects the vertical consistency between the electromagnetic inversion results and the actual logging data, which plays a key role in judging stratigraphic stratification and lithologic changes. According to relevant geological theories and exploration experience, the vertical trend consistency contributes significantly to the accuracy of the inversion results.

[0089] The lateral consistency parameter is used to assess the degree of match between the stratigraphic interfaces of the electromagnetic inversion resistivity profile and the modeling layer data. It is of great significance for understanding the lateral continuity and structural characteristics of the strata. In geological structural research, lateral stratigraphic information helps to determine the location and morphology of structures such as faults and folds. Although the lateral consistency parameter is also important, its direct impact on the accuracy of the inversion results is slightly weaker than that of the longitudinal trend fit.

[0090] The resistivity change matching degree of the marker layer can reflect the degree of conformity between the formation resistivity change in the inversion result and the actual geological conditions. In terms of sedimentary environment analysis, the resistivity change matching degree can provide important information. However, it is greatly affected by factors such as the selection of the marker layer and has a relatively high uncertainty.

[0091] The actual weight distribution needs to be determined based on the actual situation and experience of the study area, and depends on the completeness of the data. Overall, the size of the first weight, second weight and third weight is generally: first weight > second weight > third weight.

[0092] When the comprehensive evaluation result is less than the preset value, it is unqualified and the model needs to be adjusted and re-modeled for electromagnetic inversion.

[0093] According to the electromagnetic inversion results evaluation system of the embodiment of the present invention, the vertical trend fit reflects the vertical consistency of the electromagnetic inversion results and the actual logging data, which plays a key role in determining the stratigraphic stratification and lithologic changes. The lateral consistency parameter is used to evaluate the degree of match between the stratigraphic interface of the electromagnetic inversion resistivity profile and the modeled stratification data, which is of great significance for understanding the lateral continuity and structural characteristics of the stratigraphic formation. In geological structural research, lateral stratigraphic information helps to determine the location and morphology of structures such as faults and folds. The resistivity change match reflects the degree of conformity between the stratigraphic resistivity change in the inversion result and the actual geological conditions. The present invention introduces three types of data: well logging, seismic, and electromagnetic results, and has a variety of evaluation parameters. The evaluation features include vertical trend fit, lateral consistency parameter, and resistivity change match. It is not limited to single linear feature evaluation and includes horizontal and vertical multi-dimensional evaluation. By verifying and evaluating the electromagnetic inversion results from multiple dimensions, the accuracy of the inversion results is improved, providing a reliable basis for geological exploration and related decision-making.

[0094] Since the electromagnetic inversion result evaluation system adopts all the technical solutions of the electromagnetic inversion result evaluation method of the above embodiment, it has at least all the beneficial effects brought by the technical solutions of the above embodiment, and the specific process will not be repeated here.

[0095] In addition, an embodiment of the present invention further provides a control device, which includes: a memory, a processor, and a computer program stored in the memory and executable on the processor. The processor and the memory may be connected via a bus or other means.

[0096] The memory, as a non-transient computer-readable storage medium, can be used to store non-transient software programs and non-transient computer executable programs. In addition, the memory may include a high-speed random access memory and may also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage device. In some embodiments, the memory may optionally include a memory remotely arranged relative to the processor, and these remote memories may be connected to the processor via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0097] The non-transient software program and instructions required to implement the electromagnetic inversion result evaluation method of the above embodiment are stored in the memory, and when executed by the processor, the electromagnetic inversion result evaluation method of the above embodiment is executed.

[0098] The device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separate, i.e., they may be located in one place or distributed across multiple network units. Some or all of the modules may be selected based on actual needs to achieve the objectives of this embodiment.

[0099] In addition, an embodiment of the present invention also provides a computer-readable storage medium, which stores computer-executable instructions. The computer-executable instructions are executed by a processor or controller, for example, by the processor of the above embodiment, so that the above processor can execute the electromagnetic inversion results evaluation method in the above embodiment.

[0100] Those skilled in the art will appreciate that all or some of the steps and systems disclosed above can be implemented as software, firmware, hardware, or any suitable combination thereof. Some or all of the physical components may be implemented as software executed by a processor, such as a central processing unit, digital signal processor, or microprocessor, or as hardware, or as an integrated circuit, such as an application-specific integrated circuit. Such software may be distributed on computer-readable media, which may include computer storage media (or non-transitory media) and communication media (or transient media). As is well known to those skilled in the art, the term computer storage media includes volatile and non-volatile, removable and non-removable media implemented in any method or technology for storing information, such as computer-readable instructions, data structures, program modules, or other data. Computer storage media includes, but is not limited to, RAM, ROM, EEPROM, flash memory or other memory technologies, CD-ROM, digital versatile disks (DVDs) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other medium that can be used to store the desired information and can be accessed by a computer. Furthermore, as is well known to those skilled in the art, communication media typically embodies computer-readable instructions, data structures, program modules, or other data in a modulated data signal such as a carrier wave or other transport mechanism, and may include any information delivery media.

[0101] The embodiments of the present invention are described in detail above with reference to the accompanying drawings, but the present invention is not limited to the above embodiments. Various changes can be made within the scope of knowledge possessed by ordinary technicians in the relevant technical field without departing from the scope of the present invention.

Claims

1. A method for evaluating electromagnetic inversion results, characterized in that: The method comprises: Obtaining logging resistivity data and electromagnetic inversion resistivity data of calibration wells in the study area, electromagnetic inversion resistivity profile layering data and modeling layering data of the marker layer in the study area, and the average electromagnetic inversion resistivity change rate and average modeling inter-stratum resistivity change rate of the marker layer and the strata above and below the marker layer in the horizontal direction at each measuring point in the study area; Calculating a vertical trend fit based on the well logging resistivity data and the electromagnetic inversion resistivity data, wherein the vertical trend fit is used to characterize the consistency between the vertical trend of the well logging resistivity and the trend of the electromagnetic inversion resistivity; Calculating a lateral consistency parameter based on the electromagnetic inversion resistivity profile layered data and the modeling layered data, wherein the lateral consistency parameter is used to characterize the consistency of the electromagnetic inversion resistivity profile top and bottom interface layered elevation values and the modeling top and bottom interface layered elevation values of each measuring point of each measuring line in the marker layer in the lateral direction; Calculate the resistivity change matching degree based on the average electromagnetic inversion resistivity change rate of each measuring point and the average modeled inter-stratum resistivity change rate of the marker layer and the strata above and below the marker layer in the study area in the horizontal direction. The resistivity change matching degree is used to characterize the consistency between the average electromagnetic inversion resistivity change rate and the average modeled inter-stratum resistivity change rate of the marker layer and the strata above and below the marker layer in the horizontal direction; A comprehensive evaluation result is calculated based on the longitudinal trend consistency, the transverse consistency parameter and the resistivity change matching degree.

2. The electromagnetic inversion results evaluation method according to claim 1, characterized in that: Calculating the vertical trend consistency based on the well logging resistivity data and the electromagnetic inversion resistivity data includes: Normalizing the logging resistivity data according to a vertical resolution of the electromagnetic inversion resistivity data to obtain normalized resistivity data, wherein the electromagnetic inversion resistivity data includes a plurality of electromagnetic inversion resistivities at different altitudes, and the normalized resistivity data includes a plurality of normalized resistivities corresponding one-to-one to the plurality of electromagnetic inversion resistivities; calculating a difference between two electromagnetic inversion resistivities at adjacent altitudes to obtain a plurality of first differences; calculating a difference between two normalized resistivities at adjacent altitudes to obtain a plurality of second differences; The longitudinal trend consistency is calculated based on a plurality of the first differences and a plurality of the second differences.

3. The electromagnetic inversion results evaluation method according to claim 2, characterized in that: The calculating the longitudinal trend consistency according to the plurality of first differences and the plurality of second differences comprises: Recording the number of the first differences and the second differences as a total number of differences; Compare the first difference value and the corresponding second difference value of each group. If the first difference value and the corresponding second difference value are both greater than zero, equal to zero, or less than zero, then add one to the count; otherwise, do not count, and obtain the total count of all groups. The longitudinal trend consistency is calculated based on the total number of counts and the total number of differences.

4. The electromagnetic inversion results evaluation method according to claim 1, characterized in that: The electromagnetic inversion resistivity profile layering data includes the marker layer in the horizontal direction, the electromagnetic inversion resistivity profile top interface layering altitude value and the electromagnetic inversion resistivity profile bottom interface layering altitude value of each measuring point of each measuring line; the modeling layering data includes the marker layer in the horizontal direction, the modeling top interface layering altitude value and the modeling bottom interface layering altitude value of each measuring point of each measuring line; The calculating of the lateral consistency parameter according to the electromagnetic inversion resistivity profile layered data and the modeling layered data comprises: Determine the thickness value corresponding to each measuring point according to the modeling top interface layer elevation value and the modeling bottom interface layer elevation value of each measuring point on each measuring line; Calculating the difference between the top interface layer elevation value of the electromagnetic inversion resistivity profile and the top interface layer elevation value of the modeling at each measuring point on each measuring line to obtain a plurality of third differences; Calculating the difference between the bottom interface layer elevation value of the electromagnetic inversion resistivity profile and the modeling bottom interface layer elevation value at each measuring point on each measuring line to obtain a plurality of fourth differences; Calculating an average error value of the top interface of each measuring line according to the multiple thickness values and the multiple third difference values of each measuring line; Calculating an average bottom interface error value of each measuring line according to the plurality of thickness values and the plurality of fourth difference values of each measuring line; The lateral consistency parameter is calculated according to the average error value of the top interface and the average error value of the bottom interface of each survey line.

5. The electromagnetic inversion results evaluation method according to claim 4, characterized in that: The calculating the lateral consistency parameter according to the average error value of the top interface and the average error value of the bottom interface of each survey line includes: Calculate the average of the top interface average error value and the bottom interface average error value of each survey line to obtain a comprehensive error value; The lateral consistency parameter of each survey line is obtained by subtracting the comprehensive error value from 1.

6. The electromagnetic inversion result evaluation method according to claim 1, characterized in that: The average electromagnetic inversion resistivity change rate includes a first average electromagnetic inversion resistivity change rate of the marker layer and the marker layer, and a second average electromagnetic inversion resistivity change rate of the marker layer and the marker layer below the marker layer. The average modeled inter-stratum resistivity change rate includes a first average modeled inter-stratum resistivity change rate of the marker layer and the marker layer, and a second average modeled inter-stratum resistivity change rate of the marker layer and the marker layer below the marker layer. The calculating of the resistivity change matching degree according to the average electromagnetic inversion resistivity change rate of each measuring point and the average modeling inter-stratum resistivity change rate in the horizontal direction of the marker layer and the strata above and below the marker layer in the study area includes: The resistivity change matching degree of the upper part of the marker layer is calculated based on the first average electromagnetic inversion resistivity change rate of each measuring point and the first average modeling inter-stratum resistivity change rate; The resistivity change matching degree of the lower part of the marker layer is calculated based on the second average electromagnetic inversion resistivity change rate of each measuring point and the second average modeling inter-stratum resistivity change rate; The resistivity change matching degree is calculated based on the resistivity change matching degree of the upper portion of the marker layer and the resistivity change matching degree of the lower portion of the marker layer.

7. The electromagnetic inversion results evaluation method according to claim 1, characterized in that: The calculating of a comprehensive evaluation result based on the longitudinal trend consistency, the transverse consistency parameter and the resistivity change matching degree includes: The comprehensive evaluation result is calculated based on the longitudinal trend matching degree and a preset first weight, the transverse consistency parameter and a preset second weight, and the resistivity change matching degree and a preset third weight.

8. An electromagnetic inversion results evaluation system, characterized in that: The system comprises: a data acquisition unit, configured to acquire logging resistivity data and electromagnetic inversion resistivity data of calibration wells in a study area, electromagnetic inversion resistivity profile layering data and modeling layering data of a marker layer in the study area, and average electromagnetic inversion resistivity change rates of each measuring point and average modeling inter-stratum resistivity change rates of the marker layer and the strata above and below the marker layer in the horizontal direction; a vertical trend consistency calculation unit, configured to calculate a vertical trend consistency based on the well logging resistivity data and the electromagnetic inversion resistivity data, wherein the vertical trend consistency is used to characterize the consistency between the vertical trend of the well logging resistivity and the trend of the electromagnetic inversion resistivity; a lateral consistency parameter calculation unit, configured to calculate a lateral consistency parameter based on the electromagnetic inversion resistivity profile layered data and the modeling layered data, wherein the lateral consistency parameter is used to characterize the consistency between the top and bottom interface layered elevation values of the electromagnetic inversion resistivity profile and the modeling top and bottom interface layered elevation values of each measuring point of each measuring line in the horizontal direction of the marker layer; a resistivity change matching degree calculation unit, configured to calculate a resistivity change matching degree based on the average electromagnetic inversion resistivity change rate of each measuring point and the average modeled inter-stratum resistivity change rate of the marker layer and the strata above and below the marker layer in the study area in the horizontal direction, wherein the resistivity change matching degree is used to characterize the consistency between the average electromagnetic inversion resistivity change rate and the average modeled inter-stratum resistivity change rate of the marker layer and the strata above and below the marker layer in the horizontal direction; A comprehensive evaluation result calculation unit is used to calculate a comprehensive evaluation result based on the longitudinal trend matching degree, the transverse consistency parameter and the resistivity change matching degree.

9. A control device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the electromagnetic inversion results evaluation method according to any one of claims 1 to 7 is implemented.

10. A computer-readable storage medium storing computer-executable instructions, characterized in that: The computer-executable instructions are used to execute the electromagnetic inversion results evaluation method as described in any one of claims 1 to 7.

Images