CN103278866A - Evaluation method of shale oil resource potential in shale strata series - Google Patents

Abstract

The invention relates to an evaluation method of a shale oil resource potential in shale strata series. The evaluation method comprises the steps that essential data are sorted, wells are selected evenly in work areas, models are built, organic geochemistry well logging is verified, shale oil grading evaluation standard is built, a classification standard of the shale strata series is put forward, light dydrocarbon and heavy hydrocarbon in the shale oil in the shale strata series are revised, saturation adsorption oil mass of shale is calculated, shale oil resource potential is evaluated in a grading mode according to the grading standard, shale saturation adsorption oil mass is determined with application of a geochemical data statistic analysis graphic method, then movable quantity of shale oil is calculated, and shale oil resource quantity of a sandstone thin interlayer of the shale strata series is calculated with utilization of the 'hole-oil saturation' method. According to the evaluation method of the shale oil resource potential in the shale strata series, due to the fact that the concept of the shale strata series is put forward and built, the scientific space position of the shale oil is clearly determined. The evaluation method further puts forwards a calculation method of the shale saturation adsorption oil quantity and the shale oil movable quantity, and significance is provided for judging whether a work area has a shale oil development value.

Description

Oil shale fuel resource potential evaluation method in a kind of mud shale series of strata

Technical field

The present invention relates to a kind of unconventional oil-gas exploration evaluation method, particularly about oil shale fuel resource potential evaluation method in a kind of mud shale series of strata.

Background technology

The north America region of successful exploration exploitation shale oil gas about the basic standard that shale oil gas has enrichment, recoverability is in the world: abundance of organic matter is higher, degree of ripeness is higher, general organic matter type better, brittle mineral content higher (clay mineral content is lower).Yet do not propose and relevant concept and the algorithm of setting up mud shale series of strata and saturated absorption oil mass.Its major reason is, the north America region of the shale of successful exploration exploitation at present oil gas is substantially in relatively uniform marine bed, it diminishes mutually, homogenieity is stronger, utilizes limited drilling well and lab analysis data can determine the favourable growth district of shale oil gas and the locus of science substantially.

Though domestic have some scholars (Xu Yuangang etc., 2012; Wan Changmao etc., 2012; Liu Xiaoping etc., 2011; Zhu Dan etc., 2012) once putting forward the mud shale series of strata aspect the research of shale gas, but because lack the rational and effective parameter, accurate concept and the criteria for classifying of mud shale series of strata are not proposed, simultaneously main is the related difficult parameters of mud shale series of strata with the obtaining in early days of resource evaluation, so be difficult to be approved widely and apply.About the criteria for classifying of oil shale fuel mud shale series of strata, so far and do not see relevant document or research report report.Many scholars have realized that the rifted-basin in China terrestrial facies basin, especially east, and it accelerates mutually, mud shale is organic and inorganic nonuniformity is strong.Yet but almost nobody attempts to come the tax of clear and definite oil shale fuel to deposit the locus by setting up the mud shale series of strata criteria for classifying.Initial field operations personnel carry out the evaluation of oil shale fuel resource potential, be that oil shale fuel in underground all mud shales (not containing short lap) is estimated, result calculated is very surprising, and there is not grading evaluation, which resource to be effectively to exploit in the recent period as, which resource is to can be used as the distant view resource, and the wait technical progress could effectively be exploited, and which resource does not have economic worth.Although the viewpoint that oil shale fuel content is very big in the underground mud shale is correct, if do not carry out the consideration of resource classification and series of strata, the result who obtains does not have reference value for the exploitation of oil shale fuel.Thinner as mud shale thickness, all be the condition (the general minimum 15m of being of pressure break radius) that adopts the exploitation of fractured horizontal well technology and do not satisfy present oil shale fuel, the oil shale fuel stock number that calculates like this is higher.In sum, the technology and method that does not still have oil shale fuel mobility/amount evaluation aspect at present.

Summary of the invention

At the problems referred to above, the objective of the invention is to have proposed a kind of mud shale series of strata criteria for classifying, deposit the locus with the tax of clear and definite oil shale fuel, solve in the past the blindness to the evaluation of oil shale fuel.Determine the saturated absorption oil mass of mud shale, but and calculate the oil shale fuel momentum, but solve the difficult problem that people in the past can't estimate the oil shale fuel momentum, for oil shale fuel actual investigation and prospecting in oil field provides technical support.

In order to realize the foregoing invention purpose, the present invention is by the following technical solutions: oil shale fuel resource potential evaluation method in a kind of mud shale series of strata, may further comprise the steps: 1) arrangement basic data, described basic data comprises: geologic data, logging trace data, organic and inorganic geochemistry data; 2) in the work area, evenly select well, and encrypt at core space well number, according to selected yardstick with the gridding of described work area, 1～3 mouthful of described well of primary election in each grid; 3) organic geochemistry well logging modeling and checking: utilize unconventional petroleum resources to estimate platform selecting representative well wherein and carry out organic nonuniformity modeling, comprise pyrolysis S ₁, chloroform bitumen " A " and TOC, further select for use another representative well that the model of setting up is verified; 4) set up oil shale fuel grading evaluation scheme: according to " three fens property " feature of mud shale oil content, be that the III level is disperseed resource with expression by enrichment degree with the oil shale fuel resource division, the II level is with expression poor efficiency resource, and the I level is with expression enrichment resource; 5) according to U.S.'s oil shale fuel Development characteristics and China's mud shale deposition characteristics, set up the criteria for classifying of mud shale series of strata, and according to the described criteria for classifying, write mud shale series of strata identification system, automatically identify the mud shale series of strata of every mouthful of well by described mud shale series of strata identification system; 6) lighter hydrocarbons, the heavy hydrocarbon of oil shale fuel in the described mud shale series of strata are proofreaied and correct; 7) oil shale fuel resource potential grading evaluation: carry out oil shale fuel resource potential grading evaluation with the described grading evaluation scheme in volumetric method and the step 4); 8) determine the saturated absorption oil mass of described mud shale by geochemistry data statistical study graphical method, but and further use material balance principle to calculate the oil shale fuel momentum; 9) utilize " hole-oil saturation " method to calculate the oil shale fuel stock number of sandstone short lap in the described mud shale series of strata, finish oil shale fuel resource potential grading evaluation.

Described 5) in the step, the described criteria for classifying of mud shale series of strata is: mud shale series of strata The initial segment is the mud shale section that thickness surpasses 2m, and the interior sandstone proportion of any 10m is less than 1/3 downwards, and sandstone thickness is less than 2m continuously, mud shale Ro greater than 0.5%, TOC greater than 0.5%, pyrolysis S ₁Greater than 0.5mgHC/g, concluding paragraph is the mud stone section, and mud shale and sandstone cumulative thickness are greater than 30m.

Described 5) in the step, according to the described criteria for classifying, write described mud shale series of strata identification system, automatically identify the described mud shale series of strata of every mouthful of well by described mud shale series of strata identification system, specific as follows: 1. to having selected the logging data of well to put in order in the described work area; 2. utilize the organic nonuniformity model of well logging to calculate continuous described TOC, described pyrolysis S ₁, described chloroform bitumen " A ", obtain the described TOC of described mud stone section, described pyrolysis S ₁, described chloroform bitumen " A " mean value; 3. according to terrestrial heat flow, underground temperature gradient, formation temperature, reservoir pressure, the data such as described Ro, Tmax of basin/depression, make the degree of ripeness thermal evolution history of described basin/depression, determine every mouthful of described Ro of well at last greater than 0.5% degree of depth section; 4. according to described mud shale series of strata identification system require form with the above-mentioned data of handling well, the described data of handling well comprise well location, the degree of depth, lithology, described TOC average, described pyrolysis S ₁Average, described chloroform bitumen " A " average, write an identification system, with the described identification system of described mud shale series of strata the described mud shale series of strata of every mouthful of well are identified, for described mud shale series of strata internal classification stock number and movable RESERVE EVALUATION lay the foundation.

Described 8) in the step, determine the saturated absorption oil mass of described mud shale Q by described geochemistry data statistical study graphical method _{The saturated absorption oil mass of mud shale}But, and further utilize material balance principle to calculate described oil shale fuel momentum Q _{But oil shale fuel momentum}But, described oil shale fuel momentum Q _{But oil shale fuel momentum}Computing formula be:

Q _{But oil shale fuel momentum}=Q _{Original place oil shale fuel stock number}-Q _{The saturated absorption oil mass of mud shale}

=V×ρ×(A%×K _A-TOC ^*×K×K _TOC)

=Q _{Original place oil shale fuel stock number}-V * ρ * TOC ^** K * K _TOC

In the formula, ρ is mud shale density, and unit is g/cm ³V is the volume of mud shale, and unit is m ³TOC ^*Be the TOC average of three rank shale, unit is %, I level S ₁Greater than the 2.0mgHC/g rock, II level S ₁Between the 0.5-2.0mgHC/g rock, III level S ₁Less than the 0.5mgHC/g rock; K _ABe the coefficient of restitution of chloroform bitumen " A ", dimensionless; K _TOCBe the coefficient of restitution of organic carbon TOC, dimensionless; K is saturated adsorption coefficient, dimensionless.

The present invention is owing to take above technical scheme, it has the following advantages: 1, the present invention is owing to propose and set up mud shale series of strata concept and standard, the clear and definite science locus of evaluation object oil shale fuel, take into account the validity of oil shale fuel in the series of strata and the exploitation condition of oil shale fuel, abandoned the oil shale fuel resource that does not possess exploitation value at present.Therefore, make vast oil field scientific research personnel to the evaluation of the oil shale fuel evaluation from all mud shale sections blindly, be specially the evaluation in the mud shale series of strata that the proposition of this concept and standard has science and universality.2, the present invention determines the saturated absorption oil mass of mud shale owing to be according to geochemistry data statistical study graphical method, but and further calculate the oil shale fuel momentum, estimate basin/depression/hollow interior oil shale fuel resource potential easily.Therefore, but solved the geology difficult problem that people in the past can't estimate, have no way of doing it the oil shale fuel momentum effectively, oil shale fuel classification stock number and movable stock number have scientifically been calculated, the idiographic flow of oil shale fuel resource evaluation has been proposed, for oil shale fuel actual investigation and prospecting in oil field provides technical support, evaluation result has scientific and characteristics with a high credibility.3, the inventive method is significant for the exploitation value whether the judgement work area possesses oil shale fuel, and has important use value for unconventional shale oil gas resource exploration and exploitation.The present invention can be widely used in the exploratory development process of oil shale fuel resource in the various mud shale series of strata.

Description of drawings

Fig. 1 is schematic flow sheet of the present invention

Fig. 2 be the husky three times inferior sections of the hollow SHAHEJIE FORMATION of BHW basin BN select well achievement synoptic diagram

Fig. 3 is R ₂₅S is calculated in the well logging of-CNL mud shale ₁With actual measurement S ₁Concern synoptic diagram

Fig. 4 is R ₂₅-CNL mud shale amplitude difference and actual measurement S ₁Concern synoptic diagram

Fig. 5 is X1 well R ₂₅-CNLS ₁Modeling effect synoptic diagram

Fig. 6 is the checking of checking well X2(modeling well X1) modelling verification resultant effect synoptic diagram

Fig. 7 is the hollow husky three sections pyrolysis S of BN ₁Three fens property feature synoptic diagram

Fig. 8 is three fens property of the hollow husky three sections oil content chloroform bitumens of BN " A " feature synoptic diagram

Fig. 9 is mud shale series of strata criteria for classifying synoptic diagram

Figure 10 is the hollow husky three times dark 1 mud shale series of strata connection of the new Bohai Sea of the inferior section Shao sieve 52-sieve 48-sieve 69-67-well profile synoptic diagram of BHW basin BN

Figure 11 is the equivalent synoptic diagram of the hollow Es3x A of BN coefficient of restitution

Figure 12 is the hollow Es3x S of BN ₁The equivalent synoptic diagram of coefficient of restitution

Figure 13 is the equivalent synoptic diagram of chloroform bitumen resource intensity in the hollow Es3x mud shale of the BN series of strata

Figure 14 is pyrolysis S in the hollow Es3x mud shale of the BN series of strata ₁The equivalent synoptic diagram of method resource intensity

Figure 15 uses pyrolysis S in the hollow husky three times inferior section mud shale series of strata of BN ₁The movable resource intensity synoptic diagram of method

Figure 16 is the hollow husky three times inferior section short lap resource intensity synoptic diagram of BN

Embodiment

The present invention proposes mud shale series of strata concept, and propose and set up the criteria for classifying of mud shale series of strata according to U.S.'s oil shale fuel Development characteristics and China mud shale deposition characteristics.The criteria for classifying of the mud shale series of strata that the present invention proposes is: mud shale series of strata The initial segment is the mud shale section that thickness surpasses 2m, the interior sandstone ratio of 10m is no more than 1/3rd arbitrarily downwards, and sandstone thickness is no more than the different exploitation of 2m(block continuously, this numerical value should be different), Ro in the mud shale (maturity of organic matter) is greater than 0.5%, TOC (organic carbon content) is greater than 0.5%, i.e. pyrolysis S ₁Greater than 0.5mgHC/g, concluding paragraph is the mud stone section, and mud shale and sandstone cumulative thickness are not less than 30m.

Selected well is carried out the mud shale series of strata divide, can guarantee that the oil shale fuel stock number of calculating is reasonable, reliable, mud shale oil comprises present fine and close oil gas, and intension all is to compose the oil gas that exists in mud shale and the short lap thereof.The thin excessively enforcement that is unfavorable for fracturing technique of short lap, short lap is blocked up then to be the conventional gas and oil category, sandstone in the mud shale series of strata (or carbonatite) interlayer is " dessert " in the shale oil-gas exploration and development often, and be in mud shale (TOC〉0.5%, Ro 0.5%) sandstone (being full of oil) in the scope more merits attention.Therefore, need earlier underground mud shale series of strata to be divided, have only when the oil shale fuel resource abundance reaches a certain amount of in the series of strata, just have commercial value, otherwise exploitation has little significance.

The present invention is oil shale fuel resource potential evaluation methods in a kind of mud shale series of strata, after proposing and setting up the mud shale series of strata, carries out the evaluation of oil shale fuel resource potential in these series of strata, may further comprise the steps (as shown in Figure 1):

1) basic data arrangement: original geologic data, logging trace data and the organic and inorganic geochemistry data in selected work area put in order.Wherein geologic data comprises well logging, layering, well location, rock density, crude oil property etc., the logging trace data comprise interval transit time, resistivity, natural gamma, neutron, density etc., and the organic and inorganic geochemistry data comprises Ro, TOC, chloroform bitumen " A ", pyrolysis S ₁, isotope mass spectrometry, saturated hydrocarbon gas chromatography, crude oil family component, total rock analysis etc.

2) select well: gridding is carried out by certain yardstick in well to be selected work area divide, such as with 5km * 5km, 10km * 10km(only as example, but be not limited thereto) carry out gridding, 1～3 mouthful of well of primary election in each grid, the principle of selecting well is to distribute comparatively uniformly in the work area, and core space well number is suitably encrypted.

Select the requirement of well:

1. drill the zone of interest position at least 50m there is no and bore to meet tomography;

2. logging data is arranged;

3. layering data;

4. has R at least ₂₅A kind of in a kind of and AC (interval transit time), CNL (neutron) in (2.5 meters apparent resistivities), RT (dark side direction resistivity), three kinds of resistivity logging curves of CRT (true formation resistivity), three kinds of logging traces of DEN (density).

3) organic geochemistry well logging modeling and checking: utilize the unconventional petroleum resources evaluation platform in the prior art to carry out organic nonuniformity modeling (with reference to figure 5), comprise pyrolysis S ₁, chloroform bitumen " A ", TOC, select for use the many wells of geochemistry data to verify then, guarantee that model is extrapolated to the precision of other wells more than 80%.

4) set up oil shale fuel grading evaluation scheme: according to " three fens property " feature of mud shale oil content, be that III level resource disperses resource, II level resource to show that poor efficiency resource and I level resource are to show enrichment resource, S to show by enrichment degree with the oil shale fuel resource division ₁Value is the III level when 0～0.5mgHC/g, S ₁Value is the II level when 0.5～2mgHC/g, S ₁It is the I level greater than 2mgHC/g the time.

5) criteria for classifying that proposes according to the present invention is set up a mud shale series of strata identification system, automatically identifies the mud shale series of strata of every mouthful of well by the identification system, sets up the identification system and may further comprise the steps:

1. to having selected the logging data of well to put in order in the work area, all lithology are returned into two classes: sandstone and mud stone, wherein mud stone comprises oolitic limestone, calcareous algae limestone, coal and is the rock etc. of principal ingredient name with mud stone or shale.Sandstone comprises peat, grouan, bentonitic clay, pierite and is the rock etc. of principal ingredient name with siltstone, sandstone, conglomerate.Continuous sandstone or continuous mud stone are merged into one section, will select the logging data of well to be processed into sandstone mud stone alternate form at last;

2. the continuous T OC, the pyrolysis S that utilize the organic nonuniformity method of well logging to calculate ₁, chloroform bitumen " A ", the TOC, the pyrolysis S that calculate above-mentioned corresponding well mud stone section ₁, chloroform bitumen " A " average;

3. according to terrestrial heat flow, underground temperature gradient, formation temperature, reservoir pressure, Ro, the Tmax data such as (maximum pyrolysis peak temperature) of basin/depression, make the degree of ripeness thermal evolution history of basin/depression, determine at last every mouthful of well Ro exactly greater than 0.5% degree of depth section;

4. with the above-mentioned data of handling well, comprise well location, the degree of depth, lithology, TOC average, pyrolysis S ₁Data such as average, chloroform bitumen " A " average form a mud shale series of strata identification system file, and the mud shale series of strata of every mouthful of well are identified, for next step classification stock number and movable stock number of calculating in the mud shale series of strata lays the foundation.

6) lighter hydrocarbons, the heavy hydrocarbon of oil content are proofreaied and correct: the lighter hydrocarbons school is the component hydrocarbon-generating dynamics method according to prior art, recovers the lighter hydrocarbons loss percentage and proofreaies and correct out the windage of lighter hydrocarbons.The heavy hydrocarbon school is the S that proposes according to Wang Anqiao (1987) _{1 heavy hydrocarbon}Correction coefficient is recovered.

7) oil shale fuel resource potential grading evaluation: the grade scale with volumetric method and step 4) is estimated.

1. the chloroform bitumen resource potential is estimated principle

What chloroform bitumen " A " reflected is the content of solvable organic matter in the deposition rock, recently represents with the percentage that accounts for rock quality usually.As the synthesis result of living hydrocarbon and the effect of row's hydrocarbon, from essence, the abundance that is actually remaining organic matter in the hydrocarbon source rock of chloroform bitumen " A " reflection, it is the index of a resid amount.Therefore, to estimate the residual oil mass of hydrocarbon source rock comparatively suitable for the index of using chloroform bitumen " A ".

Carry out the formula that is calculated as follows of mud shale oil mass by original chloroform bitumen " A ":

$\begin{array}{cc}Q=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{S}_i\×H_i\×{\mathrm{\ρ}}_i\×k_a\×A_i& (i=\mathrm{1,2},\·\·\·n)\end{array}---\left(1\right)$

In the formula, Q is residual oil mass, and unit is 10 ⁴T/km ²S _iBe the hydrocarbon source rock area after the gridding, unit is 1km ²Because thickness and abundance of organic matter, type and the degree of ripeness of same source bed exist obvious variation in the plane, for improving evaluation precision, hydrocarbon source rock distributive province, work area is divided in the plane the grid regions of several 1km * 1km, (be mud shale oil mass intensity, unit is 10 to calculate the stock number of each grid regions respectively ⁴T/km ²), the summation that adds up then can obtain oil shale fuel stock number in the mud shale series of strata of work area; H _iBe the thickness of well logging resolution, be generally 0.125m; ρ _iBe hydrocarbon source rock density, unit is g/cm ³, according to well-log information this to get hydrocarbon source rock density be 2.53g/cm ³A _iBe that organic log response model calculates chloroform bitumen " A ", it has represented the value of average original chloroform bitumen " A " of the hydrocarbon source rock of well logging resolution thickness (0.125m); k _aBe the lighter hydrocarbons compensation correction coefficient of chloroform bitumen " A ", can obtain dimensionless by correction coefficient table; N is the number of the chloroform bitumen " A " of the individual well that calculates of organic log response model, generally along with buried depth increases, just can calculate the value of chloroform bitumen " A " every a well logging resolution (0.125m), because the zone thickness difference of various objectives layer, so also difference of n value.

In sum, can calculate the unit area (1km of work area individual well according to formula (1) ²) residual oil mass (10 ⁴T), same quadrat method can obtain the residual oil mass of other wells of work area successively, like this, makes the well of the mud shale oil mass of calculating cover the distribution plane of work area hydrocarbon source rock as far as possible.For improving evaluation precision, utilizing Surfer identification system that data are carried out gridding handles, hydrocarbon source rock distributive province, work area is divided in the plane the grid regions of several 1km * 1km, to the summation that adds up of the stock number intensity (resid amount intensity) of each grid regions, can obtain work area oil residues total amount (10 ⁴T).

2. pyrolysis S ₁The method resource potential is estimated principle

Rock pyrolysis data pyrolysis S ₁Be on free state (mgHC/g rock), be rock in the pyrolysis temperature-rise period before 300 ℃ thermal evaporation come out, be the hydrocarbon product in the source rock (rock) Already in.Pyrolysis S ₁Also can be used as the index of weighing resid amount, its numerical value has reflected remaining organic abundance in the hydrocarbon source rock.

Identical with principle, the method for using original chloroform bitumen " A " calculating mud shale oil mass, original pyrolysis S ₁The formula that calculates the mud shale oil mass is as follows:

$\begin{array}{c}Q=\underset{i=1}{\overset{n}{\mathrm{\&Sigma;}}}{10}^{-1}\&times;S_i\&times;H_i\&times;{\mathrm{\&rho;}}_i\&times;k_{{\mathrm{<figure-callout\; id="1"\; label="s"\; filenames="HDA00003315938100051.png,HDA00003315938100052.png"\; state="\{\{state\}\}">s</figure-callout>}}_1}\&times;S_{1_i}\end{array}(i=\mathrm{1,2},\&CenterDot;\&CenterDot;\&CenterDot;n)---\left(2\right)$

In the formula, 10 ^-1Be unit conversion, dimensionless; Q is residual oil mass, and unit is 10 ⁴T/km ²S _iBe the hydrocarbon source rock area after the gridding, unit is 1km ²H _iBe the thickness of well logging resolution, be generally 0.125m; ρ _iBe hydrocarbon source rock density, unit is g/cm ³, according to well-log information this to get hydrocarbon source rock density be 2.53g/cm ³

Be that organic log response model calculates

, it has represented hydrocarbon source rock on average original of well logging resolution thickness (0.125m)

Value; Be pyrolysis S ₁Correction coefficient can obtain by correction coefficient table; N is the individual well that calculates of organic log response model Number.

Use rock pyrolysis data pyrolysis S ₁The method that obtains work area mud shale oil mass is identical, repeats no more here.

8) by geochemistry data statistical study graphical method, determine the saturated absorption oil mass of mud shale, but and further calculate the oil shale fuel momentum.

Oil shale fuel is to be stored in big section to be rich in organic mud shale and the liquid hydrocarbon in the inorganic and organic hole in the interbedded sand (wherein part is adsorbed on organic matter, pyrite and rock particles surface), mode is deposited in the tax of oil shale fuel and reservoir space obviously is different from the conventional gas and oil Tibetan, be rich in organic mud shale and contain a large amount of oil generation matrix kerogens and big molecule hydrocarbon, adsorbed oil is very capable.The mud shale oil generating quantity satisfied self absorption residual hydrocarbons amount, water-soluble residual hydrocarbons amount, kapillary or storage pore shutoff residual quantity etc. various forms of residual after, unnecessary oil is discharged.Use the oil mass in the mud shale that current logging evaluation technique computes comes out, formed by ADSORPTION STATE oil mass and free state oil mass two parts.For the oil shale fuel exploitation, the free state oil mass has more actual value.Accordingly, the corresponding residual oil mass of oil extraction thresholding is saturated absorption oil mass as can be known, and we adopt saturated adsorption coefficient (ratio of saturated absorption oil mass and content of organic matter TOC) to come the height of quantitatively characterizing mud shale adsorptive power.But the oil shale fuel momentum can be definite like this:

But the oil shale fuel momentum=saturated absorption oil mass of original place oil shale fuel stock number-mud shale

Wherein oil shale fuel stock number in original place can calculate by the organic nonuniformity model of logging evaluation, and key is determining of the saturated absorption oil mass of mud shale, and it is to determine by geochemistry data statistical study graphical method, at first by chloroform bitumen " A ", S ₁+ S ₂, TOC and the degree of depth relation determine the oil extraction thresholding according to graphical method, secondly determine Q in conjunction with the oil extraction thresholding _{The saturated absorption oil mass of mud shale}, calculate Q by material balance method at last _{But oil shale fuel momentum}:

Q _{But oil shale fuel momentum}=Q _{Original place oil shale fuel stock number}-Q _{The saturated absorption oil mass of mud shale}

=V×ρ×(A%×K _A-TOC ^*×K×K _TOC) （3）

=Q _{Original place oil shale fuel stock number}-V * ρ * TOC ^** K * K _TOC

In the formula: ρ is mud shale density, and unit is _g/ cm ³V is the volume of mud shale, and unit is m ³TOC ^*Be the TOC average of three rank shale, unit is %, I level S ₁Greater than the 2.0mgHC/g rock, II level S ₁Between the 0.5-2.0mgHC/g rock, III level S ₁Less than the 0.5mgHC/g rock; K _ABe the coefficient of restitution of chloroform bitumen " A ", dimensionless; K _TOCBe the coefficient of restitution of organic carbon TOC, dimensionless; K is saturated adsorption coefficient, dimensionless.

Use existing logging evaluation method, be not difficult to calculate dispersion resource, poor efficiency resource and enrichment resource in the mud shale.But owing to be subjected to sedimentary environment, mineral composition and the wherein influence of organic abundance, type, degree of ripeness and expulsive efficiency, oleaginousness in the different mud shales has evident difference, and because mud shale has densification, low hole, especially hypotonic feature, the exploitation difficulty of oil shale fuel is bigger.For exploitation, but what of oil shale fuel momentum then are directly to determine the research work area whether to be worth investment, drilling well or carry out drilling well, pressure break in which piece, which degree of depth.But the present invention determines method according to oil shale fuel momentum evaluation method and saturated absorption oil mass that material balance principle proposes, has fundamentally solved the evaluation difficult problem of oil shale fuel mobility.

9) utilize " hole-oil saturation " method to calculate the oil shale fuel stock number of sandstone short lap in the mud shale series of strata, finish appraisal procedure.

In the sandstone short lap oil shale fuel resource potential estimate principle and mud shale resource intensity to ask for principle identical, utilize the stock number of volumetric method calculating sandstone short lap equally, principle is as follows, sees formula (4):

$Q=\underset{i=1}{\overset{n}{\mathrm{\&Sigma;}}}{h}_i\&CenterDot;S_i\&CenterDot;{\mathrm{\&phi;}}_i\&CenterDot;S_{\mathrm{oi}}\&CenterDot;{\mathrm{\&rho;}}_{\mathrm{oi}}(i=\mathrm{1,2},\&CenterDot;\&CenterDot;\&CenterDot;n)---\left(4\right)$

In the formula, Q is sandstone short lap oleaginousness, and unit is 10 ⁴T/km ²N is the number of plies of sandstone short lap; h _iSandstone short lap thickness, unit is m; S _iBe the source rock area after the gridding, unit is 1km ²φ _iBe factor of porosity, unit is %; S _OiBe oil saturation, unit is %; ρ _OiBe oil density, unit is _g/ cm ³

Wherein, oil saturation S _OiCan be obtained by following formula, see formula (5), (6)

S _Oi=1-S _w （5）

$\frac{1}{R_t}={[{V_{\mathrm{sh}}}^{1-\frac{{\mathrm{<figure-callout\; id="2"\; label="V\; sh"\; filenames="HDA00003315938100022.png,HDA00003315938100032.png"\; state="\{\{state\}\}">V</figure-callout>}}_{\mathrm{sh}}}{2}}\&CenterDot;\sqrt{\frac{1}{I_{\mathrm{sh}}}}+\frac{{\mathrm{\&phi;}}^{\frac{m}{2}}}{\sqrt{a\&CenterDot;R_W}}]}^2\&CenterDot;{S_W}^n---\left(6\right)$

And shaliness index is calculated by following formula, sees formula (7):

Shale index is obtained by following formula, sees formula (8):

$V_{\mathrm{sh}}=\frac{2^{\mathrm{GcuR}\&CenterDot;I_{\mathrm{sh}}}-1}{2^{\mathrm{GcuR}}-1}---\left(8\right)$

Calculate factor of porosity by sound wave, see formula (9), (10):

1·Δt=φ·Δt _f+V _sh·Δt _sh+(1-φ-V _sh)Δt _ma （9）

$\mathrm{\&phi;}=\frac{\mathrm{\&Delta;t}-V_{\mathrm{sh}}\&CenterDot;{\mathrm{\&Delta;t}}_{\mathrm{sh}}-{\mathrm{\&Delta;t}}_{\mathrm{ma}}+V_{\mathrm{sh}}\&CenterDot;{\mathrm{\&Delta;t}}_{\mathrm{ma}}}{{\mathrm{\&Delta;t}}_f-{\mathrm{\&Delta;t}}_{\mathrm{ma}}}---\left(10\right)$

In the formula, Sw is water saturation, and unit is %; A is the coefficient relevant with lithology, value between 0.6 to 1.2, dimensionless; M is porosity exponent, dimensionless; N is saturation exponent, dimensionless; I _ShBe shaliness index, dimensionless; V _ShBe shale index, unit is %; R _WBe resistivity of water, unit is Ω m; R _tBe formation resistivity, unit is Ω m; GR _Max, GR _Min, GR _SurveyBe maximum, minimum and the measured value of well logging GR value; GcuR is experience factor, and tertiary stratum is 3.7, and old stratum is 2; Δ t is interval transit time, and unit is μ s/ft; Δ t _fBe the fluid interval transit time, unit is μ s/ft; Δ t _ShBe the shale interval transit time, unit is μ s/ft; Δ t _MaBe the skeleton interval transit time, unit is μ s/ft.

Through type (4)～(10) can calculate the oil shale fuel stock number of sandstone short lap in the mud shale series of strata, finish whole appraisal procedure.

Enumerate specific embodiment below, the present invention there is further understanding.

Embodiment: be example with the husky three times inferior sections of the hollow SHAHEJIE FORMATION of BHW basin BN, use oil shale fuel resource potential evaluation method in the mud shale series of strata, but calculate oil shale fuel momentum and classification stock number in the hollow husky three times inferior section mud shale series of strata of BN.Concrete steps are:

1) basic data arrangement: collect hollow geologic data, logging trace data, the organic and inorganic geochemistry data of BN, put out the hollow husky three times inferior sections of BN in order and extend above-mentioned Various types of data in the 50m scope up and down, for the proposition of selecting well, the organic nonuniformity evaluation of well logging, oil shale fuel grading evaluation standard, mud shale series of strata concept and standard and determine, work such as the lighter hydrocarbons of oil content and heavy hydrocarbon correction prepare.

2) select well: in the hollow husky three times inferior sections of BN, select well, with work area 10km * 10km gridding (according to the hollow actual work area of BN situation), first 1～3 mouthful of well in each grid.

Select the requirement of well:

1. drill the hollow husky three times inferior sections of BN at least 50m there is no to bore and meet tomography;

2. logging data is arranged;

3. the layering data is arranged;

4. has R at least ₂₅, a kind of and AC, CNL and three kinds of well loggings of DEN among the RT, CRT three resistivity logging curves

A kind of in the curve.

Based on the above well condition of selecting, in the hollow husky three times inferior sections of BN, select 28 mouthfuls of wells (as shown in Figure 2), as can be seen from the figure, selected well location is according to available data control plane spread to greatest extent, and suitably encrypt at core space well number, for next step organic nonuniformity evaluation of logging well is prepared.

3) organic geochemistry well logging modeling and checking:

1. model is set up

At the research of work area oil shale fuel, the representative well many to geochemistry data, namely X1, X2 well carry out the cluster sampling analysis, and organic nonuniformity evaluation is conducive to log well.

Utilize R ₂₅Or RT and curves such as AC, CNL or DEN logging trace data make up, and set up the organic nonuniformity of well logging (TOC, chloroform bitumen " A " and S at husky three hypomeres of X1 well (2911-3140m) ₁) model.

This modeling is primarily aimed at mud shale to be carried out, and as shown in Figure 3, the modeling effect is better, as can be seen from the figure, and R ₂₅Measured value S is calculated in the well logging of-CNL mud shale ₁With calculated value S ₁The match relation is better, and the Trendline slope shows that close to 1 calculated value and measured value are comparatively approaching.As shown in Figure 4, R as can be seen from this figure ₂₅-CNL mud shale measured value and amplitude difference match relation better, related coefficient square all greater than 0.72(facies relationship number average greater than more than 0.85) and have tangible linear relationship, converge near the straight line, show that model credibility is higher.X1 well R ₂₅-CNLS ₁Modeling effect (as shown in Figure 5), as can be seen from this figure, calculated value and measured value relation are coincide fine continuously.

This adopts multiple difference (R in BHW basin BN is hollow ₂₅-AC, R ₂₅-CNL, RT-AC, RT-CNL, R ₂₅-DEN, RT-DEN) logging trace combination,, find out from institute's established model calculated value and measured value graph of a relation that it is better that measured value and calculated value match concern, can promote at local area in order to can suitably choose when using.Set up R in the hollow husky three times inferior sections of BHW basin BN ₂₅-CNL, R ₂₅The pyrolysis S of four kinds of curve combinations such as-AC, RT-AC, RT-CNL ₁The logging evaluation model, the parameter following (as shown in table 1) of different logging suite combinations and abundance of organic matter relational model:

The organic nonuniformity model parameter of table 1 logging evaluation

2. modelling verification

Be extrapolated to confidence level the unknown of other wells after model is set up, need verify model credibility.Eyeball (pyrolysis S is chosen in this modelling verification ₁, chloroform bitumen " A ", TOC) more X2 well, utilize the RT-AC model of mud shale to calculate X2 well continuous T OC value, calculated value and measured value are compared.Have only two mouthfuls of sealed coring well in the target area, wherein come modeling with the X1 well, verify with left unique a bite well X2, the checking of the modeling well of checking well X2(X1) modelling verification resultant effect fine (as shown in Figure 6), as can be seen from the figure, measured value and calculated value are very approaching, illustrate that model can be extrapolated to other wells.

3. model is used

This research is chosen 28 mouthfuls of wells according to well-log information abundance, drilling depth and well location planar distribution situation and is logged well, and estimates organic nonuniformity model and uses, and calculates continuous T OC, the pyrolysis S of every mouthful of well ₁And chloroform bitumen " A ", for next step mud shale series of strata are divided and mud shale oil resource grading evaluation is prepared.

4) set up oil shale fuel grading evaluation scheme: according to " three fens property " feature of mud shale oil content, be that the III level disperses resource, II level to represent that poor efficiency resource and I level are with expression enrichment resource with expression by enrichment degree with the oil shale fuel resource division.

The hollow husky three sections pyrolysis S of BN ₁And chloroform bitumen " A " three fens property feature such as Fig. 7, shown in Figure 8, as can be seen from the figure, chloroform bitumen " A " and pyrolysis S ₁Increase with TOC, the process of increase has three fens property, is divided into stablizing low value section, ascent stage and stable high value section.Accordingly, mark TOC and chloroform bitumen " A " and TOC and pyrolysis S ₁The scatter diagram enveloping outer enclosure, in three sections boundarys TOC boundary line of drawing, as the TOC grade scale of oil shale fuel classification resource evaluation, stablizing the low value section is that III level resource is to show dispersion (invalid) resource, ascent stage is II level resource to show the poor efficiency resource, and stable high value section is that I level resource is to show the enrichment resource.To stablize low value section, ascent stage boundary and the corresponding chloroform bitumen of ascent stage mid point " A " and pyrolysis S ₁Be the boundary, the chloroform bitumen that draws " A " and S ₁The boundary line is as chloroform bitumen " A " and pyrolysis S ₁Grade scale, this grading evaluation standard is with pyrolysis S ₁" three fens property " be as the criterion.Grading evaluation standard following (as shown in table 2):

Table 2BN-three times inferior section mud shale oil content level evaluation criterions of sand

5) criteria for classifying that proposes according to the present invention is set up a mud shale series of strata identification system, automatically identifies the mud shale series of strata of every mouthful of well by the identification system, and hollow mud shale series of strata identification may further comprise the steps to BHW basin BN:

1. to having selected the logging data of well to put in order in the hollow husky three times inferior sections of BN, all lithology are returned into two classes: sandstone and mud stone, wherein mud stone comprises oolitic limestone, calcareous algae limestone, coal and is the rock etc. of principal ingredient name with mud stone or shale; Sandstone comprises peat, grouan, bentonitic clay, pierite and is the rock etc. of principal ingredient name with siltstone, sandstone, conglomerate.Continuous sandstone or mud stone section are merged into one section, will select the logging data of well to be processed into sandstone, mud stone alternate form at last;

2. the continuous T OC, the pyrolysis S that utilize the organic nonuniformity model of well logging to calculate ₁, chloroform bitumen " A ", the TOC, the pyrolysis S that calculate above-mentioned corresponding well mud stone section ₁, chloroform bitumen " A " average;

3. according to the hollow terrestrial heat flow of BHW basin/BN, underground temperature gradient, formation temperature, reservoir pressure, data such as degree of ripeness Ro, Tmax, make the hollow degree of ripeness thermal evolution history of BHW basin/BN, determine at last every mouthful of well Ro exactly greater than 0.5% degree of depth section;

4. with the above-mentioned data of handling well, comprise well location, the degree of depth, lithology, TOC average, pyrolysis S ₁Columns certificates such as average, chloroform bitumen " A " average form an identification system by the requirement of mud shale series of strata identification system.With mud shale series of strata identification system the mud shale series of strata of every mouthful of well of the hollow husky three times inferior sections of BN are identified.Set up the dark 1 mud shale series of strata connection of the new Bohai Sea of Shao sieve 52-sieve 48-sieve 69-67-well profile (as shown in figure 10), these upright three connection well profiles (this shows a connection well profile) of in the hollow husky three times inferior section mud shale series of strata of BHW basin BN, building together.Can determine the locus of husky three times inferior section mud shale series of strata can tentatively obtain conclusion by three connection well profiles, the mud shale series of strata in the hollow husky three times inferior sections of BN are more stable enrichment positions, oil shale fuel space.

6) lighter hydrocarbons of oil content, heavy hydrocarbon are proofreaied and correct: it is to recover the windage that the lighter hydrocarbons loss percentage is proofreaied and correct out lighter hydrocarbons according to component hydrocarbon-generating dynamics method that lighter hydrocarbons are proofreaied and correct, and it is the pyrolysis S that proposes according to Wang Anqiao (1987) that heavy hydrocarbon is proofreaied and correct ₁The heavy hydrocarbon correction coefficient is recovered.

1. the lighter hydrocarbons of chloroform bitumen " A " recover

Chloroform bitumen " A " loses light hydrocarbon component (C in extractive process _14-), when utilizing chloroform bitumen to calculate the oil shale fuel stock number, need the light hydrocarbon component of loss be recovered.

It is as follows that chloroform bitumen " A " recovers formula:

" A " _Original=" A "+K _{The A-coefficient}* " A " * C _{Stable hydrocarbon}+ _{Aromatic hydrocarbon}(11)

In the formula, " A " _OriginalBe to recover back chloroform bitumen " A ", unit is %; " A " is actual measurement chloroform bitumen " A ", and unit is %; C _{Stable hydrocarbon+aromatic hydrocarbon}Be stable hydrocarbon and aromatic hydrocarbon proportion in actual measurement chloroform bitumen " A " (different oil fields are as the criterion with measured data separately, can be definite value, also can be the function of degree of ripeness), unit is %; K _{The A-coefficient}Be chloroform bitumen " A " lighter hydrocarbons correction coefficient, dimensionless is the C that calculates in the component dynamics _6-13With C ₁₃₊Ratio, C _6-13/ C ₁₃₊Value should be noted that the C of dissimilar organic matters as following table (as shown in table 3) _6-13/ C ₁₃₊Chart has than big difference, and this stock number is calculated chloroform bitumen " A " the lighter hydrocarbons correction coefficient of the chart that utilizes (1997) such as Behar.

Table 3 WenZK16Jing, Behar laboratory sample (I type kerogen) different evolution stages C _6-13/ C ₁₃₊Tables of data

Can get different depth chloroform bitumen " A " coefficient of restitution by chloroform bitumen " A " lighter hydrocarbons coefficient of restitution chart, and then draw chloroform bitumen " A " coefficient of restitution equal-value maps in husky three times inferior section mud shale series of strata, the hollow Es3x A of BN coefficient of restitution equal-value map as shown in figure 11, can recover the stock number of having calculated more accurately in conjunction with chloroform bitumen " A " coefficient of restitution equal-value map, calculate underground real resource amount, original place oil shale fuel stock number.

2. pyrolysis S ₁Recovery

Pyrolysis S ₁Drawn by the Rock-Eval analysis, the hydrocarbon that volatilizes when being no more than 300 ℃ for the rock sample heating, but represent extracting free hydrocarbon content in the rock, also be solvable hydrocarbon, be C basically _7-33Hydrocarbon.Yet carry out the used sample of pyrolysis analysis and experienced remaining sample for a long time, the lighter hydrocarbons (C in the rock core often _6-13Or C _6-14) existing more loss, the pyrolysis S that adopts pyrolysis analysis to obtain ₁Can not represent pyrolysis S in the subsurface rock ₁Original content.Simultaneously, in the past studies show that the cracking hydrocarbon S that pyrolysis analysis obtains ₂In the contribution of the solvable hydrocarbon of part is still arranged.Therefore, also need cracking hydrocarbon S ₂In solvable hydrocarbon recover.Pyrolysis S ₁The recovery formula is as follows:

S _{1 is original}=S ₁+ S _{1 heavy hydrocarbon}+ S _{1 lighter hydrocarbons}(12)

In the formula, S _{1 is original}Be the rock Central Plains molten hydrocarbon amount that begins, unit is mg/g; S ₁Be pyrolysis analysis method actual measurement pyrolysis S ₁, unit is mg/g; S _{1 heavy hydrocarbon}Be cracking hydrocarbon S ₂In solvable hydrocarbon, unit is mg/g; S _{1 lighter hydrocarbons}The lighter hydrocarbons of loss when being the core sample placement, unit is mg/g.

S _{1 heavy hydrocarbon}Restoration methods: get source rock sample, do pyrolysis test pyrolysis S ₁, cracking hydrocarbon S _2,Other gets same source rock sample and carries out chloroform bitumen " A " extracting, to rock sample after the extracting carry out pyrolysis test S _2', cracking hydrocarbon S then ₂With S _2'Difference (Δ S ₂) be and enter cracking hydrocarbon S ₂In the heavy component free hydrocarbon, dissimilar mud shale Δ S ₂/ S ₁Be worth following (as shown in table 4):

Table 4S _{1 heavy hydrocarbon}Coefficient of restitution table (according to Wang Anqiao, 1987)

S _{1 heavy hydrocarbon}The recovery formula is as follows:

S _{1 heavy hydrocarbon}=S ₁* △ S ₂/ S ₁(13)

In the formula, S _{1 heavy hydrocarbon}Be cracking hydrocarbon S ₂In solvable hydrocarbon, unit is mg/g; S ₁Be pyrolysis analysis method actual measurement pyrolysis S ₁, unit is mg/g; △ S ₂/ S ₁Be S _{1 heavy hydrocarbon}Coefficient of restitution, dimensionless.

S _{1 lighter hydrocarbons}Restoration methods: its each component is carried out component kinetic model parameter calibration according to source rock PY-GC experimental result, and use according to the geologic condition heating rate, calculate each component proportion in the hydro carbons that the hydrocarbon source rock different evolution stages generates, think that the composition of hydrocarbon is approximate identical in living hydrocarbon, row hydrocarbon and the residual hydrocarbon, can obtain that lighter hydrocarbons account for heavy hydrocarbon (S in the residual hydrocarbon ₁+ Δ S ₂) ratio.Δ S wherein _{1 lighter hydrocarbons}Can obtain according to following table (as shown in table 5).Therefore, the coefficient of restitution (lighter hydrocarbons/S of lighter hydrocarbons in the mud shale in the time of can drawing different evolution degree according to the component dynamic method ₁), i.e. K _{Lighter hydrocarbons}K when table 5 provides the organic different evolution degree of Jiyang depression I type _{Lighter hydrocarbons}In addition, it should be noted that the S of dissimilar organic matters _{1 lighter hydrocarbons}The coefficient of restitution chart has than big difference.

Table 5 WenZK16Jing (I type kerogen) different evolution stages K _{Lighter hydrocarbons}Table

S _{1 lighter hydrocarbons}The recovery formula is as follows:

S _{1 lighter hydrocarbons}=S ₁* K _{Lighter hydrocarbons}(14)

In the formula, S _{1 lighter hydrocarbons}The lighter hydrocarbons of loss when being the core sample placement, unit is mg/g; S ₁Be pyrolysis analysis method actual measurement pyrolysis S ₁, unit is mg/g; K _{Lighter hydrocarbons}Be the coefficient of restitution of lighter hydrocarbons in the mud shale, dimensionless.

By S _{1 lighter hydrocarbons}Coefficient of restitution chart and S _{1 heavy hydrocarbon}The coefficient of restitution table can get different depth pyrolysis S ₁Coefficient of restitution, and then draw pyrolysis S in husky three times inferior section mud shale series of strata ₁The coefficient of restitution equal-value map.In conjunction with pyrolysis S ₁Coefficient of restitution equal-value map (as shown in figure 12) can be used pyrolysis S more accurately ₁Method is calculated, is recovered stock number.

7) oil shale fuel resource potential grading evaluation: calculate the interior classification stock number of mud shale series of strata with volumetric method and grading evaluation scheme.

In conjunction with the computing method in the embodiment step 7), with oil shale fuel total resources (as shown in figure 13) in oil shale fuel classification stock number and the series of strata in the hollow husky three times inferior section mud shale series of strata of chloroform bitumen calculating BN.As can be seen from the figure to concentrate on BN hollow in the high value of chloroform bitumen " A " mud shale oil resource intensity districts in husky three times inferior section mud shale series of strata, and mxm. can reach 2846.12 * 10 ⁴T/km ², the resource intensity value of most of well is greater than 500 * 10 in the work area ⁴T/km ², with chloroform bitumen calculate the classification stock number of husky three times inferior sections see the following form (as shown in table 6).

Pyrolysis S in the hollow Es3x mud shale of the BN series of strata ₁The method resource intensity as shown in figure 14, as can be seen from the figure S in husky three times inferior section mud shale series of strata ₁The high value of mud shale oil resource intensity district mainly is distributed near dark 5 wells of the hollow Bohai Sea of BN, and mxm. can reach 2590.82 * 10 ⁴T/km ², the intensity level of most of area is greater than 450 * 10 in the work area ⁴T/km ², use pyrolysis S ₁Method calculate the classification stock number of husky three times inferior sections see the following form (as shown in table 5).

8) determine the saturated absorption oil mass of mud shale by geochemistry data statistical study graphical method, but and further calculate mud shale series of strata oil shale fuel momentum;

At first by being the hollow chloroform bitumen of BN " A ", S ₁+ S ₂, TOC and the degree of depth relation determine the oil extraction thresholding according to graphical method, secondly determine Q in conjunction with the oil extraction thresholding _{The saturated absorption oil mass of mud shale}, calculate Q by above-mentioned volumetric method (seeing formula 1, formula 2) _{The saturated absorption oil mass of mud shale}, calculate Q by material balance method at last _{But oil shale fuel momentum}, see formula (3), chloroform bitumen " A "/TOC in BN is hollow ^*100 be 15, S ₁/ TOC ^*1000 is 75, but quantitative calculate original place oil shale fuel stock numbers and oil shale fuel momentum (as table 6) in husky three times inferior section mud shale series of strata, use pyrolysis S ₁Method gets movable resource intensity (as shown in figure 15) in the hollow husky three times inferior section mud shale series of strata of BN.

Movable stock number and classification stock number in the hollow husky three times inferior section mud shale series of strata of table 6BN

9) utilize " hole-oil saturation " method to calculate the oil shale fuel stock number of sandstone short lap in the mud shale series of strata.

In conjunction with embodiment 9) in concrete steps, " hole-oil saturation " method of utilization calculates the hollow husky three times inferior section short lap resource intensities (as shown in figure 16) of BN, as can be seen from the figure, the high value of the resource intensity district of husky three times interior sandstone short laps of inferior section mud shale series of strata mainly concentrated near adopted 109 wells and adopted 282 wells in BN was hollow, and high value can reach 97.39 * 10 ⁴T/km ², the intensity level of nearly half area is arranged less than 10 * 10 in the work area ⁴T/km ², sandstone short lap stock number is 2.3 hundred million tons in husky three times inferior section mud shale series of strata.

The various embodiments described above only are used for explanation the present invention; wherein the structure of each parts, connected mode and manufacture craft etc. all can change to some extent; every equivalents and improvement of carrying out on the basis of technical solution of the present invention all should do not got rid of outside protection scope of the present invention.

Claims (5)

1. oil shale fuel resource potential evaluation method in the mud shale series of strata may further comprise the steps:

1) arrangement basic data, described basic data comprises: geologic data, logging trace data, organic and inorganic geochemistry data;

2) in the work area, evenly select well, and encrypt at core space well number, according to selected yardstick with the gridding of described work area, 1～3 mouthful of described well of primary election in each grid;

3) organic geochemistry well logging modeling and checking: utilize unconventional petroleum resources to estimate platform selecting representative well wherein and carry out organic nonuniformity modeling, comprise pyrolysis S ₁, chloroform bitumen " A ", TOC, further select for use another representative well that the model of setting up is verified;

4) set up oil shale fuel grading evaluation standard: according to " three fens property " feature of mud shale oil content, be that the III level is disperseed resource with expression by enrichment degree with the oil shale fuel resource division, the II level is with expression poor efficiency resource, and the I level is with expression enrichment resource;

5) according to U.S.'s oil shale fuel Development characteristics and China's mud shale deposition characteristics, set up the criteria for classifying of mud shale series of strata, and according to the described criteria for classifying, write mud shale series of strata identification system, automatically identify the mud shale series of strata of every mouthful of well by described mud shale series of strata identification system;

6) lighter hydrocarbons, the heavy hydrocarbon of described oil shale fuel in the described mud shale series of strata are proofreaied and correct;

7) oil shale fuel resource potential grading evaluation: carry out oil shale fuel resource potential grading evaluation with the described grading evaluation scheme in volumetric method and the step 4);

8) determine the saturated absorption oil mass of described mud shale by geochemistry data statistical study graphical method, but and further calculate the oil shale fuel momentum;

9) utilize " hole-oil saturation " method to calculate the oil shale fuel stock number of sandstone short lap in the described mud shale series of strata, finish oil shale fuel resource potential grading evaluation.

2. oil shale fuel resource potential evaluation method in a kind of mud shale series of strata according to claim 1, it is characterized in that: described 5) in the step, the described criteria for classifying of mud shale series of strata is: mud shale series of strata The initial segment is the mud shale section that thickness surpasses 2m, the interior sandstone proportion of any 10m is less than 1/3 downwards, and sandstone thickness is less than 2m continuously, mud shale Ro greater than 0.5%, TOC greater than 0.5%, pyrolysis S ₁Greater than 0.5mgHC/g, concluding paragraph is the mud stone section, and mud shale and sandstone cumulative thickness are greater than 30m.

3. oil shale fuel resource potential evaluation method in a kind of mud shale series of strata according to claim 1 and 2, it is characterized in that: described 5) in the step, according to the described criteria for classifying, write described mud shale series of strata identification system, automatically identify the described mud shale series of strata of every mouthful of well by described mud shale series of strata identification system, specific as follows:

1. to having selected the logging data of well to put in order in the described work area;

2. utilize the organic nonuniformity model of well logging to calculate continuous described TOC, described pyrolysis S ₁, described chloroform bitumen " A ", obtain the described TOC of described mud stone section, described pyrolysis S ₁, described chloroform bitumen " A " mean value;

3. according to terrestrial heat flow, underground temperature gradient, formation temperature, reservoir pressure, the data such as described Ro, Tmax of basin/depression, make the degree of ripeness thermal evolution history of described basin/depression, determine every mouthful of described Ro of well at last greater than 0.5% degree of depth section;

4. according to described mud shale series of strata identification system require form with the above-mentioned data of handling well, the described data of handling well comprise well location, the degree of depth, lithology, described TOC average, described pyrolysis S ₁Average, described chloroform bitumen " A " average, write an identification system, with the described identification system of described mud shale series of strata the described mud shale series of strata of every mouthful of well are identified, for described mud shale series of strata internal classification stock number and movable RESERVE EVALUATION lay the foundation.

4. oil shale fuel resource potential evaluation method in a kind of mud shale series of strata according to claim 1 and 2 is characterized in that: described 8) in the step, determine the saturated absorption oil mass of described mud shale Q by described geochemistry data statistical study graphical method _{The saturated absorption oil mass of mud shale}But, and further utilize material balance principle to calculate described oil shale fuel momentum Q _{But oil shale fuel momentum}But, described oil shale fuel momentum Q _{But oil shale fuel momentum}Computing formula be:

Q _{But oil shale fuel momentum}=Q _{Original place oil shale fuel stock number}-Q _{The saturated absorption oil mass of mud shale}

=V×ρ×(A%×K _A-TOC ^*×K×K _TOC)

=Q _{Original place oil shale fuel stock number}-V * ρ * TOC ^** K * K _TOC

In the formula, ρ is mud shale density, and unit is g/cm ³V is the volume of mud shale, and unit is m ³TOC ^*Be the TOC average of three rank shale, unit is %, I level S ₁Greater than the 2.0mgHC/g rock, II level S ₁Between the 0.5-2.0mgHC/g rock, III level S ₁Less than the 0.5mgHC/g rock; K _ABe the coefficient of restitution of chloroform bitumen " A ", dimensionless; K _TOCBe the coefficient of restitution of organic carbon TOC, dimensionless; K is saturated adsorption coefficient, dimensionless.

5. oil shale fuel resource potential evaluation method in a kind of mud shale series of strata according to claim 3 is characterized in that: described 8) in the step, determine the saturated absorption oil mass of described mud shale Q by described geochemistry data statistical study graphical method _{The saturated absorption oil mass of mud shale}But, and further utilize material balance principle to calculate described oil shale fuel momentum Q _{But oil shale fuel momentum}But, described oil shale fuel momentum Q _{But oil shale fuel momentum}Computing formula be:

Q _{But oil shale fuel momentum}=Q _{Original place oil shale fuel stock number}-Q _{The saturated absorption oil mass of mud shale}

=V×ρ×(A%×K _A-TOC ^*×K×K _TOC)

=Q _{Original place oil shale fuel stock number}-V * ρ * TOC ^** K * K _TOC

In the formula, ρ is mud shale density, and unit is g/cm ³V is the volume of mud shale, and unit is m ³TOC ^*Be the TOC average of three rank shale, unit is %, I level S ₁Greater than the 2.0mgHC/g rock, II level S ₁Between the 0.5-2.0mgHC/g rock, III level S ₁Less than the 0.5mgHC/g rock; K _ABe the coefficient of restitution of chloroform bitumen " A ", dimensionless; K _TOCBe the coefficient of restitution of organic carbon TOC, dimensionless; K is saturated adsorption coefficient, dimensionless.

Images