Fault stress plugging coefficient plane characterization rapid implementation method and device, medium and equipment

Abstract

The invention relates to a fault stress plugging coefficient plane characterization rapid implementation method and device, a medium and equipment. The method includes reading a regional stress direction angle and a fault polygon orientation data file. And calculating the fault stress plugging coefficient of each fault according to the regional stress direction angle and the fault polygonal orientation data file. And forming a plane graph by gridding the fault stress plugging coefficient of each fault. And reading the corresponding fault plane polygon, and externally tangent the plane graph by using the corresponding fault plane polygon to obtain a stress plugging coefficient plane graph reflecting the form of the fault plane. And respectively obtaining stress plugging coefficient plane graphs of different faults and combining the stress plugging coefficient plane graphs to obtain a stress plugging coefficient plane graph describing all the faults. According to the method, calculation is simple, only regional stress direction information needs to be provided, drilled well information does not need to be obtained, the fault plugging performance can be directly analyzed in a well-free and few-well area and an area with fast continental facies stratum phase change, and calculation efficiency is high.

Description

TECHNICAL FIELD

[0001] The present application relates to the offshore oil technology field, and in particular to a fault stress plugging coefficient plane characterization fast implementation method, device, medium and equipment. BACKGROUND

[0002] With the gradual deepening of oil and gas exploration, the exploration objects are more and more complex, the deep fault activity is developed, and the fault block is complex. The fault plugging evaluation is an indispensable link for the optimization of complex fault blocks and the evaluation of trap effectiveness.

[0003] The current fault plugging analysis method includes qualitative and quantitative analysis methods. The qualitative analysis method mainly includes lithology butt joint plugging analysis, acoustic travel time analysis method, physical and chemical property indication method, oil-water interface analysis method, fault zone filler analysis method, fault activity analysis method and the like. These methods can qualitatively judge the plugging of the fault, but cannot quantify the plugging capacity of the fault. The quantitative analysis method mainly includes: mudstone smearing coefficient (CSP) analysis method, mudstone smearing factor (CSSF) analysis method, fault mud ratio SGR analysis method and stratum model calculation method based on stress and cross section normal stress plugging analysis method. Among them, the stratum model calculation method is to obtain the continuous ground stress profile along the wellbore by using logging data, and then to model and calculate by using the data obtained by the experimental method to predict the ground stress. The corresponding ground stress is calculated by combining the logging data to establish the stratum model, so as to quantitatively evaluate the plugging of the fault to the oil and gas. The cross section normal stress plugging analysis method is to calculate the normal stress of the fault surface according to the fact that the hydrostatic pressure in the stratum rock does not affect the plugging of the fault, and only the gravity of the overlying stratum rock skeleton produces pressure, so as to quantitatively evaluate the plugging of the cross section.

[0004] The above-mentioned several quantitative fault plugging analysis methods have the following deficiencies: (1) the quantitative analysis methods based on lithology (CSP, CSSF, SGR) and stress (stratum model calculation method, cross section normal stress plugging analysis method) all need to obtain the drilled well data, i.e. the lithology, acoustic wave or density of the stratum, and the like. The actual production research target area is often deep and lacks drilling data, so it cannot be used to analyze the fault plugging. (2) For continental sedimentation, the stratum sedimentation characteristics change quickly laterally, and the stratum lithology (lithology, acoustic wave or density) information of a single well cannot replace the regional stratum information, so the reliability of the fault plugging analysis results of the above-mentioned methods cannot be guaranteed. (3) The calculation efficiency is low: the CSP and CSSF analysis methods need to manually pick up the fault throw, the SGR method and the stratum model calculation method need to combine the three-dimensional modeling of the seismic data, and the cross section normal stress method not only needs three-dimensional modeling, but also needs to accurately obtain multiple calculation parameters such as the depth corresponding to each point of the cross section, the average density of the overlying stratum, the stratum water density and the cross section dip angle. SUMMARY

[0005] The technical problem solved by the present application is to provide a fault stress plugging coefficient plane characterization fast implementation method, device, medium and equipment to solve the above-mentioned defects.

[0006] The technical solution adopted by the present application to solve its technical problem is: a fault stress plugging coefficient plane characterization fast implementation method, comprising the following steps: S1, reading regional stress direction angle and fault polygon orientation data file; S2, calculating the fault stress plugging coefficient of each fault according to the regional stress direction angle and the fault polygon orientation data file; the fault stress plugging coefficient is a stress plugging ability characterization of the fault on the plane; S3, forming a plane graph by gridding the fault stress plugging coefficient of each fault; S4, reading the corresponding fault plane polygon, and using the corresponding fault plane polygon to circumscribe the above-mentioned plane graph to obtain a stress plugging coefficient plane graph reflecting the plane shape of the fault; S5, obtaining different fault stress plugging coefficient plane graphs respectively and merging them to obtain a stress plugging coefficient plane graph depicting all faults.

[0007] Further, in the fault stress plugging coefficient plane characterization fast implementation method described in the present application, the fault polygon orientation data file is a.txt text format data file, and the file data contains three columns, the first column is the horizontal coordinate, the second column is the vertical coordinate, and the third column is the fault indicator, wherein the fault indicator is represented by a number, different numbers represent different faults, and the same fault indicator position is continuous in sequence, and step S2 comprises: calculating a regional stress unit vector according to the regional stress direction angle, and determining the number of faults according to the fault indicator; The fault stress plugging coefficient calculation step for each fault comprises: selecting the fault polygon to be calculated, determining the fault horizontal coordinate vector and the fault vertical coordinate vector; resampling the fault horizontal coordinate vector and the fault vertical coordinate vector to obtain the resampled fault horizontal coordinate vector and the fault vertical coordinate vector; calculating the tangent vector, slope, second derivative and curvature of each sampling point according to the resampled fault horizontal coordinate vector and the fault vertical coordinate vector; calculating the normal vector according to the tangent vector and the curvature; calculating the included angle between the normal vector and the regional stress unit vector according to the regional stress unit vector and the normal vector; calculating the fault stress plugging coefficient of each sampling point of the fault according to the second derivative, the slope, the included angle between the normal vector and the regional stress unit vector.

[0008] Further, in the fault stress sealing coefficient plane representation rapid implementation method, the fault stress sealing coefficient calculation step for each fault comprises: The regional stress unit vector is Wherein, , ; wherein, is the regional stress direction angle; Read the fault polygon disc data file fault indicator to determine the number of faults, recorded as NumFlt; and determine the position of different indicators in the fault polygon disc data file, recorded to the matrix T, wherein the first list in the ith row of the matrix T indicates the position of the first indicator in the fault data file, and the second list in the ith row of the matrix T indicates the position of the last indicator in the fault data file; Starting from the first fault, the fault horizontal coordinate vector Fx and the fault vertical coordinate vector Fy of the first fault are read from the fault file according to the position information of the first fault indicator in the matrix T; A preset resampling rate is adopted, and a cubic spline interpolation method is used to resample the fault horizontal coordinate vector Fx and the fault vertical coordinate vector Fy, to obtain the resampled fault horizontal coordinate vector Fltx and the fault vertical coordinate vector Flty.

[0009] Further, in the fault stress sealing coefficient plane representation rapid implementation method, the fault stress sealing coefficient calculation step for each fault comprises: S2-1, according to the resampled fault horizontal coordinate vector Fltx and the fault vertical coordinate vector Flty, calculate the tangent vector, slope, second derivative and curvature of each sampling point of the fault; Wherein, The difference between the adjacent horizontal coordinates of the ith sampling point is: Vx(i)=Fltx(i+1)-Fltx(i); The difference between the adjacent vertical coordinates of the ith sampling point is: Vy(i)=Flty(i+1)-Flty(i); The slope of the ith sampling point is: k(i)=Vy(i) / Vx(i); The second derivative of the ith sampling point is: dk(i)=(k(i+1)-k(i)) / Vx(i); The curvature of the ith sampling point is: c(i)=dk(i) / (1+k(i) 2 ) 3 / 2 ; The tangent vector of the ith sampling point is: ; S2-2, determining the i-th sampling point normal vector: If the curvature c(i)≥0, the i-th sampling point normal vector is: ; If the curvature c(i)<0, the i-th sampling point normal vector is: ; S2-3, calculating the angle between the regional stress and the fault point normal vector : Using the Euclidean dot product formula to express the cosine value between two vectors, the cosine value between the regional stress unit vector and the fault point normal vector is , and the angle between the regional stress and the fault point tangent is ; S2-4, calculating the fault stress sealing coefficient of each sampling point of the fault : When the second derivative dk(i)<0, if the slope k(i)≥0, ; otherwise ; When the second derivative dk(i)≥0, if the slope k(i)≥0, ; otherwise .

