CN110424956A

CN110424956A - Evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation

Info

Publication number: CN110424956A
Application number: CN201910618385.0A
Authority: CN
Inventors: 林会喜; 王圣柱; 汪誉新; 曲彦胜; 马雪洁; 于洪洲; 张奎华; 张关龙; 赵建成; 牛晓燕
Original assignee: China Petrochemical Corp; Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Current assignee: China Petrochemical Corp; Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Priority date: 2019-07-09
Filing date: 2019-07-09
Publication date: 2019-11-08
Anticipated expiration: 2039-07-09
Also published as: CN110424956B

Abstract

The present invention provides evaluation unit in a kind of shale oil Resources calculation and saves coefficient weights quantization assignment method, comprising: step 1, according to the practical geologic structure in research area and preservation condition feature, divides different evaluation units；Step 2, typical well or outcrop profile are chosen；Step 3, different evaluation unit, the abundance of organic matter of different lithology lithofacies and oil content parameters relationship plate are established；Step 4, different evaluation unit is established, point lithofacies save coefficient weights assignment quantitative model；Step 5, divide evaluation unit, divide Lithofacies Types amount of computational resources.Evaluation unit preservation coefficient weights quantization assignment method is more scientific in the shale oil Resources calculation, as a result more objective, solves the problems, such as that being directed to retrofit strategies shale oil RESERVE EVALUATION in the past is not allowed, and can provide technical support for the practical investigation and prospecting of shale oil.

Description

Evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation

Technical field

The present invention relates to oil field development technical field, evaluation unit in a kind of shale oil Resources calculation is especially related to It saves coefficient weights and quantifies assignment method.

Background technique

In recent years, the exploration and development of shale oil gas achieves huge success in north America region, has pushed global shale The upsurge of oil-gas exploration research, shale oil unconventional oil and gas resource exploration have been to be concerned by more and more people.Global shale oil gas Research be still in infancy at present.At oil-gas exploration initial stage, shale oil-gas exploration direction is selected in shale oil resource assessment It selects and is of great significance with play assessment.About conventional resource assessment, forefathers propose analogy method, statistic law, origin cause of formation method, basin A variety of research methods such as simulation, have obtained the extensive utilization of oil exploration personnel, and technical method maturation is perfect.Lu Shuanfang (2012), the scholars such as Xue Haitao (2016), Wang Xuefei (2013), Wang Min (2013) are directed to Song-liao basin Qingshankou group, Nenjiang group, Numerous studies have been carried out in the shale oil resource assessment such as Bonan Sag of Jiyang Depression SHAHEJIE FORMATION, achieve research achievement abundant, Experiment test establishes that mud shale is organic and inorganic heterogeneity Logging estimation model with conventional logging means in conjunction with；Based on rock Stone sample experiments and Organic Matter Hydrocarbon Generation dynamics research form rock oil content (adsorbed hydrocarbons S₁With chloroform bitumen " A ") heavy hydrocarbon, Lighter hydrocarbons restoration methods propose classification RESERVE EVALUATION Research Thinking, have effectively pushed domestic shale oil resource assessment to work, taken Obtained good exploration effects.

Shale oil resource assessment at present uses different meters according to the occurrence form (matrix type and sandwich type) of shale oil gas Calculation method.For matrix type shale oil stock number, rock pyrolysis adsorbed hydrocarbons S is proposed₁Method and rock extracting chloroform bitumen " A " Method；For sandwich type shale oil stock number, oil saturation method is proposed.At present about the main concentration of exploration research of shale oil The relatively stable sedimentary basin (area) for inheriting sexual development, such as Song-liao basin, Bohai gulf basin are being constructed, and it is strong to constructing The shale oil research of reconstructed basin is related to very few.The shale oil RESERVE EVALUATION proposed for stable type sedimentary basin (area) Method considers deficiency due to saving element, causes when applying to reconstructed basin (area) shale oil RESERVE EVALUATION, calculate As a result bigger than normal, there are apparent unworthiness.If carrying out exploration constituency according to current shale oil Resources calculation result to comment Valence will appear obvious deviation, or even the understanding that must be made mistake.Therefore it needs to consider tectonic reworking zone geologic structure from different places Complexity, the big geology characteristic of preservation condition superiority and inferiority difference explore the applicable shale oil Resources Assessment Method in Renewal and development type basin. Coefficient weights are saved we have invented evaluation unit in a kind of new shale oil Resources calculation thus and quantify assignment method, are solved Above technical problem.

Summary of the invention

The object of the present invention is to provide a kind of progress shale oil RESERVE EVALUATION is more scientific, the evaluation result obtained is more Evaluation unit saves coefficient weights quantization assignment method in objective reasonable shale oil Resources calculation.

The purpose of the present invention can be achieved by the following technical measures: evaluation unit saves system in shale oil Resources calculation Number weight quantifies assignment method, and evaluation unit preservation coefficient weights quantization assignment method includes: in the shale oil Resources calculation Step 1, according to the practical geologic structure in research area and preservation condition feature, different evaluation units is divided；Step 2, it chooses typical Drilling well or outcrop profile；Step 3, different evaluation unit, the abundance of organic matter of different lithology lithofacies and oil content parameter is established to close It is plate；Step 4, different evaluation unit is established, point lithofacies save coefficient weights assignment quantitative model；Step 5, divide evaluation single Member divides Lithofacies Types amount of computational resources.

The purpose of the present invention can be also achieved by the following technical measures:

In step 1, conduct a research area's tectonic interpretation first, layer distribution characteristics that have a definite purpose and main faults spread Feature, using target zone overlap of beds or pinching distribution, boundary fault as boundary condition, binding purpose layer buried depth of strata, region lid Layer, the analysis of stratum aqueous nature preservation condition, divide different evaluation units.

In step 2, basic demand is at least a bite typical case well or section in the every a kind of evaluation unit marked off, And drilling well or section have systematic abundance of organic matter TOC, pyrolysis adsorbed hydrocarbons S₁With chloroform bitumen " A " basic geochemical testing Data have representativeness of sample, at least to there is 20 or more data points.

Step 3 includes:

Step 3a carries out lithofacies division；

Step 3b, sub-unit, a point lithofacies establish abundance of organic matter and oil content parameters relationship plate.

In step 3a, based on rock core, appear, thin slice, the macro microscopic observation of scanning electron microscope and and total rock X diffraction analysis, according to Rock deposition construction, mineral constituent, pore structure characteristic, using " rock constituents-sedimentary structure-organic matter " lithofacies division side Case divides different Lithofacies Types.

In step 3b, the rock abundance of organic matter TOC of different lithofacies in each evaluation unit, pyrolysis adsorbed hydrocarbons S are counted respectively₁ Or chloroform bitumen " A " measured data, system comparative analysis is carried out to data, rejects individual abnormal points；According to research area's shale oil Lithofacies Types splitting scheme and research need, and establish the abundance of organic matter TOC and oil-containing of different evaluation unit, different Lithofacies Types Property S₁Or chloroform bitumen " A " scatterplot relationship plate.

Step 4 includes:

Step 4a, selection criteria area；

Step 4b builds different evaluation unit oil content comparison bridge；

Step 4c carries out evaluation unit and saves coefficient relative weighting assignment.

In step 4a, it is contemplated that the oil content of identical lithofacies rock caused by different evaluation unit preservation condition superiority and inferiority is different Difference, according to research the practical geological condition in area, choose preservation condition construction of stable the most superior area be standard regions, It saves coefficient relative weighting assignment and takes 1.0, other kinds of evaluation unit saves that coefficient relative weighting assignment is in contrast to be obtained.

In step 4b, it is assumed that the sample free hydrocarbon loss amount of the identical lithofacies of different evaluation unit and its initial oil rate positive Closing is that loss late is consistent, therefore, will not be generated to the loss amount before test analysis to the use of measured data later period after sample acquisition It influences；Using abundance of organic matter as bridge, the comparison of different evaluation unit rock sample oil content is realized.

In step 4c, according to abundance of organic matter TOC and oil content adsorbed hydrocarbons S₁Or chloroform bitumen " A " scatterplot data, it carries out Data quantization relationship analysis obtains the TOC and adsorbed hydrocarbons S of different Lithofacies Types₁Or the quantization of chloroform bitumen " A " between the two is public Formula:

C_{Oil base}=a*TOC^b (1)

C_{Oil folder}=c+d*TOC (2)

In formula: C_{Oil base}、C_{Oil folder}Respectively matrix type, sandwich type rock oil content parameter are pyrolyzed adsorbed hydrocarbons S₁, unit mg/ G, or extracting chloroform bitumen " A ", unit %；A, b, c, d are equation regression coefficient；TOC is rock abundance of organic matter, unit %.

Step 4c includes:

Step 1, single-point save coefficient with respect to assignment；According to the corresponding absorption of the identical TOC of the same lithofacies of different evaluation unit Hydrocarbon S₁Or chloroform bitumen " A " relative amount quantitative relation, the preservation coefficient quantization model of different evaluation unit is established, and then determine The preservation coefficient relative weighting assignment of each evaluation unit；According to the TOC of above-mentioned foundation and oil content quantitative relationship formula, by other The rock oil content quantity of the rock oil content of a certain abundance of organic matter of evaluation unit abundance of organic matter corresponding with standard regions is closed System determines the preservation coefficient relative weighting assignment of the corresponding lithofacies of corresponding evaluation unit；

Matrix type lithofacies save coefficient relative weighting assignment formula are as follows:

K_{Unit base i}=C_{Oil base i}/C_{Oil disjunction mark is quasi-}=a_i*TOC_{Base i}^b_i/a_{Disjunction mark}*TOC_{Disjunction mark}^b_{Disjunction mark} (3)

In formula: K_{Unit base i}Coefficient relative weighting assignment is saved for matrix type lithofacies evaluation unit i；C_{Oil base i}For matrix type lithofacies Evaluation unit i oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；C_{Oil disjunction mark is quasi-}For matrix type Lithofacies evaluation unit standard regions oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；a_i, b_i For matrix type lithofacies evaluation unit i equation regression coefficient；a_{Disjunction mark}, b_{Disjunction mark}For matrix type lithofacies evaluation unit standard regions equation regression system Number；TOC_{Base i}For matrix type lithofacies evaluation unit i abundance of organic matter, unit %；TOC_{Disjunction mark}For matrix type lithofacies evaluation unit standard Area's abundance of organic matter, unit %；I is the evaluation unit i marked off；

Sandwich type lithofacies save coefficient relative weighting assignment formula are as follows:

K_{Unit presss from both sides i}=C_{Oil presss from both sides i}/C_{Oil presss from both sides standard}=(C_i+d_i*TOC_{Press from both sides i})/(C_{Folder mark}+d_{Folder mark}*TOC_{Folder mark})(4)

In formula: K_{Unit presss from both sides i}Coefficient relative weighting assignment is saved for sandwich type lithofacies evaluation unit i；C_{Oil presss from both sides i}For sandwich type lithofacies Evaluation unit i oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；C_{Oil presss from both sides standard}For sandwich type Lithofacies evaluation unit standard regions oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；c_i, d_i For sandwich type lithofacies evaluation unit i equation regression coefficient；C_{Folder mark}, d_{Folder mark}For sandwich type lithofacies evaluation unit standard regions equation regression system Number；TOC_{Press from both sides i}For sandwich type lithofacies evaluation unit i abundance of organic matter, unit %；TOC_{Folder mark}For sandwich type lithofacies evaluation unit standard Area's abundance of organic matter, unit %；I is the evaluation unit i marked off；

Step 2, different evaluation unit save the comprehensive assignment of coefficient relative weighting；Had according to research area's difference Lithofacies Types Multiple single-points are saved system on the basis of determining that abundance of organic matter compares section and comparison interval by machine matter abundance distribution feature The opposite assignment of number seeks arithmetic mean of instantaneous value as different units and saves the comprehensive assignment of coefficient weights；Including following two step:

1. different lithofacies abundance of organic matter Characteristics of Distribution；Sample is surveyed in the identical lithofacies of analysis different evaluation unit On the basis of TOC distribution characteristics, reasonable TOC comparison section is chosen；Statistical sample point, point lithofacies divide organic matter section, draw Research area's abundance of organic matter distribution histogram, the identical lithofacies abundance of organic matter TOC distribution characteristics of comparative analysis different evaluation unit, Determine that different evaluation unit oil content TOC compares section；

Coefficient comprehensive weight assignment is saved 2. determining；According to the abundance of organic matter TOC characteristic distributions of different lithofacies, master is determined Want TOC distributed area (m, n]；The certain interval TOC inter is chosen, and then determines contrast sample's point and sample between evaluation unit This number Num takes each sample to count number and saves the guarantor that coefficient weights correspond to lithofacies with respect to assignment average value as corresponding evaluation unit Deposit coefficient comprehensive weight assignment；

Num=(n-m)/inter (5)

K_{Unit i base is comprehensive}=(K_{Unit i base 1}+K_{Unit i base 2}+…+K_{Unit i base j})/Num (6)

K_{Unit i folder is comprehensive}=(K_{Unit i folder 1}+K_{Unit i folder 2}+…+K_{Unit i presss from both sides j})/Num (7)

J=1,2,3 ..., Num

In formula: m left side low value between toc region, unit is %；N right side high level between toc region, unit is %；Inter is TOC comparative analysis value interval, unit is %；Num comparative analysis number of samples between different evaluation unit, unit are； K_{Unit i base j}Coefficient relative weighting assignment, decimal are saved for j-th of sample point of evaluation unit i matrix type lithofacies；K_{Unit i presss from both sides j}It is single for evaluation J-th of sample point of first i sandwich type lithofacies saves coefficient relative weighting assignment, decimal；K_{Unit i base is comprehensive}For evaluation unit i matrix type lithofacies Shale oil saves coefficient relative weighting assignment mean values, decimal；K_{Unit i folder is comprehensive}For the preservation of evaluation unit i sandwich type lithofacies shale oil Coefficient relative weighting assignment mean values, decimal.

In steps of 5, in rock abundance of organic matter TOC, pyrolysis free hydrocarbon S₁Lighter hydrocarbons compensation, heavy hydrocarbon recovery, chloroform bitumen The compensation of " A " lighter hydrocarbons and S₁, chloroform bitumen " A " saturation these key parameters of adsorption coefficient determine on the basis of, for different pages Pyrolysed hydrocarbon S is respectively adopted in shale oil occurrence form₁Method or chloroform bitumen "A" method and oil saturation method carry out the original money of shale oil The calculating of source amount according to different evaluation unit, the corresponding opposite preservation coefficient of different Lithofacies Types, determines page on this basis The stock number now of shale oil.

Evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation in the present invention, for construction Reconstructed basin geologic structure is complicated, the big feature of preservation condition difference, the Research Thinking of sub-unit evaluation is proposed, according to difference A large amount of sample actual measurements abundance of organic matter (TOC) of evaluation unit and oil content (S₁Or chloroform bitumen " A ") date comprision base On plinth, the preservation coefficient weights assignment quantization plate of different evaluation unit, different lithofacies is established, determines that the corresponding coefficient that saves is assigned It is worth weight.It is corresponding to save coefficient multiplied by corresponding evaluation unit on the basis of conventional method original shale oil Resources calculation Obtain stock number now.The evaluation method is more scientific, as a result more objective, solves and provided in the past for retrofit strategies shale oil The inaccurate problem of source amount evaluation, can provide technical support for the practical investigation and prospecting of shale oil.

It has the advantage that the 1, present invention due to considering reconstructed basin geologic structure difference big feature and preservation Influence of the condition superiority and inferiority to shale oil enrichment degree proposes the Research Thinking of sub-unit evaluation, according to actual measurement sample statistics, divides Evaluation unit, point lithofacies establish abundance of organic matter (TOC) and oiliness parameter (S respectively₁Or chloroform bitumen " A ") relationship plate, The quantization assignment for realizing the preservation coefficient relative weighting of different evaluation unit, avoids the subjectivity of artificial assignment, while gram The oil-containing sex differernce as caused by lithologic character and lithofacies difference is taken, using pyrolysis adsorbed hydrocarbons S₁It is determined respectively with chloroform bitumen "A" method The preservation coefficient of different units lithofacies, comprehensive study area measured data distribution characteristics, selection calibration comparison section are chosen synthesis and are commented Valence coefficient fundamentally solves and is transformed basin shale oil RESERVE EVALUATION using stable region shale oil Resources Assessment Method As a result inaccurate problem.2, the present invention is due to considering the oil content difference of different lithofacies, using sub-unit, divides the resources of lithofacies Thinking is evaluated, solves the problems, such as that previous methods resource assessment precision is low, evaluation precision is substantially increased, using the method for the present invention It is more scientific to carry out shale oil RESERVE EVALUATION, the evaluation result obtained is more objective rationally.

Detailed description of the invention

Fig. 1 is that evaluation unit saves the tool that coefficient weights quantify assignment method in shale oil Resources calculation of the invention The flow chart of body embodiment；

Southeast reed grass ditch group evaluation unit divides sectional view subject to Fig. 2；

Southeast reed grass ditch group evaluation unit divides plan view subject to Fig. 3；

Fig. 4 is depressed area matrix type lithofacies rock abundance of organic matter (TOC) and pyrolysis adsorbed hydrocarbons (S₁) relational graph；

Fig. 5 is structurally stable region matrix type lithofacies rock abundance of organic matter (TOC) and pyrolysis adsorbed hydrocarbons (S₁) relational graph；

Fig. 6 is to thrust transformation territorial matrix type lithofacies rock abundance of organic matter (TOC) and pyrolysis adsorbed hydrocarbons (S₁) relational graph；

Fig. 7 is outcrop area matrix type lithofacies rock abundance of organic matter (TOC) and pyrolysis adsorbed hydrocarbons (S₁) relational graph；

Fig. 8 is depressed area sandwich type lithofacies rock abundance of organic matter (TOC) and pyrolysis adsorbed hydrocarbons (S₁) relational graph；

Fig. 9 is structurally stable region sandwich type lithofacies rock abundance of organic matter (TOC) and pyrolysis adsorbed hydrocarbons (S₁) relational graph；

Figure 10 is to thrust transformation area's sandwich type lithofacies rock abundance of organic matter (TOC) and pyrolysis adsorbed hydrocarbons (S₁) relational graph；

Figure 11 is outcrop area sandwich type lithofacies rock abundance of organic matter (TOC) and pyrolysis adsorbed hydrocarbons (S₁) relational graph；

Figure 12 is different evaluation unit matrix type lithofacies rock abundance of organic matter (TOC) and chloroform bitumen " A " relational graph；

Figure 13 is different evaluation unit sandwich type lithofacies rock abundance of organic matter (TOC) and chloroform bitumen " A " relational graph；

Southeast reed grass ditch group matrix type lithofacies TOC distribution histogram subject to Figure 14；

Southeast reed grass ditch group sandwich type lithofacies TOC distribution histogram subject to Figure 15.

Specific embodiment

To enable above and other objects, features and advantages of the invention to be clearer and more comprehensible, preferably implementation is cited below particularly out Example, and cooperate shown in attached drawing, it is described in detail below.

Coefficient weights, which are saved, as shown in FIG. 1, FIG. 1 is evaluation unit in shale oil Resources calculation of the invention quantifies assignment The flow chart of method.

1) evaluation unit divides: according to the practical geologic structure in research area and preservation condition feature, dividing different evaluation lists Member；Evaluation unit classification, conduct a research area's tectonic interpretation first, layer distribution characteristics that have a definite purpose and main faults spread Feature, using target zone overlap of beds or pinching distribution, boundary fault as boundary condition, binding purpose layer buried depth of strata, region lid Layer, the analysis of stratum aqueous nature (water type, salinity) preservation condition, divide different evaluation units；

2) typical well (or outcrop profile) is chosen: basic demand be in the every a kind of evaluation unit marked off at least A bite typical case well (or section), and drilling well (or section) has systematic abundance of organic matter (TOC), pyrolysis adsorbed hydrocarbons (S₁) and chlorine Imitative pitch " A " basic geochemistry test data, has preferable representativeness of sample, generally at least to there is 20 or more data Point；

3) establish different evaluation unit, different lithology lithofacies (matrix type and sandwich type) abundance of organic matter (TOC) with contain Oil cut rate parameter (S₁Or chloroform bitumen " A ") relationship plate: according to research area's shale oil Lithofacies Types splitting scheme (shale phase, mud Lithofacies, sandstone phase, carbonate facies) and research need, establish the abundances of organic matter (TOC) of different evaluation unit point lithofacies and contain Oiliness (S₁Or chloroform bitumen " A ") relationship plate；

Evaluation unit, lithologic character and lithofacies type is divided to establish TOC and oil content parameter (adsorbed hydrocarbons S respectively₁Or chloroform bitumen " A ") Relationship plate.It is assumed that the rock sample of the comparable identical lithofacies of abundance of organic matter (TOC), in the original not influenced by tectonic reworking Its oil content should substantially quite under beginning state, and preservation condition superiority and inferiority caused by later structural is transformed makes different evaluation unit now Abundance of organic matter is suitable, the identical rock oil content of lithofacies has differences.In order to exclude the rock pore structures of different lithofacies because The influence of oil content caused by element and oil and gas loss otherness can divide the type of different lithofacies, respectively construct different evaluation unit, The TOC and (S of different lithofacies₁Or chloroform bitumen " A ") quantitative relationship model.Lithofacies Types are specifically included to divide and plate foundation:

1. lithofacies divide

In order to make the oil content of different evaluation unit sample that there is comparability, need according to research area's lithology development feelings Condition carries out lithofacies division, and on this basis, identical lithofacies compare.Method particularly includes: based on rock core, appears, thin slice, sweeps Retouch the macro microscopic observation of Electronic Speculum and and total rock X diffraction analysis, according to rock deposition construct (lamellar, stratiform, bulk), mineral constituent (quartz, feldspar, carbonate rock, clay mineral), pore structure characteristic, using " rock constituents-sedimentary structure-organic matter " lithofacies Splitting scheme divides different Lithofacies Types.

2. sub-unit, a point lithofacies establish abundance of organic matter and oil content parameters relationship plate

The rock abundance of organic matter (TOC) of different lithofacies in each evaluation unit, pyrolysis adsorbed hydrocarbons (S are counted respectively₁Or chloroform Pitch " A ") measured data, system comparative analysis is carried out to data, rejects individual abnormal points.According to research area's shale oil lithofacies class Type splitting scheme (shale phase, mud stone phase, sandstone phase, carbonate facies) and research need, and establish different evaluation unit, different rocks The abundance of organic matter (TOC) and oiliness (S of facies type₁Or chloroform bitumen " A ") scatterplot relationship plate.

4) different evaluation unit is established, point lithofacies save coefficient weights assignment quantitative model: with data such as originLab It analyzes software and carries out different evaluation unit, the abundance of organic matter (TOC) of lithofacies and pyrolysis adsorbed hydrocarbons (S₁Or chloroform bitumen " A ") survey The network analysis for trying data, obtains corresponding quantitative formula；

Determine the preservation coefficient relative weighting assignment of different evaluation unit, different lithofacies, comprising the following steps:

1. standard regions select

It, can in view of the difference of the oil content of identical lithofacies rock caused by different evaluation unit preservation condition superiority and inferiority difference According to the research practical geological condition in area, the area for choosing preservation condition construction of stable the most superior is standard regions, saves system Number relative weighting assignment take 1.0, and other kinds of evaluation unit saves that coefficient relative weighting assignment is in contrast to be obtained.

2. building different evaluation unit oil content comparison bridge

Rock abundance of organic matter (TOC) by preservation condition and later period slacking although also influenced, relative to pyrolysis Adsorbed hydrocarbons S₁For (or chloroform bitumen " A "), degree of susceptibility is relatively weak.It is assumed that the sample of the identical lithofacies of different evaluation unit Product free hydrocarbon loss amount and its initial oil rate are positively correlated (i.e. loss late is consistent), therefore, to before test analysis after sample acquisition Loss amount the use of measured data later period will not be had an impact.It can be realized different using abundance of organic matter as bridge in research The comparison of evaluation unit rock sample oil content.

3. evaluation unit saves coefficient relative weighting assignment method

According to abundance of organic matter (TOC) and oil content (adsorbed hydrocarbons S₁Or chloroform bitumen " A ") scatterplot data, it uses OriginLab (or excel) Data Analysis Software carries out data quantization relationship analysis, obtains TOC and the suction of different Lithofacies Types Attached hydrocarbon S₁The quantitative formula of (or chloroform bitumen " A ") between the two:

C_{Oil base}=a*TOC^b (1)

C_{Oil folder}=c+d*TOC (2)

It includes two steps that evaluation unit, which saves coefficient relative weighting assignment:

Step 1, single-point save coefficient with respect to assignment.According to the corresponding absorption of the identical TOC of the same lithofacies of different evaluation unit Hydrocarbon S₁(or chloroform bitumen " A ") relative amount quantitative relation establishes the preservation coefficient quantization model of different evaluation unit, and then really The preservation coefficient relative weighting assignment of fixed each evaluation unit.Specific method is closed according to the TOC of above-mentioned foundation and oil content quantization It is formula, by the rock of the rock oil content of a certain abundance of organic matter of other evaluation units abundance of organic matter corresponding with standard regions Oil content quantitative relation determines the preservation coefficient relative weighting assignment of the corresponding lithofacies of corresponding evaluation unit.

In formula: K_{Unit base i}Coefficient relative weighting assignment is saved for matrix type lithofacies evaluation unit i；C_{Oil base i}For matrix type lithofacies Evaluation unit i oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；C_{Oil disjunction mark is quasi-}For matrix type Lithofacies evaluation unit standard regions oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；a_i, b_i For matrix type lithofacies evaluation unit i equation regression coefficient；a_{Disjunction mark}, b_{Disjunction mark}For matrix type lithofacies evaluation unit standard regions equation regression system Number；TOC_{Base i}For matrix type lithofacies evaluation unit i abundance of organic matter, unit %；TOC_{Disjunction mark}For matrix type lithofacies evaluation unit standard Area's abundance of organic matter, unit %；I is the evaluation unit i marked off.

In formula: K_{Unit presss from both sides i}Coefficient relative weighting assignment is saved for sandwich type lithofacies evaluation unit i；C_{Oil presss from both sides i}For sandwich type lithofacies Evaluation unit i oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；C_{Oil presss from both sides standard}For sandwich type Lithofacies evaluation unit standard regions oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；c_i, d_i For sandwich type lithofacies evaluation unit i equation regression coefficient；C_{Folder mark}, d_{Folder mark}For sandwich type lithofacies evaluation unit standard regions equation regression system Number；TOC_{Press from both sides i}For sandwich type lithofacies evaluation unit i abundance of organic matter, unit %；TOC_{Folder mark}For sandwich type lithofacies evaluation unit standard Area's abundance of organic matter, unit %；I is the evaluation unit i marked off.

Step 2, different evaluation unit save the comprehensive assignment of coefficient relative weighting.Mainly according to research area's difference lithofacies Type abundance of organic matter distribution characteristics, on the basis of determining that abundance of organic matter comparison section and comparison are spaced, by multiple single-points It saves coefficient and seeks arithmetic mean of instantaneous value as the comprehensive assignment of different units preservation coefficient weights with respect to assignment.Including following two step:

1. different lithofacies abundance of organic matter Characteristics of Distribution.In order to guarantee that the difference of the identical lithofacies of different evaluation unit has The oil content Relative quantification relationship that machine matter abundance section obtains is representative, needs lithofacies identical in analysis different evaluation unit On the basis of surveying sample TOC distribution characteristics, reasonable TOC comparison section is chosen.Statistical sample point, point lithofacies divide organic matter Research area's abundance of organic matter distribution histogram, the identical lithofacies abundance of organic matter of comparative analysis different evaluation unit are drawn in section (TOC) distribution characteristics determines that different evaluation unit oil content TOC compares section.

Coefficient comprehensive weight assignment is saved 2. determining.According to abundance of organic matter (TOC) characteristic distributions of different lithofacies, determine Main TOC distributed area (m, n].Comprehensive analysis chooses certain interval TOC (inter), such as 0.1%, 0.2% or 0.5%, and then determine contrast sample's point and sample number (Num) between evaluation unit, it takes each sample to count number and saves coefficient weights phase The preservation coefficient comprehensive weight assignment of lithofacies is corresponded to as corresponding evaluation unit to assignment average value.

Num=(n-m)/inter (5)

J=1,2,3 ..., Num

5) divide evaluation unit, divide Lithofacies Types Resources calculation: in rock abundance of organic matter (TOC), pyrolysis free hydrocarbon (S1) lighter hydrocarbons compensation, heavy hydrocarbon recovery, the compensation of chloroform bitumen " A " lighter hydrocarbons and S₁, the crucial ginseng of chloroform bitumen " A " saturation adsorption coefficient etc. On the basis of number determines, for different shale oil occurrence forms (matrix type, sandwich type), pyrolysed hydrocarbon S is respectively adopted₁Method, (or Chloroform bitumen " A ") and oil saturation method carry out shale oil virgin abundance, virgin biomass calculating.On this basis, according to different evaluation The corresponding opposite preservation coefficient of unit, different Lithofacies Types, determines the stock number now of shale oil.

Rock pyrolysis free hydrocarbon (S₁) lighter hydrocarbons compensation, heavy hydrocarbon recovery, the compensation of chloroform bitumen " A " lighter hydrocarbons and saturation adsorption coefficient Equal key parameters, according to the correlation technique that Wang Min (2014), Xue Haitao (2016) propose, binding area rock organic matter is rich Degree, pyrolysis parameters of rock (S₁), chloroform bitumen " A " extracting, maturity Ro, race composition measured data specifically determine, using pyrolysis adsorb Hydrocarbon S₁Method (or chloroform bitumen "A" method) calculates matrix type shale oil virgin abundance, virgin biomass, calculates sandwich type using oil saturation method Shale oil virgin abundance, virgin biomass, in conjunction with different evaluation unit, different lithofacies (shale oil occurrence form), accordingly preservation coefficient weights are comprehensive Assignment is closed, shale oil stock number now is obtained.

By taking the Middle P ermian reed grass ditch group shale oil resource assessment of Junggar Basin Southeastern Margin ring Bogda Mountain as an example, propose Sub-unit, the shale oil Resources calculation Research Thinking for dividing lithofacies, the measured data analysis based on a large amount of samples, point evaluation are single Member, a point lithofacies establish abundance of organic matter (TOC) and oil content parameter (adsorbed hydrocarbons S₁Or chloroform bitumen " A ") quantitative relationship model, Realize the quantization assignment for saving coefficient relative weighting of different evaluation unit, different lithofacies shale oil.It is original in conventional method On the basis of shale oil Resources calculation, respectively multiplied by different evaluation unit, the corresponding preservation coefficient of different lithofacies to get to now Modern shale oil stock number.Shale oil stock number storage unit coefficient assignment method proposed by the present invention is more scientific, obtains As a result more objective reasonable, it fundamentally solves and type basin is transformed using currently a popular shale oil RESERVE EVALUATION method When ground shale oil Resources calculation there is error in result, can provide technical support for the practical investigation and prospecting of shale oil.

The present invention is that evaluation unit saves coefficient weights quantization assignment method in a kind of shale oil Resources calculation, specific to wrap Include following steps:

1) evaluation unit divides: conduct a research area's tectonic interpretation first, layer distribution characteristics that have a definite purpose and main disconnected Distribution Characteristics are split, using target zone overlap of beds or pinching distribution, boundary fault as boundary condition, binding purpose layer buried depth of strata, Regional cap rock, the analysis of stratum aqueous nature (water type, salinity) preservation condition, divide different evaluation units.Since fracture belt is One complicated three-dimensional geologic has certain width, field section fracture belt structure observation statistics, in tomography two sides 2km model Interior diaclase development is enclosed, preservation condition is poor, should give consideration when using tomography as evaluation unit boundary.According to quasi- Southeastern Margin The practical geology design feature of ring Bogda Mountain, ten thousand geologic maps of boundary faults and 1:20 such as the distribution of reed grass ditch group, Fukang tomography Reed grass ditch group earth's surface exposure situation, in combination with up to 1, Miquan 1, surprise 1,1 ground Chai Can drilling stratum aqueous nature, strata pressure feelings Condition, comprehensive analysis determine evaluation unit boundary, mark off depressed area, structurally stable region, Thrust tectonics transformation area and earth's surface exposure 4 class evaluation unit of area.

2) typical well (or outcrop profile) is chosen: basic demand be in the every a kind of evaluation unit marked off at least A bite typical case well (or section), and drilling well (or section) has systematic abundance of organic matter (TOC), pyrolysis adsorbed hydrocarbons (S₁) and chlorine Imitative pitch " A " basic geochemistry test data, has preferable representativeness of sample, generally at least to there is 20 or more data Point.Well (section) requirement is selected based on above-mentioned, Ji174Jing, Ji251Jing and Ji22Jing is chosen in depressed area, is selected in structurally stable region 1 well of new lucky ginseng and odd 1 well are taken, 1 well of rich ginseng, 4 well of quasi- page, drilling 1 is chosen in Thrust tectonics transformation area, chooses Jing Jing in outcrop area Cunette, small Longkou, the big Longkou in west, Dahuangshan, Yao Moshan, three work river sections carry out congestion system sampling analysis.Each evaluation unit Inside there are different lithology measured data point 29~172, sample has preferable representative.

3) in the step, evaluation unit, lithologic character and lithofacies type is divided to establish TOC and oil content parameter (adsorbed hydrocarbons S respectively₁ Or chloroform bitumen " A ") relationship plate.It is assumed that the rock sample of the comparable identical lithofacies of abundance of organic matter (TOC), not by structure Making its oil content under the reset condition that transformation influences should substantially quite, and preservation condition superiority and inferiority caused by later structural is transformed makes existing Modern different evaluation unit abundance of organic matter is suitable, the identical rock oil content of lithofacies has differences.In order to exclude different lithofacies The influence of oil content caused by rock pore structure factor and oil and gas loss otherness can divide the type of different lithofacies, respectively structure Build the TOC and (S of different evaluation unit, different lithofacies₁Or chloroform bitumen " A ") quantitative relationship model.Specifically include Lithofacies Types It divides and plate is established:

1. lithofacies divide

Using above-mentioned Lithofacies Identification means and splitting scheme, ring Bogda Mountain is marked off into matrix type lithofacies (shale, mud Rock, Sandy Silt) and sandwich type lithofacies (sandstone, carbonate rock) two types.

(recess is divided according to ring Bogda Mountain reed grass ditch group matrix type and sandwich type lithofacies development characteristics and geologic element Transformation area and earth's surface exposure area are thrust in area, structurally stable region), sub-unit, a point lithofacies establish abundance of organic matter and oil content (adsorbed hydrocarbons S₁, chloroform bitumen " A ") corresponding plate (as shown in Fig. 4~Figure 13).

4) in the step, the preservation coefficient relative weighting assignment of different evaluation unit, different lithofacies is determined, including following Step:

1. standard regions select

Ring Bogda Mountain is influenced by Himalayan is tectogenetic, and geologic structure is complicated, different units preservation condition Difference is big.Jimusaer Sag has certain distance apart from BogdaM ountain, and tectonically active basins obviously weaken, reed grass ditch group deposition with To show as Continuous Settlement evolution Feature.In order to reduce influence of the tectonism to rock oil content, the lucky wooden Sa of selection in research The standard regions that your recess is compared as reference.

2. building different evaluation unit oil content comparison bridge

Rock abundance of organic matter (TOC) by preservation condition and later period slacking although also influenced, relative to pyrolysis For adsorbed hydrocarbons S1 (or chloroform bitumen " A "), degree of susceptibility is relatively weak.It is assumed that the sample of the identical lithofacies of different evaluation unit Product free hydrocarbon loss amount and its initial oil rate are positively correlated (i.e. loss late is consistent), therefore, to before test analysis after sample acquisition Loss amount the use of measured data later period will not be had an impact.It can be realized different using abundance of organic matter as bridge in research The comparison of evaluation unit rock sample oil content.

3. evaluation unit saves coefficient relative weighting assignment method

C_{Oil base}=a*TOC^b (1)

C_{Oil folder}=c+d*TOC (2)

According to the organic matter of 2 kinds of Lithofacies Types samples of reed grass ditch group matrix type and sandwich type of Bogda Mountain difference zone Abundance (TOC) and oil content (S₁, chloroform bitumen " A ") relationship analysis, establish its quantitative formula respectively:

S_{1 depressed area base}=0.1748*TOC^1.3565 (3)

S_{1 stable region base}=0.1254*TOC^1.3122 (4)

S_{1 thrusts Qu Ji}=0.0416*TOC^1.4753 (5)

S_{1 outcrop area base}=0.0289*TOC^1.3364 (6)

“A”_{Depressed area base}=0.0650*TOC^1.2298 (7)

“A”_{Stable region base}=0.0489*TOC^1.1257 (8)

“A”_{Thrust Qu Ji}=0.0241*TOC^0.9237 (9)

“A”_{Outcrop area base}=0.0137*TOC^1.0231 (10)

S_{1 depressed area folder}=-0.5123+1.1577*TOC (11)

S_{1 stable region folder}=-0.2804+0.5098*TOC (12)

S_{1 thrusts area's folder}=-0.1272+0.2101*TOC (13)

S_{1 outcrop area folder}=-0.1285+0.1824*TOC (14)

“A”_{Depressed area folder}=-0.4648+1.3015*TOC (15)

“A”_{Stable region folder}=-0.3038+0.7316*TOC (16)

“A”_{Thrust area's folder}=-0.2021+0.3824*TOC (17)

“A”_{Outcrop area folder}=-0.1520+0.2592*TOC (18)

In formula: S_{1 depressed area base}Adsorbed hydrocarbons (S is pyrolyzed for depressed area evaluation unit matrix type lithofacies₁) content, unit mg/g； S_{1 stable region base}Adsorbed hydrocarbons (S is pyrolyzed for structurally stable region evaluation unit matrix type lithofacies₁) content, unit mg/g；S_{1 thrusts Qu Ji}For construction punching Disconnected area's evaluation unit matrix type lithofacies are pyrolyzed adsorbed hydrocarbons (S₁) content, unit mg/g；S_{1 outcrop area base}For surface occurence area evaluation unit base Matter type lithofacies are pyrolyzed adsorbed hydrocarbons (S₁) content, unit mg/g；S_{1 depressed area folder}Adsorbed hydrocarbons are pyrolyzed for depressed area evaluation unit sandwich type lithofacies (S₁) content, unit mg/g；S_{1 stable region folder}Adsorbed hydrocarbons (S is pyrolyzed for structurally stable region evaluation unit sandwich type lithofacies₁) content, unit mg/g；S_{1 thrusts area's folder}Area's evaluation unit sandwich type lithofacies, which are thrust, for construction is pyrolyzed adsorbed hydrocarbons (S₁) content, unit mg/g；S_{1 outcrop area folder}For Surface occurence area evaluation unit sandwich type lithofacies are pyrolyzed adsorbed hydrocarbons (S₁) content, unit mg/g；"A"_{Depressed area base}It is evaluated for depressed area single First matrix type lithofacies rock extracts chloroform bitumen " A " content, unit %；"A"_{Stable region base}For structurally stable region evaluation unit matrix type Lithofacies rock extracts chloroform bitumen " A " content, unit %；"A"_{Thrust Qu Ji}Area's evaluation unit matrix type lithofacies rock is thrust for construction Extract chloroform bitumen " A " content, unit %；"A"_{Outcrop area base}Chloroform is extracted for surface occurence area evaluation unit matrix type lithofacies rock Pitch " A " content, unit %；"A"_{Depressed area folder}Chloroform bitumen " A " content is extracted for depressed area evaluation unit sandwich type lithofacies rock, Unit %；"A"_{Stable region folder}Chloroform bitumen " A " content, unit % are extracted for structurally stable region evaluation unit sandwich type lithofacies rock； “A”_{Thrust area's folder}Area's evaluation unit sandwich type lithofacies rock, which is thrust, for construction extracts chloroform bitumen " A " content, unit %；"A"_{Outcrop area folder}For Surface occurence area evaluation unit sandwich type lithofacies rock extracts chloroform bitumen " A " content, unit %；TOC is that rock organic matter is rich Angle value, unit %.

K_{Unit i base}=C_{Oil base i}/C_{Oil disjunction mark is quasi-}=a_i*TOC_{Base i}^b_i/a_Mark*TOC_{Disjunction mark}^b_Mark (19)

In formula: K_{Unit i base}Coefficient relative weighting assignment is saved for matrix type lithofacies evaluation unit i；C_{Oil base i}For matrix type lithofacies Evaluation unit i oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；C_{Oil disjunction mark is quasi-}For matrix type Lithofacies evaluation unit standard regions oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；a_i, b_i For matrix type lithofacies evaluation unit i equation regression coefficient；a_Mark, b_MarkFor matrix type lithofacies evaluation unit standard regions equation regression system Number；TOC_{Base i}For matrix type lithofacies evaluation unit i abundance of organic matter, unit %；TOC_{Disjunction mark}For matrix type lithofacies evaluation unit standard Area's abundance of organic matter, unit %；I is the evaluation unit i marked off.

According to the different evaluation unit quantitative formula that above-mentioned Bogda Mountain reed grass ditch group matrix type lithofacies are established, by not With the comparison of evaluation unit and Jimusaer Sag (standard regions), determines that it is corresponded to and save coefficient relative weighting assignment, it is specific public Formula are as follows:

K_{The stable region S1 base}=S_{1 stable region base}/S_{1 depressed area base}=(0.1254*TOC^1.3122)/(0.1748*TOC^1.3565) (20)

K_{S1 thrusts Qu Ji}=S_{1 thrusts Qu Ji}/S_{1 depressed area base}=(0.0416*TOC^1.4753)/(0.1748*TOC^1.3565) (21)

K_{The outcrop area S1 base}=S_{1 outcrop area base}/S_{1 depressed area base}=(0.0289*TOC^1.3364)/(0.1748*TOC^1.3565) (22)

In formula: K_{The stable region S1 base}Coefficient relative weighting assignment, decimal are saved for structurally stable region evaluation unit matrix type lithofacies； K_{S1 thrusts Qu Ji}Area's evaluation unit matrix type lithofacies, which are thrust, for construction saves coefficient relative weighting assignment, decimal；K_{The outcrop area S1 base}For ground expression Head Section evaluation unit matrix type lithofacies save coefficient relative weighting assignment, decimal；S_{1 depressed area base}For depressed area evaluation unit matrix type Lithofacies are pyrolyzed adsorbed hydrocarbons (S₁) content, unit mg/g；S_{1 stable region base}Adsorbed hydrocarbons are pyrolyzed for structurally stable region evaluation unit matrix type lithofacies (S₁) content, unit mg/g；S_{1 thrusts Qu Ji}Area's evaluation unit matrix type lithofacies, which are thrust, for construction is pyrolyzed adsorbed hydrocarbons (S₁) content, unit mg/g；S_{1 outcrop area base}Adsorbed hydrocarbons (S is pyrolyzed for surface occurence area evaluation unit matrix type lithofacies₁) content, unit mg/g；TOC is rock Abundance of organic matter value, unit %.

K_{Unit presss from both sides i}=C_{Oil presss from both sides i}/C_{Oil presss from both sides standard}=(C_i+d_i*TOC_{Press from both sides i})/(C_Mark+d_Mark*TOC_{Folder mark}) (23)

In formula: K_{Unit presss from both sides i}Coefficient relative weighting assignment is saved for sandwich type lithofacies evaluation unit i；C_{Oil presss from both sides i}For sandwich type lithofacies Evaluation unit i oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；C_{Oil presss from both sides standard}For sandwich type Lithofacies evaluation unit standard regions oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；c_i, d_i For sandwich type lithofacies evaluation unit i equation regression coefficient；c_Mark, d_MarkFor sandwich type lithofacies evaluation unit standard regions equation regression system Number；TOC_{Press from both sides i}For sandwich type lithofacies evaluation unit i abundance of organic matter, unit %；TOC_{Folder mark}For sandwich type lithofacies evaluation unit standard Area's abundance of organic matter, unit %；I is the evaluation unit i marked off.

According to the different evaluation unit quantitative formula that above-mentioned Bogda Mountain reed grass ditch group sandwich type lithofacies are established, by not With the comparison of evaluation unit and Jimusaer Sag (standard regions), determines that it is corresponded to and save coefficient relative weighting assignment, it is specific public Formula are as follows:

K_{The stable region S1 folder}=S_{1 stable region folder}/S_{1 depressed area folder}=(- 0.2804+0.5098*TOC)/(- 0.5123+1.1577* TOC) (24)

K_{S1 thrusts area's folder}=S_{1 thrusts area's folder}/S_{1 depressed area folder}=(- 0.1272+0.2101*TOC)/(- 0.5123+1.1577* TOC) (25)

K_{The outcrop area S1 folder}=S_{1 outcrop area folder}/S_{1 depressed area folder}=(- 0.1285+0.1824*TOC)/(- 0.5123+1.1577* TOC) (26)

In formula: K_{The stable region S1 folder}Coefficient relative weighting assignment, decimal are saved for structurally stable region evaluation unit sandwich type lithofacies； K_{S1 thrusts area's folder}Area's evaluation unit sandwich type lithofacies, which are thrust, for construction saves coefficient relative weighting assignment, decimal；K_{The outcrop area S1 folder}For ground expression Head Section evaluation unit sandwich type lithofacies save coefficient relative weighting assignment, decimal；S_{1 depressed area folder}For depressed area evaluation unit sandwich type Lithofacies adsorbed hydrocarbons S₁Content, unit mg/g；S_{1 stable region folder}For structurally stable region evaluation unit sandwich type lithofacies adsorbed hydrocarbons S₁Content, it is single Position mg/g；S_{1 thrusts area's folder}Area evaluation unit sandwich type lithofacies adsorbed hydrocarbons S is thrust for construction₁Content, unit mg/g；S_{1 outcrop area folder}For earth's surface Outcrop area evaluation unit sandwich type lithofacies adsorbed hydrocarbons S₁Content, unit mg/g；TOC is rock abundance of organic matter value, unit %.

Abundance of organic matter is surveyed according to Bogda Mountain different evaluation unit matrix type and sandwich type lithofacies rock sample (TOC) it statisticallys analyze, matrix type lithofacies abundance of organic matter (TOC) distributed area is 0.35%~30.8%, and sandwich type lithofacies have Machine matter abundance (TOC) distributed area Relatively centralized is 0.14%~5.01% (as shown in Figure 14,15).Different evaluation unit TOC Distribution characteristics comprehensive analysis, matrix type lithofacies select TOC comparison section be (1.0,12.0], sandwich type lithofacies select TOC compare Section be (0.5,4.0].

Coefficient comprehensive weight assignment is saved 2. determining.According to abundance of organic matter (TOC) characteristic distributions of different lithofacies, determine Main TOC distributed area (m, n].Comprehensive analysis chooses certain interval TOC (inter), such as 0.1%, 0.2% or 0.5%, and then determine contrast sample's point and sample number (Num) between evaluation unit, taking each sample point to save, coefficient weights are opposite to be assigned Value average value corresponds to the preservation coefficient comprehensive weight assignment of lithofacies as corresponding evaluation unit.

Num=(n-m)/inter (27)

K_{Unit i base is comprehensive}=(K_{Unit i base 1}+K_{Unit i base 2}+…+K_{Unit i base j})/Num (28)

K_{Unit i folder is comprehensive}=(K_{Unit i folder 1}+K_{Unit i folder 2}+…+K_{Unit i presss from both sides j})/Num (29)

J=1,2,3 ..., Num

According to the comprehensive assignment formula of above-mentioned preservation coefficient relative weighting, 4 classes that Bogda Mountain marks off are evaluated respectively The reed grass ditch group matrix type and sandwich type shale oil of unit save coefficient weights, according to averaged between corresponding toc region, make To save the comprehensive assignment of coefficient weights, wherein pyrolysis adsorbed hydrocarbons (S₁) matrix type lithofacies TOC comparison interval take 0.2%, sandwich type Lithofacies TOC comparison section takes 0.1%；Chloroform bitumen " A " matrix type lithofacies TOC comparison interval takes 0.2%, sandwich type lithofacies TOC Comparison section takes 0.1%.

Matrix type lithofacies adsorbed hydrocarbons S₁Save the comprehensive assignment calculation formula of coefficient weights are as follows:

Num=(n-m)/inter=(12-1)/0.2=55 (30)

K_{The stable region S1 base is comprehensive}=(K_{The stable region S1 base 1}+K_{The stable region S1 base 2}+…+K_{The stable region S1 base j})/Num=(0.7127+0.7078+ ...+0.6440)/ 55=0.67 (31)

K_{It is comprehensive that S1 thrusts Qu Ji}=(K_{S1 thrusts Qu Ji 1}+K_{S1 thrusts Qu Ji 2}+…+K_{S1 thrusts Qu Ji j})/Num=(0.2453+0.2492+ ...+ 0.31670)/55=0.29 (32)

K_{The outcrop area S1 base is comprehensive}=(K_{The outcrop area S1 base 1}+K_{The outcrop area S1 base 2}+…+K_{The outcrop area S1 base j})/Num=(0.1598+0.1594+ ...+0.1531)/ 55=0.16 (33)

J=1,2,3 ..., 55

In formula: m left side low value between toc region, unit is %；N right side high level between toc region, unit is %；Inter is TOC comparative analysis value interval, unit is %；Num comparative analysis number of samples between different evaluation unit, unit are； K_{The stable region S1 base is comprehensive}For structurally stable region matrix type lithofacies adsorbed hydrocarbons (S₁) save the comprehensive assignment of coefficient weights, decimal；K_{It is comprehensive that S1 thrusts Qu Ji}For structure It makes and thrusts territorial matrix type lithofacies adsorbed hydrocarbons (S₁) save the comprehensive assignment of coefficient weights, decimal；K_{The outcrop area S1 base is comprehensive}For surface occurence territorial matrix Type lithofacies adsorbed hydrocarbons (S₁) save the comprehensive assignment of coefficient weights, decimal.

Sandwich type lithofacies adsorbed hydrocarbons S₁Save the comprehensive assignment calculation formula of coefficient weights are as follows:

Num=(n-m)/inter=(4-0.5)/0.1=30 (34)

K_{The stable region S1 folder is comprehensive}=(K_{The stable region S1 folder 1}+K_{The stable region S1 folder 2}+…+K_{The stable region S1 presss from both sides j})/Num=(0.3684+0.3779+ ...+0.4270)/ 30=0.41 (35)

K_{It is comprehensive that S1 thrusts area's folder}=(K_{S1 thrusts area's folder 1}+K_{S1 thrusts area's folder 2}+…+K_{S1 thrusts area folder j})/Num=(0.1365+0.0.1425+ ...+ 0.1732)/30=0.16 (36)

K_{The outcrop area S1 folder is comprehensive}=(K_{The outcrop area S1 folder 1}+K_{The outcrop area S1 folder 2}+…+K_{The outcrop area S1 presss from both sides j})/Num=(0.0949+0.1032+ ...+0.1460)/ 30=0.13 (37)

J=1,2,3 ..., 30

In formula: m left side low value between toc region, unit is %；N right side high level between toc region, unit is %；Inter is TOC comparative analysis value interval, unit is %；Num comparative analysis number of samples between different evaluation unit, unit are； K_{The stable region S1 folder is comprehensive}The comprehensive assignment of coefficient weights, decimal are saved for structurally stable region matrix type lithofacies chloroform bitumen " A "；K_{It is comprehensive that S1 thrusts area's folder}For Construction thrusts territorial matrix type lithofacies chloroform bitumen " A " and saves the comprehensive assignment of coefficient weights, decimal；K_{The outcrop area S1 folder is comprehensive}For surface occurence area Matrix type lithofacies chloroform bitumen " A " saves the comprehensive assignment of coefficient weights, decimal.

Utilize rock pyrolysis adsorbed hydrocarbons (S₁) parametric method determine each evaluation unit matrix type lithofacies save coefficient it is opposite It is 1.0 that the comprehensive assignment of weight, which is respectively as follows: depressed area, and structurally stable region 0.68, thrusting transformation area is 0.27, and outcrop area is 0.16；The comprehensive assignment of sandwich type lithofacies preservation coefficient relative weighting of each evaluation unit determined using this method is respectively as follows: recessed Falling into area is 1.0, structurally stable region 0.41, and thrusting transformation area is 0.16, outcrop area 0.13.

Matrix type lithofacies chloroform bitumen " A " saves the comprehensive assignment calculation formula of coefficient weights are as follows:

Num=(n-m)/inter=(12-1)/0.2=55 (38)

K_{" stable region A " base is comprehensive}=(K_{" stable region A " base 1}+K_{" stable region A " base 2}+…+K_{" stable region A " base j})/Num=(0.7382+0.7264+ ...+ 0.5808)/55=0.63 (39)

K_{" to thrust Qu Ji comprehensive by A "}=(K_{" A " thrusts Qu Ji 1}+K_{" A " thrusts Qu Ji 2}+…+K_{" A " thrusts Qu Ji j})/Num=(0.3506+0.3345+ ...+ 0.1733)/55=0.22 (40)

K_{" outcrop area A " base is comprehensive}=(K_{" outcrop area A " base 1}+K_{" outcrop area A " base 2}+…+K_{" outcrop area A " base j})/Num=(0.2030+0.1966+ ...+ 0.1261)/55=0.15 (41)

J=1,2,3 ..., 55

In formula: m left side low value between toc region, unit is %；N right side high level between toc region, unit is %；Inter is TOC comparative analysis value interval, unit is %；Num comparative analysis number of samples between different evaluation unit, unit are； K_{" stable region A " base is comprehensive}The comprehensive assignment of coefficient weights, decimal are saved for structurally stable region matrix type lithofacies chloroform bitumen " A "；K_{" to thrust Qu Ji comprehensive by A "}For Construction thrusts territorial matrix type lithofacies chloroform bitumen " A " and saves the comprehensive assignment of coefficient weights, decimal；K_{" outcrop area A " base is comprehensive}For surface occurence area Matrix type lithofacies chloroform bitumen " A " saves the comprehensive assignment of coefficient weights, decimal.

Sandwich type lithofacies chloroform bitumen " A " saves the comprehensive assignment calculation formula of coefficient weights are as follows:

Num=(n-m)/inter=(4-0.5)/0.1=30 (42)

K_{" stable region A " folder it is comprehensive}=(K_{" stable region A " folder 1}+K_{" stable region A " folder 2}+…+K_{" stable region A " folder j})/Num=(0.4276+0.4668+ ...+ 0.5529)/30=0.53 (43)

K_{" A " thrust area folder it is comprehensive}=(K_{" A " thrust area folder 1}+K_{" A " thrust area folder 2}+…+K_{" A " thrust area folder j})/Num=(0.0865+0.1470+ ...+ 0.2796)/30=0.25 (44)

K_{" outcrop area A " folder it is comprehensive}=(K_{" outcrop area A " folder 1}+K_{" outcrop area A " folder 2}+…+K_{" outcrop area A " folder j})/Num=(0.0111+0.0660+ ...+ 0.1863)/30=0.16 (45)

J=1,2,3 ..., 30

In formula: m left side low value between toc region, unit is %；N right side high level between toc region, unit is %；Inter is TOC comparative analysis value interval, unit is %；Num comparative analysis number of samples between different evaluation unit, unit are； K_{" stable region A " folder it is comprehensive}The comprehensive assignment of coefficient weights, decimal are saved for structurally stable region sandwich type lithofacies chloroform bitumen " A "；K_{" A " thrust area folder it is comprehensive}For Construction thrusts area's sandwich type lithofacies chloroform bitumen " A " and saves the comprehensive assignment of coefficient weights, decimal；K_{" outcrop area A " folder it is comprehensive}For surface occurence area Sandwich type lithofacies chloroform bitumen " A " saves the comprehensive assignment of coefficient weights, decimal.

With rock pyrolysis adsorbed hydrocarbons (S₁) method principle is identical, each evaluation for utilizing rock extracting chloroform bitumen "A" method to determine It is 1.0 that the matrix type lithofacies of unit, which save the comprehensive assignment of coefficient relative weighting to be respectively as follows: depressed area, structurally stable region 0.67, Thrusting transformation area is 0.26, outcrop area 0.17；The sandwich type lithofacies of each evaluation unit determined using this method save coefficient It is 1.0 that the comprehensive assignment of relative weighting, which is respectively as follows: depressed area, and structurally stable region 0.53, thrusting transformation area is 0.25, and outcrop area is 0.16。

5) in the step, rock pyrolysis free hydrocarbon (S₁) lighter hydrocarbons compensation, heavy hydrocarbon recovery, the compensation of chloroform bitumen " A " lighter hydrocarbons And the key parameters such as saturation adsorption coefficient, according to the correlation technique that Wang Min (2014), Xue Haitao (2016) propose, binding area Rock abundance of organic matter, pyrolysis parameters of rock (S₁), chloroform bitumen " A " extracting, maturity Ro, race composition measured data specifically determine, Using pyrolysis adsorbed hydrocarbons S₁Method (or chloroform bitumen "A" method) calculates matrix type shale oil virgin abundance, virgin biomass, using oil saturation Method calculates sandwich type shale oil virgin abundance, virgin biomass, accordingly protects in conjunction with different evaluation unit, different lithofacies (shale oil occurrence form) The comprehensive assignment of coefficient weights is deposited, shale oil stock number now is obtained.

Petrochemical industry exploratory area in ring BogdaM ountain area, using " stable type " the regional shale oil Resources calculation side proposed early period Method, original shale oil total resources is up to 50.2~50.6 hundred million tons, wherein adsorbed hydrocarbons (S₁) the original oil shale money of method calculating matrix type Source amount is 25.2 hundred million tons, and it is 25.6 hundred million tons that chloroform bitumen "A" method, which calculates the original oil shale stock number of matrix type, oil saturation method Calculating the original shale oil stock number of sandwich type is 25.0 hundred million tons.Using " reworked " proposed by the present invention regional shale oil stock number Calculation method, on the basis of original shale oil Resources calculation, different evaluation unit and lithofacies are obtained multiplied by the corresponding coefficient that saves The total resources of shale oil now arrived is 15.4~16.8 hundred million tons, the one third or so of only original shale oil stock number, In, adsorbed hydrocarbons (S₁) shale oil stock number is 7.2 hundred million tons to method calculating matrix type now, shale oil stock number is sandwich type now 8.2 hundred million tons；Chloroform bitumen "A" method calculates matrix type, and shale oil stock number is 6.3 hundred million tons now, sandwich type shale oil resource now Amount is 10.5 hundred million tons.

Above-described embodiment for explaining only the invention, rather than limits the invention, for those skilled in the art For member, on the basis of application principle method disclosed by the invention, the equivalents or improvement of progress are each fallen within of the invention Protection scope.

Claims

1. evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation, which is characterized in that the shale oil Evaluation unit preservation coefficient weights quantization assignment method includes: in Resources calculation

Step 1, according to the practical geologic structure in research area and preservation condition feature, different evaluation units is divided；

Step 2, typical well or outcrop profile are chosen；

Step 3, different evaluation unit, the abundance of organic matter of different lithology lithofacies and oil content parameters relationship plate are established；

Step 4, different evaluation unit is established, point lithofacies save coefficient weights assignment quantitative model；

Step 5, divide evaluation unit, divide Lithofacies Types amount of computational resources.

2. evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation according to claim 1, It is characterized in that, in step 1, conduct a research area's tectonic interpretation first, have a definite purpose layer distribution characteristics and main faults Distribution Characteristics, using target zone overlap of beds or pinching distribution, boundary fault as boundary condition, binding purpose floor buried depth of strata, area Domain cap rock, the analysis of stratum aqueous nature preservation condition, divide different evaluation units.

3. evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation according to claim 1, It is characterized in that, in step 2, basic demand at least a bite typical case well or is cutd open in the every a kind of evaluation unit marked off Face, and drilling well or section have systematic abundance of organic matter TOC, pyrolysis adsorbed hydrocarbons S₁With chloroform bitumen " A " basic geochemistry Test data has representativeness of sample, at least to there is 20 or more data points.

4. evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation according to claim 1, It is characterized in that, step 3 includes:

Step 3a carries out lithofacies division；

5. evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation according to claim 4, It is characterized in that, in step 3a, based on rock core, appear, thin slice, the macro microscopic observation of scanning electron microscope and total rock X diffraction analysis, according to Rock deposition construction, mineral constituent, pore structure characteristic, using " rock constituents-sedimentary structure-organic matter " lithofacies division side Case divides different Lithofacies Types.

6. evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation according to claim 4, It is characterized in that, counting the rock abundance of organic matter TOC of different lithofacies in each evaluation unit, pyrolysis absorption respectively in step 3b Hydrocarbon S₁Or chloroform bitumen " A " measured data, system comparative analysis is carried out to data, rejects individual abnormal points；According to research area page Shale oil Lithofacies Types splitting scheme and research need, establish different evaluation unit, different Lithofacies Types abundance of organic matter TOC with Oiliness S₁Or chloroform bitumen " A " scatterplot relationship plate.

7. evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation according to claim 1, It is characterized in that, step 4 includes:

Step 4a, selection criteria area；

Step 4b builds different evaluation unit oil content comparison bridge；

8. evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation according to claim 7, It is characterized in that, in step 4a, it is contemplated that identical lithofacies rock caused by different evaluation unit preservation condition superiority and inferiority is different contains The difference of oil cut rate, according to the research practical geological condition in area, the area for choosing preservation condition construction of stable the most superior is standard Area saves coefficient relative weighting assignment and takes 1.0, and it is in contrast that other kinds of evaluation unit saves coefficient relative weighting assignment It obtains.

9. evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation according to claim 7, It is characterized in that, in step 4b, it is assumed that the sample free hydrocarbon loss amount of the identical lithofacies of different evaluation unit and its initial oil rate Being positively correlated is that loss late is consistent, therefore, will not to the use of measured data later period to the loss amount before test analysis after sample acquisition It has an impact；Using abundance of organic matter as bridge, the comparison of different evaluation unit rock sample oil content is realized.

10. evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation according to claim 7, It is characterized in that, in step 4c, according to abundance of organic matter TOC and oil content adsorbed hydrocarbons S₁Or chloroform bitumen " A " scatterplot data, into The analysis of row data quantitative relationship, obtains the TOC and adsorbed hydrocarbons S of different Lithofacies Types₁Or the quantization of chloroform bitumen " A " between the two Formula:

C_{Oil base}=a*TOC^b (1)

C_{Oil folder}=c+d*TOC (2)

In formula: C_{Oil base}、C_{Oil folder}Respectively matrix type, sandwich type rock oil content parameter are pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or take out Mention chloroform bitumen " A ", unit %；A, b, c, d are equation regression coefficient；TOC is rock abundance of organic matter, unit %.

11. evaluation unit saves coefficient weights quantization assignment side in shale oil Resources calculation according to claim 10 Method, which is characterized in that step 4c includes:

Step 1, single-point save coefficient with respect to assignment；According to the corresponding adsorbed hydrocarbons S of the identical TOC of the same lithofacies of different evaluation unit₁ Or chloroform bitumen " A " relative amount quantitative relation, the preservation coefficient quantization model of different evaluation unit is established, and then determination is respectively commented The preservation coefficient relative weighting assignment of valence unit；According to the TOC of above-mentioned foundation and oil content quantitative relationship formula, other are evaluated The rock oil content quantitative relation of the rock oil content of a certain abundance of organic matter of unit abundance of organic matter corresponding with standard regions is true Surely the preservation coefficient relative weighting assignment of the corresponding lithofacies of evaluation unit is corresponded to；

In formula: K_{Unit base i}Coefficient relative weighting assignment is saved for matrix type lithofacies evaluation unit i；C_{Oil base i}For the evaluation of matrix type lithofacies Unit i oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；C_{Oil disjunction mark is quasi-}For matrix type lithofacies Evaluation unit standard regions oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；a_i, b_iFor base Matter type lithofacies evaluation unit i equation regression coefficient；a_{Disjunction mark}, b_{Disjunction mark}For matrix type lithofacies evaluation unit standard regions equation regression coefficient； TOC_{Base i}For matrix type lithofacies evaluation unit i abundance of organic matter, unit %；TOC_{Disjunction mark}Have for matrix type lithofacies evaluation unit standard regions Machine matter abundance, unit %；I is the evaluation unit i marked off；

K_{Unit presss from both sides i}=C_{Oil presss from both sides i}/C_{Oil presss from both sides standard}=(C_i+d_i*TOC_{Press from both sides i})/(C_{Folder mark}+d_{Folder mark}*TOC_{Folder mark}) (4)

In formula: K_{Unit presss from both sides i}Coefficient relative weighting assignment is saved for sandwich type lithofacies evaluation unit i；C_{Oil presss from both sides i}For the evaluation of sandwich type lithofacies Unit i oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；C_{Oil presss from both sides standard}For sandwich type lithofacies Evaluation unit standard regions oil content is pyrolyzed adsorbed hydrocarbons S₁, unit mg/g, or extracting chloroform bitumen " A ", unit %；c_i, d_iFor folder Stratotype lithofacies evaluation unit i equation regression coefficient；C_{Folder mark}, d_{Folder mark}For sandwich type lithofacies evaluation unit standard regions equation regression coefficient； TOC_{Press from both sides i}For sandwich type lithofacies evaluation unit i abundance of organic matter, unit %；TOC_{Folder mark}Have for sandwich type lithofacies evaluation unit standard regions Machine matter abundance, unit %；I is the evaluation unit i marked off；

Step 2, different evaluation unit save the comprehensive assignment of coefficient relative weighting；According to research area's difference Lithofacies Types organic matter Multiple single-points are saved coefficient phase on the basis of determining that abundance of organic matter comparison section and comparison are spaced by abundance distribution feature Arithmetic mean of instantaneous value is sought as different units to assignment and saves the comprehensive assignment of coefficient weights；Including following two step:

1. different lithofacies abundance of organic matter Characteristics of Distribution；Identical sample TOC points of lithofacies actual measurement in analysis different evaluation unit On the basis of cloth feature, reasonable TOC comparison section is chosen；Statistical sample point, point lithofacies divide organic matter section, draw research Area's abundance of organic matter distribution histogram, the identical lithofacies abundance of organic matter TOC distribution characteristics of comparative analysis different evaluation unit determine Different evaluation unit oil content TOC compares section；

Coefficient comprehensive weight assignment is saved 2. determining；According to the abundance of organic matter TOC characteristic distributions of different lithofacies, determine main TOC distributed area (m, n]；The certain interval TOC inter is chosen, and then determines contrast sample's point and sample number between evaluation unit Num takes each sample to count number and saves the preservation system that coefficient weights correspond to lithofacies with respect to assignment average value as corresponding evaluation unit Number synthesis weight assignment；

Num=(n-m)/inter (5)

J=1,2,3 ..., Num

In formula: m left side low value between toc region, unit is %；N right side high level between toc region, unit is %；Inter is TOC pairs Score extract value interval, unit is %；Num comparative analysis number of samples between different evaluation unit, unit are；K_{Unit i base j}For J-th of sample point of evaluation unit i matrix type lithofacies saves coefficient relative weighting assignment, decimal；K_{Unit i presss from both sides j}For evaluation unit i interlayer J-th of sample point of type lithofacies saves coefficient relative weighting assignment, decimal；K_{Unit i base is comprehensive}For evaluation unit i matrix type lithofacies shale oil Save coefficient relative weighting assignment mean values, decimal；K_{Unit i folder is comprehensive}Coefficient phase is saved for evaluation unit i sandwich type lithofacies shale oil To weight assignment mean values, decimal.

12. evaluation unit saves coefficient weights quantization assignment method in shale oil Resources calculation according to claim 1, It is characterized in that, in steps of 5, in rock abundance of organic matter TOC, the compensation of pyrolysis free hydrocarbon S1 lighter hydrocarbons, heavy hydrocarbon restores, chloroform drips Green " A " lighter hydrocarbons compensation and S₁, chloroform bitumen " A " saturation these key parameters of adsorption coefficient determine on the basis of, for different Pyrolysed hydrocarbon S is respectively adopted in shale oil occurrence form₁It is original that method or chloroform bitumen "A" method and oil saturation method carry out shale oil The calculating of stock number according to different evaluation unit, the corresponding opposite preservation coefficient of different Lithofacies Types, determines on this basis The stock number now of shale oil.