The present invention relates to the technical field of zone evaluation method in petroleum exploration, in particular relates to an oil and gas zone effectiveness evaluation method and apparatus.
The oil and gas zone effectiveness evaluation refers to the evaluation of whether or not an industrial oil and gas reservoir exists in a zone, and the oil and gas zone effectiveness evaluation is the key to determining the success rate of drilling, the oil and gas discovery and the growth rate of the reserves.
The oil and gas zone effectiveness evaluation methods in the prior arts include: R factor principle component analysis method, multi-factor weighted evaluation method, multi-map superimposition evaluation method, semi-quantitative evaluation method, and evaluation methods based on resource quantity and economical efficiency.
The R factor principle component analysis method cannot evaluate a single zone. Since the conditions of the zones vary greatly, the evaluation result has a low coincidence rate. The weight coefficient of the multi-factor weighted evaluation method is determined artificially, thus the indeterminacy is large, and the determinative effect of controlling parameters is not taken into consideration. The multi-map superimposition evaluation method cannot realize a quantitative evaluation, thus for lower exploration areas, it is unlikely to accurately obtain a single factor map of the evaluated region, and accordingly, the evaluation coincidence rate is low. For the semi-quantitative evaluation method, the controlling factor parameter values of different oil and gas zones differ greatly from each other, the applicable range of the evaluation parameter standard is limited, and accurate parameter values cannot be obtained for lower exploration areas, and thus the evaluation coincidence rate is low. For the evaluation methods based on resource quantity and economical efficiency, the resource quantity cannot determine whether an industrial oil and gas reservoir can be obtained in a zone, and it is hard to evaluate the economical efficiency of the zone before drilling, and thus the evaluation coincidence rate is low.
The current oil and gas zone evaluation technology mainly includes the geological condition analogy, mathematic model calculation, resource quantity evaluation and multi-factor superposition, etc., which are influenced greatly by human experience and level of data acquisition, and cannot establish a unified standard of zone effective evaluation, i.e., cannot achieve a unified evaluation of zone effectiveness, and thus the zone evaluation coincidence rate is low.
The purpose of the present invention is to provide an oil and gas zone effectiveness evaluation method and apparatus which can improve the evaluation coincidence rate.
The above purpose of the present invention can be realized by adopting the following technical solutions: an oil and gas zone effectiveness evaluation method, comprising: obtaining a controlling factor parameter distribution of a target interval, wherein the controlling factor parameter includes: stratum dip angle, fault, unconformity within an effective hydrocarbon generation range, reservoir thickness, cap rock thickness, oil generation abundance and gas generation abundance; obtaining an evaluation parameter lower limit value of an industrial oil and gas reservoir, wherein the evaluation parameter includes: reservoir thickness, porosity, cap rock thickness, oil generation abundance and gas generation abundance; calculating an accumulation parameter value of each grid coordinate point according to the controlling factor parameter distribution; calculating an evaluation parameter processing value of each grid coordinate point according to the evaluation parameter, the evaluation parameter lower limit value and the accumulation parameter value of the grid coordinate point; and calculating a zone evaluation value of the target interval according to the evaluation parameter processing value of the grid coordinate point.
As a preferred embodiment, said calculating an accumulation parameter value of each grid coordinate point according to the controlling factor parameter distribution specifically includes: obtaining a dimensionless flag value of each grid coordinate point according to the controlling factor parameter; obtaining the accumulation parameter value of each grid coordinate point according to the dimensionless flag value of each grid coordinate point.
As a preferred embodiment, the accumulation parameter value of each grid coordinate point is obtained according to the following formula:
wherein, AR—the accumulation parameter value, dimensionless, when AR≥1, it represents an effective accumulation region, and when AR<1, it represents a non-effective accumulation region; Parai—dimensionless parameter values Sths, Tss, Das, Rds and Fs; n—the number of Parai; Sths—a cap rock thickness flag value, dimensionless; Tss—a conduction system flag value, dimensionless; Das—a stratum dip angle flag value, dimensionless; Rds—an oil and gas migration distance flag value, dimensionless; Fs—a fault flag value.
As a preferred embodiment, the cap rock thickness flag value is determined according to the following formula:
wherein, Sths—the cap rock thickness flag value, dimensionless; Sth—the cap rock thickness of the target interval, m; Sthlim t—a cap rock thickness lower limit value of the industrial oil and gas reservoir, m.
As a preferred embodiment, the conduction system flag value is determined according to the following formula:
wherein, Tss—the conduction system flag value, dimensionless; H—the reservoir thickness of the target interval, m; Hlim t—a reservoir thickness lower limit value of the industrial oil and gas reservoir, m; UF—unconformity of direct contact with the source rock within the effective hydrocarbon generation range and/or a fault; when there is unconformity of direct contact with the source rock within the effective hydrocarbon generation range and/or a fault, the UF is 1, and when there is no unconformity of direct contact with the source rock within the effective hydrocarbon generation range and/or fault, the UF is 0.
As a preferred embodiment, the effective hydrocarbon generation range of the source rock is determined according to the following steps: obtaining an oil generation abundance lower limit value and a gas generation abundance lower limit value of the industrial oil and gas reservoir; determining the effective hydrocarbon generation range of the source rock according to the oil generation abundance lower limit value and the gas generation abundance lower limit value.
As a preferred embodiment, the stratum dip angle flag value is determined according to the following formula:
Das=sin(α),
wherein, Das—the stratum dip angle flag value, dimensionless; α—the stratum dip angle, degree.
As a preferred embodiment, the oil and gas migration distance flag value is determined according to the following formula:
wherein, Rds the oil and gas migration distance flag value, dimensionless; Rd—an equivalent radius of the effective hydrocarbon generation range, km; Lhg—a distance to a boundary of the effective hydrocarbon generation range from outside the effective hydrocarbon generation range, km, which is 0 when within the effective hydrocarbon generation range; a—an empirical coefficient, which takes 3 when the target interval contains conventional oil and gas, takes 1.2 when the target interval contains dense oil and gas, and takes 1 when the target interval contains shale oil and gas.
As a preferred embodiment, when the fault that cuts through the cap rock of the target interval is an open fault, the fault flag value is −1; and when the fault that cuts through the cap rock of the target interval is a sealed fault, the fault flag value is 1.
As a preferred embodiment, the evaluation parameter processing value of each grid coordinate point is calculated according to the following formula:
wherein, Pi—the evaluation parameter processing value of the ith grid coordinate point; Si—the ith evaluation parameter value; Si_lim t—the ith evaluation parameter lower limit value; n—the number of the evaluation parameter.
As a preferred embodiment, the zone evaluation value of the target interval is calculated according to the following formula:
wherein, Vplay— the zone evaluation value, −1˜1, when Vplay≥0, it represents an effective zone distribution area, and when Vplay≤0, it represents a non-effective zone distribution area; Pi—the evaluation parameter processing value of the ith grid coordinate point; n—the number of the grid coordinate point.
As a preferred embodiment, the f(Pi) is determined according to the following formula:
wherein, Pi— the evaluation parameter processing value of the ith grid coordinate point.
As a preferred embodiment, obtaining a stratum dip angle distribution of the target interval specifically includes: obtaining structure of the target interval according to geophysical data of the target interval; obtaining the stratum dip angle distribution of the target interval according to the structure of the target interval.
As a preferred embodiment, the stratum dip angle of the target interval is determined according to the following steps: obtaining a vertical line between an effective hydrocarbon generation center of the source rock and a top surface structural contour of the target interval; obtaining an inclined angle of the vertical line; obtaining an inclined direction of the inclined angle according to a positional relationship between the vertical line and a horizontal line; obtaining the stratum dip angle according to the inclined angle and the inclined direction.
As a preferred embodiment, the fault includes: a first fault that cuts through the cap rock of the target interval, a second fault that communicates the source rock and the reservoir within the effective hydrocarbon generation range of the source rock of the target interval, and a third fault that communicates the unconformity of the target interval and the reservoir.
As a preferred embodiment, obtaining the oil generation abundance distribution and the gas generation abundance distribution of the target interval specifically includes: obtaining an organic carbon content distribution of the target interval according to well logging data of the target interval and the organic carbon content obtained from a core analysis of the source rock; obtaining a vitrinite reflectance distribution according to the vitrinite reflectance obtained from the core analysis of the source rock of the target interval; performing a basin simulation operation to the organic carbon content distribution and the vitrinite reflectance distribution to obtain the oil generation abundance distribution and the gas generation abundance distribution of the target interval.
As a preferred embodiment, obtaining the reservoir thickness lower limit value of the industrial oil and gas reservoir specifically includes: obtaining the reservoir thickness lower limit value according to the reservoir thickness and oil and gas testing data of the industrial oil and gas reservoir, wherein, a reservoir type of the industrial oil and gas reservoir is the same as that of the target interval.
As a preferred embodiment, obtaining a porosity lower limit value of the industrial oil and gas reservoir specifically includes: obtaining the porosity lower limit value according to porosity analysis data and gas testing data of the industrial oil and gas reservoir.
As a preferred embodiment, obtaining a cap rock thickness lower limit value of the industrial oil and gas reservoir specifically includes: obtaining the cap rock thickness lower limit value according to the cap rock thickness of the industrial oil and gas reservoir, wherein lithology of the cap rock of the industrial oil and gas reservoir is consistent with that of the target interval.
An oil and gas zone effectiveness evaluation apparatus, comprising: a controlling factor parameter obtaining module for obtaining a controlling factor parameter distribution of a target interval, wherein the controlling factor parameter includes: stratum dip angle, fault, unconformity within an effective hydrocarbon generation range, reservoir thickness, cap rock thickness, oil generation abundance and gas generation abundance; an evaluation parameter lower limit value obtaining module for obtaining an evaluation parameter lower limit value of an industrial oil and gas reservoir, wherein the evaluation parameter includes: reservoir thickness, porosity, cap rock thickness, oil generation abundance and gas generation abundance; an accumulation parameter value calculating module for calculating an accumulation parameter value of each grid coordinate point according to the controlling factor parameter distribution; an evaluation parameter processing value calculating module for calculating an evaluation parameter processing value of each grid coordinate point according to the evaluation parameter, the evaluation parameter lower limit value and the accumulation parameter value of the grid coordinate point; and a zone evaluation value calculating module for calculating a zone evaluation value of the target interval according to the evaluation parameter processing value of the grid coordinate point.
The oil and gas zone effectiveness evaluation method provided by the present application has the following advantageous effect: the oil and gas zone effectiveness evaluation method according to the embodiments of the present application calculates the evaluation parameter processing value of each grid coordinate point according to the evaluation parameter, the evaluation parameter lower limit value and the accumulation parameter value of the grid coordinate point, thereby determining the accumulation range of the target interval, and then calculates the zone evaluation value of the target interval according to the evaluation parameter processing value of the grid coordinate point. In this way, the oil and gas zone effectiveness evaluation is conducted quantitatively. Therefore, compared to the prior arts, the oil and gas zone effectiveness evaluation method of the present application can have a unified standard of oil and gas zone effectiveness evaluation. Accordingly, the influence of artificial factors can be eliminated. Besides, since the evaluation parameter lower limit value of the oil and gas reservoir which has got industrial value is employed, the identification rate of the oil and gas zone effectiveness is improved, and thereby the coincidence rate of the oil and gas zone effectiveness evaluation is improved. As such, the success rate of drilling and the speed of finding an effective zone can be improved, so that the costs for oil and gas exploration can be reduced.
In order to explain the technical solutions in the embodiments in the present invention more clearly, the following will introduce the figures needed in the description of the embodiments briefly. Obviously, the figures in the following description are only some embodiments of the present invention. For ordinary persons skilled in the art, other figures may also be obtained based on these figures without paying creative efforts.
The following will make a clear and comprehensive description to the technical solutions in the embodiments of the present invention in combination with the figures in the embodiments of the present invention. Obviously, the embodiments described herein are only a part of, rather than all of the embodiments of the present application. Based on the embodiments of the present invention, all other embodiments obtained by the ordinary skilled persons in the field without paying creative efforts should pertain to the scope of protection of the present invention.
Referring to
In can be seen from the above technical solution that: the oil and gas zone effectiveness evaluation method according to the embodiment of the present application obtains the evaluation parameter processing value of each grid coordinate point according to the evaluation parameter, the evaluation parameter lower limit value and the accumulation parameter value of each grid coordinate point, thereby determining the accumulation range of the target interval, and then obtains the zone evaluation value of the target interval according to the evaluation parameter processing value of the grid coordinate point. In this way the oil and gas zone effectiveness evaluation is conducted quantitatively. Therefore, compared to the prior arts, the oil and gas zone effectiveness evaluation method of the present application can have a unified standard of oil and gas zone effectiveness evaluation. Accordingly, the influence of artificial factors can be eliminated. Besides, since the evaluation parameter lower limit value of the oil and gas reservoir which gas got industrial value is employed, the identification rate of the oil and gas zone effectiveness is improved, and thereby the coincidence rate of the oil-gas effectiveness evaluation is improved. As such, the success rate of drilling and the speed of finding an effective zone can be improved, so that the costs for oil-gas exploration can be reduced.
As shown in
As shown in
S101: obtaining structure of the target interval according to geophysical data of the target interval;
S103: obtaining the stratum dip angle distribution of the target interval according to the structure of the target interval.
As shown in
S1010: obtaining a vertical line between an effective hydrocarbon generation center of the source rock and a top surface structural contour of the target interval;
S1030: obtaining an inclined angle of the vertical line;
S1050: obtaining an inclined direction of the inclined angle according to a positional relationship between the vertical line and a horizontal line;
to be specific, when the vertical line between the effective hydrocarbon generation center of the source rock and the top surface structural contour of the target interval is above the horizontal line of the effective hydrocarbon generation center of the source rock, the stratum dip angle is positive; when the vertical line between the effective hydrocarbon generation center of the source rock and the top surface structural contour of the target interval is below the horizontal line of the effective hydrocarbon generation center of the source rock, the stratum dip angle is negative; and when the vertical line between the effective hydrocarbon generation center of the source rock and the top surface structural contour of the target interval coincides with the horizontal line of the effective hydrocarbon generation center of the source rock, the stratum dip angle is 0;
S1070: obtaining the stratum dip angle according to the inclined angle and the inclined direction.
In a specific embodiment, the fault includes: a first fault that cuts through the cap rock of the target interval, a second fault that communicates the source rock and the reservoir within the effective hydrocarbon generation range of the source rock of the target interval, and a third fault that communicates the unconformity of the target interval and the reservoir.
In a specific embodiment, the unconformity within the effective hydrocarbon generation range can be an unconformity between the source rock and the cap rock.
In a specific embodiment, obtaining a reservoir thickness distribution of the target interval can be realized by collecting geophysical data of the target interval.
In a specific embodiment, obtaining a porosity distribution of the target interval can be realized by collecting porosity analysis data and geophysical data of the target interval.
In a specific embodiment, obtaining a cap rock thickness distribution of the target interval can be realized by collecting geophysical data of the target interval.
As shown in
S105: collecting well logging data of the target interval and an organic carbon content obtained from core analysis of the source rock to obtain the organic carbon content distribution of the target interval;
S107: collecting a vitrinite reflectance obtained from core analysis of the source rock of the target interval to obtain the vitrinite reflectance distribution;
S109: performing a basin simulation operation to the organic carbon content distribution and the vitrinite reflectance distribution to obtain the oil generation abundance distribution and the gas generation abundance distribution.
As shown in
In a specific embodiment, obtaining a reservoir thickness lower limit value of the industrial oil and gas reservoir specifically includes:
obtaining the reservoir thickness lower limit value according to the reservoir thickness and oil and gas testing data of the industrial oil and gas reservoir, wherein, a reservoir type of the industrial oil and gas reservoir is the same as that of the target interval.
In a specific embodiment, obtaining a porosity lower limit value of the industrial oil and gas reservoir specifically includes: obtaining the porosity lower limit value according to porosity analysis data and gas testing data of the industrial oil and gas reservoir.
In a specific embodiment, obtaining a cap rock thickness lower limit value of the industrial oil and gas reservoir specifically includes: obtaining the cap rock thickness lower limit value according to the cap rock thickness of the industrial oil and gas reservoir, wherein lithology of the cap rock of the industrial oil and gas reservoir is consistent with that of the target interval.
In a specific embodiment, the oil generation abundance lower limit value and the gas generation abundance lower limit value of the industrial oil and gas reservoir are used to determine the effective hydrocarbon generation range of the source rock.
The oil and gas zone effectiveness evaluation method according to the present application is employed to evaluate the oil and gas zone effectiveness of a target interval of the mountain foreland district of Kuqa in Tarim Basin. The target interval is an oil and gas zone formed by source rocks of the Triassic and Jurassic systems and clastic rock reservoirs of the Bashijiqike formation of the Cretaceous system, and the trap of the target interval mainly accumulates natural gas.
The geophysical data of the target interval is obtained, and the structure, the stratum dip angle distribution, the fault and the unconformity within the effective hydrocarbon generation range of the target interval are obtained.
As shown in
As shown in
As shown in
As shown in
As shown in
As shown in
The well logging data of the target interval and the organic carbon content (TOC) obtained from core analysis of the source rock are collected to obtain the organic carbon content distribution of the target interval; As shown in
As shown in
Data of the organic carbon content (TOC) and the vitrinite reflectance (Ro) obtained from core analysis of the source rock of the target interval are collected, and the natural gas of the target interval come from the Jurassic and Triassic systems. The well logging data and seismic data of the target interval are collected, the thickness of the source rock is determined by seismic interpretation, and the TOC distribution is determined by well logging interpretation. The gas generation abundance distribution (
As shown in
As shown in
S51: obtaining a dimensionless flag value of each reference point according to the controlling factor parameter;
S53: obtaining the accumulation parameter value of each reference value according to the dimensionless flag value.
In this embodiment, the accumulation parameter value of each reference point is obtained according to the following formula:
wherein, AR—the accumulation parameter value, dimensionless, when AR≥1, it represents an effective accumulation region, and when AR<1, it represents a non-effective accumulation region; Parai—dimensionless parameter values Sths, Tss, Das, Rds and Fs; n—the number of Parai; Sths—a cap rock thickness flag value, dimensionless; Tss—a conduction system flag value, dimensionless; Das—a stratum dip angle flag value, dimensionless; Rds—an oil and gas migration distance flag value, dimensionless; Fs— a fault flag value.
In this embodiment, the flag value of the cap rock is determined according to the following formula:
wherein, Sths—the cap rock thickness flag value, dimensionless; Sth—the cap rock thickness of the target interval, m; Sthlim t—a cap rock thickness lower limit value of the industrial oil and gas reservoir, m.
In this embodiment, the conduction system flag value is determined according to the following formula:
wherein, Tss—the conduction system flag value, dimensionless; H—the reservoir thickness of the target interval, m; Hlim t—a reservoir thickness lower limit value of the industrial oil and gas reservoir, m; UF—unconformity of direct contact with the source rock within the effective hydrocarbon generation range and/or a fault; when there is unconformity of direct contact with the source rock within the effective hydrocarbon generation range and/or a fault, the UF is 1, and when there is no unconformity of direct contact with the source rock within the effective hydrocarbon generation range and/or fault, the UF is 0.
As shown in
S55: collecting an oil generation abundance lower limit value and a gas generation abundance lower limit value of the industrial oil and gas reservoir;
S57: determining the effective hydrocarbon generation range of the source rock according to the oil generation abundance lower limit value and the gas generation abundance lower limit value.
In this embodiment, the stratum dip angle flag value is determined according to the following formula:
Das=sin(α),
wherein, Das—the stratum dip angle flag value, dimensionless; α—the stratum dip angle, degree.
In this embodiment, the oil and gas migration distance flag value is determined according to the following formula:
wherein, Rds—the oil and gas migration distance flag value, dimensionless; Rd—an equivalent radius of the effective hydrocarbon generation range, km; Lhg—a distance to a boundary of the effective hydrocarbon generation range from outside the effective hydrocarbon generation range, km, which is 0 when within the effective hydrocarbon generation range; a—an empirical coefficient, which takes 3 when the target interval contains conventional oil and gas, takes 1.2 when the target interval contains dense oil and gas, and takes 1 when the target interval contains shale oil and gas.
In this embodiment, when the fault that cuts through the cap rock of the target interval is an open fault, the fault flag value is −1; and when the fault that cuts through the cap rock of the target interval is a sealed fault, the fault flag value is 1.
As shown in
In this embodiment, the evaluation parameter processing value of each grid coordinate point is determined according to the following formula:
wherein, Pi—the evaluation parameter processing value of the ith grid coordinate point; Si—the ith evaluation parameter value; Si_lim t—the ith evaluation parameter lower limit value; n—the number of the evaluation parameter.
As shown in
wherein, Vplay— the zone evaluation value, −1˜1, when Vplay≥0, it represents an effective zone distribution area, and when Vplay<0, it represents a non-effective zone distribution area; Pi—the evaluation parameter processing value of the ith grid coordinate point; n—the number of the grid coordinate point.
In this embodiment, the f(Pi) is determined according to the following formula:
wherein, Pi—the evaluation parameter processing value of the ith grid coordinate point.
As shown in
The oil and gas zone effectiveness evaluation apparatus according to the embodiment of the present application calculates the evaluation parameter processing value of each grid coordinate point according to the evaluation parameter, the evaluation parameter lower limit value and the accumulation parameter value of the grid coordinate point, thereby determining the accumulation range of the target interval, and then calculates the zone evaluation value of the target interval according to the evaluation parameter processing value of the grid coordinate point. In this way the oil and gas zone effectiveness evaluation is conducted quantitatively. Therefore, compared to the prior arts, the oil and gas zone effectiveness evaluation apparatus of the present application can have a unified standard of oil and gas zone effectiveness evaluation. Accordingly, the influence of artificial factors can be eliminated. Besides, since the evaluation parameter lower limit value of the oil and gas reservoir which gas got industrial value is employed, the identification rate of the oil and gas zone effectiveness is improved, and thereby the coincidence rate of the oil and gas effectiveness evaluation is improved. As such, the success rate of drilling and the speed of finding an effective zone can be improved, so that the costs for oil-gas exploration can be reduced.
In this embodiment, the controlling factor parameter obtaining module 11 is used for obtaining the controlling factor parameter distribution of the target interval, wherein the controlling factor parameter includes: stratum dip angle, fault, unconformity within an effective hydrocarbon generation range, reservoir thickness, cap rock thickness, oil generation abundance and gas generation abundance.
In a specific embodiment, the controlling factor parameter obtaining module can obtain the stratum dip angle distribution of the target interval according to the following steps:
1) obtaining structure of the target interval according to geophysical data of the target interval;
2) obtaining the stratum dip angle distribution of the target interval according to the structure of the target interval.
Further, the stratum dip angle of the target interval is determined according to the following steps:
1) obtaining a vertical line between an effective hydrocarbon generation center of the source rock and a top surface structural contour of the target interval;
2) obtaining an inclined angle of the vertical line:
3) obtaining an inclined direction of the inclined angle according to a positional relationship between the vertical line and a horizontal line;
to be specific, when the vertical line between the effective hydrocarbon generation center of the source rock and the top surface structural contour of the target interval is above the horizontal line of the effective hydrocarbon generation center of the source rock, the stratum dip angle is positive; when the vertical line between the effective hydrocarbon generation center of the source rock and the top surface structural contour of the target interval is below the horizontal line of the effective hydrocarbon generation center of the source rock, the stratum dip angle is negative; and when the vertical line between the effective hydrocarbon generation center of the source rock and the top surface structural contour of the target interval coincides with the horizontal line of the effective hydrocarbon generation center of the source rock, the stratum dip angle is 0; 4) obtaining the stratum dip angle according to the inclined angle and the inclined direction.
In a specific embodiment, the fault includes: a first fault that cuts through the cap rock of the target interval, a second fault that communicates the source rock and the reservoir within the effective hydrocarbon generation range of the source rock of the target interval, and a third fault that communicates the unconformity of the target interval and the reservoir.
In a specific embodiment, the unconformity within the effective hydrocarbon generation range can be an unconformity between the source rock and the cap rock.
In a specific embodiment, obtaining, by the controlling factor parameter obtaining module 11, the reservoir thickness distribution of the target interval can be realized by collecting the geophysical data of the target interval.
In a specific embodiment, obtaining, by the controlling factor parameter obtaining module 11, the porosity distribution of the target interval can be realized by collecting the porosity analysis data and geophysical data of the target interval.
In a specific embodiment, obtaining, by the controlling factor parameter obtaining module 11, the cap rock thickness distribution of the target interval can be realized by collecting the geophysical data of the target interval.
In a specific embodiment, obtaining, by the controlling factor parameter obtaining module 11, the oil generation abundance distribution and the gas generation abundance distribution of the target interval specifically includes:
1) collecting well logging data of the target interval and an organic carbon content obtained from core analysis of the source rock to obtain the organic carbon content distribution of the target interval;
2) collecting a vitrinite reflectance obtained from core analysis of the source rock of the target interval to obtain the vitrinite reflectance distribution;
3) performing a basin simulation operation to the organic carbon content distribution and the vitrinite reflectance distribution to obtain the oil generation abundance distribution and the gas generation abundance distribution.
In a specific embodiment, the evaluation parameter lower limit value obtaining module 13 can obtain the lower limit value of the reservoir thickness according to the reservoir thickness and oil and gas testing data of the industrial oil and gas reservoir, wherein, a reservoir type of the industrial oil and gas reservoir is the same as that of the target interval.
In a specific embodiment, the evaluation parameter lower limit value obtaining module 13 can obtain the porosity lower limit value according to porosity analysis data and gas testing data of the industrial oil and gas reservoir.
In a specific embodiment, the evaluation parameter lower limit value obtaining module 13 can obtain the lower limit value of the cap rock thickness according to the cap rock thickness of the industrial oil and gas reservoir, wherein, the lithology of the cap rock of the industrial oil and gas reservoir is consistent with that of the cap rock of the target interval.
In a specific embodiment, the oil generation abundance lower limit value and the gas generation abundance lower limit value of the industrial oil and gas reservoir are used to determine the effective hydrocarbon generation range of the source rock.
The oil and gas zone effectiveness evaluation apparatus according to the present application is employed to evaluate the oil and gas zone effectiveness of a target interval of the mountain foreland district of Kuqa in Tarim basin. The target interval is an oil and gas zone formed by source rocks of the Triassic and Jurassic systems and clastic rock reservoirs of the Bashijiqike formation of the Cretaceous system, and the trap of the target interval mainly accumulates natural gas.
The geophysical data of the target interval is obtained, and the structure, the stratum dip angle distribution, the fault and the unconformity within an effective hydrocarbon generation range of the target interval are obtained.
As shown in
As shown in
As shown in
As shown in
As shown in
As shown in
The well logging data of the target interval and the organic carbon content (TOC) obtained from core analysis of the source rock are collected to obtain the organic carbon content distribution of the target interval;
As shown in
As shown in
Data of the organic carbon content (TOC) and the vitrinite reflectance (Ro) obtained from core analysis of the source rock of the target interval are collected, and the natural gas of the target interval come from the Jurassic and Triassic systems. The well logging data and seismic data of the target interval are collected, the thickness of the source rock is determined by seismic interpretation, and the TOC distribution is determined by well logging interpretation. The gas generation abundance distribution (
The accumulation parameter value calculating module 15 can calculate the accumulation parameter value of each reference point according to the following steps:
1) obtaining a dimensionless flag value of each reference point according to the controlling factor parameter;
2) obtaining the accumulation parameter value of each reference value according to the dimensionless flag value.
In this embodiment, the accumulation parameter value of each reference point is obtained according to the following formula:
wherein, AR—the accumulation parameter value, dimensionless, when AR≥1, it represents an effective accumulation region, and when AR<1, it represents a non-effective accumulation region; Parai dimensionless parameter values Sths, Tss, Das, Rds and Fs; n—the number of Parai; Sths—a cap rock thickness flag value, dimensionless; Tss—a conduction system flag value, dimensionless; Das—a stratum dip angle flag value, dimensionless; Rds—an oil and gas migration distance flag value, dimensionless; Fs—a fault flag value.
In this embodiment, the cap rock thickness flag value is determined according to the following formula:
wherein, Sths—the cap rock thickness flag value, dimensionless; Sth—the cap rock thickness of the target interval, m; Sthlim t—a cap rock thickness lower limit value of the industrial oil and gas reservoir, m.
In this embodiment, the conduction system flag value is determined according to the following formula:
wherein, Tss—the conduction system flag value, dimensionless; H—the reservoir thickness of the target interval, m; Hlim t—a reservoir thickness lower limit value of the industrial oil and gas reservoir, m; UF—unconformity of direct contact with the source rock within the effective hydrocarbon generation range and/or a fault; when there is unconformity of direct contact with the source rock within the effective hydrocarbon generation range and/or a fault, the UF is 1, and when there is no unconformity of direct contact with the source rock within the effective hydrocarbon generation range and/or fault, the UF is 0.
As shown in
1) collecting an oil generation abundance lower limit value and a gas generation abundance lower limit value of the industrial oil and gas reservoir;
2) determining the effective hydrocarbon generation range of the source rock according to the oil generation abundance lower limit value and the gas generation abundance lower limit value.
In this embodiment, the stratum dip angle flag value is determined according to the following formula:
Das=sin(α),
wherein, Das—the stratum dip angle flag value, dimensionless; α—the stratum dip angle, degree.
In this embodiment, the oil and gas migration distance flag value is determined according to the following formula:
wherein, Rds—the oil and gas migration distance ag value, dimensionless; Rd—an equivalent radius of the effective hydrocarbon generation range, km; Lhg—a distance to a boundary of the effective hydrocarbon generation range from outside the effective hydrocarbon generation range, km, which is 0 when within the effective hydrocarbon generation range; a—an empirical coefficient, which takes 3 when the target interval contains conventional oil and gas, takes 1.2 when the target interval contains dense oil and gas, and takes 1 when the target interval contains shale oil and gas.
In this embodiment, when the fault that cuts through the cap rock of the target interval is an open fault, the flag value of the fault is −1; and when the fault that cuts through the cap rock of the target interval is a sealed fault, the flag value of the fault is 1.
As shown in
wherein, Pi—the evaluation parameter processing value of the ith grid coordinate point; Si—the ith evaluation parameter value; Si_lim t—the ith evaluation parameter lower limit value; n—the number of the evaluation parameter.
In this embodiment, the zone evaluation value calculating module 19 can determine the zone evaluation value of the target interval according to the following formula:
wherein, Vplay— the zone evaluation value, −1˜1, when Vplay≥0, it represents an effective zone distribution area, and when Vplay<0, it represents a non-effective zone distribution area; Pi—the evaluation parameter processing value of the ith grid coordinate point; n—the number of the grid coordinate point.
In this embodiment, f((Pi) is determined according to the following formula:
wherein, Pi—the evaluation parameter processing value of the ith grid coordinate point.
The above embodiments are only several embodiments of the present invention. Those skilled in the art can make various modifications or transformations to the embodiments of the present invention according to the contents disclosed in the present application document without deviating from the scope of the present invention.
Number | Date | Country | Kind |
---|---|---|---|
201710919263.6 | Sep 2017 | CN | national |
Number | Name | Date | Kind |
---|---|---|---|
10190998 | Zhou | Jan 2019 | B1 |
20130046476 | Schultz et al. | Feb 2013 | A1 |
20170051605 | Jackson | Feb 2017 | A1 |
20190080273 | Hou | Mar 2019 | A1 |
Number | Date | Country |
---|---|---|
103777245 | May 2014 | CN |
104166796 | Nov 2014 | CN |
104453873 | Mar 2015 | CN |
104965979 | Oct 2015 | CN |
105068144 | Nov 2015 | CN |
105469159 | Apr 2016 | CN |
106501854 | Dec 2016 | CN |
2012154912 | Jan 2013 | WO |
Entry |
---|
Chinese Search Report dated Nov. 14, 2018 for Chinese application No. 201710919263.6. |
Chinese Office Action dated Nov. 22, 2018 for Chinese application No. 201710919263.6. |
Weiyu et al; “Assessment for Petroleum Resource in East Greenland Basin”; Marine Geology Frontiers vol. 32, No. 8; Aug. 2016; pp. 1-11. |
Number | Date | Country | |
---|---|---|---|
20190100997 A1 | Apr 2019 | US |