[0010] Further, in the fault stress sealing coefficient plane representation rapid implementation method described in the application, step S2 further comprises: Outputting a fault stress sealing coefficient text file; the data in the file contains three columns, the first column is the horizontal coordinate of each sampling point, the second column is the vertical coordinate of each sampling point, and the third column is the corresponding fault stress sealing coefficient of each sampling point.

[0011] Further, in the fault stress sealing coefficient plane representation rapid implementation method described in the application, step S3 comprises: Reading the fault stress sealing coefficient text file, gridding the data in the file, and forming a plane graph.

[0012] Further, in the fault stress sealing coefficient plane representation rapid implementation method described in the application, step S1 further comprises: setting the regional stress angle of the north direction as 0°, and rotating counterclockwise, the stress angle range is 0-360°.

[0013] In addition, the application also provides a fault stress sealing coefficient plane representation rapid implementation device, comprising: A data reading module for reading the regional stress direction angle and the fault polygon disk data file; The fault stress plugging coefficient calculation module is configured to calculate the fault stress plugging coefficient of each fault according to the regional stress direction angle and the fault polygon data file; the fault stress plugging coefficient is used to represent the stress plugging capability of the fault in a plane; The fault stress plugging coefficient plane fast mapping module is configured to form a plane graph by gridding the fault stress plugging coefficient of each fault; read the corresponding fault plane polygon, and use the corresponding fault plane polygon to circumscribe the above plane graph to obtain a stress plugging coefficient plane graph reflecting the plane shape of the fault; obtain different fault stress plugging coefficient plane graphs respectively and combine them to obtain a stress plugging coefficient plane graph depicting all faults.

[0014] In addition, the present application also provides a computer readable storage medium, which stores a computer program, and the computer program is suitable for being loaded by a processor to execute the steps of the fault stress plugging coefficient plane representation fast implementation method.

[0015] In addition, the present application also provides a computer device, which comprises a memory and a processor, and the memory stores a computer program; the processor executes the steps of the fault stress plugging coefficient plane representation fast implementation method by calling the computer program stored in the memory.

[0016] The fault stress plugging coefficient plane representation fast implementation method, device, medium and equipment have the following beneficial effects: the method makes the calculation of the fault stress plugging coefficient simple, can analyze the influence of the fault trend and shape on the fault plugging, only needs to provide the regional stress direction information, does not need to obtain the drilled well information, can directly analyze the fault plugging in the area without wells or with few wells and the area with fast facies change in continental strata, has high calculation efficiency, can quickly and efficiently calculate the stress plugging coefficients of multiple faults, and quickly realizes the plane representation and mapping of the stress plugging coefficients. BRIEF DESCRIPTION OF DRAWINGS

[0017] The present application will be further described below in combination with the drawings and embodiments, and the drawings are as follows: Figure 1 FIG. 1 is a flowchart of the fault stress plugging coefficient plane representation fast implementation method provided by the embodiments of the present application; Figure 2 FIG. 2 is a plane schematic diagram of the regional stress acting on the fault; Figure 3 FIG. 3 is a structure schematic diagram of the fault stress plugging coefficient plane representation fast implementation device provided by the embodiments of the present application; Figure 4 FIG. 4 is a fault stress plugging coefficient plane graph of a certain area in the east of the South China Sea provided by some embodiments of the present application. DETAILED DESCRIPTION

[0018] In order to make the technical features, objectives and effects of the present application more clearly understood, the specific embodiments of the present application will be described in detail with reference to the drawings. In the following description, it should be understood that the "front", "back", "upper", "lower", "left", "right", "vertical", "horizontal", "vertical", "horizontal", "top", "bottom", "inner", "outer", "head", "tail" and other indications of orientation or positional relationship are based on the orientation or positional relationship shown in the drawings, constructed and operated in a particular orientation, and are only for the convenience of describing the technical solutions, and do not indicate that the indicated device or element must have a particular orientation, therefore, it cannot be understood as a limitation on the present application.

[0019] It should also be noted that unless otherwise explicitly specified and limited, the terms "mounting", "connecting", "connecting", "fixing", "setting" and the like should be broadly understood, for example, it can be fixedly connected, or it can be detachably connected, or it can be integrated; it can be mechanically connected, or it can be electrically connected; it can be directly connected, or it can be indirectly connected through an intermediate medium, or it can be the internal communication of two elements or the interaction relationship between two elements. When an element is referred to as "on" or "below" another element, the element can be "directly" or "indirectly" above the other element, or there can be one or more intervening elements. The terms "first", "second", "third" and the like are only for the convenience of describing the technical solutions, and cannot be understood as indicating or implying relative importance or implicitly indicating the number of indicated technical features, therefore, the features with "first", "second", "third" and the like can be explicitly or implicitly included one or more of the features. For those skilled in the art, the specific meaning of the above terms in the present application can be understood according to the specific circumstances.

[0020] In the following description, specific details such as specific system structures, technologies, etc. are presented for the purpose of illustration, not for the purpose of limitation, so as to thoroughly understand the embodiments of the present application. However, it should be clear to those skilled in the art that the present application can also be implemented in other embodiments without these specific details. In other cases, detailed description of well-known systems, devices, circuits and methods is omitted to avoid unnecessary details that hinder the description of the present application.

[0021] As shown in Figure 1 In a preferred embodiment, the fault stress plugging coefficient plane characterization rapid implementation method of the present embodiment includes the following steps: S1. Read the regional stress direction angle and fault polygon disk orientation data files. It can be understood that disk orientation data represents the directional variation characteristics of the disk surface, i.e., the trend of data variation in a certain direction. In this application, the fault polygon disk orientation data is either fault polygon descending disk data or fault polygon ascending disk data. That is, both ascending and descending disk data can represent the fault's shape on the plane; only one needs to be selected, and it is not necessary to include both in the calculation. It should be noted that the regional stress angle in the due north direction is set to 0°, and rotating counterclockwise, the stress angle range is 0-360°.

[0022] It can be understood that when regional stress F acts on a point of the fault, it can be decomposed into a force F1 parallel to the fault strike (parallel to the tangent) and a force F2 perpendicular to the fault strike (parallel to the normal), such as... Figure 2 As shown. The force perpendicular to the strike of a fault exerts either a tensile or compressive effect on the strata. The sealing capacity of the fault can be determined based on the magnitude and direction of the vertical stress. If the vertical stress exhibits a compressive effect, the larger the value, the better the sealing performance of the fault. Definition The fault stress sealing coefficient is used to describe the stress sealing capacity of a fault in a plane, and its mathematical expression is: in When it exhibits a tensile effect, When it manifests as a squeezing effect, . It is the angle between stress F and force F2 perpendicular to the fault strike.

[0023] Before calculating the fault stress sealing coefficient, it is necessary to determine the regional stress direction angle θ and the fault plane file (fault polygon ascending or descending disk).

[0024] S2. Calculate the fault stress sealing coefficient for each fault based on the regional stress direction angle and the fault polygon disk orientation data file.

[0025] S3. The fault stress sealing coefficient of each fault is gridded to form a planar diagram.

[0026] S4. Read the corresponding fault plane polygon and use the corresponding fault plane polygon to circumferentially tangent the above plane diagram to obtain a stress sealing coefficient plane diagram that reflects the plane morphology of the fault.

[0027] S5. Obtain the stress sealing coefficient plane diagrams of different faults and merge them to obtain the stress sealing coefficient plane diagrams that characterize all faults.

[0028] This embodiment simplifies the calculation of fault stress sealing coefficients, allowing analysis of the impact of fault strike and morphology on fault sealing performance. It only requires regional stress direction information and does not need to obtain well drilling information. It can directly analyze fault sealing performance in areas with few or no wells and in regions with rapid facies changes in continental strata. It has high calculation efficiency, can quickly and efficiently calculate the stress sealing coefficients of multiple faults, and can quickly perform planar characterization and mapping.

[0029] In some embodiments, the fault polygon disk orientation data file is a .txt text file containing three columns: the first column is the x-axis, the second column is the y-axis, and the third column is the fault indicator. The fault indicator is represented by numbers, with different numbers representing different faults, and the positions of the same fault indicator are consecutive. It is understood that the text format of the fault polygon disk orientation data file can also be other text formats. This is merely an exemplary embodiment, and other text formats should also be within the scope of protection of this invention.

[0030] Step S2 includes: The regional stress unit vector is calculated based on the regional stress direction angle, and the number of faults is determined based on the fault indicator.

[0031] The calculation steps for the fault stress sealing coefficient of each fault include: selecting the fault polygon to be calculated and determining the fault's abscissa and ordinate vectors. Resampling the fault's abscissa and ordinate vectors to obtain resampled abscissa and ordinate vectors. Calculating the tangent vector, slope, second derivative, and curvature of each sampling point based on the resampled abscissa and ordinate vectors. Calculating the normal vector of each sampling point using the orthogonality between the tangent vector and the normal vector and the curvature value. Calculating the angle α between the normal vector and the regional stress unit vector based on the cosine similarity between the regional stress unit vector and the normal vector. Calculating the fault stress sealing coefficient for each sampling point of the fault based on the second derivative, slope, and the cosine value of the angle α between the normal vector and the regional stress unit vector.

[0032] In some embodiments, the calculation steps for the fault stress sealing coefficient of each fault include: The unit vector of regional stress is ,in , .in, The angle represents the direction of regional stress.

[0033] The number of faults, denoted as NumFlt, is determined by reading the fault indicators from the fault polygon disk-oriented data file. The positions of different indicators within the fault polygon disk-oriented data file are then recorded in matrix T. In matrix T, the first column of row i represents the first occurrence of the i-th indicator in the fault data file, and the second column of row i represents the last occurrence of the i-th indicator in the fault data file.

[0034] Starting from the first fault, based on the position information of the first fault indicator in matrix T, the fault horizontal coordinate vector Fx and fault vertical coordinate vector Fy of the first fault are read from the fault file.

[0035] A preset resampling rate is used, and cubic spline interpolation is employed to resample the fault abscissa vector Fx and fault ordinate vector Fy, resulting in the resampled fault abscissa vector Fltx and fault ordinate vector Flty.

[0036] The steps for calculating the fault stress sealing coefficient for each fault include: S2-1. Calculate the tangent vector, slope, second derivative, and curvature of each sampling point of the fault based on the resampled fault abscissa vector Fltx and fault ordinate vector Flty. in, The difference between the adjacent x-coordinates of the i-th sampling point is: Vx(i) = Fltx(i+1) - Fltx(i); The difference between the adjacent ordinates of the i-th sampling point is: Vy(i) = Flty(i+1) - Flty(i); The slope of the i-th sampling point is: k(i) = Vy(i) / Vx(i); The second derivative of the i-th sampling point is: dk(i) = (k(i+1) - k(i)) / Vx(i); The curvature of the i-th sampling point is: c(i) = dk(i) / (1 + k(i)) 2 ) 3 / 2 ; The cutting vector for the i-th sampling point is: ; S2-2, Determine the normal vector of the i-th sampling point: If curvature c(i) ≥ 0, then the normal vector of the i-th sampling point is: ; If curvature c(i) < 0, then the normal vector of the i-th sampling point is: ; S2-3, Calculate the angle between the stress in the calculation area and the normal vector at each point of the fault. : Using the Euclidean dot product formula to express the cosine value between two vectors, the cosine value of the unit vector of regional stress and the normal vector at each point of the fault is... The angle between the regional stress and the tangent at each point of the fault is... ; S2-4. Calculate the fault stress sealing coefficient at each sampling point of the fault. : When the second derivative dk(i) < 0, if the slope k(i) ≥ 0, ;otherwise ; When the second derivative dk(i)≥0, if the slope k(i)≥0, ;otherwise .

[0037] It can be understood that the fault stress sealing coefficient is calculated iteratively from the first fault to the NumFlt fault. When the resampling rate is set to 1000, the fault stress sealing coefficient for each fault is calculated iteratively from the first sampling point to the 1000th sampling point. The fault stress sealing coefficient at each sampling point is the fault stress sealing coefficient for that fault.

[0038] In some embodiments, step S2 further includes: Output a text file containing the fault stress sealing coefficients. This file contains three columns: the first column represents the x-coordinate of each sampling point, the second column represents the y-coordinate of each sampling point, and the third column represents the fault stress sealing coefficient for each sampling point.

[0039] In some embodiments, step S3 includes: The fault stress sealing coefficient text file is read in, the data in the file is gridded, and a planar diagram is generated.

[0040] refer to Figure 3 In another preferred embodiment, the device for rapid realization of plane characterization of fault stress sealing coefficient in this embodiment includes: The data reading module is used to read regional stress direction angles and fault polygon disk orientation data files.

[0041] The fault stress sealing coefficient calculation module is used to calculate the fault stress sealing coefficient of each fault based on the regional stress direction angle and the fault polygon disk orientation data file. The fault stress sealing coefficient characterizes the stress sealing capacity of a fault in a plane.

[0042] The fault stress sealing coefficient planar mapping module is used to generate a planar map by meshing the fault stress sealing coefficient of each fault. It reads the corresponding fault plane polygon and uses the corresponding fault plane polygon to circumferentially tangent to the above planar map to obtain a stress sealing coefficient planar map reflecting the planar morphology of that fault. Different fault stress sealing coefficient planar maps are obtained separately and then merged to obtain a stress sealing coefficient planar map characterizing all faults.

[0043] This embodiment simplifies the calculation of fault stress sealing coefficients, allowing analysis of the impact of fault strike and morphology on fault sealing performance. It only requires regional stress direction information and does not need to obtain well drilling information. It can directly analyze fault sealing performance in areas with few or no wells and in regions with rapid facies changes in continental strata. It has high calculation efficiency, can quickly and efficiently calculate the stress sealing coefficients of multiple faults, and can quickly perform planar characterization and mapping.

[0044] In another preferred embodiment, the computer-readable storage medium of this embodiment stores a computer program adapted for loading by a processor to execute the steps of the rapid implementation method for planar characterization of fault stress sealing coefficient as described in the above embodiment. This embodiment simplifies the calculation of fault stress sealing coefficient, allows analysis of the influence of fault strike and morphology on fault sealing performance, requires only regional stress direction information and does not need to obtain drilled well information, and can directly analyze fault sealing performance in areas with few or no wells and in regions with rapid facies changes in continental strata. It has high computational efficiency, can quickly and efficiently calculate the stress sealing coefficients of multiple faults, and rapidly achieve planar characterization and mapping of them.

[0045] In another preferred embodiment, the computer device of this embodiment includes a memory and a processor. The memory stores a computer program, and the processor executes the steps of the rapid implementation method for planar characterization of fault stress sealing coefficient as described in the above embodiment by calling the computer program stored in the memory. This embodiment simplifies the calculation of fault stress sealing coefficient, allows analysis of the influence of fault strike and morphology on fault sealing performance, requires only regional stress direction information and does not need to obtain drilled well information, and can directly analyze fault sealing performance in areas with few or no wells and in areas with rapid facies changes in continental strata. It has high computational efficiency, can quickly and efficiently calculate the stress sealing coefficients of multiple faults, and rapidly perform planar characterization mapping on them.

[0046] Figure 4This paper presents a planar map of the fault stress sealing coefficient in a region of the eastern South China Sea. This embodiment selects a fault-developed area in the Pearl River Estuary Basin and uses the proposed rapid plane characterization method for fault stress sealing analysis to analyze its fault sealing properties, thereby evaluating the exploration potential of the region. Fifty faults were selected for analysis. Fault plane polygons were obtained through seismic data interpretation, and these polygons were divided into ascending and descending disk data. The descending disk data was used in this calculation. The data was edited into a .txt text file containing three columns: the first column is the horizontal axis, the second column is the vertical axis, and the third column is the fault indicator (represented by numbers, with different numbers indicating different faults, and the positions of the same fault indicator being consecutive). Furthermore, statistical analysis revealed that the regional stress field direction in this area is 330°. First, the regional stress angle and the descending disk text file were read using the data reading module. Then, the stress sealing coefficient calculation module was used to calculate the stress sealing coefficient for each fault at the corresponding sampling point. Then, the fault stress sealing coefficient planar mapping module was used to complete the planar mapping of the fault stress sealing coefficients of 50 faults in this region, such as... Figure 4 As shown in the figure, the stress in the NE-trending fault exhibits compression characteristics, which is conducive to oil and gas plugging. Blocks A5 (0.4-0.7), C3-1 / 2 (0.4-0.7), D2 (0.4-0.7), D3-1 (0.5-0.65), and D3-2 / 3 (0.65-0.99) show better plugging effects.

[0047] The computer-readable storage medium of the present invention can be any computer-readable storage medium capable of storing program code, such as a USB flash drive, a portable hard drive, a read-only memory (ROM), a magnetic disk, or an optical disk.

[0048] The processor of this invention provides computing and control capabilities to support the operation of the entire device. It should be understood that, in the embodiments of this application, the processor may be a Central Processing Unit (CPU), or it may be other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or any conventional processor.

[0049] Those skilled in the art will further recognize that the units and algorithm steps of the various examples described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of both. To clearly illustrate the interchangeability of hardware and software, the components and steps of the various examples have been generally described in terms of functionality in the foregoing description. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementations should not be considered beyond the scope of this invention.

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

[0051] It is understood that the above embodiments only illustrate preferred embodiments of the present invention, and their descriptions are relatively specific and detailed, but they should not be construed as limiting the scope of the present invention. It should be noted that those skilled in the art can freely combine the above technical features without departing from the concept of the present invention, and can also make several modifications and improvements, all of which fall within the protection scope of the present invention. Therefore, all equivalent transformations and modifications made with respect to the scope of the claims of the present invention should fall within the scope of the claims of the present invention.

Claims

1. A method for quickly realizing a planar representation of fault stress sealing coefficient, characterized in that, The method comprises the following steps: S1, reading regional stress direction angle and fault polygon data file; S2, calculating fault stress sealing coefficient of each fault according to regional stress direction angle and fault polygon data file; The fault stress sealing coefficient is a stress sealing ability of fault on a plane; S3, forming a plan by gridding the fault stress sealing coefficient of each fault; S4, reading corresponding fault plane polygon, and using the corresponding fault plane polygon to circumscribe the above plan to obtain stress sealing coefficient plan reflecting the plane shape of the fault; S5, obtaining different fault stress sealing coefficient plan respectively and merging them to obtain stress sealing coefficient plan of all faults.

2. The method according to claim 1, wherein, The fault polygon data file is a.txt text format data file, and the file data comprises three columns, the first column is the horizontal coordinate, the second column is the vertical coordinate, and the third column is the fault indicator, wherein the fault indicator is represented by numbers, different numbers represent different faults, and the same fault indicator position is continuous in sequence, and step S2 comprises: According to the regional stress direction angle, a regional stress unit vector is calculated, and the number of faults is determined according to the fault indicator; The fault stress sealing coefficient calculation step for each fault comprises: Selecting the fault polygon to be calculated, determining the fault horizontal coordinate vector and the fault vertical coordinate vector; Resampling the fault horizontal coordinate vector and the fault vertical coordinate vector to obtain the resampled fault horizontal coordinate vector and the resampled fault vertical coordinate vector; Calculating the tangent vector, the slope, the second derivative and the curvature of each sampling point according to the resampled fault horizontal coordinate vector and the resampled fault vertical coordinate vector; Calculating the normal vector according to the tangent vector and the curvature; Calculating the included angle between the normal vector and the regional stress unit vector according to the regional stress unit vector and the normal vector; Calculating the fault stress sealing coefficient of each sampling point of the fault according to the second derivative, the slope, the included angle between the normal vector and the regional stress unit vector.

3. The method of claim 2, wherein the method further comprises: The fault stress sealing coefficient calculation step for each fault comprises: The regional stress unit vector is wherein , ; wherein is the regional stress direction angle; Reading the fault indicator of the fault polygon data file to determine the number of faults, recorded as NumFlt, and determining the position of different indicators in the fault polygon data file, recorded in matrix T, wherein the first column of the i-th row of matrix T indicates the position of the i-th indicator in the fault data file for the first time, and the second column of the i-th row of matrix T indicates the position of the i-th indicator in the fault data file for the last time; Starting from the first fault, reading the fault horizontal coordinate vector Fx and the fault vertical coordinate vector Fy of the first fault in the fault file according to the position information of the first fault indicator in the matrix T; Pre-setting a resampling rate, and resampling the fault horizontal coordinate vector Fx and the fault vertical coordinate vector Fy by using a cubic spline interpolation method to obtain the resampled fault horizontal coordinate vector Fltx and the resampled fault vertical coordinate vector Flty.

4. The method according to claim 3, wherein, The fault stress sealing coefficient calculation step for each fault comprises: S2-1, calculating the tangent vector, slope, second derivative and curvature of each sampling point of the fault according to the fault horizontal coordinate vector Fltx and the fault vertical coordinate vector Flty after resampling; Wherein, The difference between the adjacent horizontal coordinates of the i th sampling point is: Vx(i)=Fltx(i+1)-Fltx(i); The difference between the adjacent vertical coordinates of the i th sampling point is: Vy(i)=Flty(i+1)-Flty(i); The slope of the i th sampling point is: k(i)=Vy(i) / Vx(i); The second derivative of the i th sampling point is: dk(i)=(k(i+1)-k(i)) / Vx(i); The curvature of the ith sample point is: c(i) = dk(i) / (1 + k(i)) 2 ) 3 / 2 ; The tangent vector of the i-th sampling point is: ; S2-2, determining the normal vector of the i th sampling point: If the curvature c(i) > 0, then the i-th sample point normal vector is: ; If the curvature c(i) < 0, then the i-th sample point normal vector is: ; S2-3, calculating an angle between the regional stress and the fault point normal vector : Using the Euclidean dot product formula to express the cosine value between two vectors, the cosine value between the regional stress unit vector and the normal vector of each point of the fault is , and the angle between the regional stress and the tangent of each point of the fault is . S2-4, calculate the fault stress sealing coefficient of each sampling point of the fault : When the second derivative dk(i) < 0, if the slope k(i) ≥ 0, ; otherwise ; When the second derivative dk(i) ≥ 0, if the slope k(i) ≥ 0, ; otherwise .

5. The method of claim 2, wherein the method further comprises: Step S2 further comprises: Outputting a fault stress sealing coefficient text file; the data in the file contains three columns, the first column is the horizontal coordinate of each sampling point, the second column is the vertical coordinate of each sampling point, and the third column is the corresponding fault stress sealing coefficient of each sampling point.

6. The method according to claim 5, wherein, Step S3 comprises: Reading the fault stress sealing coefficient text file, gridding the data in the file, and forming a plan view.

7. The method of claim 1, wherein, Step S1 further comprises: setting the regional stress angle of the north direction as 0°, and rotating counterclockwise, with the stress angle range being 0-360°.

8. A device for rapid implementation of planar representation of fault stress seal coefficient, characterized in that, Comprise: A data reading module for reading regional stress direction angle and fault polygon orientation data file; A fault stress sealing coefficient calculation module for calculating the fault stress sealing coefficient of each fault according to the regional stress direction angle and the fault polygon orientation data file; The fault stress sealing coefficient is a stress sealing capacity of the fault on the plane; A fault stress sealing coefficient planar fast mapping module for forming a plan view by gridding the fault stress sealing coefficient of each fault; Reading the corresponding fault planar polygon and using the corresponding fault planar polygon to circumscribe the above plan view to obtain a stress sealing coefficient plan view reflecting the planar morphology of the fault; Different fault stress sealing coefficient plan views are obtained respectively and combined to obtain a stress sealing coefficient plan view depicting all faults.

9. A computer-readable storage medium, characterized in that, The computer readable storage medium stores a computer program, and the computer program is suitable for being loaded by a processor to execute the steps of the fault stress sealing coefficient planar characterization fast implementation method according to any one of claims 1-7.

10. A computer device, comprising: Comprise a memory and a processor, the memory stores a computer program, and the processor executes the steps of the fault stress sealing coefficient planar characterization fast implementation method according to any one of claims 1-7 by calling the computer program stored in the memory. Comprise a memory and a processor, the memory stores a computer program, and the processor executes the steps of the fault stress sealing coefficient planar characterization fast implementation method according to any one of claims 1-7 by calling the computer program stored in the memory.

Images