CN111123357B - Method and device for determining evaluation index of unconventional oil and gas reservoir - Google Patents

Method and device for determining evaluation index of unconventional oil and gas reservoir Download PDF

Info

Publication number
CN111123357B
CN111123357B CN201811290163.2A CN201811290163A CN111123357B CN 111123357 B CN111123357 B CN 111123357B CN 201811290163 A CN201811290163 A CN 201811290163A CN 111123357 B CN111123357 B CN 111123357B
Authority
CN
China
Prior art keywords
reservoir
evaluation index
seismic
organic carbon
carbon content
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811290163.2A
Other languages
Chinese (zh)
Other versions
CN111123357A (en
Inventor
黄平
肖富森
张福宏
安辉
邹定永
马波
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201811290163.2A priority Critical patent/CN111123357B/en
Publication of CN111123357A publication Critical patent/CN111123357A/en
Application granted granted Critical
Publication of CN111123357B publication Critical patent/CN111123357B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48Processing data
    • G01V1/50Analysing data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/624Reservoir parameters

Abstract

The invention discloses a method and a device for determining an unconventional oil and gas reservoir evaluation index, and belongs to the technical field of oil and gas field geophysical prospecting engineering. The method comprises the following steps: acquiring an evaluation index of residual organic carbon content of the earthquake, an evaluation index of a relatively good reservoir, a relatively good reservoir thickness and a fracture density evaluation index; and obtaining an unconventional hydrocarbon reservoir evaluation index according to the evaluation index of the residual organic carbon content of the earthquake, the evaluation index of the relatively high-quality reservoir, the thickness of the relatively high-quality reservoir and the fracture density evaluation index. The invention provides comprehensive evaluation factors suitable for the unconventional oil and gas reservoir formation aiming at the unconventional oil and gas reservoir formation control factors, quantifies the comprehensive evaluation factors, achieves the aims of accurately evaluating the unconventional oil and gas reservoir and predicting the distribution of an oil and gas enrichment region, and improves the drilling success rate.

Description

Method and device for determining evaluation index of unconventional oil and gas reservoir
Technical Field
The invention belongs to the technical field of oil and gas field geophysical prospecting engineering, and particularly relates to a method and a device for determining an unconventional oil and gas reservoir evaluation index.
Background
The unconventional oil and gas reservoirs generally have the characteristics of complex geological conditions, poor physical properties of reservoirs and strong heterogeneity, and have high oil and gas exploitation difficulty and complex yield decreasing rule, so that the medium-and-long-term productivity prediction is difficult, and further the exploitation degree, the recoverable reserve, the final recovery ratio, the development age limit and the like in undevelopable periods cannot be accurately predicted. Therefore, it is necessary to provide an unconventional hydrocarbon reservoir evaluation index, provide comprehensive evaluation factors suitable for unconventional hydrocarbon reservoir formation aiming at the control factors of unconventional hydrocarbon reservoir formation, quantify the comprehensive evaluation factors, achieve accurate evaluation of unconventional hydrocarbon reservoirs, predict the distribution conditions of unconventional hydrocarbon reservoirs, improve the drilling success rate and improve the development efficiency of unconventional hydrocarbon reservoirs.
The method for determining the evaluation index of the oil and gas reservoir provided by the related technology mainly comprises the following steps: acquiring geological conditions and logging information of the oil and gas reservoir, and determining a reservoir space type and an interpretation threshold value of the reservoir of the oil and gas reservoir; obtaining a wave impedance data body of a reservoir through seismic inversion, obtaining a functional relation between the wave impedance data body and the reservoir pore space according to an interpretation threshold value and a reservoir space type of the reservoir, calculating the porosity and the reservoir thickness of the reservoir through the functional relation, and further drawing a plane distribution diagram of the porosity and the thickness of the reservoir; and drawing a reservoir fracture plane distribution diagram through seismic fracture prediction. And completing the analysis evaluation and prediction of the reservoir through the planar distribution diagram of the porosity and the thickness of the reservoir and the planar distribution diagram of the fracture of the reservoir.
The inventors found that the related art has at least the following technical problems:
the evaluation of the oil and gas reservoir only aims at evaluating the porosity size, the reservoir thickness, the crack development condition, the structure position and the like of the reservoir, and the evaluation index result obtained by evaluating and analyzing the method is not accurate and complete aiming at the unique geological condition and the reservoir physical property of the unconventional oil and gas reservoir.
Disclosure of Invention
The invention discloses a method and a device for determining an unconventional oil and gas reservoir evaluation index, which can solve the technical problem.
In one aspect, an embodiment of the present invention provides a method for determining an unconventional hydrocarbon reservoir evaluation index, where the method includes:
obtaining an evaluation index of residual organic carbon content of an earthquake, obtaining an evaluation index of a relatively high-quality reservoir, obtaining the thickness of the relatively high-quality reservoir and obtaining an evaluation index of fracture density;
obtaining an unconventional hydrocarbon reservoir evaluation index according to the evaluation index of the residual organic carbon content of the earthquake, the evaluation index of the relatively high-quality reservoir, the thickness of the relatively high-quality reservoir and the fracture density evaluation index;
wherein the unconventional hydrocarbon reservoir evaluation index is as follows:
Figure GDA0003422099730000021
wherein a, b, c, d and e are constants, i is the number of sample points, SUPTOC (i) is the evaluation index of the residual organic carbon content of the earthquake, SUPVEL (i) is the evaluation index of the relatively good reservoir,
Figure GDA0003422099730000022
supf (i) is the fracture density evaluation index for the relative premium reservoir thickness.
In an alternative embodiment, the obtaining an evaluation index of residual organic carbon content of the earthquake comprises:
obtaining a relation function of logging residual organic carbon content data and seismic data;
obtaining the residual organic carbon content of the earthquake according to the relation function of the logging residual organic carbon content data and the seismic data, and obtaining the residual organic carbon content evaluation index of the earthquake according to the residual organic carbon content of the earthquake, wherein the residual organic carbon content evaluation index of the earthquake is as follows:
SUPTOC(i)=(TOCseismic(i)-TOCmin)/(TOCmax-TOCmin);
wherein, TOCseismicAs a function of the relation between the logging residual organic carbon content data and the seismic data, TOCminAs the TOCseismicLower limit value of, TOCmaxAs the TOCseismicIs measured.
In an alternative embodiment, the obtaining a relation function of the logging residual organic carbon content data and the seismic data comprises:
acquiring logging residual organic carbon content data and seismic data;
obtaining a relation function between the logging residual organic carbon content data and the seismic data according to the logging residual organic carbon content data and the seismic data
The relation function between the logging residual organic carbon content data and the seismic data is as follows:
TOCseismic(i)=F[S(i)&TOCwell(i)];
wherein, TOCwell(i) (ii) is the logging residual organic carbon content data, S (i) is the seismic data, F]And obtaining a relation function of the logging residual organic carbon content data and the seismic data.
In an alternative embodiment, the obtaining the relatively good reservoir evaluation index includes:
obtaining the relatively good reservoir evaluation index according to the reservoir seismic inversion speed, the reservoir seismic inversion speed lower limit value and the reservoir seismic inversion speed minimum value:
SUPVEL(i)=(Vmax-V(i))/(Vmax-Vmin);
wherein V (i) is the reservoir seismic inversion velocity, VmaxIs the reservoir seismic inversion velocity lower limit, VminAnd the minimum value of the reservoir seismic inversion velocity is obtained.
In an alternative embodiment, obtaining the relatively good reservoir thickness comprises:
obtaining the thickness of the relatively high-quality reservoir according to the time thickness of the relatively high-quality reservoir and the seismic inversion speed of the relatively high-quality reservoir:
Figure GDA0003422099730000031
wherein HtFor said relatively good reservoir time thickness, VgoodAnd the seismic inversion speed of the relatively good reservoir is obtained.
In an alternative embodiment, obtaining the fracture density evaluation index comprises:
obtaining the fracture density evaluation index according to the fracture density, the fracture density threshold value and the maximum fracture density value:
SUPF(i)=(FM(i)-FMmin)/(FMmax-FMmin);
wherein i is the number of sample points, FM (i) is the fracture density, FMminIs the threshold value of the crack density, FMmaxIs the maximum fracture density.
In an alternative embodiment, obtaining the log residual organic carbon content data comprises:
acquiring the data of the content of the logging residual organic carbon according to the acoustic time difference:
TOCwell(i)=a′*AC(i)2+b′*AC(i)+c′;
wherein AC is the acoustic time difference.
In an alternative embodiment, acquiring the seismic data comprises:
acquiring the seismic data according to the seismic wavelets and the seismic reflection coefficients:
S(i)=R(i)*W(i);
wherein W (i) is seismic wavelet, and R (i) is seismic reflection coefficient.
In an alternative embodiment, the value of the constant a ranges from 0.001 to 0.01, the value of the constant b ranges from 0.01 to 0.06, and the absolute value of the constant c is less than 3.5.
In another aspect, the embodiment of the present invention further provides a device for determining an unconventional reservoir evaluation index, wherein the device is used in any one of the methods described above.
The technical scheme provided by the embodiment of the invention has the beneficial effects that at least:
the comprehensive evaluation factors suitable for the reservoir of the unconventional oil and gas reservoir are provided aiming at the reservoir formation control factors of the unconventional oil and gas reservoir and quantized, so that the aims of accurately evaluating the unconventional oil and gas reservoir and predicting the distribution of enriched oil and gas are fulfilled, and the drilling success rate is improved.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present application, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
FIG. 1 is a schematic flow chart of a method for determining an unconventional reservoir evaluation index provided by an embodiment of the invention;
FIG. 2 is a diagram of the TOC of the residual organic carbon content of the well logging provided by the embodiment of the inventionwell(i) A prediction plot of a relationship function between the data and seismic data S (i);
FIG. 3 is a schematic diagram of an unconventional reservoir evaluation index determination device provided by an embodiment of the invention;
fig. 4 is a schematic structural diagram of a first obtaining module in the device for determining the unconventional reservoir evaluation index according to the embodiment of the invention.
Detailed Description
To make the objects, technical solutions and advantages of the present disclosure more apparent, embodiments of the present disclosure will be described in detail with reference to the accompanying drawings.
The embodiment of the invention provides a method for determining an unconventional hydrocarbon reservoir evaluation index, which comprises the following steps of:
in step S101, obtaining an evaluation index of residual organic carbon content of an earthquake;
in step S102, a relatively good reservoir evaluation index is obtained;
in step S103, a relatively good reservoir thickness is obtained;
in step S104, a crack density evaluation index is acquired;
in step S105, obtaining an unconventional oil and gas reservoir evaluation index according to the residual organic carbon content evaluation index, the relatively good reservoir thickness and the fracture density evaluation index of the earthquake.
Wherein, the unconventional oil and gas reservoir evaluation index is as follows:
Figure GDA0003422099730000041
wherein a, b, c, d and e are constants, i is the number of sample points, SUPTOC (i) is the evaluation index of the residual organic carbon content of the earthquake, SUPVEL (i) is the evaluation index of the relatively good reservoir,
Figure GDA0003422099730000051
supf (i) is the fracture density evaluation index for relatively good reservoir thickness.
Next, steps S101 to S105 shown in fig. 1 will be explained.
Optionally, in step S101, acquiring an evaluation index of residual organic carbon content of the earthquake includes step S1011 and step S1012:
step S1011, obtaining a relation function of logging residual organic carbon content data and seismic data to obtain the seismic residual organic carbon content;
and step S1012, obtaining the residual organic carbon content of the earthquake according to the relation function of the logging residual organic carbon content data and the seismic data, and obtaining the evaluation index of the residual organic carbon content of the earthquake according to the residual organic carbon content of the earthquake.
Wherein the evaluation index of the residual organic carbon content of the earthquake obtained according to the residual organic carbon content of the earthquake is as follows:
SUPTOC(i)=(TOCseismic(i)-TOCmin)/(TOCmax-TOCmin);
wherein, TOCseismicAs a function of the relation between the logging residual organic carbon content data and the seismic data, TOCminAs TOCseismicLower limit value of, TOCmaxAs TOCseismicIs measured.
Analysis combined with geological understanding and oil and gas wells to determine TOC contributing to oil and gas productionseismicThe lower limit of the value, and then the evaluation index suptoc (i) of the residual organic carbon content of the earthquake is calculated according to the above formula.
According to the size of the SUPTOC (i), the reservoir can be evaluated by combining the actual well production condition. The larger the SUPTOC (i) value, the better the storage capacity of unconventional reservoirs.
Further, the TOC of the seismic residual organic carbon content is obtained by obtaining a relation function between the logging residual organic carbon content data and the seismic dataseismicIn time, the method provided by the embodiment of the application is used for measuring the TOC of the residual organic carbon content data of the wellwellThe distribution rule of the seismic data is analyzed, and the distribution rule of the seismic data is analyzed. Based on the analysis results, in general, the TOC data of the residual organic carbon content of the well logwellIf the seismic data also meet the positive-over distribution, searching the TOC of the logging residual organic carbon content datawellEstablishing a corresponding relation with the seismic data S (i) as follows:
TOCseismic(i)=F[S(i)&TOCwell(i)](ii) a Where i is 1, 2, 3, … … n, n is the number of spots, TOCwell(i) For logging residual organic carbon content data, F [ [ alpha ] ]]TOC data for established log residual organic carbon contentwell(i) And seismic data S (i). For example, FIG. 3 is a representation of a correspondence function between the two, wherein velocity is the primary parameter of the seismic data S (i).
Further, the TOC is inverted by stochastic seismic inversion according to the above formulaseismic
The seismic inversion technology is a process of imaging or solving the spatial structure and physical properties of an underground rock stratum by utilizing earth surface observation seismic data and taking known geological rules and well drilling and logging data as constraints. In the embodiment of the invention, the seismic data S (i) and TOC are fitted according to the sampling datawell(i) Approximate functional relationship, inverting TOC using the above relationshipseismic(i) The numerical value of (c).
Optionally, for step S1011, obtaining a relation function between logging residual organic carbon content data and seismic data, including steps S10111 to S10113:
step S10111, obtaining the TOC of the logging residual organic carbon content datawell(i);
In an alternative embodiment, logging residual organic carbon content data T is obtainedOCwellThe method comprises the following steps:
obtaining the data of the content of the residual organic carbon of the logging according to the acoustic time difference:
testing the residual organic carbon content TOC of the rock core sample in a laboratory, and fitting the residual organic carbon content TOC tested in the laboratory with the logging acoustic time difference to obtain the logging residual organic carbon content data TOCwell
The fitting relation is: TOCwell(i)=a′*AC(i)2+ b 'ac (i) + c'; where i is 1, 2, 3, … … n, n is the number of samples, a ', b ', c ' are constants, and AC is the acoustic time difference.
Alternatively, the value of the constant a may range from 0.001 to 0.01, for example, 0.001, 0.002, 0.003, 0.004, 0.005, 0.006, 0.007, 0.008, 0.009, 0.01, etc. The value of the constant b may range from 0.01 to 0.06, e.g., 0.01, 0.02, 0.03, 0.04, 0.05, 0.06, etc. The absolute value of the constant c is less than 3.5, e.g., 3.4, 3.3, 3.2, 3.1, 3.0, 2.0, 2.4, 1.4, etc.
Step S10112, acquiring seismic data S (i);
in an alternative embodiment, acquiring seismic data s (i) comprises:
acquiring seismic data according to the seismic wavelets and the seismic reflection coefficients: s (i) ═ r (i) × w (i); where i is 1, 2, 3, … … n, n is the number of sample points, w (i) is the seismic wavelet, and r (i) is the seismic reflection coefficient.
Seismic wavelets w (i) refer to a piece of signal with a defined start time, limited energy and a certain duration, which is the basic unit in a seismic recording. It is generally considered that the seismic wave generated by the seismic source excitation is only a sharp pulse with extremely short duration, the high-frequency component of the sharp pulse is attenuated rapidly as the sharp pulse propagates in the viscoelastic medium, the waveform is increased, and a seismic wavelet is formed, wherein the seismic wavelet generally has a duration of 2 to 3 phases, about 90 milliseconds, and then propagates underground in the form of the seismic wavelet.
The seismic wavelet is generally generated by being excited by a seismic source in the field, and the acquired seismic wavelet data is summarized and sorted by acquiring seismic wavelet data to obtain the seismic wavelet S (i) in the embodiment of the invention.
The seismic reflection coefficient R (i) is the ratio of the reflection of the seismic reflection through the formation to the incident.
And S10113, acquiring a relation function between the logging residual organic carbon content data and the seismic data according to the logging residual organic carbon content data.
Wherein, the relation function between the logging residual organic carbon content data and the seismic data is as follows:
TOCseismic(i)=F[S(i)&TOCwell(i)];
wherein, TOCwell(i) For logging residual organic carbon content data, S (i) is seismic data, F [, ]]And obtaining a relation function of the obtained logging residual organic carbon content data and the seismic data.
For step S102, obtaining a relative goodness reservoir evaluation index supvel (i) includes: and obtaining a relatively high-quality reservoir stratum evaluation index SUPVEL according to the reservoir stratum seismic inversion speed, the reservoir stratum seismic inversion speed lower limit value and the reservoir stratum seismic inversion speed minimum value.
SUPVEL(i)=(Vmax-V(i))/(Vmax-Vmin)。
Wherein i is the number of sample points, V (i) is the reservoir seismic inversion velocity, VmaxIs the reservoir seismic inversion velocity lower limit, VminIs the minimum value of the reservoir seismic inversion velocity.
And (3) inverting the seismic velocity data volume by using a pre-stack or post-stack seismic inversion technology, wherein the higher the velocity is, the denser the rock stratum is, and when the velocity is less than a certain threshold value, the property of the reservoir stratum can be reflected. Generally speaking, reservoirs with different lithologies have different threshold values, the reservoir threshold value data can be determined according to whether the existing well can form industrial airflow, and the quality of the reservoir can be reflected through speed change. The higher the speed is, the more compact the reservoir is, the smaller the pore space is, the oil and gas are to be accumulated in the compact reservoir, the capillary resistance of oil-gas-water multiphase fluid in the capillary needs to be overcome, and the more compact the reservoir is, the larger the capillary resistance is. And (3) explaining the reservoir data by combining logging, determining a speed lower limit value corresponding to the relatively high-quality reservoir of the unconventional oil and gas reservoir through analysis, and calculating the relatively high-quality reservoir evaluation index of the unconventional oil and gas reservoir by using the following formula according to the analysis result:
it is understood that the above-mentioned well logging interpretation of reservoir data is to perform a test by a logging personnel and interpret the test result, and the final data after interpretation is the well logging interpretation reservoir data required by the embodiment of the present invention.
For step S103, a relatively good reservoir thickness is obtained
Figure GDA0003422099730000071
The method comprises the following steps: obtaining the thickness of the relatively high-quality reservoir according to the time thickness of the relatively high-quality reservoir and the seismic inversion speed of the relatively high-quality reservoir:
Figure GDA0003422099730000072
wherein HtFor relatively good reservoir time thickness, VgoodThe seismic inversion speed of the relatively good reservoir is obtained.
According to the relative quality reservoir evaluation index, combining the completed well for analysis, determining the evaluation index quality reservoir threshold value SUPVEL (i)minOr determining a higher threshold value according to the well yield as the threshold value of the high-quality reservoir, wherein the evaluation index value of the relatively high-quality reservoir to be counted needs to be greater than or equal to the threshold value, the number of the evaluation index values which are greater than or equal to the threshold value is counted, and the sum of the numbers is the time thickness H of the relatively high-quality reservoirtTime thickness of relatively good reservoir HtMultiplying by the average of the velocities corresponding to the points, i.e. the seismic inversion velocity V of the relatively good reservoirgoodTo obtain relatively good reservoir thickness
Figure GDA0003422099730000081
For step S104, obtaining the fracture density evaluation index supf (i) includes: acquiring a fracture density evaluation index according to the fracture density, the fracture density threshold value and the maximum fracture density value:
SUPF(i)=(FM(i)-FMmin)/(FMmax-FMmin)。
wherein i is the number of sample points, FM (i) is the crack density, FMminIs the threshold value of crack density, FMmaxIs the maximum fracture density.
And calculating the fracture density by using a prestack fracture inversion technology, combining imaging logging information or other logging interpretation fracture results, determining a fracture density interpretation threshold value, and calculating a fracture density evaluation index according to the formula. Wherein, FM (i) is obtained by calculation of prestack crack prediction software, and the maximum value FM of crack density is selected by screening the obtained FM (i) numerical valuemaxAnd crack density threshold FMmin
For step S105, obtaining the unconventional reservoir evaluation index p (i) includes:
according to the reservoir forming control factors of unconventional oil and gas reservoirs, the seismic data obtained in the steps, the evaluation index SUPTOC (i) of the residual organic carbon content of the earthquake, the evaluation index SUPVEL (i) of the relatively good reservoir, and the thickness of the relatively good reservoir
Figure GDA0003422099730000082
And fracture density evaluation index SUPF (i), wherein the obtained unconventional oil and gas reservoir evaluation index is as follows:
Figure GDA0003422099730000083
wherein a, b, c, d and e are constants, i is 1, 2, 3, … … n, and n is the number of sample points; in general, a, b and d are less than 1, c is less than 0.1 and e is less than 3 absolute.
In another aspect, an embodiment of the present invention further provides an apparatus for determining an unconventional hydrocarbon evaluation index of a hydrocarbon reservoir, as shown in fig. 4, the apparatus including:
the first obtaining module 401 is configured to obtain an evaluation index of residual organic carbon content of an earthquake;
a second obtaining module 402, configured to obtain a relatively high-quality reservoir evaluation index;
a third obtaining module 403, configured to obtain a relatively high-quality reservoir thickness;
a fourth obtaining module 404, configured to obtain a fracture density evaluation index;
a fifth obtaining module 405, configured to obtain an unconventional hydrocarbon reservoir evaluation index according to the residual organic carbon content evaluation index, the relatively high-quality reservoir thickness, and the fracture density evaluation index of the earthquake;
wherein, the unconventional oil and gas reservoir evaluation index is as follows:
Figure GDA0003422099730000091
wherein a, b, c, d and e are constants, i is the number of sample points, SUPTOC (i) is the evaluation index of the residual organic carbon content of the earthquake, SUPVEL (i) is the evaluation index of the relatively good reservoir,
Figure GDA0003422099730000092
supf (i) is the fracture density evaluation index for relatively good reservoir thickness.
In an optional implementation, the first obtaining module 401 includes:
the first obtaining unit 4011 is configured to obtain a relation function between logging residual organic carbon content data and seismic data;
the second obtaining unit 4012 is used for obtaining the residual organic carbon content of the earthquake;
and the third obtaining unit 4013 is configured to obtain an evaluation index of the residual organic carbon content of the earthquake according to the residual organic carbon content of the earthquake. The evaluation index of the residual organic carbon content of the earthquake is as follows:
SUPTOC(i)=(TOCseismic(i)-TOCmin)/(TOCmax-TOCmin);
wherein, TOCseismicAs a function of the relation between the logging residual organic carbon content data and the seismic data, TOCminAs TOCseismicLower limit value of, TOCmaxAs TOCseismicIs measured.
In an alternative embodiment, the first obtaining unit 4011 includes:
the first acquisition subunit is used for acquiring the residual organic carbon content data of the logging;
the second acquisition subunit is used for acquiring seismic data;
and the third obtaining subunit is used for obtaining a relation function between the logging residual organic carbon content data and the seismic data according to the logging residual organic carbon content data and the seismic data. The relationship function is:
TOCseismic(i)=F[S(i)&TOCwell(i)];
wherein, TOCwell(i) For logging residual organic carbon content data, S (i) is seismic data, F [, ]]And obtaining a relation function of the obtained logging residual organic carbon content data and the seismic data.
In an optional implementation manner, the second obtaining module 402 is configured to obtain a relatively good reservoir evaluation index according to the reservoir seismic inversion velocity, the reservoir seismic inversion velocity lower limit value, and the reservoir seismic inversion velocity minimum value:
SUPVEL(i)=(Vmax-V(i))/(Vmax-Vmin);
wherein i is the number of sample points, V (i) is the reservoir seismic inversion velocity, VmaxIs the reservoir seismic inversion velocity lower limit, VminIs the minimum value of the reservoir seismic inversion velocity.
In an alternative embodiment, the third obtaining module 403 is configured to obtain a relatively good reservoir thickness, and includes:
obtaining the thickness of the relatively high-quality reservoir according to the time thickness of the relatively high-quality reservoir and the seismic inversion speed of the relatively high-quality reservoir:
Figure GDA0003422099730000101
wherein HtFor relatively good reservoir time thickness, VgoodThe seismic inversion speed of the relatively good reservoir is obtained.
In an alternative embodiment, the fourth obtaining module 404 is configured to obtain a fracture density evaluation index, and includes:
acquiring a fracture density evaluation index according to the fracture density, the fracture density threshold value and the maximum fracture density value:
SUPF(i)=(FM(i)-FMmin)/(FMmax-FMmin);
wherein i is the number of sample points, FM (i) is the crack density, FMminIs the threshold value of crack density, FMmaxIs the maximum fracture density.
In an alternative embodiment, the first acquiring subunit is configured to acquire logging residual organic carbon content data according to acoustic moveout: TOCwell(i)=a′*AC(i)2+b′*AC(i)+c′;
Wherein i is the number of sample points, a ', b ', c ' are constants, and AC is the acoustic time difference.
In an alternative embodiment, the second acquiring subunit is configured to acquire the seismic data according to the seismic wavelet and the seismic reflection coefficient: s (i) ═ r (i) × w (i);
wherein i is the number of sample points, W (i) is the seismic wavelet, and R (i) is the seismic reflection coefficient.
All the above optional technical solutions may be combined arbitrarily to form the optional embodiments of the present disclosure, and are not described herein again.
The above description is only an illustrative embodiment of the present invention, and should not be taken as limiting the scope of the invention, and any modifications, equivalents, improvements and the like that are within the spirit and principle of the present invention should be included in the scope of the present invention.

Claims (7)

1. A method of determining an unconventional reservoir evaluation index, the method comprising:
obtaining an evaluation index of residual organic carbon content of an earthquake, obtaining an evaluation index of a relatively high-quality reservoir, obtaining the thickness of the relatively high-quality reservoir and obtaining an evaluation index of fracture density; the step of obtaining the evaluation index of the residual organic carbon content of the earthquake comprises the following steps: obtaining a relation function of logging residual organic carbon content data and seismic data; according to the relation function of the logging residual organic carbon content data and the seismic dataObtaining the residual organic carbon content of the earthquake, and obtaining the evaluation index of the residual organic carbon content of the earthquake according to the residual organic carbon content of the earthquake, wherein the evaluation index of the residual organic carbon content of the earthquake is as follows: SUPTOC (i) ═ TOCseismic(i)-TOCmin)/(TOCmax-TOCmin) (ii) a Wherein, TOCseismicAs a function of the relation between the logging residual organic carbon content data and the seismic data, TOCminAs the TOCseismicLower limit value of, TOCmaxAs the TOCseismicMaximum value of (d); the acquiring of the relatively good reservoir evaluation index comprises the following steps: obtaining the relatively good reservoir evaluation index according to the reservoir seismic inversion speed, the reservoir seismic inversion speed lower limit value and the reservoir seismic inversion speed minimum value: supvel (i) ═ Vmax-V(i))/(Vmax-Vmin) (ii) a Wherein V (i) is the reservoir seismic inversion velocity, VmaxIs the reservoir seismic inversion velocity lower limit, VminIs the reservoir seismic inversion velocity minimum; the reservoir seismic inversion speed is used for reflecting the properties of a reservoir, and is positively correlated with the compactness of the reservoir; obtaining the fracture density evaluation index comprises: obtaining the fracture density evaluation index according to the fracture density, the fracture density threshold value and the maximum fracture density value: supf (i) ═ FM (i) — FMmin)/(FMmax-FMmin) (ii) a Wherein FM (i) is the fracture density, FMminIs the threshold value of the crack density, FMmaxIs the fracture density maximum;
obtaining an unconventional hydrocarbon reservoir evaluation index according to the evaluation index of the residual organic carbon content of the earthquake, the evaluation index of the relatively high-quality reservoir, the thickness of the relatively high-quality reservoir and the fracture density evaluation index;
wherein the unconventional hydrocarbon reservoir evaluation index is as follows:
Figure FDA0003422099720000011
wherein a, b, c, d and e are constant(ii) the number of samples, SUPTOC (i) the evaluation index for the residual organic carbon content of the earthquake, SUPVEL (i) the evaluation index for the relatively good reservoir,
Figure FDA0003422099720000012
supf (i) is the fracture density evaluation index for the relative premium reservoir thickness.
2. The method for determining the unconventional hydrocarbon reservoir evaluation index of claim 1, wherein the obtaining the relation function of the logging residual organic carbon content data and the seismic data comprises:
acquiring logging residual organic carbon content data and seismic data;
acquiring a relation function between the logging residual organic carbon content data and the seismic data according to the logging residual organic carbon content data and the seismic data;
the relation function between the logging residual organic carbon content data and the seismic data is as follows:
TOCseismic(i)=F[S(i)&TOCwell(i)];
wherein, TOCwell(i) (ii) is the logging residual organic carbon content data, S (i) is the seismic data, F]And obtaining a relation function of the logging residual organic carbon content data and the seismic data.
3. The method of claim 1, wherein obtaining the relatively good reservoir thickness comprises:
obtaining the thickness of the relatively high-quality reservoir according to the time thickness of the relatively high-quality reservoir and the seismic inversion speed of the relatively high-quality reservoir:
Figure FDA0003422099720000021
wherein HtFor said relatively good reservoir time thickness, VgoodAnd the seismic inversion speed of the relatively good reservoir is obtained.
4. The method of claim 2, wherein obtaining the log residual organic carbon content data comprises:
acquiring the data of the content of the logging residual organic carbon according to the acoustic time difference:
TOCwell(i)=a′*AC(i)2+b′*AC(i)+c′;
wherein AC is the acoustic time difference.
5. The method of claim 2, wherein obtaining the seismic data comprises:
acquiring the seismic data according to the seismic wavelets and the seismic reflection coefficients:
S(i)=R(i)*W(i);
wherein W (i) is seismic wavelet, and R (i) is seismic reflection coefficient.
6. The method of claim 1, wherein the constant a is in the range of 0.001-0.01, the constant b is in the range of 0.01-0.06, and the constant c is less than 3.5 in absolute value.
7. An apparatus for determining an unconventional reservoir evaluation index, wherein the apparatus is used in the method of any one of claims 1 to 6.
CN201811290163.2A 2018-10-31 2018-10-31 Method and device for determining evaluation index of unconventional oil and gas reservoir Active CN111123357B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811290163.2A CN111123357B (en) 2018-10-31 2018-10-31 Method and device for determining evaluation index of unconventional oil and gas reservoir

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811290163.2A CN111123357B (en) 2018-10-31 2018-10-31 Method and device for determining evaluation index of unconventional oil and gas reservoir

Publications (2)

Publication Number Publication Date
CN111123357A CN111123357A (en) 2020-05-08
CN111123357B true CN111123357B (en) 2022-03-29

Family

ID=70494286

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811290163.2A Active CN111123357B (en) 2018-10-31 2018-10-31 Method and device for determining evaluation index of unconventional oil and gas reservoir

Country Status (1)

Country Link
CN (1) CN111123357B (en)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4360359A (en) * 1981-03-13 1982-11-23 Conoco Inc. Method for relating shallow electrical anomalies to the presence of deeper hydrocarbon reservoirs
CN103256046A (en) * 2013-04-28 2013-08-21 北京大学 Unconventional oil and gas reservoir horizontal well section full-fracture-length fracturing parameter analog method and device
CN105257252A (en) * 2015-06-08 2016-01-20 中国石油集团川庆钻探工程有限公司 Method for optimally selecting shale gas horizontal well clustering perforation well section by utilizing logging data
CN105426650A (en) * 2015-06-08 2016-03-23 中国石油集团川庆钻探工程有限公司 Method for selecting out shale gas fracturing modification well section by using logging data
CN105842753A (en) * 2016-03-24 2016-08-10 中国石油大学(北京) Method and device for recovering residual carbon content in hydrocarbon source rocks
CN106526669A (en) * 2016-09-19 2017-03-22 中国石油化工股份有限公司 Shale oil and gas reservoir seismic reservoir prediction method
CN107807407A (en) * 2017-09-30 2018-03-16 中国石油天然气股份有限公司 A kind of petroleum zone efficiency evaluation method and apparatus

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2895887B1 (en) * 2012-09-12 2019-11-06 BP Exploration Operating Company Limited System and method for determining retained hydrocarbon fluid
CN104977618B (en) * 2014-04-09 2017-09-22 中国石油集团东方地球物理勘探有限责任公司 A kind of method evaluated shale gas reservoir and find dessert area

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4360359A (en) * 1981-03-13 1982-11-23 Conoco Inc. Method for relating shallow electrical anomalies to the presence of deeper hydrocarbon reservoirs
CN103256046A (en) * 2013-04-28 2013-08-21 北京大学 Unconventional oil and gas reservoir horizontal well section full-fracture-length fracturing parameter analog method and device
CN105257252A (en) * 2015-06-08 2016-01-20 中国石油集团川庆钻探工程有限公司 Method for optimally selecting shale gas horizontal well clustering perforation well section by utilizing logging data
CN105426650A (en) * 2015-06-08 2016-03-23 中国石油集团川庆钻探工程有限公司 Method for selecting out shale gas fracturing modification well section by using logging data
CN105842753A (en) * 2016-03-24 2016-08-10 中国石油大学(北京) Method and device for recovering residual carbon content in hydrocarbon source rocks
CN106526669A (en) * 2016-09-19 2017-03-22 中国石油化工股份有限公司 Shale oil and gas reservoir seismic reservoir prediction method
CN107807407A (en) * 2017-09-30 2018-03-16 中国石油天然气股份有限公司 A kind of petroleum zone efficiency evaluation method and apparatus

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
川中侏罗系致密油藏烃源预测技术研究;黄平,等;《天然气勘探与开发》;20131031;第54-58页 *
页岩储层可压性评价关键指标体系;侯冰,等;《石油化工高等学校学报》;20141231;第42-49页 *
页岩气储层研究现状及存在问题探讨;胡明毅,等;《特种油气藏》;20150430;第1-8页 *

Also Published As

Publication number Publication date
CN111123357A (en) 2020-05-08

Similar Documents

Publication Publication Date Title
CN1040364C (en) Method for analising seismic traces using a statistical calibration technique for deduction of geological characteristics
CN105277982B (en) A kind of mud shale total content of organic carbon earthquake prediction method
US6789018B1 (en) Mapping reservoir rocks using frequency spectral broadening and the presence of the slow-wave
CN1734288A (en) Method for using data expressed in differing time scales to construct a model of a heterogeneous environment described by several parameters
CN112255688B (en) Method for inverting formation pressure by three-dimensional earthquake based on rock physics theory
CN111897011B (en) Reservoir pore characteristic determination method, device and equipment
CN104278991A (en) Multivariate well logging computing method for total organic carbon and hydrocarbon generating potential of hydrocarbon source rocks in salt-lake facies
CN112922591B (en) Shale reservoir lithofacies dessert prediction method and system
CN103630939A (en) Air layer identification and evaluation method
CN111381292B (en) Logging interpretation method and device for predicting sandstone hydrocarbon-bearing reservoir
CN112363226A (en) Geophysical prediction method for unconventional oil and gas favorable area
CN113216945B (en) Quantitative evaluation method for permeability of tight sandstone reservoir
CN102253414A (en) Reservoir detecting method based on analysis of earthquake lines
CN111123357B (en) Method and device for determining evaluation index of unconventional oil and gas reservoir
CN106568918B (en) Shale organic carbon content TOC prediction method
CN103197348A (en) Method using internal samples at reservoirs to carry out weighting and compile logging crossplot
CN112711068A (en) Method and device for predicting effective oil and gas reservoir in sandstone
CN112835124B (en) Crack effectiveness evaluation method based on imaging logging and array acoustic logging data
CN112505764B (en) High-porosity hydrocarbon-containing sandstone reservoir prediction method and device
CN115099014A (en) Natural gas well geological exploration reserves calculation method based on logging while drilling
CN112211628A (en) Stratum dividing method based on multiple logging curves
RU2210094C1 (en) Method of geophysical prospecting to establish filtration capacitive properties of oil-and-gas bearing deposits in interwell space
CN108931813B (en) Method for identifying reef cover, reef core and reef base
CN114428049B (en) Method for calculating asphalt content of ancient carbonate reservoir
CN113050168B (en) Crack effectiveness evaluation method based on array acoustic logging and acoustic remote detection logging data

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant