Coalfield exploration borehole stratum data processing method

Abstract

The present application relates to the technical field of data processing, and discloses a coalfield exploration borehole stratum data processing method, which comprises the following steps: spatial interpolation is performed on measured borehole data to generate initial virtual borehole stratum data; conflict detection is performed on multiple check criteria that do not meet the requirements, and the check criteria that have a conflict relationship are marked as conflict criterion pairs; the current highest priority criterion is taken as the target criterion to adjust the stratum interface elevation value; the adjusted stratum interface elevation value is re-substituted into the check, and if another criterion still does not meet the requirements, another criterion is rotated to continue the correction, and when only a single check criterion does not meet the requirements, the stratum interface elevation value is corrected; all stratum interface elevation values are summarized to generate corrected virtual borehole stratum data, which is combined with the initial borehole stratum data to generate an enhanced borehole stratum data set; and the present application can improve the efficiency of the cooperative correction of the virtual borehole stratum interface elevation under the condition of multiple geological constraint conflicts.

Description

Technical Field

[0001] This invention relates to the field of data processing technology, and in particular to a method for processing formation data from coalfield exploration boreholes. Background Technology

[0002] In coalfield exploration, measured boreholes are typically sparsely distributed, making it difficult to directly support a continuous characterization of stratigraphic interfaces across the entire area. Current processing procedures generally involve first spatially interpolating the stratigraphic interface elevations of the measured boreholes to generate virtual borehole stratigraphic data. Then, various geological constraints, such as thickness continuity, dip gradient, and fault-fold compatibility, are introduced to verify the stratigraphic data at the virtual borehole locations. Based on this, stratigraphic interface elevations that do not meet the constraints are adjusted, aiming to increase borehole control density and improve the geological rationality of the stratigraphic data.

[0003] In exploration areas with complex geological structures, the same virtual borehole location often fails to meet multiple verification criteria simultaneously, and the correction requirements for the same stratigraphic interface elevation are inconsistent among these criteria. Existing methods typically verify and correct each constraint individually, which can easily lead to situations where adjusting the elevation according to one criterion results in another criterion still not being met, or even previously met criteria being violated again. This makes it difficult for elevation correction under multiple constraints to converge collaboratively, thus reducing the reliability of virtual borehole stratigraphic data. Therefore, improving the efficiency of collaborative elevation correction of virtual borehole stratigraphic interfaces under multiple geological constraint conflicts has become an urgent problem to be solved. Summary of the Invention

[0004] This invention provides a method for processing formation data from coalfield exploration boreholes to solve the problems mentioned in the background art.

[0005] To achieve the above objectives, the present invention provides a method for processing formation data from coalfield exploration boreholes, comprising: Z1. Spatial interpolation is performed on the measured borehole data in the initial borehole formation dataset to generate initial virtual borehole formation data, and a multi-geological constraint verification system is established based on the initial borehole formation data. Z2. Based on the multi-geological constraint verification system, the formation data at each virtual borehole location is verified, the verification criteria that are not met are determined, and conflict detection is performed on multiple verification criteria that are not met. Verification criteria that have conflict relationships are marked as conflict criterion pairs. Z3. Based on the fault density data and fold amplitude data in the initial borehole formation data, determine the priority order of each verification criterion, take the current highest priority criterion as the target criterion, and adjust the formation interface elevation value until the target criterion is met. Z4. Substitute the adjusted formation interface elevation value back into the verification. If the other criterion in the conflict criterion pair is still not satisfied, switch to another criterion to continue the correction until all verification criteria are satisfied or the preset maximum number of iterations is reached. When only a single verification criterion is not satisfied, refer to the measured data corresponding to the single verification criterion to correct the formation interface elevation value until the single verification criterion is satisfied. Z5. Summarize the corrected formation interface elevation values ​​of all virtual borehole points to generate corrected virtual borehole formation data, and merge it with the initial borehole formation dataset to generate an enhanced borehole formation dataset.

[0006] Compared with the prior art, the present invention has the following beneficial effects: 1. This invention addresses the problem of conflicting correction requirements for virtual borehole formation data under the simultaneous action of multiple geological constraints. Based on the identification of conflict criteria pairs, it can determine the target criteria according to priority and complete the directional adjustment of the formation interface elevation. When necessary, it can rotate and correct among conflict criteria, making it easier to coordinate and satisfy multiple constraints at the same virtual borehole location, thereby stably achieving elevation correction under conflict conditions.

[0007] 2. This invention enables the merging of corrected virtual borehole formation data with initial measured borehole data to form an enhanced borehole formation dataset with higher reliability and usability, thereby improving the efficiency of collaborative correction of virtual borehole formation interface elevations under multiple geological constraint conflicts. Attached Figure Description

[0008] Figure 1 This is a flowchart illustrating a method for processing formation data from coalfield exploration boreholes, as provided in an embodiment of the present invention.

[0009] The realization of the objective, functional features and advantages of the present invention will be further explained in conjunction with the embodiments and with reference to the accompanying drawings. Detailed Implementation

[0010] It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.

[0011] This application provides a method for processing formation data from coalfield exploration boreholes. The executing entity of this method includes, but is not limited to, at least one of the following electronic devices that can be configured to execute the method provided in this application: a server, a terminal, etc. In other words, the method for processing formation data from coalfield exploration boreholes can be executed by software or hardware installed on a terminal device or a server device. The server includes, but is not limited to, a single server, a server cluster, a cloud server, or a cloud server cluster. The server can be an independent server or a cloud server that provides basic cloud computing services such as cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communication, middleware services, domain name services, security services, content delivery networks (CDNs), and big data and artificial intelligence platforms.

[0012] Reference Figure 1 The diagram shown is a flowchart illustrating a method for processing formation data from coalfield exploration boreholes according to an embodiment of the present invention. In this embodiment, the method for processing formation data from coalfield exploration boreholes includes: Z1. Spatial interpolation is performed on the measured borehole data in the initial borehole formation dataset to generate initial virtual borehole formation data, and a multi-geological constraint verification system is established based on the initial borehole formation data. In practice, the processor reads the initial borehole formation dataset from the memory, extracts the planar coordinate data field of each measured borehole and the elevation value field of each measured borehole at each target formation interface, and forms a measured sample table.

[0013] The initial borehole formation dataset is a structured dataset, which includes, for example, borehole identifiers, plane coordinates, target formation interface identifiers, measured elevations, fault attitude and displacement related fields, fold axis traces and limb attitude related fields, etc.

[0014] The target stratigraphic interface identifier can be specified by configuration and is used to filter interface records participating in this batch of processing.

[0015] Based on the distribution density of the measured borehole plane coordinates within the target interpolation region, the processor generates a virtual borehole location table within the target interpolation region and writes the plane coordinates and location identifiers of each virtual borehole location.

[0016] In some embodiments, this application preferably uses a regular grid to generate virtual points: first, the measured average point distance is calculated, and then the grid spacing is determined according to one-half to one-third of the average point distance.

[0017] For example, when the average point spacing is about 500 meters, the grid spacing can be between 200 and 250 meters to achieve a balance between increasing data density and controlling interpolation uncertainty.

[0018] For each virtual borehole location, the processor uses the plane coordinates of that location as the center and selects the actual boreholes that are within the search radius from the actual boreholes to participate in the interpolation, according to the preset search radius.

[0019] The preset search radius is a planar distance threshold used for spatial interpolation neighborhood retrieval. It is provided by the configuration parameter table before entering the interpolation module. The preferred data type is floating point, and the unit is consistent with the planar coordinates.

[0020] The data processing meaning is as follows: taking the current virtual borehole point plane coordinates as the center, the measured borehole record is written into the interpolation borehole set only when the measured borehole plane coordinates fall into a circular neighborhood with a radius not greater than the threshold.

[0021] Therefore, this parameter directly determines the number of input records involved in the weighted overlay, which in turn affects the calculation result of the initial elevation field and the computational complexity of operations such as distance calculation and sorting.

[0022] In some embodiments of this application, the processor first scans all measured borehole plane coordinates within the target interpolation area, calculates the average nearest neighbor distance, and then takes the preset search radius as the product of the average nearest neighbor distance and a multiplier coefficient, wherein the multiplier coefficient is preferably between 1.5 and 2.5.

[0023] For example, when the average nearest neighbor distance is about 500 meters, the preset search radius can be between 800 meters and 1200 meters.

[0024] The rationale for selecting this interval is as follows: if the search radius is too small, a large number of virtual point neighborhood sets will be empty or there will be too few samples, forcing the interpolation module to frequently expand the circle or output null values, reducing the integrity of batch processing; if the search radius is too large, remote measured records will be included in the weighted calculation, reducing the spatial resolution of local data and significantly increasing the overhead of distance calculation and sorting.

[0025] Furthermore, the minimum number of measured boreholes and the maximum search radius can be configured simultaneously within the neighborhood.

[0026] When the number of records retrieved with the current preset search radius is less than the lower limit of the number of records, the processor expands the temporary search radius by a fixed increment and searches again until the number of records reaches the lower limit of the number of records or reaches the maximum search radius.

[0027] It should be noted that the temporary radius used for expanding the search area is a temporary variable during the operation; the preset search radius refers to the baseline search radius that is written into the configuration as the default value for the first retrieval and used as the main threshold for neighborhood filtering.

[0028] Based on the spatial distance between the measured borehole and the virtual borehole points involved in the interpolation, the influence weight of each measured elevation on the virtual point elevation is calculated.

[0029] In some embodiments, this application preferably employs inverse distance weighting: the weight is inversely proportional to the power of the distance, the power exponent is preferably two, and the weights of the same virtual point are normalized so that the sum of all weights is one.

[0030] Each measured elevation is multiplied by its corresponding weight to obtain a weighted elevation component, which is then superimposed to obtain the initial elevation value of the virtual point. The initial elevation values ​​of all virtual points are collected to generate an initial virtual borehole formation data table.

[0031] Through the interpolation process of controlling the neighborhood retrieval with a preset search radius, the processor expands the sparse measured records into a regularized virtual point elevation field, providing a unified data granularity for subsequent constraint verification.

[0032] Among them, a multi-geological constraint verification system is established based on the initial borehole formation data, which in terms of data processing is to generate a set of executable verification rule objects and their allowable range parameters.

[0033] In some embodiments, this application preferably includes: Based on the statistical results of the thickness and dip angle changes of adjacent measured boreholes, the allowable range of thickness change and the allowable interval of dip angle change rate are generated. Based on this, the judgment function of formation thickness continuity constraint and formation dip angle gradual change constraint is constructed. For example, the allowable zone is determined by the 10th to 90th percentile of the relative change rate or by the mean plus or minus 1.5 times the standard deviation, so that the verification function can output the deviation direction and deviation magnitude.

[0034] The fault spatial distribution range is generated from the fault-related fields and expanded into a fault influence zone polygon or grid mask. The fold influence zone mask is generated from the fold-related fields. The intersection of the two masks is obtained to obtain the fault-fold composite influence zone. Based on this, a decision function for the spatial compatibility constraint of fault-fold is constructed.

[0035] The external width can be used as an auxiliary configuration item, such as setting it as a multiple of the break distance or a fixed buffer distance.

[0036] The above constraints are assigned data type identifiers and integrated into a multi-geological constraint verification system, which is a set of rules with criterion identifiers as keys and decision functions and allowable parameters as values, for Z2 to call one by one.

[0037] Z2. Based on the multi-geological constraint verification system, the formation data at each virtual borehole location is verified, the verification criteria that are not met are determined, and conflict detection is performed on multiple verification criteria that are not met. Verification criteria that have conflict relationships are marked as conflict criterion pairs. In practice, the processor traverses the virtual borehole location table, substitutes the formation data field of each location into the judgment function of each criterion in the verification system, obtains the deviation direction data and deviation magnitude level, and writes them into the verification result table.

[0038] In some embodiments, this application preferably quantifies the deviation into five levels, from level zero to level four: level zero indicates that it falls within the allowable data range, and levels one to four indicate slight, moderate, significant, and severe deviation, respectively.

[0039] Mapping rules can be determined by the relative proportion of the excess beyond the allowable boundary. For example, if the excess is less than 10% of the allowable bandwidth, it is classified as Level 1; 10% to 30% is classified as Level 2; 30% to 60% is classified as Level 3; and greater than 60% is classified as Level 4.

[0040] Based on the deviation direction data and the deviation magnitude level, the verification criteria that exceed the preset deviation tolerance level are marked as unmet verification criteria, and the point identifier corresponding to the unmet verification criteria and the criterion identifier of the unmet verification criteria are recorded.

[0041] The preset deviation tolerance level is a level comparison threshold used in the constraint verification module to determine whether a certain criterion is "not satisfied". It is provided by the configuration parameter table before entering the verification judgment, and the data type is preferably an integer enumeration type.

[0042] The processor substitutes the formation data of each virtual borehole point into each verification criterion to obtain the deviation level; if and only if the deviation level is greater than the preset allowable deviation level, the criterion is marked as unsatisfied and written to the point identifier and the criterion identifier to form an unsatisfied record table.

[0043] From a data processing perspective, this parameter compresses continuous or semi-continuous deviations into binary decisions, thereby controlling the scale of records entering the subsequent conflict detection module.

[0044] In some embodiments of this application, preferably, under the above-mentioned zero-level to four-level quantization system, the preset deviation tolerance level is set to zero, that is, any level greater than zero is included in the non-satisfaction set, so as to achieve strict verification; when the input data noise is large and needs to be suppressed or corrected, the preset value can be relaxed to level one, so that only moderate and above deviations trigger the subsequent conflict detection and correction module.

[0045] Setting the value too low will expand the candidate correction set and increase the burden of subsequent iterations; setting the value too high will miss real conflicts and reduce the constraint consistency of the augmented dataset.

[0046] When two or more verification criteria are not met at the same virtual borehole location, the correction direction data of each unmet verification criterion are paired up, and the directional difference between the two correction direction data in each pair is determined. When the directional difference exceeds the preset directional difference tolerance threshold, the two verification criteria in the pair are marked as having a conflict relationship. All pairings with conflict relationships are extracted, and the two verification criteria in the same pair are recorded as a conflict criterion pair.

[0047] The corrected direction data includes a direction pointing component and a direction correction degree component.

[0048] The direction component can be represented by two discrete symbols, "rise" and "fall", or mapped to positive one and negative one; the direction correction degree component can be represented by the deviation amplitude level, or converted into the elevation adjustment requirement.

[0049] In some embodiments, the pairing and direction difference are preferably calculated as follows: the correction direction data of each non-satisfied verification criterion are paired to generate each pairing group; for each pairing group, the first direction pointing component and the second direction pointing component are compared; when the two are consistent, the first direction correction degree component is subtracted from the second direction correction degree component to generate the difference correction degree component, which is used as the direction difference of the pairing group; when the two are inconsistent, the first direction correction degree component and the second direction correction degree component are added to generate the sum correction degree component, which is used as the direction difference of the pairing group.

[0050] The preset directional difference tolerance threshold is a difference comparison threshold used in the conflict detection module to determine whether two non-compliant criteria constitute a conflict criterion pair. It is provided by the configuration parameter table before entering the conflict detection phase.

[0051] Its input is the directional difference of the above pairing groups, and the output is a logical flag indicating whether there is a conflict: when the directional difference is greater than the threshold, the pairing group is marked as having a conflict relationship and written into the conflict criterion table; if a single item that does not form a conflict does not meet the record, it enters the single item correction branch of Z4.

[0052] In some embodiments of this application, when the correction degree component is represented by a level, the preset allowable threshold for directional difference is two; when the correction degree component is represented by an elevation in meters, the threshold may be twice the preset step size described below.

[0053] The basis for the value is as follows: When the two criteria are in the same direction and the degree of correction is similar, the update requirements of the same elevation field are compatible and should not enter the conflict rotation process. Otherwise, the situation that can be met in one go will be misjudged as a conflict, wasting iteration resources. When the two criteria are in opposite directions or in the same direction but the degree of correction is too different, if the elevation field is still updated independently according to the criteria, it is very easy to have oscillations of "satisfying one while destroying the other". They must be marked as conflict criteria pairs and uniformly scheduled by priority and rotation mechanism.

[0054] The smaller the threshold, the more conflicting pairs there are, and the more frequently the subsequent steps three and four are triggered; the larger the threshold, the more pairs are considered to be cooperative, and the system is more inclined to single-path correction.

[0055] Those skilled in the art can adjust and solidify the oscillation rate and average number of iterations in the configuration during batch trial calculations.

[0056] By first filtering with a preset deviation tolerance level and then determining conflicts with a preset directional difference tolerance threshold, the original verification output can be compressed into a schedulable set of conflict criteria pairs and a set of individual items to be corrected, thus avoiding data oscillations caused by indiscriminate item-by-item correction.

[0057] Z3. Based on the fault density data and fold amplitude data in the initial borehole formation data, determine the priority order of each verification criterion, take the current highest priority criterion as the target criterion, and adjust the formation interface elevation value until the target criterion is met. In practice, the processor calculates the fault density value and fold amplitude value based on fields such as the number of faults aggregated in the evaluation block, the range of fault influence, the traces of fold axes, and the attitude of the two limbs.

[0058] Summing the two yields the priority score for fault-fold spatial coordination constraints. Then, based on the spatial matching relationship between fault density distribution, thickness distribution, and dip angle distribution, the thickness constraint correction coefficient and dip angle constraint correction coefficient are determined. The sum of the variance of fault density distribution and the variance of fold amplitude distribution are multiplied by the above coefficients to obtain the priority scores for thickness continuity constraints and dip angle gradual change constraints.

[0059] A priority sequence is generated by sorting the scores from highest to lowest, and the highest score is taken as the target criterion.

[0060] Based on the target criteria, the formation interface elevation values ​​at virtual borehole locations that have been marked with conflict criteria but do not meet the target criteria are adjusted until the target criteria are met.

[0061] Specifically, from the virtual borehole locations involved in the conflict criteria, the formation interface elevation values ​​at locations that do not meet the target criteria are selected as the elevation values ​​to be adjusted; the correction direction data corresponding to the target criteria are obtained, and a single incremental adjustment is performed on the elevation value to be adjusted according to the adjustment direction indicated by its direction component to generate a first intermediate elevation value; the first intermediate elevation value is back-substituted into the target criteria for verification, and if it passes, it is used as the adjusted formation interface elevation value; if it fails, the first intermediate elevation value is used as the new elevation value to be adjusted, and the single incremental adjustment and back-substitution verification are repeated until the verification passes.

[0062] In some embodiments of this application, the incremental magnitude used in the single incremental adjustment is the same as the preset step size in Z4, so that the directional adjustment and subsequent rotation correction are consistent in terms of update granularity; the data type, preferred value and the effects of the preset step size being too small or too large are described in detail in Z4 in conjunction with the rotation correction process.

[0063] The directional adjustment process is also constrained by the upper limit of the number of iterations. This upper limit is the same as the preset maximum number of iterations in Z4, and its meaning and preferred range are also explained in detail in Z4.

[0064] The data processing in this step is essentially as follows: under the constraint of conflict pairs, first perform a search with a step size and an upper limit on whether the judgment condition is met for the highest priority criterion, so as to provide a stable intermediate state for subsequent rotations.

[0065] Z4. Substitute the adjusted formation interface elevation value back into the verification. If the other criterion in the conflict criterion pair is still not satisfied, switch to another criterion to continue the correction until all verification criteria are satisfied or the preset maximum number of iterations is reached. When only a single verification criterion is not satisfied, refer to the measured data corresponding to the single verification criterion to correct the formation interface elevation value until the single verification criterion is satisfied. In practice, the adjusted stratigraphic interface elevation values ​​are substituted into the first and second criteria of the conflict criteria for verification. If both criteria are met, the correction is stopped.

[0066] When the first criterion meets the verification requirements and the second criterion does not meet the verification requirements, the second criterion is taken as the current correction object; according to the adjustment direction indicated by the correction direction data and the preset step size, the formation interface elevation value at the virtual borehole point that does not meet the current correction object is incrementally adjusted to generate a new adjusted formation interface elevation value, and then the first criterion and the second criterion are substituted back for re-verification.

[0067] The preset step size is the single adjustment amount of the elevation field incremental update module. It is provided by the configuration parameter table before entering the rotation correction or single approximation correction cycle. The data type is preferably floating point type, and the unit is consistent with the formation interface elevation field.

[0068] The processor updates the elevation value to be adjusted according to the direction component of the correction direction data as follows: it adds or subtracts a preset step size from the original elevation in the adjustment direction to obtain a new elevation value, which is then used as an intermediate elevation value for back-substitution and verification.

[0069] Therefore, this parameter controls the granularity of each write-back of the elevation field.

[0070] In some embodiments, the preset step size is preferably between 0.2 meters and 1 meter.

[0071] For data batches with high accuracy requirements or narrow allowable bandwidth, a value of 0.2 to 0.5 meters can be used; for data batches with large elevation variations, a value of 0.5 to 1 meter can be used.

[0072] Too small a step size will significantly increase the number of iterations required to reach the allowable boundary, thus increasing processor time; too large a step size will easily cause overshoot near the allowable boundary, resulting in the criterion being violated again in the next round, increasing the number of iterations.

[0073] Optionally, the configuration parameter table can support a variable step size strategy: initially use a larger step size for rapid approximation, and halve the current step size when the deviation level drops to level one; however, the preset step size refers to the baseline step size used as the default single incremental write configuration.

[0074] If, after back substitution, the first criterion is met but the second criterion is still not met, then the second criterion remains the current correction target, and incremental adjustment and back substitution verification are repeated. If, after back substitution, the first criterion is not met, then the current correction target is rotated to the first criterion, and incremental adjustment and back substitution verification are repeated. If both are met, then the correction is stopped.

[0075] If the cumulative number of rotations and incremental adjustments does not simultaneously meet the preset maximum number of iterations, then the correction will stop, and the formation interface elevation value after the last adjustment will be used as the final formation interface elevation value.

[0076] The preset maximum number of iterations is the forced termination threshold of the iterative correction module, which is provided by the configuration parameter table before entering the conflict rotation correction or single approximation correction loop. The data type is preferably a positive integer.

[0077] The processor maintains an iteration counter for each virtual borehole point to be corrected or each conflict criterion for the processing instance; the counter is incremented by one each time an incremental adjustment and back-substitution check are performed.

[0078] When both criteria in the conflict criterion pair pass the verification, or when only one criterion passes the verification, the process terminates normally and the current elevation is output. When the counter reaches the preset maximum number of iterations and the conditions are not met simultaneously, the correction is forcibly stopped, and the last adjusted formation interface elevation value is written back as the final result.

[0079] Its data processing function is to ensure that batch processing can be terminated, thus avoiding infinite loops that consume computing resources.

[0080] In some embodiments of this application, the maximum number of preset iterations is preferably between fifty and two hundred.

[0081] This range is matched with the preset step size and allowable bandwidth: when the step size is about 0.5 meters and the common deviation is on the order of several meters, dozens of iterations are usually sufficient to cover directional approximation and a limited number of rotations; for data partitions with more dense conflicts, a larger value can be taken.

[0082] If the value is too small, a large number of points will be truncated before convergence, leaving many quality marker points in the enhanced dataset; if the value is too large, computational power will be wasted on points with no objective intersection solution.

[0083] The approximation correction and conflict rotation correction of a single verification criterion can share the same preset maximum number of iterations, or two upper limits can be set separately in the configuration, but both belong to the termination threshold determined before the start of the loop.

[0084] When only a single verification criterion is not met, the formation interface elevation value is corrected by referring to the measured data corresponding to the single verification criterion until the single verification criterion is met.

[0085] Specifically, when the verification result indicates that only a single verification criterion is not met and no conflicting criterion pair is formed, the virtual borehole location that does not meet the single verification criterion is taken as the location to be corrected, and its formation interface elevation value is taken as the current elevation value to be corrected. Within the evaluation block to which the location to be corrected belongs, the nearest measured borehole is selected according to the planar distance, and the measured elevation value of the measured borehole at the same target formation interface is taken as the reference elevation value. The adjustment direction is determined according to the height relationship between the current elevation value to be corrected and the reference elevation value. Along the adjustment direction, a single incremental approximation adjustment is performed on the current elevation value to be corrected according to the preset step size to generate the formation interface elevation value after approximation adjustment and substitute it back into the single verification criterion. If it passes, it is taken as the corrected formation interface elevation value. If it fails, the value after approximation adjustment is taken as the new current elevation value to be corrected, and the single incremental approximation adjustment and substitution verification are repeated until the verification passes, or until the number of approximation adjustments reaches the preset maximum number of iterations.

[0086] Therefore, in Z4, the preset step size controls the granularity of writing back the elevation field each time, and the preset maximum number of iterations provides a deterministic loop exit condition, so that the correction process of the entire batch of virtual points ends within a finite number of calculation steps, which meets the requirements of computer implementation methods for termination and resource controllability.

[0087] Z5. Summarize the corrected formation interface elevation values ​​of all virtual borehole points, generate corrected virtual borehole formation data, and merge it with the initial borehole formation dataset to generate an enhanced borehole formation dataset. In practice, for the same virtual borehole location, the corrected elevations of each target interface are arranged in order of stratigraphic age from oldest to newest, generating a stratigraphic columnar sequence record; the corrected virtual borehole stratigraphic data table is obtained by merging according to the location identifier.

[0088] Optionally, perform layer sequence and positive thickness checks on adjacent interfaces within the same location, and write any abnormalities to the review flag field.

[0089] The columnar sequence of virtual points is added as a new entry to the logical dataset that is the same as the measured columnar sequence in the initial borehole formation dataset. Then, all entries are spatially sorted according to planar coordinates and the data source identifier is written to generate an enhanced borehole formation dataset. This dataset is then written back to memory or exported to a database for use in subsequent downstream data processing tasks such as mapping, profile interpretation, or 3D modeling.

[0090] In one feasible implementation, various preset parameters can be written by a graphical configuration interface or configuration file before the task starts; different evaluation blocks of the same exploration project can be equipped with different parameter groups, and the processor switches to read according to the block identifier to which the point belongs, thereby realizing a data processing strategy with differentiated partitioning without modifying the main process code.

[0091] In one feasible implementation, in the Z1 of the above embodiments, spatial interpolation can also use other numerical interpolation kernels such as ordinary kriging or radial basis functions; however, no matter how the kernel is replaced, the initial selection of neighborhood samples is still constrained by a preset search radius to keep the computational complexity controllable.

[0092] In one feasible implementation, in the above embodiment Z2, the deviation magnitude can also be kept as a continuous floating-point type deviation amount. At this time, the preset deviation tolerance level can be replaced with an equivalent preset deviation tolerance value, and the directional difference tolerance threshold in the conflict detection is changed to a floating-point type threshold accordingly. Its comparison logic and the level-type threshold in the claims are essentially the same as the "mark if the threshold is exceeded" determination in data processing.

[0093] In one feasible implementation, Z3 and Z4 in the above embodiments can adopt variable step size: the preset step size is used as the initial value, and it is halved when it approaches the allowable boundary; the iteration count is still constrained by the preset maximum number of iterations.

[0094] In one feasible implementation, when the preset maximum number of iterations is not simultaneously met in Z4 of the above embodiment, in addition to writing back the last elevation, a "not fully converged" quality flag and the final deviation fields of the two criteria can be added to the result table for use by the downstream quality screening module.

[0095] In one feasible implementation, Z5 in the above embodiment can generate a correction log table while merging the output, recording fields such as point identifier, initial elevation, final elevation, conflict criterion pair identifier, actual iteration count, and whether it is truncated due to reaching the preset maximum iteration count, so as to audit the data processing process and evaluate whether the parameter configuration is reasonable.

[0096] In some embodiments, in order to verify the effectiveness of the coalfield exploration borehole formation data processing method proposed in this application, a series of comparative tests will be conducted.

[0097] The test subjects are two types of computer-implemented processes based on the same initial borehole formation dataset: one is the traditional item-by-item verification-item-by-item correction process without conflict detection and priority rotation; the other is the conflict detection-priority orientation-rotation correction process described in this application, which is jointly controlled by a preset search radius, a preset deviation tolerance level, a preset directional difference tolerance threshold, a preset step size, and a preset maximum number of iterations.

[0098] The details are as follows: Experimental platform and testing conditions The experimental platform consists of the following components: a data processing host, configured with a multi-core central processing unit and no less than 16 gigabytes of memory, used to perform spatial interpolation, multi-geological constraint verification, conflict detection, priority sorting, iterative correction, and dataset merging; a relational or file-based database used to store the initial borehole formation dataset, configuration parameter table, intermediate result table, and enhanced borehole formation dataset; and an optional visualization terminal used to verify the distribution of virtual borehole points and correction logs.

[0099] The test data were selected from the existing borehole logging and structural interpretation digitization results of a coalfield exploration area.

[0100] The exploration area covers approximately 45 square kilometers, with 82 boreholes measured. The target strata interface was selected from the bottom interface of the main coal seam and the interface of the overlying marker layer.

[0101] In addition to the plane coordinates and interface elevation fields, the initial borehole formation dataset also includes structured attribute fields such as fault attitude, fault displacement, fold axis traces, and attitude of the limbs.

[0102] The experiment was conducted in a unified software environment. The two comparison processes read the same input dataset and the same set of basic tolerance interval statistical results. The only differences were whether conflict detection, priority sorting and rotation correction were enabled, and whether the above-mentioned preset parameters were used in accordance with the method of this application.

[0103] The preset parameters of the method in this application are configured as described in the previous steps: the preset search radius is twice the measured average point distance; the preset deviation tolerance level is zero; the preset directional difference tolerance threshold is set to two according to the level system; the preset step size is 0.5 meters; and the preset maximum number of iterations is one hundred.

[0104] Traditional comparison methods correct each non-compliant criterion in a fixed order, with each correction using a 0.5-meter increment. However, they do not rotate between conflicting criteria and markers, and the iteration limit is set to one hundred times to ensure a fair comparison of time consumption and truncation ratio.

[0105] Test Sample Design To comprehensively evaluate the performance of this data processing method under different data characteristics, three sets of differentiated test samples were designed.

[0106] Sample Group 1: The evaluation blocks are relatively simple to construct, with low fault density and small fold amplitude. At the virtual borehole locations, most of the single criteria are not met.

[0107] This sample set aims to verify the correction efficiency and elevation field stability of the single-item measured anchor approximation branch.

[0108] Sample Group 2: Evaluation block of fault-fold composite influence zone development, where the same virtual borehole location has a high probability of triggering two or three verification criteria simultaneously.

[0109] This sample set aims to verify the collaborative convergence capability of the conflict criterion for identification, priority-oriented adjustment, and rotation correction.

[0110] Sample group 3: Mixed samples of sparse and dense areas of measured boreholes, used to evaluate the integrity and usability of the enhanced borehole formation dataset after Z1 interpolation and Z5 merging under different control densities.

[0111] Experimental steps Step 1: Data preparation and system establishment.

[0112] Press Z1 to perform spatial interpolation based on the preset search radius, generate initial virtual borehole formation data, and establish a multi-geological constraint verification system; record the total number of virtual borehole points and the success rate of initial elevation field generation.

[0113] Step 2: Verification and Conflict Detection.

[0114] Substitute each virtual point into the verification system according to Z2, generate a non-compliance record table based on the preset deviation tolerance level, and generate a conflict criterion check table based on the preset directional difference tolerance threshold; count the percentage of non-compliance points, the percentage of conflict points, and the number of conflict criterion check tables.

[0115] Step 3: Priority adjustment and iterative correction.

[0116] For the method in this application, the target criterion is determined and adjusted according to Z3, and then the rotation correction or single-item approximation correction is performed according to Z4 with a preset step size, which is constrained by the preset maximum number of iterations; for the traditional comparison method, the fixed criteria are incrementally corrected one by one until the iteration limit is met or reached.

[0117] Record the average number of iterations, the proportion of points that simultaneously meet all relevant criteria, the proportion of points that are truncated due to reaching the preset maximum number of iterations, and the processing time per batch.

[0118] Step 4: Merge and compare the datasets.

[0119] Generate an enhanced borehole formation dataset using Z5; summarize the distribution of elevation differences before and after correction for both methods and the pass rate of sequence validity checks on the same set of points, and output correction logs for verification.

[0120] Experimental results Sample Group 1 In sample group one, both methods were able to restore the majority of points that were not initially satisfied to the desired condition.

[0121] The method in this application directly enters the measured anchor approximation branch, resulting in a lower average number of iterations than the traditional method that "tries according to multiple criteria in sequence"; the proportion of truncation due to reaching the preset maximum number of iterations is close to zero.

[0122] The results show that, under conditions of no or weak conflict data, the method of this application is at least no worse than the traditional step-by-step correction, and has an advantage in terms of the number of iterations.

[0123] Sample Group 2 In sample group 2, the traditional item-by-item correction method often exhibits the phenomenon of repeated updates where "after satisfying the thickness continuity constraint, the tilt angle gradual change constraint is broken, and then the tilt angle is corrected again, which in turn breaks the thickness or coordination constraint." The average number of iterations is too high, the simultaneous satisfaction rate is too low, and a certain proportion of points still do not satisfy the constraints simultaneously when the iteration limit is reached.

[0124] The method of this application first identifies conflicting criteria pairs by setting a preset deviation tolerance level and a preset directional difference tolerance threshold, then determines the target criteria by priority scoring, and rotates and corrects between conflicting criteria with a preset step size. This can improve the simultaneous satisfaction rate and reduce the number of invalid back-to-the-elevation fields. Even if a small number of points are truncated due to reaching the preset maximum number of iterations, quality marks can still be retained in the results table to facilitate downstream screening.

[0125] Sample Group 3 and Enhanced Dataset In sample group three, both methods can complete Z5 merging and generate an enhanced borehole formation dataset.

[0126] The method in this application still ensures that the interpolation neighborhood is searchable in sparse regions with a preset search radius, and improves the convergence ratio in conflict regions with a rotation mechanism. After merging, the proportion of entries marked with "not fully converged" in the dataset is lower than that of traditional methods.

[0127] The pass rate of the stratigraphic sequence validity check and the consistency of the measured borehole elevations with those of neighboring boreholes are both within an acceptable range, indicating that the enhanced dataset is usable for subsequent mapping and 3D modeling.

[0128] Comparative analysis A comparison of the method in this application with the traditional item-by-item verification-item-by-item correction method from a data processing perspective yields the following conclusions: In terms of handling multiple constraint conflicts, this application first generates conflict criterion pairs based on the preset deviation tolerance level and the preset direction difference tolerance threshold, and then rotates them according to priority and preset step size. Compared with the traditional item-by-item correction of conflict-free management, it is easier to achieve coordinated convergence, while having a higher satisfaction rate and a lower average number of iterations.

[0129] Regarding batch processing termination and resource controllability, the preset maximum number of iterations provides a unified forced exit condition for conflict rotation and single-item approximation; combined with the preset step size to constrain the update granularity, it can avoid infinite loops and reduce the time consumption of a single batch, which is more in line with the requirements of controllable computing resources for electronic digital data processing methods.

[0130] In terms of interpolation input scale control, the preset search radius restricts the measured records participating in the weighting to a local neighborhood, ensuring that the initial virtual borehole formation data can be generated while controlling the distance calculation complexity, thus providing a stable data foundation for subsequent verification and correction.

[0131] In terms of the usability of the results, the corrected virtual borehole formation data is merged with the initial measured data to form an enhanced borehole formation dataset, which can be accompanied by correction logs and quality tags, facilitating quality audits and differentiated use by downstream modules, thereby improving the efficiency of collaborative correction of virtual borehole formation interface elevations under multiple geological constraint conflicts.

[0132] It will be apparent to those skilled in the art that the present invention is not limited to the details of the exemplary embodiments described above, and that the present invention can be implemented in other specific forms without departing from the spirit or essential characteristics of the present invention.

[0133] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit it. Although the present invention has been described in detail with reference to preferred embodiments, those skilled in the art should understand that modifications or equivalent substitutions can be made to the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention.

Claims

1. A method for processing stratigraphic data from coalfield exploration boreholes, characterized in that, The method includes: Z1. Spatial interpolation is performed on the measured borehole data in the initial borehole formation dataset to generate initial virtual borehole formation data, and a multi-geological constraint verification system is established based on the initial borehole formation data. Z2. Based on the multi-geological constraint verification system, the formation data at each virtual borehole location is verified, the verification criteria that are not met are determined, and conflict detection is performed on multiple verification criteria that are not met. Verification criteria that have conflict relationships are marked as conflict criterion pairs. Z3. Based on the fault density data and fold amplitude data in the initial borehole formation data, determine the priority order of each verification criterion, take the current highest priority criterion as the target criterion, and adjust the formation interface elevation value until the target criterion is met. Z4. Substitute the adjusted formation interface elevation value back into the verification. If the other criterion in the conflict criterion pair is still not satisfied, switch to another criterion to continue the correction until all verification criteria are satisfied or the preset maximum number of iterations is reached. When only a single verification criterion is not satisfied, refer to the measured data corresponding to the single verification criterion to correct the formation interface elevation value until the single verification criterion is satisfied. Z5. Summarize the corrected formation interface elevation values ​​of all virtual borehole points to generate corrected virtual borehole formation data, and merge it with the initial borehole formation dataset to generate an enhanced borehole formation dataset.

2. The method for processing formation data from coalfield exploration boreholes as described in claim 1, characterized in that, The process of spatially interpolating the measured borehole data in the initial borehole formation dataset to generate initial virtual borehole formation data, and establishing a multi-geological constraint verification system based on the initial borehole formation data, includes: Extract the plane coordinate data of each measured borehole and the elevation value of each measured borehole at each target stratum interface from the initial borehole stratigraphic dataset of coalfield exploration boreholes; Based on the distribution density of the plane coordinate data of each measured borehole in the target interpolation area, virtual borehole points are set up in the target interpolation area, and the plane coordinates of each virtual borehole point are determined. For each virtual borehole location, using the planar coordinates of the virtual borehole location as the center, and according to a preset search radius, select the actual boreholes that are within the search radius and participate in the interpolation from each actual borehole. Based on the spatial distance between the measured boreholes involved in the interpolation and the virtual borehole locations, the influence weight of the elevation values ​​corresponding to the measured boreholes involved in the interpolation on the formation interface elevation values ​​at the virtual borehole locations is determined. The measured elevation values ​​of each borehole involved in the interpolation are multiplied by their respective influence weights to generate weighted elevation components. All the weighted elevation components are spatially superimposed according to the planar coordinates of the corresponding virtual borehole points to generate an initial elevation value; All initial elevation values ​​are collected to generate initial virtual borehole formation data; A multi-geological constraint verification system was established based on the initial borehole formation data.

3. The method for processing formation data from coalfield exploration boreholes as described in claim 2, characterized in that, The establishment of a multi-geological constraint verification system based on the initial borehole formation data includes: The allowable range of thickness variation in the dip and strike directions is determined based on the thickness variation range between adjacent measured boreholes, and the allowable range of dip angle variation rate in the dip and strike directions is determined based on the dip angle difference between adjacent measured boreholes. The allowable range of thickness variation and the allowable range of dip angle variation rate are used together to construct the formation thickness continuity constraint and the formation dip angle gradual change constraint; Based on the initial borehole stratigraphic data, fault occurrence data, and fault displacement data of the two sides of the fault, the spatial distribution range of the fault is defined, and the fault influence zone is generated by extending outward from the spatial distribution range of the fault. The spatial distribution range of the folds is defined based on the fold axis trace data of the initial borehole formation data, and the fold influence zone is generated by expanding outward from the spatial distribution range of the folds. Acquire spatial location data of the fault influence zone and the fold influence zone, determine whether there is a spatial overlap area between them, and if so, mark the spatial overlap area as a fault-fold composite influence zone. The fault-fold spatial coordination constraint is constructed by the fault influence zone, the fold influence zone, and the fault-fold composite influence zone. The stratigraphic thickness continuity constraint, the stratigraphic dip angle gradual change constraint, and the fault-fold spatial coordination constraint are associated and integrated according to their respective data type identifiers to generate a multi-geological constraint verification system.

4. The method for processing formation data from coalfield exploration boreholes as described in claim 1, characterized in that, The aforementioned multi-geological constraint verification system verifies the formation data at each virtual borehole location, identifies verification criteria that are not met, performs conflict detection on multiple non-metrization criteria, and marks conflicting verification criteria pairs as conflict criterion pairs, including: Substitute the formation data at each virtual borehole location into the various geological constraint verification criteria in the multi-geological constraint verification system to obtain deviation direction data and deviation magnitude level; Based on the deviation direction data and the deviation magnitude level, the verification criteria that exceed the preset deviation tolerance level are marked as unmet verification criteria, and the point identifier corresponding to the unmet verification criteria and the criterion identifier of the unmet verification criteria are recorded. When two or more verification criteria are not met at the same virtual borehole location, the correction direction data of each unmet verification criterion are paired up, and the directional difference between the two correction direction data in each pair is determined. When the directional difference exceeds the preset directional difference tolerance threshold, the two verification criteria in the pairing group are marked as having a conflict relationship; Extract all pairs of pairs that have conflicting relationships, and record the two verification criteria in the same pair as conflicting criterion pairs.

5. The method for processing formation data from coalfield exploration boreholes as described in claim 4, characterized in that, When two or more verification criteria are not met at the same virtual borehole location, the correction direction data of each unmet verification criterion are paired up, and the directional difference between the two correction direction data in each pair is determined, including: When two or more verification criteria are not met at the same virtual borehole location, the correction direction data of each unmet verification criterion are paired and combined to generate each pair group. The correction direction data includes a direction pointing component and a direction correction degree component. For each pairing group, the first direction pointing component of the pairing group is compared with the second direction pointing component of the pairing group. When the first direction pointing component and the second direction pointing component are the same, the first direction correction degree component of the pairing group is subtracted from the second direction correction degree component of the pairing group to generate the difference correction degree component. When the first direction pointing component is inconsistent with the second direction pointing component, the first direction correction degree component and the second direction correction degree component are added together to generate a sum correction degree component. When the first direction pointing component is consistent with the second direction pointing component, the difference correction degree component is taken as the direction difference of the pairing group; When the first direction pointing component is inconsistent with the second direction pointing component, the sum correction degree component is taken as the direction difference of the pairing group.

6. The method for processing formation data from coalfield exploration boreholes as described in claim 1, characterized in that, The step of determining the priority ranking of each verification criterion based on the fault density data and fold amplitude data in the initial borehole formation data, and adjusting the formation interface elevation value until the target criterion is met, using the current highest priority criterion as the target criterion, includes: The fault number, fault influence range, fold axis trace, and fold limb strata attitude data of each evaluation block in the coalfield exploration area are extracted from the initial borehole stratigraphic data to generate fault density and fold amplitude values. The fault density value and the fold amplitude value are summed to generate a priority score value for the fault-fold spatial coordination constraint. The thickness constraint correction coefficient and the dip angle constraint correction coefficient are determined based on the spatial matching relationship between fault density distribution data, formation thickness distribution data, and formation dip angle distribution data. The sum of the distribution variance of the fault density value and the distribution variance of the fold amplitude value is multiplied by the thickness constraint correction coefficient and the dip angle constraint correction coefficient, respectively, to generate priority score values ​​for the formation thickness continuity constraint and the formation dip angle gradual change constraint. Sort the priority scores from high to low to generate a priority sequence, take the geological constraint verification criterion with the highest score as the current highest priority criterion, and determine the current highest priority criterion as the target criterion; Based on the target criteria, the formation interface elevation values ​​at virtual borehole locations that have been marked with conflict criteria but do not meet the target criteria are adjusted until the target criteria are met.

7. A method for processing formation data from coalfield exploration boreholes as described in claim 6, characterized in that, The adjustment of formation interface elevation values ​​at virtual borehole locations marked with conflict criterion pairs but not meeting the target criterion, based on the target criterion, until the target criterion is met, includes: From the virtual borehole locations involved in the conflict criteria, the formation interface elevation values ​​at virtual borehole locations that do not meet the target criteria are selected as the elevation values ​​to be adjusted. Obtain the correction direction data corresponding to the target criterion, and perform a single incremental adjustment on the elevation value to be adjusted according to the adjustment direction indicated by the direction component of the correction direction data to generate a first intermediate elevation value; The first intermediate elevation value is substituted back into the target criterion for verification. If the verification passes, the first intermediate elevation value is used as the adjusted stratigraphic interface elevation value. If the verification fails, the first intermediate elevation value is used as the new elevation value to be adjusted, and the single incremental adjustment and back-substitution verification are repeated until the verification passes.

8. The method for processing formation data from coalfield exploration boreholes as described in claim 7, characterized in that, The adjusted formation interface elevation value is resubmitted for verification. If another criterion in the conflict criterion pair is still not met, the other criterion is used to continue the correction until all verification criteria are met or the preset maximum number of iterations is reached. When only a single verification criterion is not met, the formation interface elevation value is corrected with reference to the measured data corresponding to the single verification criterion until the single verification criterion is met, including: The adjusted stratigraphic interface elevation values ​​are substituted into the first and second criteria in the conflict criterion pair for verification. If both criteria are met, the correction is stopped. When the first criterion meets the verification requirements but the second criterion does not meet the verification requirements, the second criterion will be used as the current correction target. According to the adjustment direction indicated by the correction direction data and the preset step size, the formation interface elevation value at the virtual borehole point that does not meet the current correction object is incrementally adjusted to generate a new adjusted formation interface elevation value, and then the first criterion and the second criterion are substituted back for re-verification. If, after back substitution, the first criterion is satisfied but the second criterion is still not satisfied, then the second criterion remains the current correction target, and incremental adjustment and back substitution verification are repeated. If the first criterion is not met after back substitution, the current correction object is rotated to meet the first criterion and the incremental adjustment and back substitution verification are repeated. If both are met, the correction is stopped. If the cumulative number of rotations and incremental adjustments does not simultaneously meet the preset maximum number of iterations, then the correction will stop, and the formation interface elevation value after the last adjustment will be used as the final formation interface elevation value. When only a single verification criterion is not met, the formation interface elevation value is corrected by referring to the measured data corresponding to the single verification criterion until the single verification criterion is met.

9. A method for processing formation data from coalfield exploration boreholes as described in claim 8, characterized in that, When only a single verification criterion is not met, the formation interface elevation value is corrected by referring to the measured data corresponding to the single verification criterion until the single verification criterion is met, including: When the verification result is that only a single verification criterion is not met and no conflicting criterion pair is formed, the virtual borehole point that does not meet the single verification criterion is taken as the point to be corrected, and the formation interface elevation value at the point to be corrected is taken as the current elevation value to be corrected. Within the evaluation block to which the point to be corrected belongs, select the measured borehole closest to the point to be corrected according to the planar distance, and use the measured elevation value of the measured borehole at the same target stratum interface as the reference elevation value. The adjustment direction is determined based on the height relationship between the current elevation value to be corrected and the reference elevation value; Along the adjustment direction, a single incremental approximation adjustment is performed on the current elevation value to be corrected to generate the adjusted formation interface elevation value. The approximation adjusted formation interface elevation value is substituted back into the single verification criterion for verification. If the verification passes, the approximation adjusted formation interface elevation value is used as the corrected formation interface elevation value. If the verification fails, the adjusted formation interface elevation value will be used as the new current elevation value to be corrected. Repeat the single incremental approximation adjustment and back-substitution verification until the verification passes, or until the number of approximation adjustments reaches the preset maximum number of iterations.

10. A method for processing formation data from coalfield exploration boreholes as described in claim 1, characterized in that, The process involves summarizing all corrected formation interface elevation values ​​from virtual borehole locations to generate corrected virtual borehole formation data, which is then merged with the initial borehole formation dataset to generate an enhanced borehole formation dataset, including: For the same virtual borehole location, the corrected stratigraphic interface elevation values ​​of each target stratigraphic interface are arranged in order of stratigraphic age from oldest to newest, generating stratigraphic columnar sequence data. All stratigraphic column sequence data are merged according to the location identifiers of each virtual borehole point to generate corrected virtual borehole stratigraphic data. The stratigraphic column sequence data of each virtual borehole point in the corrected virtual borehole stratigraphic data is added as a new data entry to the same dataset containing the measured stratigraphic column sequence data of each measured borehole in the initial borehole stratigraphic dataset, thereby generating a merged stratigraphic dataset. Based on the spatial distribution of the planar coordinates of each measured borehole and each virtual borehole point within the target area of ​​the coalfield exploration, all data entries in the merged stratigraphic dataset are spatially sorted to generate an enhanced borehole stratigraphic dataset.

Images