A data acquisition method and system applied to petroleum geological exploration
By calculating the potential coefficient, actual availability value, oil potential rate, and final potential rate, the coordinates of exploration points are optimized, solving the problems of complex exploration point layout and high cost in gravity exploration methods, and achieving efficient and accurate oil resource exploration.
Patent Information
- Application Number
- CN202510480567.1
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2025-04-17
- Publication Date
- 2025-12-09
- Estimated Expiration
- 2045-04-17
AI Technical Summary
In gravity exploration, the large number and lack of specificity of exploration points lead to complex equipment layout and operation processes, increasing costs.
By calculating the potential coefficient, true availability value, oil potential rate, correction degree, and final potential rate, the coordinate selection of exploration points is optimized, focusing on potential oil-rich areas, reducing the number of exploration points, improving data acquisition accuracy, and reducing costs.
It significantly improved the accuracy of data acquisition, reduced exploration costs, focused on potential oil-rich areas, and improved exploration efficiency.
Smart Images

Figure CN120428348B_ABST
Abstract
Description
TECHNICAL FIELD
[0001] The present application relates to the technical field of petroleum geological exploration, and particularly relates to a data collection method and system applied to petroleum geological exploration. BACKGROUND
[0002] In order to accurately reveal the distribution pattern, reserve scale and quality condition of underground petroleum resources, and provide a solid scientific basis for petroleum exploitation activities, petroleum geological exploration work is particularly crucial. Among numerous exploration technologies, gravity exploration method stands out due to its wide application and high efficiency in large-area regions, especially in the preliminary exploration stage. As a non-invasive exploration method, it analyzes geological information by precisely measuring the subtle changes of the underground gravity field, without the need for drilling operations, which not only has relatively low cost, but also has minimal disturbance to the environment.
[0003] In the process of petroleum geological exploration by using the gravity exploration method, the systematic collection of gravity data of multiple exploration points is a core link. Since low gravity value areas often indicate the existence of large underground cavities, cracks or low-density rock layers and other geological structures, these features are common in petroleum or natural gas reservoirs, therefore, exploration points with low gravity values are usually considered as important indicators of potential petroleum reservoirs. By accurately identifying these low gravity anomaly areas, valuable clues and basis can be provided for the potential location of petroleum reservoirs.
[0004] Currently, in the implementation process of the gravity exploration method, gravity measurement instruments need to operate on widely distributed exploration points, and the distance between exploration points is usually set to be small. Therefore, this method leads to the intensification of the complexity of equipment arrangement and operation process, and the data collection work faces great challenges, accompanied by significant cost increase. SUMMARY
[0005] In order to solve the technical problems of complex equipment arrangement and operation process and cost increase caused by the large number of exploration points and lack of pertinence in the prior art, the purpose of the present application is to provide a data collection method applied to petroleum geological exploration, and the technical scheme used is as follows:
[0006] An initial exploration point is determined in a target area of petroleum geological exploration, and surrounding exploration points are determined at a preset distance from the initial exploration point, respectively;
[0007] A potential coefficient is calculated according to the altitude values of the initial exploration point and the surrounding exploration points, a real usability value is calculated according to the altitude values and gravity values of the initial exploration point and the surrounding exploration points, and a petroleum potential rate is calculated based on the potential coefficient, the real usability value, the gravity values of the initial exploration point and the surrounding exploration points;
[0008] According to the altitude values of the two initial exploration points and the oil potential, a correction degree is calculated, a correction value is obtained based on the correction degree and the oil potential of the two initial exploration points corresponding to the correction degree, and an oil final potential is obtained based on the correction value and the distance between the two initial exploration points corresponding to the correction value;
[0009] According to the oil final potential and the coordinates of the initial exploration points, the coordinates of an optimal exploration point are obtained, the optimal exploration point is taken as the initial exploration point, and the process of obtaining the optimal exploration point is repeated to obtain the coordinates of a next optimal exploration point.
[0010] Further, the process of obtaining the potential coefficient includes:
[0011] The altitude value of the initial exploration point and the altitude values of the surrounding exploration points corresponding to the initial exploration point are sequentially added and then normalized to obtain a first altitude normalized value;
[0012] The potential coefficient of the initial exploration point is obtained by subtracting the first altitude normalized value from a preset value.
[0013] Further, the process of obtaining the real usability value includes:
[0014] The altitude value of the jth surrounding exploration point corresponding to the ith initial exploration point is subtracted from the altitude value of the ith initial exploration point, and then normalized to obtain a second altitude normalized value;
[0015] The gravity value of the jth surrounding exploration point corresponding to the ith initial exploration point is subtracted from the gravity value of the ith initial exploration point, and then normalized to obtain a gravity normalized value;
[0016] The second altitude normalized value and the gravity normalized value are subtracted, and then the value of the inverse number is taken as an input value of an exponential function with e as the base, and the output value of the exponential function is taken as the real usability value of the oil potential of the jth surrounding exploration point to the ith initial exploration point.
[0017] Further, the process of obtaining the oil potential includes:
[0018] The gravity value of the jth surrounding exploration point corresponding to the ith initial exploration point is subtracted from the gravity value of the ith initial exploration point, and then multiplied by the real usability value of the oil potential of the jth surrounding exploration point to the ith initial exploration point to obtain a real usable gravity value difference between the ith initial exploration point and the jth surrounding exploration point.
[0019] repeating the acquisition process of the real available gravity value difference component to obtain the real available gravity value difference component of the ith initial exploration point compared with each of the surrounding exploration points corresponding thereto;
[0020] adding the real available gravity value difference components in sequence to obtain a real available gravity value difference, and multiplying a normalized value of the real available gravity value difference by the potential coefficient of the ith initial exploration point to obtain the petroleum potential rate of the ith initial exploration point.
[0021] Further, the acquisition process of the correction degree comprises:
[0022] normalizing the square value of the difference between the altitude value of the ith initial exploration point and the altitude value of the kth initial exploration point to obtain a third altitude normalized value;
[0023] normalizing the absolute value of the difference between the petroleum potential rate of the ith initial exploration point and the petroleum potential rate of the kth initial exploration point to obtain a petroleum potential rate normalized value;
[0024] taking the absolute value of the difference between the third altitude normalized value and the petroleum potential rate normalized value as the correction degree of the petroleum potential rate of the kth initial exploration point to the ith initial exploration point.
[0025] Further, the acquisition process of the correction value comprises:
[0026] when the petroleum potential rate of the ith initial exploration point is less than the petroleum potential rate of the kth initial exploration point, multiplying the correction degree of the petroleum potential rate of the kth initial exploration point to the ith initial exploration point by the value between the petroleum potential rate of the kth initial exploration point minus the petroleum potential rate of the ith initial exploration point, and then adding the petroleum potential rate of the ith initial exploration point to obtain the correction value of the petroleum potential rate of the kth initial exploration point to the ith initial exploration point after correction.
[0027] Further, the acquisition process of the correction value further comprises:
[0028] when the petroleum potential rate of the ith initial exploration point is greater than or equal to the petroleum potential rate of the kth initial exploration point, taking the petroleum potential rate of the ith initial exploration point as the correction value of the petroleum potential rate of the kth initial exploration point to the ith initial exploration point after correction.
[0029] Further, the acquisition process of the final petroleum potential rate comprises:
[0030] an inverse of a distance between the ith initial exploration point and the kth initial exploration point as an input value of an exponential function with base e, an output value of the exponential function as a distance component of the kth initial exploration point relative to the ith initial exploration point;
[0031] repeating the process of obtaining the distance component to obtain the distance component of each initial exploration point relative to the ith initial exploration point, and sequentially adding each distance component to obtain a sum of distance components;
[0032] multiplying the distance component of the kth initial exploration point relative to the ith initial exploration point by the modified value of the kth initial exploration point to the ith initial exploration point after the distance component is divided by the sum of distance components to obtain a final petroleum potential component of the kth initial exploration point to the ith initial exploration point;
[0033] repeating the process of obtaining the final petroleum potential component to obtain the final petroleum potential component of each initial exploration point to the ith initial exploration point, and sequentially adding the final petroleum potential components to obtain the final petroleum potential of the ith initial exploration point.
[0034] Further, the process of obtaining the coordinates of the optimal exploration point comprises:
[0035] multiplying the final petroleum potential of the ith initial exploration point by the horizontal coordinate of the ith initial exploration point after the final petroleum potential of the ith initial exploration point is divided by a sum of the final petroleum potentials of each initial exploration point to obtain a horizontal coordinate component corresponding to the ith initial exploration point;
[0036] repeating the process of obtaining the horizontal coordinate component to obtain the horizontal coordinate component corresponding to each initial exploration point, and sequentially adding the horizontal coordinate components to obtain a horizontal coordinate of the optimal exploration point;
[0037] multiplying the final petroleum potential of the ith initial exploration point by the vertical coordinate of the ith initial exploration point after the final petroleum potential of the ith initial exploration point is divided by a sum of the final petroleum potentials of each initial exploration point to obtain a vertical coordinate component corresponding to the ith initial exploration point;
[0038] repeating the process of obtaining the vertical coordinate component to obtain the vertical coordinate component corresponding to each initial exploration point, and sequentially adding the vertical coordinate components to obtain a vertical coordinate of the optimal exploration point.
[0039] The embodiment of the present application also provides a data acquisition system applied to oil geological exploration, which comprises a memory, a processor and a computer program stored in the memory and running on the processor, and the processor implements the steps of the method when executing the computer program.
[0040] The present application has the following advantages:
[0041] Firstly, an initial exploration point is determined in a target area of oil geological exploration, and surrounding exploration points are determined at a preset distance from the initial exploration point. The initial exploration point and the surrounding exploration points are both basic exploration points, which are essential in oil geological exploration.
[0042] Secondly, a potential coefficient is calculated according to the altitude values of the initial exploration point and the surrounding exploration points, a real usability value is calculated according to the altitude values and gravity values of the initial exploration point and the surrounding exploration points, and an oil potential rate is calculated based on the potential coefficient, the real usability value, the gravity values of the initial exploration point and the surrounding exploration points. The oil potential rate is used to represent the probability of the existence of oil near an exploration point, and the greater the oil potential rate is, the greater the probability of the existence of oil near the exploration point is.
[0043] Then, a correction degree is calculated according to the altitude values of two optional initial exploration points and the oil potential rate, a correction value is obtained based on the correction degree and the oil potential rates of the two initial exploration points corresponding to the correction degree, and an oil final potential rate is obtained based on the correction value and the distance between the two initial exploration points corresponding to the correction value. The oil final potential rate is a final correction result of the oil potential rate, and represents the probability of the existence of oil near an exploration point more accurately.
[0044] Finally, the coordinates of an optimal exploration point are obtained according to the oil final potential rate and the coordinates of the initial exploration point, the optimal exploration point is taken as the initial exploration point, and the coordinates of a next optimal exploration point are obtained by repeating the obtaining process of the optimal exploration point. After obtaining the coordinates of the optimal exploration point, it is added to the initial exploration points, and then the coordinates of the next optimal exploration point are obtained, and the process is repeated. Through the targeted arrangement of exploration points and the focus on potential oil enrichment areas, the accuracy of data acquisition is improved and the exploration cost is effectively reduced. BRIEF DESCRIPTION OF DRAWINGS
[0045] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, and the advantages thereof, the following will briefly introduce the drawings needed in the embodiments or prior art description. Obviously, the drawings in the following description are only some embodiments of the present application, and for those skilled in the art, other drawings can be obtained based on these drawings without creative effort.
[0046] Figure 1 A flow chart of a data acquisition method for petroleum geological exploration provided by the first embodiment of the present application;
[0047] Figure 2 A flow chart of a potential coefficient acquisition process provided by the second embodiment of the present application;
[0048] Figure 3 A flow chart of a real usability value acquisition process provided by the third embodiment of the present application;
[0049] Figure 4 A flow chart of a petroleum potential rate acquisition process provided by the fourth embodiment of the present application;
[0050] Figure 5 A flow chart of a correction degree acquisition process provided by the fifth embodiment of the present application;
[0051] Figure 6 A flow chart of a petroleum final potential rate acquisition process provided by the sixth embodiment of the present application;
[0052] Figure 7 A flow chart of a coordinate of an optimal exploration point acquisition process provided by the seventh embodiment of the present application. DETAILED DESCRIPTION
[0053] In order to further illustrate the technical means and effects taken by the present application to achieve the predetermined purposes, the following will combine the drawings and the preferred embodiments to specifically describe the application, the specific embodiments, structures, features and effects of the data acquisition method and system for petroleum geological exploration according to the present application. In the following description, different "one embodiment" or "another embodiment" do not necessarily refer to the same embodiment. In addition, the specific features, structures or characteristics in one or more embodiments can be combined in any suitable form.
[0054] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which the present application belongs.
[0055] It should be noted that, in order to ensure that the calculation result is meaningful, when performing fractional operation, if the denominator is 0, a parameter adjustment factor greater than 0 is added to the denominator to prevent the denominator from being 0, and the value of the parameter adjustment factor is set by the implementer according to the actual situation, and the present application does not make special limitations.
[0056] The application provides a data acquisition method applied to petroleum geological exploration.
[0057] Please refer to Figure 1 , which shows a flowchart of a data acquisition method applied to petroleum geological exploration provided by the first embodiment of the application, and the method comprises the following steps.
[0058] S101. Determine an initial exploration point in a target area of petroleum geological exploration, and determine surrounding exploration points at a preset distance from the initial exploration point, respectively.
[0059] The initial exploration point can be set autonomously along the outer edge of the target area of petroleum geological exploration, and preferably, one initial exploration point is set every 5 km on the outer edge of the target area of petroleum geological exploration.
[0060] The preset distance can be set autonomously, and is preferably 500 m.
[0061] The number of the surrounding exploration points can be set autonomously, and is preferably 36, and the surrounding gravity exploration points can be set every 10 degrees of the central angle on the circumference with the initial exploration point as the center and the preset distance as the radius.
[0062] S102. Calculate a potential coefficient according to the altitude values of the initial exploration point and the surrounding exploration points, calculate a real availability value according to the altitude values and gravity values of the initial exploration point and the surrounding exploration points, and calculate a petroleum potential rate based on the potential coefficient, the real availability value, and the gravity values of the initial exploration point and the surrounding exploration points.
[0063] The acquisition process of the potential coefficient will be described in detail in the second embodiment, and will not be repeated here.
[0064] The acquisition process of the real availability value will be described in detail in the third embodiment, and will not be repeated here.
[0065] The acquisition process of the petroleum potential rate will be described in detail in the fourth embodiment, and will not be repeated here.
[0066] S103. Calculate a correction degree according to the elevation values of the two initial exploration points and the oil potential, obtain a correction value based on the correction degree and the oil potential of the two initial exploration points corresponding to the correction degree, and obtain a final oil potential based on the correction value, the distance between the two initial exploration points corresponding to the correction value.
[0067] The obtaining process of the correction degree will be described in detail in the fifth embodiment, and will not be described here.
[0068] The obtaining process of the correction value includes:
[0069] When the oil potential of the i-th initial exploration point is less than the oil potential of the k-th initial exploration point, the oil potential of the k-th initial exploration point minus the value between the oil potential of the i-th initial exploration point times the correction degree of the oil potential of the k-th initial exploration point to the i-th initial exploration point, and then added to the oil potential of the i-th initial exploration point to obtain the corrected correction value of the oil potential of the k-th initial exploration point to the i-th initial exploration point.
[0070] When the oil potential of the i-th initial exploration point is greater than or equal to the oil potential of the k-th initial exploration point, the oil potential of the i-th initial exploration point is taken as the corrected correction value of the oil potential of the k-th initial exploration point to the i-th initial exploration point.
[0071] Taking the corrected correction value of the oil potential of the k-th initial exploration point to the i-th initial exploration point as an example, if the oil potential of the k-th initial exploration point is smaller than the oil potential of the i-th initial exploration point, no correction is needed, if the oil potential of the k-th initial exploration point is larger than the oil potential of the i-th initial exploration point, the influence of the gravity monitoring of the i-th initial exploration point caused by sedimentation needs to be corrected, at this time the correction value of the i-th initial exploration point is not only the oil potential of the i-th initial exploration point W i , the greater the correction degree obtained above is, the closer it should be to the oil potential of the k-th initial exploration point W k , then the following formula can be constructed to represent the corrected correction value of the oil potential of the k-th initial exploration point to the i-th initial exploration point.
[0072] The correction value can be represented as:
[0073]
[0074] wherein, the P i-k represents the modified value of the petroleum potential of the kth initial exploration point to the ith initial exploration point, the W i represents the petroleum potential of the ith initial exploration point, the W k represents the petroleum potential of the kth initial exploration point, the H i-k represents the modification degree of the petroleum potential of the kth initial exploration point to the ith initial exploration point.
[0075] The petroleum final potential will be described in detail in the sixth embodiment, and will not be repeated here.
[0076] S104. Obtain the coordinate of an optimal exploration point according to the petroleum final potential and the coordinate of the initial exploration point, take the optimal exploration point as the initial exploration point, and repeat the obtaining process of the optimal exploration point to obtain the coordinate of a next optimal exploration point.
[0077] The coordinate obtaining process of the optimal exploration point will be described in detail in the seventh embodiment, and will not be repeated here.
[0078] The positions of all exploration points (initial exploration points and surrounding exploration points) can be used for gravity acquisition by using gravity measuring instruments (such as gravimeters), and the three-dimensional coordinates (x, y, z) thereof can be accurately acquired by using GPS equipment.
[0079] Low-lying terrain is often the primary place for sediment accumulation, where organic matter is rich in sediment, providing an abundant material basis for the formation of oil. Thanks to its deep sedimentary layer, these areas are more likely to be transformed into oil resources after a long geological change. In addition, low-lying areas are often associated with geological structure movements such as crustal subsidence and fault activity, which not only help the accumulation of oil, but also provide a closed storage space for oil. At the same time, the hydrological environment of the region also plays a role in promoting the accumulation of organic matter, further accelerating the oil generation process. In summary, low-lying terrain, with its unique environmental conditions, has become an area with great potential for oil resources.
[0080] Figure 2 The flow chart of the potential coefficient obtaining process provided by the second embodiment of the present application, the potential coefficient obtaining process comprising:
[0081] S201. Add the altitude value of the initial exploration point and the altitude values of each surrounding exploration point corresponding to the initial exploration point in turn, and then perform normalization processing to obtain a first altitude normalized value.
[0082] The first altitude normalization value can be expressed as:
[0083]
[0084] wherein the z i represents the altitude value of the i-th initial exploration point, the z ij represents the altitude value of the j-th surrounding exploration point corresponding to the i-th initial exploration point, the J represents the number of the surrounding exploration points, the f() represents a normalization function, preferably a range normalization function, the value range of the i is 1 to the number of the initial exploration points, and the value range of the j is 1 to the number of the surrounding exploration points.
[0085] S202. The potential coefficient of the initial exploration point is obtained by subtracting the first altitude normalization value from a preset value.
[0086] The potential coefficient can be expressed as:
[0087]
[0088] wherein the H i represents the potential coefficient of the i-th initial exploration point, and the R represents the preset value, preferably 1.
[0089] For the initial exploration point, when the sum of the altitude values of the corresponding surrounding exploration points is smaller, that is, the possibility of the existence of oil nearby is greater.
[0090] Since the density of oil is generally lower than that of most rocks, the average density of the underground around the oil reservoir is correspondingly lower, which directly leads to a relatively small gravity value (the gravity value is proportional to the density of the material) in this area. As the exploration proceeds deeper, the change in underground gravity anomaly gradually weakens, making the gravity signal of the deep oil and gas reservoirs increasingly weak. Therefore, near the top of the oil reservoir (i.e., at a lower altitude of the exploration point), the detected gravity value tends to be low; on the contrary, far from the top of the oil reservoir (i.e., at a higher altitude of the exploration point), the detected gravity value is relatively high. According to this principle, if the j-th surrounding exploration point corresponding to the i-th initial exploration point can effectively indicate the oil potential of the i-th initial exploration point, compared with the i-th initial exploration point, the higher the altitude of the j-th surrounding exploration point, the greater the gravity value measured by the j-th surrounding exploration point relative to the gravity value of the i-th initial exploration point.
[0091] Figure 3 The flowchart of the process for obtaining the real usability value provided by the third embodiment of the present application, the process for obtaining the real usability value comprises:
[0092] S301. Subtract the elevation value of the jth surrounding exploration point corresponding to the ith initial exploration point from the elevation value of the ith initial exploration point, and then normalize to obtain a second elevation normalized value.
[0093] The second elevation normalized value can be expressed as:
[0094] f(z ij -z i );
[0095] wherein z i represents the elevation value of the ith initial exploration point, z ij represents the elevation value of the jth surrounding exploration point corresponding to the ith initial exploration point, and f() represents a normalization function, preferably a range normalization function.
[0096] f(z ij -z i ) represents the height standardized difference of the jth surrounding exploration point relative to the ith initial exploration point.
[0097] S302. Subtract the gravity value of the jth surrounding exploration point corresponding to the ith initial exploration point from the gravity value of the ith initial exploration point, and then normalize to obtain a gravity normalized value.
[0098] The gravity normalized value can be expressed as:
[0099] f(G ij -G i );
[0100] wherein G i represents the gravity value of the ith initial exploration point, G ij represents the gravity value of the jth surrounding exploration point corresponding to the ith initial exploration point, and f() represents a normalization function, preferably a range normalization function.
[0101] f(G ij -G i ) represents the gravity value standardized difference of the jth surrounding exploration point relative to the ith initial exploration point.
[0102] S303. Subtract the second elevation normalized value from the gravity normalized value, and then take the value of the opposite number as the input value of the exponential function with base e, and the output value of the exponential function as the real availability value of the jth surrounding exploration point for the oil potential of the ith initial exploration point.
[0103] The real availability value can be expressed as:
[0104] S ij = exp[-|f(z ij -z i )-f(G ij -G i )|];
[0105] Wherein, the S ij represents the real availability value of the oil potential rate obtained by the jth surrounding exploration point to the ith initial exploration point, the exp represents an exponential function with e as the base, and the || represents an absolute value function.
[0106] Taking the ith initial exploration point as an example, generally, the larger the potential coefficient is, the higher the oil potential rate is, however, the potential coefficient is only effective for a larger area and cannot be accurate to a single point, therefore, on the basis of the potential coefficient of the ith initial exploration point, the oil potential rate is further determined in combination with the data of surrounding exploration points.
[0107] Considering that the oil reservoir usually exhibits a rock layer with a lower density, this characteristic can cause the gravity value of a local area to be lower than that of a surrounding area. Therefore, generally, the gravity value measured by the ith initial exploration point needs to be compared with the gravity value of the surrounding exploration point corresponding to the ith initial exploration point: the lower the gravity value of the ith initial exploration point is, the better, that is, the gravity value of the surrounding exploration point should be greater than that of the ith initial exploration point. It is known that the real availability of the surrounding exploration point is different, therefore, the oil potential rate of the ith initial exploration point can be obtained in a weighted manner.
[0108] Figure 4 A flowchart of an oil potential rate obtaining process provided by the fourth embodiment of the present application, the oil potential rate obtaining process comprising:
[0109] S401. The real available gravity value difference quantity of the ith initial exploration point compared with the jth surrounding exploration point is obtained by subtracting the gravity value of the ith initial exploration point from the gravity value of the jth surrounding exploration point corresponding to the ith initial exploration point and then multiplying the real availability value of the oil potential rate obtained by the jth surrounding exploration point to the ith initial exploration point.
[0110] The real available gravity value difference quantity can be expressed as:
[0111] S ij ×(G ij -G i );
[0112] Wherein, the S ijGrepresents the real availability value of the oil potential of the jth surrounding exploration point to the ith initial exploration point, the G i Grepresents the gravity value of the ith initial exploration point, the G ij Grepresents the gravity value of the jth surrounding exploration point corresponding to the ith initial exploration point.
[0113] S402. Repeat the process of obtaining the real available gravity value difference to obtain the real available gravity value difference of the ith initial exploration point compared with each surrounding exploration point corresponding thereto.
[0114] S403. Add the real available gravity value differences in sequence to obtain the real available gravity value difference, and multiply the normalized value of the real available gravity value difference by the potential coefficient of the ith initial exploration point to obtain the oil potential of the ith initial exploration point.
[0115] The oil potential can be represented as:
[0116]
[0117] wherein the H i Grepresents the potential coefficient of the ith initial exploration point, the J represents the number of surrounding exploration points, the f() represents a normalization function, preferably a range normalization function, the J represents the number of surrounding exploration points, and the W i Grepresents the oil potential of the ith initial exploration point.
[0118] According to the formation and evolution law of rock strata, sedimentation and erosion alternately occur, resulting in differences in rock strata thickness and ground height in different regions. Sedimentation can form thicker sedimentary layers in some regions, while erosion can weaken or remove part of the rock strata. Such differences in geological evolution make the underground oil reservoirs in a certain depth range relatively uniform, but the height of the ground rock strata varies due to different sedimentation and erosion in different regions. Therefore, for different initial exploration points, even if the depths of their oil layers are the same, the differences in the height of the ground rock strata will still result in changes in altitude. According to Newton's law of universal gravitation, the change in the gravitational field is inversely proportional to the square of the distance of an object to the observation point. Therefore, if two initial gravity exploration points are in the same area where the oil layer is located, the gravity difference should be inversely proportional to the square of the distance (after normalization, the greater the difference, the more conducive to the detection of gravity difference). Therefore, considering the effect of sedimentation on the ground rock strata, it is necessary to correct the case where the gravity monitoring is not significant due to sedimentation.
[0119] Figure 5A flow chart of the process for obtaining the correction degree provided by the fifth embodiment of the present application, the process for obtaining the correction degree comprises:
[0120] S501. Normalizing the square value of the difference between the altitude value of the i-th initial exploration point and the altitude value of the k-th initial exploration point to obtain a third altitude normalized value.
[0121] The third altitude normalized value can be expressed as:
[0122] f[(z i -z k ) 2 ];
[0123] Wherein, the z i represents the altitude value of the i-th initial exploration point, the z k represents the altitude value of the k-th initial exploration point, and the f() represents a normalization function, preferably a range normalization function.
[0124] The value range of the i is 1 to the number of the initial exploration points, the value range of the k is 1 to the number of the initial exploration points, and i is not equal to k.
[0125] S502. Normalizing the absolute value of the difference between the oil potential of the i-th initial exploration point and the oil potential of the k-th initial exploration point to obtain an oil potential normalized value.
[0126] The oil potential normalized value can be expressed as:
[0127] f[|W i -W k |);
[0128] Wherein, the W i represents the oil potential of the i-th initial exploration point, the W k represents the oil potential of the k-th initial exploration point, the f() represents a normalization function, preferably a range normalization function, and the || represents an absolute value function.
[0129] S503. Taking the absolute value of the difference between the third altitude normalized value and the oil potential normalized value as the correction degree of the oil potential of the k-th initial exploration point to the i-th initial exploration point.
[0130] The correction degree of the oil potential of the k-th initial exploration point to the i-th initial exploration point can be expressed as:
[0131] H i-k|f[(z i -z k ) 2 ]-f(|W i -W k |)|;
[0132] wherein, the H i-k represents the correction degree of the oil potential of the kth initial exploration point to the ith initial exploration point.
[0133] Since the gravity anomaly is usually local, and its influence decreases with the increase of distance, the initial key exploration points close to the ith initial exploration point are more affected by the gravity anomaly and are more sensitive to monitoring, and the correction value of the initial key exploration points has a more significant influence on the ith initial exploration point; the initial key exploration points far away from the ith initial exploration point are less affected by the anomaly, and thus the correction value of the initial key exploration points has a weaker influence on the ith initial exploration point, so the distance is combined for weighting to construct the acquisition formula of the final oil potential.
[0134] Figure 6 The flowchart of the acquisition process of the final oil potential provided by the sixth embodiment of the present application, the acquisition process of the final oil potential comprises:
[0135] S601. The reciprocal of the distance between the ith initial exploration point and the kth initial exploration point is taken as the input value of the exponential function with e as the base, and the output value of the exponential function is taken as the distance component of the kth initial exploration point relative to the ith initial exploration point.
[0136] The distance component can be represented as:
[0137] exp[-d(i,k)];
[0138] wherein, the d(i,k) represents the distance between the ith initial exploration point and the kth initial exploration point, and the exp represents the exponential function with e as the base.
[0139] S602. The acquisition process of the distance component is repeated to acquire the distance component of each initial exploration point relative to the ith initial exploration point, and each distance component is sequentially added to obtain the sum of distance components.
[0140] The sum of distance components can be represented as:
[0141]
[0142] wherein, the N represents the number of initial exploration points.
[0143] S603. Multiplying the distance component of the kth initial exploration point relative to the ith initial exploration point by the modified value of the oil potential correction of the kth initial exploration point to the ith initial exploration point after dividing the distance component by the sum of the distance components to obtain the oil final potential component of the kth initial exploration point to the ith initial exploration point.
[0144] The oil final potential component can be represented as:
[0145]
[0146] Wherein, the P i-k represents the modified value of the oil potential correction of the kth initial exploration point to the ith initial exploration point.
[0147] S604. Repeating the process of obtaining the oil final potential component to obtain the oil final potential component of each initial exploration point to the ith initial exploration point, and sequentially adding the oil final potential components to obtain the oil final potential of the ith initial exploration point.
[0148] The oil final potential can be represented as:
[0149]
[0150] Wherein, the P i represents the oil final potential of the ith initial exploration point.
[0151] By each initial exploration point, the optimal exploration point for the next oil geological exploration is determined. The above can determine the oil final potential of each initial exploration point, which reflects the richness of the oil resources that may exist near each initial exploration point. Therefore, in order to more effectively discover oil in the next exploration process, the next oil exploration point needs to be selected near the initial exploration point with a higher oil final potential, so when calculating the position of the next optimal exploration point, the position of the initial exploration point with a higher oil final potential can be given a higher weight, and finally a new exploration point position is calculated by weighted summation.
[0152] Using this method, not only can the efficiency of oil exploration be improved, but also unnecessary resource waste can be effectively reduced. With the advancement of exploration, the data collected by the initial exploration points is used to continuously optimize the selection of exploration points, forming a dynamically adjusted exploration strategy. After each exploration, the strategy can be adjusted according to new exploration data to ensure that the exploration point continuously approaches the area richest in oil resources, further improving the efficiency of exploration.
[0153] Figure 7A flow chart of a process for obtaining coordinates of an optimal exploration point is provided in the seventh embodiment of the present application. The process for obtaining coordinates of the optimal exploration point comprises:
[0154] S701. The petroleum final potential of the ith initial exploration point is divided by the sum of the petroleum final potentials of all the initial exploration points, and then multiplied by the horizontal coordinate of the ith initial exploration point to obtain a horizontal coordinate component corresponding to the ith initial exploration point.
[0155] The horizontal coordinate component corresponding to the initial exploration point can be expressed as:
[0156]
[0157] wherein the P i represents the petroleum final potential of the ith initial exploration point, the M represents the number of the initial exploration points, the x i represents the horizontal coordinate of the ith initial exploration point, the P b represents the petroleum final potential of the bth initial exploration point.
[0158] S702. The process for obtaining the horizontal coordinate component is repeated to obtain the horizontal coordinate components corresponding to all the initial exploration points, and the horizontal coordinate components are sequentially added to obtain the horizontal coordinate of the optimal exploration point.
[0159] The horizontal coordinate of the optimal exploration point can be expressed as:
[0160]
[0161] S703. The petroleum final potential of the ith initial exploration point is divided by the sum of the petroleum final potentials of all the initial exploration points, and then multiplied by the vertical coordinate of the ith initial exploration point to obtain a vertical coordinate component corresponding to the ith initial exploration point.
[0162] The vertical coordinate component corresponding to the initial exploration point can be expressed as:
[0163]
[0164] wherein the P i represents the petroleum final potential of the ith initial exploration point, the M represents the number of the initial exploration points, the y i represents the vertical coordinate of the ith initial exploration point.
[0165] S704. Repeat the process of obtaining the ordinate component to obtain the ordinate component corresponding to each of the initial exploration points, and sequentially add the ordinate components to obtain the ordinate of the optimal exploration point.
[0166] The ordinate of the optimal exploration point can be expressed as:
[0167]
[0168] The embodiment of the present application also provides a data acquisition system applied to oil geological exploration, which comprises a memory, a processor and a computer program stored in the memory and running on the processor, and the processor implements the steps of the above method when executing the computer program.
[0169] The technical features and technical effects of the data acquisition system applied to oil geological exploration provided by the embodiment of the present application are the same as those of the method provided by the embodiment of the present application, and will not be repeated here.
[0170] The present application has the following beneficial effects:
[0171] First, an initial exploration point of a first number is determined in a target area of oil geological exploration, and a surrounding exploration point of a second number is determined at a preset distance from the initial exploration point. The initial exploration point and the surrounding exploration point are both basic exploration points, which are essential exploration points in oil geological exploration.
[0172] Second, a potential coefficient is calculated according to the altitude values of the initial exploration point and the surrounding exploration point, a real usability value is calculated according to the altitude values and gravity values of the initial exploration point and the surrounding exploration point, and an oil potential rate is calculated based on the potential coefficient, the real usability value, the gravity values of the initial exploration point and the surrounding exploration point. The oil potential rate is used to represent the probability of the existence of oil near the exploration point, and the greater the oil potential rate, the greater the probability of the existence of oil near the exploration point.
[0173] Then, a correction degree is calculated according to the altitude values of the two optional initial exploration points and the oil potential rate, a correction value is obtained based on the correction degree and the oil potential rates of the two initial exploration points corresponding to the correction degree, and an oil final potential rate is obtained based on the correction value and the distance between the two initial exploration points corresponding to the correction value. The oil final potential rate is the final correction result of the oil potential rate, and represents the probability of the existence of oil near the exploration point more accurately.
[0174] Finally, the coordinates of the optimal exploration point are obtained according to the ultimate potential rate of the oil and the coordinates of the initial exploration point, the optimal exploration point is taken as the initial exploration point, the process of obtaining the coordinates of the optimal exploration point is repeated to obtain the coordinates of the next optimal exploration point. After obtaining the coordinates of the optimal exploration point, it is added to the initial exploration point, and then the coordinates of the next optimal exploration point are obtained, and the process is repeated, by laying out exploration points in a targeted manner, focusing on the potential oil enrichment area, the accuracy of data collection is significantly improved, and the exploration cost is effectively reduced.
[0175] It should be noted that the above-mentioned embodiment sequence of the present application is only for description, and does not represent the advantages and disadvantages of the embodiments. The processes depicted in the drawings do not necessarily require the specific order or continuous order shown to achieve the desired results. In some embodiments, multi-task processing and parallel processing are also possible or may be advantageous.
[0176] Each embodiment in the specification is described in a progressive manner, and the same or similar parts between each embodiment can be referred to each other, and each embodiment mainly explains the difference from other embodiments.
Claims
1. A data acquisition method applied to petroleum geological exploration, characterized in that, The method comprises: determining an initial exploration point in a petroleum geological exploration target area, and determining surrounding exploration points at a preset distance from the initial exploration point respectively; calculating a potential coefficient according to the elevation values of the initial exploration point and the surrounding exploration points, calculating a real usability value according to the elevation values and gravity values of the initial exploration point and the surrounding exploration points, and calculating a petroleum potential rate based on the potential coefficient, the real usability value, the gravity values of the initial exploration point and the surrounding exploration points; calculating a correction degree according to the elevation values of the two initial exploration points and the petroleum potential rate, obtaining a correction value based on the correction degree and the petroleum potential rates of the two initial exploration points corresponding to the correction degree, and obtaining a final petroleum potential rate based on the correction value and the distance between the two initial exploration points corresponding to the correction value; obtaining the coordinates of an optimal exploration point according to the final petroleum potential rate and the coordinates of the initial exploration point, taking the optimal exploration point as the initial exploration point, and repeating the obtaining process of the optimal exploration point to obtain the coordinates of a next optimal exploration point.
2. The data acquisition method for petroleum geological exploration according to claim 1, characterized in that, The obtaining process of the potential coefficient comprises: adding the elevation value of the initial exploration point and the elevation values of the surrounding exploration points corresponding to the initial exploration point in sequence, and then performing normalization processing to obtain a first elevation normalization value; subtracting the first elevation normalization value from a preset value to obtain the potential coefficient of the initial exploration point.
3. The data acquisition method for petroleum geological exploration according to claim 1, characterized in that, The obtaining process of the real usability value comprises: subtracting the elevation value of the i-th initial exploration point from the elevation value of the j-th surrounding exploration point corresponding to the i-th initial exploration point, and then performing normalization processing to obtain a second elevation normalization value; subtracting the gravity value of the i-th initial exploration point from the gravity value of the j-th surrounding exploration point corresponding to the i-th initial exploration point, and then performing normalization processing to obtain a gravity normalization value; subtracting the second elevation normalization value from the gravity normalization value, and then taking the inverse value as an input value of an exponential function with e as the base, and taking the output value of the exponential function as the real usability value of the j-th surrounding exploration point to the i-th initial exploration point.
4. The data acquisition method for petroleum geological exploration according to claim 1, characterized in that, The obtaining process of the petroleum potential rate comprises: subtracting the gravity value of the i-th initial exploration point from the gravity value of the j-th surrounding exploration point corresponding to the i-th initial exploration point, and then multiplying the real usability value of the j-th surrounding exploration point to the i-th initial exploration point to obtain a real usable gravity value difference of the i-th initial exploration point compared with the j-th surrounding exploration point; repeating the obtaining process of the real usable gravity value difference to obtain the real usable gravity value differences of the i-th initial exploration point compared with the surrounding exploration points corresponding to the i-th initial exploration point; and The real available gravity value difference components are sequentially added to obtain a real available gravity value difference, and a normalized value of the real available gravity value difference is multiplied by the potential coefficient of the i th initial exploration point to obtain the oil potential rate of the i th initial exploration point.
5. The data acquisition method for petroleum geological exploration according to claim 1, characterized in that, The correction degree acquisition process comprises: a square value obtained by subtracting the altitude value of the k th initial exploration point from the altitude value of the i th initial exploration point is normalized to obtain a third altitude normalized value; an absolute value obtained by subtracting the oil potential rate of the k th initial exploration point from the oil potential rate of the i th initial exploration point is normalized to obtain an oil potential rate normalized value; an absolute value of a difference between the third altitude normalized value and the oil potential rate normalized value is taken as the correction degree of the oil potential rate of the i th initial exploration point with respect to the i th initial exploration point.
6. The data acquisition method for petroleum geological exploration according to claim 1, characterized in that, The correction value acquisition process comprises: when the oil potential rate of the i th initial exploration point is less than the oil potential rate of the k th initial exploration point, the oil potential rate of the k th initial exploration point is subtracted by the oil potential rate of the i th initial exploration point, multiplied by the correction degree of the oil potential rate of the i th initial exploration point with respect to the i th initial exploration point, and then added to the oil potential rate of the i th initial exploration point to obtain the corrected correction value of the oil potential rate of the k th initial exploration point with respect to the i th initial exploration point.
7. The data acquisition method for petroleum geological exploration according to claim 1, characterized in that, The correction value acquisition process further comprises: when the oil potential rate of the i th initial exploration point is greater than or equal to the oil potential rate of the k th initial exploration point, the oil potential rate of the i th initial exploration point is taken as the corrected correction value of the oil potential rate of the k th initial exploration point with respect to the i th initial exploration point.
8. The data acquisition method for petroleum geological exploration according to claim 1, characterized in that, The final oil potential rate acquisition process comprises: an inverse of a distance between the i th initial exploration point and the k th initial exploration point is taken as an input value of an exponential function with e as a base, and an output value of the exponential function is taken as a distance component of the k th initial exploration point with respect to the i th initial exploration point; the distance component acquisition process is repeated to obtain the distance component of each initial exploration point with respect to the i th initial exploration point, and the distance components are sequentially added to obtain a distance component sum; the distance component of the k th initial exploration point with respect to the i th initial exploration point is divided by the distance component sum, and then multiplied by the corrected correction value of the oil potential rate of the k th initial exploration point with respect to the i th initial exploration point to obtain an oil final potential rate component of the k th initial exploration point with respect to the i th initial exploration point. The process of obtaining the petroleum ultimate potential component of each initial exploration point is repeated to obtain the petroleum ultimate potential component of the ith initial exploration point, and the petroleum ultimate potential components are sequentially added to obtain the petroleum ultimate potential of the ith initial exploration point.
9. The data acquisition method for petroleum geological exploration according to claim 1, characterized in that, The process of obtaining the coordinates of the optimal exploration point comprises: The petroleum ultimate potential of the ith initial exploration point is divided by the sum of the petroleum ultimate potentials of all the initial exploration points, and then multiplied by the abscissa of the ith initial exploration point to obtain the abscissa component corresponding to the ith initial exploration point; The process of obtaining the abscissa component is repeated to obtain the abscissa components corresponding to all the initial exploration points, and the abscissa components are sequentially added to obtain the abscissa of the optimal exploration point; The petroleum ultimate potential of the ith initial exploration point is divided by the sum of the petroleum ultimate potentials of all the initial exploration points, and then multiplied by the ordinate of the ith initial exploration point to obtain the ordinate component corresponding to the ith initial exploration point; The process of obtaining the ordinate component is repeated to obtain the ordinate components corresponding to all the initial exploration points, and the ordinate components are sequentially added to obtain the ordinate of the optimal exploration point.
10. A data acquisition system for use in petroleum geology exploration, the system comprising a memory, a processor, and a computer program stored in the memory and running on the processor, characterized in that, The processor implements the steps of the method of any one of claims 1-9 when executing the computer program.
Citation Information
Patent Citations
High-precision time shifting microgravity acquiring method
CN104570140A
Method and device for acquiring oil-reservoir geologic reserves with gravitational differentiation phenomenon
CN107795319A