CN115034151B - Fault vertical dominant migration channel analysis method and system - Google Patents
Fault vertical dominant migration channel analysis method and system Download PDFInfo
- Publication number
- CN115034151B CN115034151B CN202210515846.3A CN202210515846A CN115034151B CN 115034151 B CN115034151 B CN 115034151B CN 202210515846 A CN202210515846 A CN 202210515846A CN 115034151 B CN115034151 B CN 115034151B
- Authority
- CN
- China
- Prior art keywords
- point
- sampling point
- sampling
- fault
- cset
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Software Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Remote Sensing (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computer Graphics (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a fault vertical dominant migration channel analysis method and system, wherein a three-dimensional model of a fault is obtained, random sampling is carried out in the three-dimensional model of the fault to obtain a plurality of different sampling points, a fluid potential value and a fracture activity rate value at each sampling point are obtained, dominant migration search is carried out on each sampling point according to the fluid potential value and the fracture activity rate value at each sampling point to obtain a dominant migration point, a prediction region of the fault vertical dominant migration channel is obtained through calculation, and the beneficial effect of carrying out rapid and targeted fixed point on the fault vertical dominant migration channel is achieved.
Description
Technical Field
The invention belongs to the field of oil-gas exploration and data processing, and particularly relates to a fault vertical dominant migration channel analysis method and system.
Background
At present, 2 methods for researching oil and gas vertical dominant migration channels generally exist, firstly, a fault convex ridge is described through a section ridge method, namely, a fault surface buried depth contour map or a three-dimensional modeling method, the method is mainly based on fluid potential, secondly, a fault layer section with high activity rate is used as a dominant channel through an activity rate method, and the permeability is improved if the activity rate is high and the opening space is large. However, the problem exists that the section ridge is determined by the form of the fault plane, the fault can be developed at the positions of fault connection, bending, turning, tail end and the like, the positions of the fault connection and the tail end are all positions with low activity rate, the permeability is low, and the fault section with high activity rate calculated by an activity rate method cannot become a vertical dominant migration channel, so that the existing method has one-sidedness, and the fault section with high activity rate is mainly the middle section of one fault, so that the range is overlarge. In patent document CN114139330a, a method and an apparatus for determining a crude oil dominant migration channel are provided, although the crude oil dominant migration channel in the research area can be determined according to the effective hydrocarbon source rock distribution of each layer, the conductivity of the effective conductive sand body of each layer, the spatial distribution of the effective fault conductive body, and the conductivity effectiveness of each unconformity surface conductive body, the determined range is too dispersed and not accurate enough, which is not favorable for actual exploration and exploitation.
Disclosure of Invention
The present invention is directed to a method and system for analyzing a vertical dominant migration channel of a fault, which solves one or more of the problems of the prior art and provides at least one of the advantages of the present invention.
The invention provides a fault vertical dominant migration channel analysis method and system, which are used for obtaining a three-dimensional model of a fault, randomly sampling in the three-dimensional model of the fault to obtain a plurality of different sampling points, obtaining a fluid potential value and a fracture activity rate value at each sampling point, searching dominant migration of each sampling point according to the fluid potential value and the fracture activity rate value at each sampling point to obtain a dominant migration point, and calculating to obtain a prediction region of the fault vertical dominant migration channel.
In order to achieve the above object, according to an aspect of the present invention, there is provided a fault vertical dominant migration channel analysis method, including the steps of:
s100, acquiring a three-dimensional model of a fault, wherein data of each point in the three-dimensional model of the fault comprise a three-dimensional coordinate value, a fluid potential value and a fracture activity rate value of the point;
s200, randomly sampling in a three-dimensional model of a fault to obtain a plurality of different sampling points;
s300, acquiring a fluid potential value at each sampling point, and acquiring a fracture activity rate value at each sampling point;
s400, carrying out dominant migration search on each sampling machine point according to the fluid potential value and the fracture activity rate value at each sampling point to obtain a dominant migration point;
and S500, calculating to obtain a prediction area of the fault vertical dominant migration channel through the dominant migration point.
Further, in S100, the three-dimensional model of the fault is a three-dimensional model obtained by three-dimensional modeling of a geological structure of the region to be detected where the fault is located, where the source of the data for three-dimensional modeling is data obtained by a geological instrument or obtained by manual measurement, and the data includes a three-dimensional coordinate value, a fluid potential value, and a fracture activity rate value, and the types of the geological structure of the fault include a normal fault, a reverse fault, and a translational fault (where the geological instrument includes a multispectral camera, an infrared scanner, a multispectral scanner, and a microwave scanner).
The method for calculating the fluid potential value is described in document [1 ]: [1] luo Qun, pang Xiongji, jiang Zhenxue, a new method for effectively tracking oil and gas migration trajectory-the proposal and application of a section dominant migration channel [ J ]. Geological review, 2005, 51 (2): 7.;
the fracture activity rate value is calculated as described in document [2 ]: [2] zhou Xinhuai, niu Chengmin, teng Changyu the relation between the fracture activity and oil gas reservoir formation in the new tectonic movement period in the middle region of Bohai and Bohai [ J ]. Geology of oil and gas, 2009, 30 (4): 8;
the method of performing three-dimensional modeling can be referred to as: [3] yuan Ye, wang Li, xie Ruijie. Application of INPEFA technology in stratigraphic division-Subei basin Qitong recessed Nanhua block Qin three-phase example [ J ] Petroleum Experimental geology, 2018, 40 (6): 6., [4] Yuan Ye. Qitong recessed Nanhua block Qin three-phase reservoir geological modeling [ D ] Changjiang university.
Further, in S200, in the three-dimensional model of the fault, a method for randomly sampling to obtain a plurality of different sampling points includes: and randomly extracting k × 100 different sampling points in the three-dimensional model of the fault according to uniform distribution, wherein the actual geographic coordinates and the coordinates in the three-dimensional model of the fault correspond to each other one by one, k represents a positive integer, and k × 100 represents that k hundreds of different sampling points are extracted.
Further, in S300, the method for obtaining the fluid potential value at each sampling point and obtaining the fracture activity rate value at each sampling point specifically includes:
the set of sampling points is Cset, n represents the number of the sampling points, the numerical value of n is equal to k 100, the serial number of the sampling points is recorded as i, i belongs to [1,n ], and the sampling point with the serial number of i in the Cset is recorded as Cset (i);
the three-dimensional coordinates of the sampling point with the sequence number i in the Cset are loc (i), the X-axis coordinates in the loc (i) are X (i), the Y-axis coordinates in the loc (i) are Y (i), the Z-axis coordinates in the loc (i) are Z (i), and loc (i) = [ X (i), Y (i), Z (i) ];
the fluid potential value of the sampling point with the serial number i in the Cset is recorded as h (i), and the fracture activity rate value of the sampling point with the serial number i in the Cset is recorded as v (i);
Cset(i)=[loc(i), h(i), v(i)]。
further, in S400, according to the fluid potential value and the fracture activity rate value at each sampling point, a dominant migration search is performed on each sampling point, and the method for obtaining the dominant migration point specifically includes:
s401, searching for the dominant migration, wherein the process of searching for the dominant migration specifically comprises the following steps: acquiring each sampling point Cset (i) in Cset; go to S402;
s402, calculating the open square root value of the product of x (i), y (i) and z (i) of each sampling point Cset (i), selecting the sampling point with the minimum open square root value of the product of x (i), y (i) and z (i) in Cset, and acquiring the serial number of the sampling point as i1, i1 belongs to [1,n ], thereby recording the sampling point as Cset (i 1) by the serial number i 1; (the purpose of setting such a screening condition to obtain the sampling point Cset (i 1) is to select Cset (i 1) as a reference point for comparison with other sampling points;)
Further, cset (i 1) = [ loc (i 1), h (i 1), v (i 1) ], loc (i 1) = [ x (i 1), y (i 1), z (i 1) ]; go to S403;
and S403, calculating the fluid potential rate ratio of each sampling point Cset (i) according to the respective fluid potential value and the fracture activity rate value of each sampling point Cset (i), wherein the fluid potential rate ratio of the sampling point Cset (i) is Hv (i), and the calculation formula of Hv (i) is as follows:
alternatively, the first and second liquid crystal display panels may be,
in the alternative to this, either,
an exponential function is denoted by exp; the method has the advantages that the fluid potential rate ratio is calculated by combining the static element, namely the fluid potential, of the section ridge and the dynamic element, namely the fracture activity rate, firstly determining the section with high activity rate and simultaneously considering the fluid potential data characteristics of the fault, so that the position of the developed section ridge in the section with high activity rate is determined, and the method is more efficiently used for accurately positioning the position of the vertical dominant migration channel of the fault;
s404, calculating the convexity ratio of each sampling point Cset (i) according to the respective three-dimensional coordinates, wherein the convexity ratio represents a numerical value of the convex degree calculated according to the three-dimensional coordinates, the convexity ratio of the Cset (i) is recorded as cov (i), the convexity ratio cov (i) is calculated, and the calculation formula of the convexity ratio is as follows:
a cosine function is represented by a function cos (), and then, the process goes to S405;
s405, comparing the sampling points Cset (i) according to the magnitude of the value of the cov (i) corresponding to the sampling points Cset (i), selecting the sampling point with the smallest absolute value of the difference value between the value of the cov (i) corresponding to the sampling point Cset (i) and 1, and taking the sampling point as an initial sampling point;
s406, starting from the initial sampling point, respectively calculating Euclidean distances according to the three-dimensional coordinates of the initial sampling point and the three-dimensional coordinates of other sampling points except the initial sampling point in Cset, and selecting one sampling point with the largest Euclidean distance with the initial sampling point from the other sampling points as a first sampling point;
s407, connecting a line between the initial sampling point and the first sampling point, and taking the middle point of the connecting line between the initial sampling point and the first sampling point as a calibration point;
s408, calculating Euclidean distances between the Cset and the calibration point according to three-dimensional coordinates of the sampling points from other sampling points except the initial sampling point and the first sampling point, sequencing the sampling points according to the Euclidean distances between the sampling points and the calibration point from small to large, obtaining a sequence of sequenced sampling points (screening the Euclidean distances has the advantage of maximizing the search range of the initial sampling point and the other sampling points so as to accelerate the speed of searching and sampling), and turning to S409;
s409, judging whether the numerical value of k is greater than 1, if so, turning to S410, and if not, turning to S411;
s410, selecting the sampling points arranged in the front k as the superior moving points in the sequenced sampling point sequence, and obtaining a set formed by the k superior moving points as a set of the superior moving points; go to S412;
s411, taking a set consisting of the initial sampling point, the first sampling point and the 1 st sampling point in the sequenced sampling point sequence as an advantageous moving point set; go to S412;
s412, saving and outputting the dominant migration point set; finishing the dominant migration search;
the advantages of the advantage migration search are as follows: the method provided by the invention can be used for rapidly positioning a fixed point with high probability of the dominant migration channel as a dominant migration point by rapidly sampling and searching a three-dimensional model, calculating the data characteristics of a plurality of sampling points in a large-range checking manner and accurately positioning the dominant migration point with the probability of the dominant migration channel, and saving the calculation time cost while ensuring the accuracy.
Further, in S500, the method for obtaining the prediction region of the fault vertical dominant migration channel through calculation by using the dominant migration point includes:
according to each advantageous movement point in the advantageous movement point set, according to the mutual correspondence of the actual geographic coordinate and the coordinate in the three-dimensional model of the fault, the three-dimensional coordinate of each advantageous movement point has the actual geographic coordinate corresponding to the three-dimensional coordinate on the fault, the actual positioning coordinate is obtained on the to-be-detected region where the fault is located, the region where the actual positioning coordinate is located is the prediction region of the vertical advantageous movement channel of the fault, and the underground excavation in the prediction region can be used for the exploitation including the petroleum and the natural gas.
The invention also provides a fault vertical dominant migration channel analysis system, which comprises: the fault vertical dominant migration channel analysis system can be operated in computing equipment such as a desktop computer, a notebook computer, a palm computer, a cloud data center and the like, and can include, but is not limited to, a processor, a memory and a computer program stored in the memory and capable of being operated on the processor, wherein the processor executes the computer program and operates in units of the following systems:
the three-dimensional model acquisition unit is used for acquiring a three-dimensional model of a fault, and data of each point in the three-dimensional model of the fault comprise a three-dimensional coordinate value, a fluid potential value and a fracture activity rate value of the point;
the random sampling unit is used for carrying out random sampling in a three-dimensional model of a fault to obtain a plurality of different sampling points;
the numerical value acquisition unit is used for acquiring fluid potential numerical values at all sampling points and acquiring fracture activity rate numerical values at all sampling points;
the advantage migration searching unit is used for searching advantage migration of each sampling machine point according to the fluid potential value and the fracture activity rate value of each sampling point to obtain an advantage migration point;
and the prediction region acquisition unit is used for calculating and obtaining the prediction region of the fault vertical dominant migration channel through the dominant migration point.
The beneficial effects of the invention are as follows: the invention provides a fault vertical dominant migration channel analysis method and system, wherein a three-dimensional model of a fault is obtained, random sampling is carried out in the three-dimensional model of the fault to obtain a plurality of different sampling points, a fluid potential value and a fracture activity rate value at each sampling point are obtained, dominant migration search is carried out on each sampling point according to the fluid potential value and the fracture activity rate value at each sampling point to obtain a dominant migration point, a prediction region of the fault vertical dominant migration channel is obtained through calculation, and the beneficial effect of carrying out rapid and targeted fixed point on the fault vertical dominant migration channel is achieved.
Drawings
The above and other features of the present invention will become more apparent by describing in detail embodiments thereof with reference to the attached drawings in which like reference numerals designate the same or similar elements, it being apparent that the drawings in the following description are merely exemplary of the present invention and other drawings can be obtained by those skilled in the art without inventive effort, wherein:
FIG. 1 is a flow chart of a method for analyzing a vertically dominant migration path of faults;
FIG. 2 is a system diagram of a fault vertical dominant migration channel analysis system.
Detailed Description
The conception, the specific structure and the technical effects of the present invention will be clearly and completely described in conjunction with the embodiments and the accompanying drawings to fully understand the objects, the schemes and the effects of the present invention. It should be noted that the embodiments and features of the embodiments in the present application may be combined with each other without conflict.
In the description of the present invention, the meaning of a plurality of means is one or more, the meaning of a plurality of means is two or more, and larger, smaller, larger, etc. are understood as excluding the number, and larger, smaller, inner, etc. are understood as including the number. If the first and second are described for the purpose of distinguishing technical features, they are not to be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated or implicitly indicating the precedence of the technical features indicated.
Fig. 1 is a flow chart of a fault vertical dominant migration channel analysis method according to the present invention, and a fault vertical dominant migration channel analysis method and system according to an embodiment of the present invention are described below with reference to fig. 1.
The invention provides a fault vertical dominant migration channel analysis method, which specifically comprises the following steps:
s100, acquiring a three-dimensional model of the fault, wherein data of each point in the three-dimensional model of the fault comprise a three-dimensional coordinate value, a fluid potential value and a fracture activity rate value of the point;
s200, randomly sampling in a three-dimensional model of a fault to obtain a plurality of different sampling points;
s300, acquiring a fluid potential value at each sampling point, and acquiring a fracture activity rate value at each sampling point;
s400, carrying out dominant migration search on each sampling machine point according to the fluid potential value and the fracture activity rate value at each sampling point to obtain a dominant migration point;
and S500, calculating to obtain a prediction region of the fault vertical dominant migration channel through the dominant migration point.
Further, in S100, the three-dimensional model of the fault is a three-dimensional model obtained by performing three-dimensional modeling on a geological structure of the to-be-detected region where the fault is located, wherein the source of the data subjected to the three-dimensional modeling is data generated by a detection instrument or data acquired by manual measurement;
the method for calculating the fluid potential value is described in document [1 ]: [1] luo Qun, pang Xiongji, jiang Zhenxue, a new method for effectively tracking oil and gas migration trajectory-the proposal and application of a section dominant migration channel [ J ]. Geological review, 2005, 51 (2): 7.;
the fracture activity rate value is calculated as described in document [2 ]: [2] zhou Xinhuai, niu Chengmin, teng Changyu the relation between the fracture activity and oil gas reservoir formation in the new tectonic movement period in the middle region of Bohai and Bohai [ J ]. Geology of oil and gas, 2009, 30 (4): 8;
the method for performing three-dimensional modeling can be referred to as: [3] yuan Ye, wang Li, xie Ruijie. Application of INPEFA technology in stratigraphic division-Subei basin Qitong recessed Nanhua block Qin three-phase example [ J ] Petroleum Experimental geology, 2018, 40 (6): 6., [4] Yuan Ye. Qitong recessed Nanhua block Qin three-phase reservoir geological modeling [ D ] Changjiang university.
Further, in S200, in the three-dimensional model of the fault, the method of randomly sampling to obtain a plurality of different sampling points includes: the MCMC method is used for random sampling, and k x 100 different sampling points are randomly extracted according to uniform distribution in the three-dimensional model of the fault, wherein k represents a positive integer, k x 100 represents that k hundreds of different sampling points are extracted, and preferably k x 100 e (100, 1000) and k can be 3.
Further, in S300, the method for obtaining the fluid potential value at each sampling point and obtaining the fracture activity rate value at each sampling point specifically includes:
the set of sampling points is Cset, n represents the number of the sampling points, the numerical value of n is equal to k 100, the serial number of the sampling points is recorded as i, i belongs to [1,n ], and the sampling point with the serial number of i in the Cset is recorded as Cset (i);
the three-dimensional coordinates of the sampling point with the sequence number i in the Cset are loc (i), the X-axis coordinates in the loc (i) are X (i), the Y-axis coordinates in the loc (i) are Y (i), the Z-axis coordinates in the loc (i) are Z (i), and loc (i) = [ X (i), Y (i), Z (i) ];
the fluid potential value of the sampling point with the serial number i in the Cset is recorded as h (i), and the fracture activity rate value of the sampling point with the serial number i in the Cset is recorded as v (i);
Cset(i)=[loc(i), h(i), v(i)]。
further, in S400, according to the fluid potential value and the fracture rate value at each sampling point, the dominant migration search is performed on each sampling point, and the method for obtaining the dominant migration point specifically includes:
s401, starting to search for advantageous migration: calculating the convex surface rate of each sampling point Cset (i) in Cset, wherein the convex surface rate of Cset (i) is cov (i); go to S402;
s402, calculating the open square root value of the product of x (i) and y (i) and z (i) of each sampling point Cset (i), selecting the sampling point with the minimum open square root value of the product of x (i), y (i) and z (i) in Cset and obtaining the serial number of the sampling point as i1, i1 belongs to [1,n ], thereby marking the sampling point as Cset (i 1) by the serial number i 1;
further, cset (i 1) = [ loc (i 1), h (i 1), v (i 1) ], loc (i 1) = [ x (i 1), y (i 1), z (i 1) ]; go to S403;
and S403, respectively calculating the fluid potential rate ratio of each sampling point Cset (i) according to the respective fluid potential value and the fracture activity rate value of each sampling point Cset (i), recording the fluid potential rate ratio of the sampling point Cset (i) as Hv (i), wherein the calculation formula of Hv (i) is as follows:
an exponential function is denoted by exp;
s404, calculating cov (i) according to the respective three-dimensional coordinates of each sampling point Cset (i), wherein the calculation formula of cov (i) is as follows:
a cosine function is represented by a function cos (), and then, the process goes to S405;
s405, comparing the sampling points Cset (i) according to the numerical value of the cov (i) corresponding to the sampling points Cset (i), selecting a sampling point with the smallest numerical value distance 1 absolute value of the cov (i) corresponding to the sampling points Cset (i), and taking the sampling point as an initial sampling point;
s406, starting from the initial sampling point, respectively calculating Euclidean distances according to the three-dimensional coordinates of the initial sampling point and the three-dimensional coordinates of the other sampling points in the Cset, and selecting one sampling point with the largest Euclidean distance from the initial sampling point in the other sampling points in the Cset as a first sampling point;
s407, connecting a line between the initial sampling point and the first sampling point, and taking the middle point of the connection line between the initial sampling point and the first sampling point as a calibration point;
s408, calculating Euclidean distances between the Cset and the calibration point according to the three-dimensional coordinates of the sampling points from the other sampling points except the initial sampling point and the first sampling point, sequencing the sampling points according to the Euclidean distances between the sampling points and the calibration point from small to large, obtaining a sequence of sequenced sampling points, and turning to S409;
s409, judging whether the numerical value of k is greater than 1, if so, turning to S410, and if not, turning to S411;
s410, selecting k sampling points at the front of a row as dominant migration points in the sequenced sampling point sequence to obtain a dominant migration point set consisting of k dominant migration points; go to S412;
s411, selecting the 1 st sampling point, the initial sampling point and the first sampling point which are sequenced from the sequenced sampling point sequence as the dominant migration points, and obtaining a dominant migration point set consisting of k dominant migration points; go to S412;
s412, saving and outputting the dominant migration point set; and finishing the dominant migration search.
Further, in S500, the method for obtaining the prediction region of the fault vertical dominant migration channel through calculation by using the dominant migration point includes:
and acquiring actual positioning coordinates on the to-be-detected region where the fault is located according to each dominant migration point, wherein the region where the actual positioning coordinates are located is a prediction region of the vertical dominant migration channel of the fault.
The fault vertical dominant migration channel analysis system comprises: the fault vertical dominant migration channel analysis system can be operated in computing equipment such as a desktop computer, a notebook computer, a palm computer, a cloud data center and the like, and the operable system can include, but is not limited to, a processor, a memory and a computer program which is stored in the memory and can be operated on the processor.
As shown in fig. 2, the fault vertical dominant migration channel analysis system according to an embodiment of the present invention includes: a processor, a memory and a computer program stored in the memory and executable on the processor, the processor implementing the steps of one of the above embodiments of the fault vertical dominant migration channel analysis method when executing the computer program, the processor executing the computer program to run in the units of the following system:
the three-dimensional model acquisition unit is used for acquiring a three-dimensional model of a fault, and data of each point in the three-dimensional model of the fault comprise a three-dimensional coordinate value, a fluid potential value and a fracture activity rate value of the point;
the random sampling unit is used for carrying out random sampling in a three-dimensional model of a fault to obtain a plurality of different sampling points;
the numerical value acquisition unit is used for acquiring fluid potential numerical values at all sampling points and acquiring fracture activity rate numerical values at all sampling points;
the advantage migration searching unit is used for carrying out advantage migration searching on each sampling machine point according to the fluid potential value and the fracture activity rate value at each sampling point to obtain an advantage migration point;
and the prediction region acquisition unit is used for calculating and obtaining the prediction region of the fault vertical dominant migration channel through the dominant migration point.
The fault vertical dominant migration channel analysis system can be operated in computing equipment such as a desktop computer, a notebook computer, a palm computer and a cloud data center. The fault vertical dominant migration channel analysis system comprises, but is not limited to, a processor and a memory. It will be understood by those skilled in the art that the example is merely an example of a fault vertical dominant migration channel analysis method and system, and does not constitute a limitation on a fault vertical dominant migration channel analysis method and system, and may include more or less components than the fault vertical dominant migration channel analysis method and system, or some components in combination, or different components, for example, the fault vertical dominant migration channel analysis system may further include an input/output device, a network access device, a bus, and the like.
The Processor may be a Central Processing Unit (CPU), other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other Programmable logic device, discrete component Gate or transistor logic, discrete hardware components, etc. The general processor can be a microprocessor or the processor can be any conventional processor and the like, the processor is a control center of the fault vertical dominant migration channel analysis system, and various interfaces and lines are used for connecting various subareas of the whole fault vertical dominant migration channel analysis system.
The memory can be used for storing the computer program and/or module, and the processor can realize various functions of the fault vertical dominant migration channel analysis method and system by running or executing the computer program and/or module stored in the memory and calling data stored in the memory. The memory may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required by at least one function (such as a sound playing function, an image playing function, etc.), and the like; the storage data area may store data (such as audio data, a phonebook, etc.) created according to the use of the cellular phone, and the like. In addition, the memory may include high speed random access memory, and may also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), at least one magnetic disk storage device, a Flash memory device, or other volatile solid state storage device.
The invention provides a fault vertical dominant migration channel analysis method and system, which are used for acquiring a three-dimensional model of a fault, randomly sampling in the three-dimensional model of the fault to obtain a plurality of different sampling points, acquiring a fluid potential value and a fracture activity rate value at each sampling point, searching for dominant migration of each sampling point according to the fluid potential value and the fracture activity rate value at each sampling point to obtain a dominant migration point, and calculating to obtain a prediction area of a fault vertical dominant migration channel, so that the beneficial effect of performing rapid and targeted positioning on the fault vertical dominant migration channel is realized.
Although the present invention has been described in considerable detail and with reference to certain illustrated embodiments, it is not intended to be limited to any such details or embodiments or any particular embodiment, so as to effectively encompass the intended scope of the invention. Furthermore, the foregoing describes the invention in terms of embodiments foreseen by the inventor for which an enabling description was available, notwithstanding that insubstantial modifications of the invention, not presently foreseen, may nonetheless represent equivalent modifications thereto.
Claims (4)
1. A fault vertical dominant migration channel analysis method is characterized by comprising the following steps:
s100, acquiring a three-dimensional model of the fault, wherein data of each point in the three-dimensional model of the fault comprise a three-dimensional coordinate value, a fluid potential value and a fracture activity rate value of the point;
s200, randomly sampling in a three-dimensional model of a fault to obtain a plurality of different sampling points;
s300, acquiring a fluid potential value at each sampling point, and acquiring a fracture activity rate value at each sampling point;
s400, carrying out advantage migration search on each sampling machine point according to the fluid potential value and the fracture activity rate value at each sampling point to obtain an advantage migration point;
s500, calculating to obtain a prediction area of a fault vertical dominant migration channel through the dominant migration point;
in S200, in the three-dimensional model of the fault, the method of randomly sampling to obtain a plurality of different sampling points includes: randomly extracting k × 100 different sampling points in the three-dimensional model of the fault according to uniform distribution, wherein the actual geographic coordinates and the coordinates in the three-dimensional model of the fault correspond to each other one by one, k represents a positive integer, and k × 100 represents that k hundreds of different sampling points are extracted;
in S300, the method for obtaining the fluid potential value at each sampling point and obtaining the fracture activity rate value at each sampling point specifically includes:
the set of sampling points is Cset, n represents the number of the sampling points, the numerical value of n is equal to k 100, the serial number of the sampling points is recorded as i, i belongs to [1,n ], and the sampling point with the serial number of i in Cset is recorded as Cset (i);
the three-dimensional coordinates of the sampling point with the sequence number i in the Cset are loc (i), the X-axis coordinates in the loc (i) are X (i), the Y-axis coordinates in the loc (i) are Y (i), the Z-axis coordinates in the loc (i) are Z (i), and loc (i) = [ X (i), Y (i), Z (i) ];
the fluid potential value of the sampling point with the serial number i in the Cset is recorded as h (i), and the fracture activity rate value of the sampling point with the serial number i in the Cset is recorded as v (i);
Cset(i)=[loc(i), h(i), v(i)]；
in S400, according to the fluid potential value and the fracture activity rate value at each sampling point, performing dominant migration search on each sampling point, and obtaining a dominant migration point specifically includes:
s401, starting to perform dominant migration search, specifically including: acquiring each sampling point Cset (i) in Cset; go to S402;
s402, calculating the open square root value of the product of x (i) and y (i) and z (i) of each sampling point Cset (i), selecting the sampling point with the minimum open square root value of the product of x (i), y (i) and z (i) in Cset and obtaining the serial number of the sampling point as i1, i1 belongs to [1,n ], thereby marking the sampling point as Cset (i 1) by the serial number i 1;
further, cset (i 1) = [ loc (i 1), h (i 1), v (i 1) ], loc (i 1) = [ x (i 1), y (i 1), z (i 1) ]; go to S403;
and S403, calculating the fluid potential rate ratio of each sampling point Cset (i) according to the respective fluid potential value and the fracture activity rate value of each sampling point Cset (i), wherein the fluid potential rate ratio of the sampling point Cset (i) is Hv (i), and the calculation formula of Hv (i) is as follows:
an exponential function is denoted by exp;
s404, calculating the convexity ratio of each sampling point Cset (i) according to the respective three-dimensional coordinates, wherein the convexity ratio represents a numerical value of the convex degree calculated according to the three-dimensional coordinates, the convexity ratio of the Cset (i) is recorded as cov (i), the convexity ratio cov (i) is calculated, and the calculation formula of the convexity ratio is as follows:
a cosine function is represented by a function cos (), and then, the process goes to S405;
s405, comparing the sampling points Cset (i) according to the magnitude of the value of the cov (i) corresponding to the sampling points Cset (i), selecting the sampling point with the smallest absolute value of the difference value between the value of the cov (i) corresponding to the sampling point Cset (i) and 1, and taking the sampling point as an initial sampling point;
s406, starting from the initial sampling point, respectively calculating Euclidean distances according to the three-dimensional coordinates of the initial sampling point and the three-dimensional coordinates of the other sampling points in Cset, and selecting one of the other sampling points with the largest Euclidean distance from the initial sampling point as a first sampling point;
s407, connecting a line between the initial sampling point and the first sampling point, and taking the middle point of the connection line between the initial sampling point and the first sampling point as a calibration point;
s408, calculating Euclidean distances between the Cset and the calibration point according to the three-dimensional coordinates of the sampling points from the other sampling points except the initial sampling point and the first sampling point, sequencing the sampling points according to the Euclidean distances between the sampling points and the calibration point from small to large, obtaining a sequence of sequenced sampling points, and turning to S409;
s409, judging whether the numerical value of k is greater than 1, if so, turning to S410, and if not, turning to S411;
s410, selecting k sampling points in front of the sequence as dominant migration points in the sequenced sampling point sequence to obtain a dominant migration point set consisting of k dominant migration points; go to S412;
s411, taking a set consisting of the initial sampling point, the first sampling point and the 1 st sampling point in the sequenced sampling point sequence as an advantageous moving point set; go to S412;
s412, saving and outputting the dominant migration point set; and finishing the dominant migration search.
2. The fault vertical dominant migration channel analysis method according to claim 1, wherein in S100, the three-dimensional model of the fault is a three-dimensional model obtained by three-dimensional modeling of a geological structure of a region to be detected where the fault is located, and the data subjected to three-dimensional modeling is obtained by a geological instrument or by manual measurement and comprises a three-dimensional coordinate value, a fluid potential value and a fracture activity rate value.
3. The method for analyzing the vertical dominant migration channel of the fault as claimed in claim 2, wherein in S500, the method for calculating the predicted region of the vertical dominant migration channel of the fault from the dominant migration point comprises:
and acquiring actual positioning coordinates on the to-be-detected region where the fault is located according to each dominant migration point, wherein the region where the actual positioning coordinates are located is a prediction region of the vertical dominant migration channel of the fault.
4. A fault vertical dominant migration channel analysis system, comprising: a processor, a memory and a computer program stored in the memory and running on the processor, the processor implementing the steps in a fault vertical dominant migration corridor analysis method according to any one of claims 1 to 3 when executing the computer program, the fault vertical dominant migration corridor analysis system running in a computing device of a desktop computer, a notebook computer, a palm computer and a cloud data center.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210515846.3A CN115034151B (en) | 2022-05-12 | 2022-05-12 | Fault vertical dominant migration channel analysis method and system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210515846.3A CN115034151B (en) | 2022-05-12 | 2022-05-12 | Fault vertical dominant migration channel analysis method and system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115034151A CN115034151A (en) | 2022-09-09 |
| CN115034151B true CN115034151B (en) | 2023-02-03 |
Family
ID=83121184
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210515846.3A Active CN115034151B (en) | 2022-05-12 | 2022-05-12 | Fault vertical dominant migration channel analysis method and system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115034151B (en) |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9542511B2 (en) * | 2013-12-27 | 2017-01-10 | Halliburton Energy Services, Inc. | Synthetic gas-oil-ratio determination for gas dominant fluids |
| CN112489208A (en) * | 2019-09-11 | 2021-03-12 | 中国石油化工股份有限公司 | Crack sheet extraction method based on ant algorithm and three-dimensional geological model construction method |
| CN114139330B (en) * | 2020-09-04 | 2022-11-04 | 中国石油天然气股份有限公司 | Method and device for determining dominant migration channel of crude oil |
-
2022
- 2022-05-12 CN CN202210515846.3A patent/CN115034151B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN115034151A (en) | 2022-09-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Huisman et al. | Principles of geographic information systems | |
| Filzmoser et al. | Identification of local multivariate outliers | |
| Liu et al. | Multivariate analysis by data depth: descriptive statistics, graphics and inference,(with discussion and a rejoinder by liu and singh) | |
| CN110414723A (en) | The method, apparatus and system of fractured hydrocarbon reservoir history matching based on microseismic event | |
| Liu et al. | Petrophysical characteristics and log identification of lacustrine shale lithofacies: A case study of the first member of Qingshankou Formation in the Songliao Basin, Northeast China | |
| CN113420768A (en) | Core category determination method and device, electronic equipment and storage medium | |
| WO2021036780A1 (en) | Three-dimensional collection method and apparatus for magnetotelluric data, and terminal device | |
| Carvalho et al. | Soil classification system from cone penetration test data applying distance-based machine learning algorithms | |
| WO2012026923A1 (en) | Angle invariant visual analysis | |
| CN112711646B (en) | Ore finding method and device based on ground science information, electronic equipment and storage medium | |
| CN110675498A (en) | Long and large strip-shaped three-dimensional geological environment modeling method and device and terminal equipment | |
| CN115034151B (en) | Fault vertical dominant migration channel analysis method and system | |
| CN110555352A (en) | interest point identification method, device, server and storage medium | |
| CA2965867C (en) | Determining non-linear petrofacies using cross-plot partitioning | |
| CN111221038A (en) | Method and device for quantitatively predicting thickness of thin reservoir | |
| CN113838078B (en) | Method and device for identifying and extracting coal mining subsidence ground cracks and storage medium | |
| CN115358327A (en) | Ocean thermocline data visualization method, device, equipment and medium based on PCA-SVM | |
| Raith et al. | Visual Eddy Analysis of the Agulhas Current. | |
| CN115186074A (en) | Meta analysis-based method for simulating spatial distribution pattern of pH value of soil | |
| CN114743395A (en) | Signal lamp detection method, device, equipment and medium | |
| CN108053411B (en) | Remote sensing detection positioning method for river bottom tunnel under boundary heat exchange constraint | |
| CN107133360B (en) | Construction method of large-scale remote sensing image feature point library | |
| CN108108512B (en) | Reservoir lithofacies characterization method and device | |
| CN114998560B (en) | Positive fault section three-dimensional attribute-based related wrinkle processing method and system | |
| Wu et al. | Quantification of Connectivity in Images |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |