CN107808019A - Oil shale in-situ conversion temperature, time, conversion ratio plate method for building up and system - Google Patents

Oil shale in-situ conversion temperature, time, conversion ratio plate method for building up and system Download PDF

Info

Publication number
CN107808019A
CN107808019A CN201610811514.4A CN201610811514A CN107808019A CN 107808019 A CN107808019 A CN 107808019A CN 201610811514 A CN201610811514 A CN 201610811514A CN 107808019 A CN107808019 A CN 107808019A
Authority
CN
China
Prior art keywords
oil shale
conversion ratio
pyrolysis
temperature
time
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610811514.4A
Other languages
Chinese (zh)
Other versions
CN107808019B (en
Inventor
马中良
郑伦举
余晓露
徐勤
葛颖
赵中熙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Exploration and Production Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Exploration and Production Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201610811514.4A priority Critical patent/CN107808019B/en
Publication of CN107808019A publication Critical patent/CN107808019A/en
Application granted granted Critical
Publication of CN107808019B publication Critical patent/CN107808019B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/16Enhanced recovery methods for obtaining hydrocarbons
    • E21B43/24Enhanced recovery methods for obtaining hydrocarbons using heat, e.g. steam injection
    • E21B43/241Enhanced recovery methods for obtaining hydrocarbons using heat, e.g. steam injection combined with solution mining of non-hydrocarbon minerals, e.g. solvent pyrolysis of oil shale
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Mining

Abstract

The invention discloses a kind of oil shale in-situ conversion temperature, time, conversion ratio plate method for building up and system, this method includes intending production zone selection oil shale sample from oil shale, and the oil shale sample is tested using rock pyrolysis analysis instrument;Experimental result obtains the relation data between pyrolysis activation energy, conversion ratio and frequency factor;Established based on the relation data when conversion ratio is different numerical value, the Core analysis between the inverse and time logarithm of oil shale in-situ conversion temperature.Using this method can quickly know the conversion of mining area oil shale in-situ in what temperature, when reach the conversion ratio of which kind of degree, provide important parameter for engineering decision and cost accounting.

Description

Oil shale in-situ conversion temperature, time, conversion ratio plate method for building up and system
Technical field
The invention belongs to oil shale recovery and research field, more particularly to a kind of oil shale in-situ conversion temperature, time, turn Rate plate method for building up and system.
Background technology
Oil shale is a kind of high ash-deposition rock containing flammable organic matter, and its organic matter is mainly kerogen/kerogen. In the case of completely cutting off air or oxygen, the organic matter (kerogen/kerogen) in oil shale can be pyrolyzed, and produce shale oil, destructive distillation Gas and solid carbonaceous residues and a small amount of pyrolysis water.Global oil shale resources very abundant, according to incompletely statistics its contain stock number 10,000,000,000,000 tons are there are about, than conventional petroleum stock number (271,000,000,000 tons) more more than 50%.
Oil shale exploitation at present mainly has two ways:Converted in-situ is exploited and ground destructive distillation.Ground destructive distillation refers to oily page Rock delivers to ground through strip mining transformation or underground mining, after crushing and screening to required granularity or lumpiness, into gas retort in heating it is dry Evaporate, produce shale oil gas and shale semicockes or shale lime-ash.Converted in-situ exploitation refers to bury oil shale in underground without opening Adopt, directly try to heat destructive distillation in underground, the shale oil gas of production is exported to ground again.Compared with the destructive distillation of ground, original position is opened Adopting has the advantages that to save strip mining transformation expense, reduces the destroyed degree of ground vegetation and floor space less.
For oil shale in-situ conversion exploitation, final temperature, heat time and the final oil gas conversion ratio of heating, The quality of oil gas and production cost are closely bound up.In the prior art obtain oil shale heating-up temperature-when m- conversion ratio between relation The time of method and cost are higher, such as Shell Co. Ltd passes through showing in 160 acres of section's roller state block in 2003~2005 years Field experiment, obtains relations of the block mahogany with oil shale heating-up temperature-heat time-conversion ratio.Except the scene of carrying out The technology of experiment, other belong to the in-situ retorting technology of conceptual phase, and also trial establishes oil shale with laboratory physical simulation experiment The quantitative relationship of Temperature-time-conversion ratio, although the information of some Temperature-time-conversion ratios can be obtained, indoor thing The problem of managing simulated experiment maximum is simulated time yardstick, and different from oil shale destructive distillation, destructive distillation conversion in situ will be by prolonging The long heat time realizes conversion at relatively low temperature, could reduce production cost, and the experiment of indoor short time can only pass through high temperature Oil shale conversion oil gas could be realized.Therefore, how fast quantification judge oil shale in-situ conversion production temperature-when m- conversion Relation between rate, establish oil shale in-situ conversion production temperature-when m- conversion ratio plate seem most important.
The content of the invention
One of technical problems to be solved by the invention are to need to provide a kind of fast quantification to judge that oil shale in-situ converts Production temperature-when m- conversion ratio between relation, establish oil shale in-situ conversion production temperature-when m- conversion ratio plate side Method.
In order to solve the above-mentioned technical problem, embodiments herein provide firstly a kind of oil shale in-situ conversion temperature, Time, conversion ratio plate method for building up, including intend production zone from oil shale and choose oil shale sample, and using rock pyrolysis point Analyzer is tested to the oil shale sample;Experimental result obtain pyrolysis activation energy, conversion ratio and frequency factor it Between relation data;Established based on the relation data when conversion ratio is different numerical value, oil shale in-situ conversion temperature is fallen Core analysis between number and time logarithm.
Preferably, the experimental result obtains the relation number between pyrolysis activation energy, conversion ratio and frequency factor According to, including:The delta data that yield is pyrolyzed under different heating rates is recorded using rock pyrolysis analysis instrument;Produced according to the pyrolysis The delta data of rate obtains the corresponding relation between pyrolysis Auditory steady-state responses and conversion ratio;The heat is based on using rock pyrolysis analysis instrument The delta data of solution yield is simulated the corresponding relation to determine pyrolysis activation energy and be pyrolyzed between Auditory steady-state responses;According to pyrolysed hydrocarbon The corresponding relation between corresponding relation and pyrolysis activation energy and pyrolysis Auditory steady-state responses between yield and conversion ratio establishes pyrolysis activation Can be with the relation curve of conversion ratio.
Preferably, simulated using delta data of the rock pyrolysis analysis instrument based on the pyrolysis yield to determine heat Also include determining between pyrolysis activation energy and frequency factor in the step of solving activation energy and the corresponding relation being pyrolyzed between Auditory steady-state responses Corresponding relation.
Preferably, it is described to be established based on the relation data when conversion ratio is different numerical value, oil shale in-situ conversion temperature Core analysis between the inverse and time logarithm of degree, including:According to the characteristic of the heat scission reaction with first order reaction feature Establish under the differential responses for the fixed extent of reaction or for same reaction degree, inverse and the logarithm of time of temperature Linear model;The relation data is substituted into the linear model, obtain when conversion ratio is different numerical value the inverse of temperature and when Between logarithm between Core analysis.
Preferably, the characteristic of the heat scission reaction with first order reaction feature is as shown in following formula:
Wherein, E represents pyrolysis activation energy, and T represents temperature, and t represents the time, and R is gas constant, and A is constant to be determined.
Preferably, when being tested using rock pyrolysis analysis instrument to the oil shale sample, including:By the oily page Rock sample comminution, make its particle diameter between 0.07mm~0.15mm;A certain amount of sample is weighed every time, is recorded respectively several times not With the measurement data under the rate of heat addition.
Preferably, when it is different numerical value that foundation, which obtains conversion ratio, inverse and the time of oil shale in-situ conversion temperature After Core analysis between logarithm, also using the Core analysis prediction oil shale in-situ conversion temperature, the time with And conversion ratio.
Embodiments herein additionally provides a kind of oil shale in-situ conversion temperature, time, conversion ratio plate establish system, Including:Sample analysis module, it is intended production zone from oil shale and chooses oil shale sample, and using rock pyrolysis analysis instrument to institute Oil shale sample is stated to be tested;Data acquisition module, its experimental result obtain pyrolysis activation energy, conversion ratio and frequency Relation data between the factor;Plate establishes module, and it is based on the relation data and established when conversion ratio is different numerical value, oil Core analysis between the inverse and time logarithm of shale in-situ conversion temperature.
Preferably, the plate is established module and established according to following steps when conversion ratio is different numerical value, and oil shale is former Core analysis between the inverse and the logarithm of time of position conversion temperature, according to the heat scission reaction with first order reaction feature Characteristic is established under the differential responses for the fixed extent of reaction or for same reaction degree, the inverse and pair of time of temperature Several linear models;The relation data is substituted into the linear model, obtains the inverse of temperature when conversion ratio is different numerical value Core analysis between the logarithm of time.
Preferably, the characteristic of the heat scission reaction with first order reaction feature is as shown in following formula:
Wherein, E represents pyrolysis activation energy, and T represents temperature, and t represents the time, and R is gas constant, and A is constant to be determined.
Compared with prior art, one or more of such scheme embodiment can have the following advantages that or beneficial to effect Fruit:
Kinetic model is obtained using the oil shale pyrolysis curve under different heating rates, and it is former according to chemical kinetics Reason, establish oil shale in-situ conversion temperature-when m- conversion ratio plate, can quickly know the conversion of mining area oil shale in-situ assorted Temperature, the conversion ratio which kind of degree when reached, important parameter is provided for engineering decision and cost accounting.
Other advantages, the target of the present invention, and feature will be illustrated in the following description to a certain extent, and And to a certain extent, based on will be apparent to those skilled in the art to investigating hereafter, Huo Zheke To be instructed from the practice of the present invention.The target and other advantages of the present invention can be wanted by following specification, right Specifically noted structure is sought in book, and accompanying drawing to realize and obtain.
Brief description of the drawings
Accompanying drawing is used for providing to the technical scheme of the application or further understanding for prior art, and constitution instruction A part.Wherein, the accompanying drawing for expressing the embodiment of the present application is used for the technical side for explaining the application together with embodiments herein Case, but do not form the limitation to technical scheme.
Fig. 1 be according to the oil shale in-situ conversion temperature of the embodiment of the present invention, the time, conversion ratio plate method for building up stream Journey schematic diagram;
Fig. 2 is according to typical the oil shale pyrolysis Activated factor Ⅶ and conversion ratio relation curve of an of the invention example;
Fig. 3 according to the typical oil shale in-situ of an of the invention example conversion temperature-when m- conversion ratio Core analysis;
Fig. 4 is the knot that system is established according to the oil shale in-situ conversion temperature of the embodiment of the present invention, time, conversion ratio plate Structure schematic diagram.
Embodiment
Embodiments of the present invention are described in detail below with reference to drawings and Examples, and how the present invention is applied whereby Technological means solves technical problem, and the implementation process for reaching relevant art effect can fully understand and implement according to this.This Shen Each feature that please be in embodiment and embodiment, can be combined with each other under the premise of not colliding, the technical scheme formed Within protection scope of the present invention.
Fig. 1 be according to the oil shale in-situ conversion temperature of the embodiment of the present invention, the time, conversion ratio plate method for building up stream Journey schematic diagram, as illustrated, this method comprises the following steps:
Step S110, intend production zone from oil shale and choose oil shale sample, and using rock pyrolysis analysis instrument to oily page Rock sample product are tested.
Step S120, experimental result obtains the relation data between pyrolysis activation energy, conversion ratio and frequency factor.
Step S130, established based on above-mentioned relation data when conversion ratio is different numerical value, oil shale in-situ conversion temperature Inverse and the logarithm of time between Core analysis.
Specifically, in step s 110, the typical oil shale sample comminution that production zone is intended in selection is standby, chooses sample Quantity and the degree of crushing and processing determine according to requirement of experiment.For example, at least 20g oil shale sample is chosen, is crushed The powder for being particle diameter between 0.07mm~0.15mm is standby.During experiment, 5g samples are weighed every time, utilize the types of Rock-Eval 6 Rock pyrolysis analysis instrument records to experimental result.
In the step s 120, the change number that yield is pyrolyzed under different heating rates is recorded first with rock pyrolysis analysis instrument According to.For example, it is that hydrocarbon products are pyrolyzed yield delta data under 5,15 and 25 DEG C/min to record the rate of heat addition respectively.Then, according to heat The delta data for solving yield obtains the corresponding relation being pyrolyzed between Auditory steady-state responses and conversion ratio.Rock pyrolysis analysis instrument is recycled to be based on The delta data of pyrolysis yield is simulated the corresponding relation to determine pyrolysis activation energy and be pyrolyzed between Auditory steady-state responses.For example, will The result of the same rate of heat addition synthesizes a pyrolysis yield curve, utilizes the type rock pyrolysis analysis instrument dynamics of Rock-Eval 6 Simulation softwardSimulated, by the pyrolysis data input of 3 different heating rates into dynamics simulation software, led to Cross the relation that simulation obtains the pyrolysis Activated factor Ⅶ of the oil shale and is pyrolyzed between Auditory steady-state responses.Finally, according to pyrolysis Auditory steady-state responses Corresponding relation between conversion ratio and pyrolysis activation energy and the corresponding relation being pyrolyzed between Auditory steady-state responses establish pyrolysis activation energy with The relation curve of conversion ratio.
In addition, simulated using delta data of the rock pyrolysis analysis instrument based on pyrolysis yield to determine pyrolysis activation When energy and the corresponding relation being pyrolyzed between Auditory steady-state responses, pair being pyrolyzed between activation energy and frequency factor is also determined simultaneously by simulation It should be related to.
Next in step s 130, established according to the characteristic of the heat scission reaction with first order reaction feature for fixation The extent of reaction or differential responses for same reaction degree under, the linear model of the inverse and the logarithm of time of temperature.
Specifically, it is approximately tool first order reaction feature by the process of oil shale organic matter (kerogen) evolution generation oil gas Heat scission reaction, that is, the speed reacted is directly proportional to the first power of reactant concentration, as shown in expression formula (1):
In formula, t is the reaction time, and unit s, C are the concentration of reactant, and k is rate constant.
Reaction rate is adapted to Arrhenius equation with temperature and the relation of activation energy, as shown in expression formula (2):
K=k0e-E/RT (2)
In formula, k0For frequency molecule, E is pyrolysis activation energy, and R is gas constant, and T is absolute temperature
Expression formula (1) is integrated, and expression formula (2) is substituted into and can obtained:
In formula, C0It is the concentration of (t=0) reactant when reaction starts, C is the concentration in moment t reactant.
Expression formula (3) both sides are taken the logarithm, obtained:
For the extent of reaction of a certain fixation, or the differential responses for same reaction degree, C are also constant, therefore, table It can be expressed as up to formula (4):
In formula, A is constant to be determined, i.e., the inverse (1/T) of temperature has linear relationship with the logarithm (lnt) of time.
Further, by the relation data between the pyrolysis activation energy, conversion ratio and the frequency factor that are obtained in step S120 The linear model as shown in expression formula (5) is substituted into, obtains the inverse (1/T) of temperature and pair of time when conversion ratio is different numerical value Core analysis between number (lnt).
The relational expression of the inverse (1/T) of temperature under differentiated yields and the logarithm (lnt) of time is put into same Among grid board, that is, establish oil shale in-situ conversion temperature-when m- conversion ratio plate.
If prediction judges what temperature the oil shale in-situ converts in, when reaches great conversion ratio, only need true Just cultellation would know that fixed two factors in plate.
The plate method for building up of the embodiment of the present invention, obtained using not plate with the oil shale pyrolysis curve under heating rate Kinetic model, oil shale pyrolysis activation energy (E) distribution and frequency factor (A) are determined, according to principles of chemical kinetics, establish oil Shale in-situ conversion temperature-when m- conversion ratio plate, can quickly know the conversion of mining area oil shale in-situ in what temperature, assorted Time reaches the conversion ratio of which kind of degree, and important parameter is provided for engineering decision and cost accounting.
The embodiment of the present invention can effectively predict institute under differentiated yields in view of the shortcomings of the prior art, in order to establish one kind The temperature and time needed, exploit decision-making for oil shale in-situ thermal conversion and science parameter is provided, the present invention is according to chemomotive force Learn principle, there is provided a kind of oil shale in-situ conversion temperature-when m- conversion ratio plate method for building up, can quickly obtain not With temperature-time matching relationship required during conversion ratio.
The embodiment of the present invention is further described below by an example.
Present example provides resources in Maoming Basin phyllanthus emblica nest group oil shale in-situ conversion temperature-when m- conversion ratio plate established Journey.
Weigh the resources in Maoming Basin phyllanthus emblica preferable oil shale sample 20g of nest group lithologic homogeneity, crush, particle diameter 0.07mm~ Between 0.15mm.
5g samples are weighed every time, different heating rates (5,15 Hes are recorded using the type rock pyrolysis analysis instrument of Rock-Eval 6 25 DEG C/min) under be pyrolyzed yield change, be as shown in table 15 DEG C/min of the rate of heat addition pyrolysed hydrocarbon yield data, similarly, obtain The rate of heat addition 15 and pyrolysed hydrocarbon yield data during 25 DEG C/min, as shown in table 1:
Table 1 is pyrolyzed Auditory steady-state responses change (5 DEG C/min of the rate of heat addition)
The pyrolysis data of 3 different heating rates are inputted, utilize the type rock pyrolysis analysis instrument kinetic simulations of Rock-Eval 6 Intend softwarePyrolysis activation energy (E) distribution and the frequency factor (A) of the oil shale are obtained, as shown in table 2:
Table 2 is pyrolyzed Activated factor Ⅶ and frequency factor
The relation of Auditory steady-state responses and conversion ratio is known that by table 1, the relation of Auditory steady-state responses and activation energy is known that by table 2, Using Auditory steady-state responses as bridge, the corresponding data of activation energy and conversion ratio is drawn, establishes activation energy and conversion rate curve, as shown in Figure 2.
Process by oil shale organic matter (kerogen) evolution generation oil gas is approximately the thermal cracking for having first order reaction feature Reaction,
For the extent of reaction of a certain fixation, or the differential responses for same reaction degree, the inverse (1/T) of temperature with The logarithm (lnt) of time has the linear relationship as shown in expression formula (5).
Expression is substituted into using the conversion ratio in table 1, table 2, activation energy, frequency factor data, R (8.3144J/ (mol*K)) Formula (5), obtain the inverse (1/T) of temperature and the relational expression of the logarithm (lnt) of time under differentiated yields;Such as individually below When to be conversion ratio be 10%, 60% and 90%, the relational expression of the inverse (1/T) of temperature and the logarithm (lnt) of time:
Conversion ratio 10%, required activation energy=238kj/mol, frequency factor=2.565E+16moll-1s-1, Then lnt=28625/T-37.783.
Conversion ratio 60%, required activation energy=264kj/mol, frequency factor=2.565E+166moll-1s-1, Then lnt=31752/T-37.783.
Conversion ratio 90%, required activation energy=280kj/mol, frequency factor=2.565E+166moll-1s-1, Then lnt=33677/T-37.783.
The relational expression of the inverse (1/T) of temperature under differentiated yields and the logarithm (lnt) of time is put into same Among grid board, that is, establish oil shale in-situ conversion temperature-when m- conversion ratio plate, as shown in figure 3, in figure T unit For K, t unit is s.
If prediction judges what temperature the oil shale in-situ converts in, when reaches great conversion ratio, only need true Just cultellation would know that (as shown in table 3) fixed two factors in plate.
The converted in-situ Temperature-time of table 3-conversion ratio relation
Such as understand that if being heated to 350 DEG C exploits area's oil shale from Fig. 3, table 3, conversion ratio needs 0.97 small up to 10% When, conversion ratio needs 6 days up to 60%, and conversion ratio needs 4.5 months up to 90%.If plan is heated to 200 DEG C and exploits area's oil page Rock, it can be seen that be not-so-practical thing, even if only conversion 10% is also required to 240 years.On the basis of these data parameters, so that it may To calculate the output and input of intending input engineering, rational production temperature is selected, is that scientific and reasonable decision-making and economic benefit are maximum Change and important evidence is provided.
Fig. 4 is the knot that system is established according to the oil shale in-situ conversion temperature of the embodiment of the present invention, time, conversion ratio plate Structure schematic diagram, as illustrated, the system includes:
Sample analysis module 41, it performs step S110 operation, and data acquisition module 42, it performs step S120 behaviour Make, plate establishes module 43, and it performs step S130 operation.Here is omitted.
Although disclosed herein embodiment as above, described content only to facilitate understand the present invention and adopt Embodiment, it is not limited to the present invention.Any those skilled in the art to which this invention pertains, this is not being departed from On the premise of the disclosed spirit and scope of invention, any modification and change can be made in the implementing form and in details, But the scope of patent protection of the present invention, still should be subject to the scope of the claims as defined in the appended claims.

Claims (10)

1. a kind of oil shale in-situ conversion temperature, time, conversion ratio plate method for building up, including:
Intend production zone from oil shale and choose oil shale sample, and the oil shale sample is carried out using rock pyrolysis analysis instrument Experiment;
Experimental result obtains the relation data between pyrolysis activation energy, conversion ratio and frequency factor;
Established based on the relation data when conversion ratio is different numerical value, inverse and the time pair of oil shale in-situ conversion temperature Core analysis between number.
2. according to the method for claim 1, it is characterised in that the experimental result obtains pyrolysis activation energy, turned Relation data between rate and frequency factor, including:
The delta data that yield is pyrolyzed under different heating rates is recorded using rock pyrolysis analysis instrument;
The corresponding relation between pyrolysis Auditory steady-state responses and conversion ratio is obtained according to the delta data of the pyrolysis yield;
Simulated using delta data of the rock pyrolysis analysis instrument based on the pyrolysis yield to determine pyrolysis activation energy and heat Solve the corresponding relation between Auditory steady-state responses;
According to the corresponding relation between pyrolysis Auditory steady-state responses and conversion ratio and pyrolysis activation energy and the corresponding pass being pyrolyzed between Auditory steady-state responses System establishes the relation curve of pyrolysis activation energy and conversion ratio.
3. according to the method for claim 2, it is characterised in that be based on the pyrolysis yield using rock pyrolysis analysis instrument Delta data the step of being simulated to determine pyrolysis activation energy and the corresponding relation being pyrolyzed between Auditory steady-state responses in also include it is true Surely the corresponding relation being pyrolyzed between activation energy and frequency factor.
4. according to the method in any one of claims 1 to 3, it is characterised in that described to be established based on the relation data When conversion ratio is different numerical value, the Core analysis between the inverse and time logarithm of oil shale in-situ conversion temperature, including:
Established according to the characteristic of the heat scission reaction with first order reaction feature for the fixed extent of reaction or for identical anti- Under the differential responses for answering degree, the linear model of the inverse and the logarithm of time of temperature;
The relation data is substituted into the linear model, obtains the inverse and pair of time of temperature when conversion ratio is different numerical value Core analysis between number.
5. according to the method for claim 4, it is characterised in that there is the characteristic of the heat scission reaction of first order reaction feature such as Shown in following formula:
<mrow> <mi>ln</mi> <mi>t</mi> <mo>=</mo> <mfrac> <mi>E</mi> <mrow> <mi>R</mi> <mi>T</mi> </mrow> </mfrac> <mo>-</mo> <mi>A</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mi>T</mi> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mi>E</mi> <mi>R</mi> </mfrac> <mo>-</mo> <mi>A</mi> </mrow>
Wherein, E represents pyrolysis activation energy, and T represents temperature, and t represents the time, and R is gas constant, and A is constant to be determined.
6. method according to any one of claim 1 to 5, it is characterised in that using rock pyrolysis analysis instrument to institute When stating oil shale sample and being tested, including:
By the oil shale sample comminution, make its particle diameter between 0.07mm~0.15mm;
A certain amount of sample is weighed every time, records the measurement data under different heating rates respectively several times.
7. according to the method for claim 1, it is characterised in that when it is different numerical value that foundation, which obtains conversion ratio, oil shale After Core analysis between the inverse and the logarithm of time of converted in-situ temperature, also oil is predicted using the Core analysis Temperature, time and the conversion ratio of shale in-situ conversion.
8. a kind of oil shale in-situ conversion temperature, time, conversion ratio plate establish system, including:
Sample analysis module, it is intended production zone from oil shale and chooses oil shale sample, and using rock pyrolysis analysis instrument to institute Oil shale sample is stated to be tested;
Data acquisition module, its experimental result obtain the relation number between pyrolysis activation energy, conversion ratio and frequency factor According to;
Plate establishes module, and it is based on the relation data and established when conversion ratio is different numerical value, oil shale in-situ conversion temperature Core analysis between the inverse and time logarithm of degree.
9. system according to claim 8, it is characterised in that the plate, which establishes module and established to work as according to following steps, to be turned When rate is different numerical value, the Core analysis between the inverse and the logarithm of time of oil shale in-situ conversion temperature,
Established according to the characteristic of the heat scission reaction with first order reaction feature for the fixed extent of reaction or for identical anti- Under the differential responses for answering degree, the linear model of the inverse and the logarithm of time of temperature;
The relation data is substituted into the linear model, obtains the inverse and pair of time of temperature when conversion ratio is different numerical value Core analysis between number.
10. system according to claim 9, it is characterised in that there is the characteristic of the heat scission reaction of first order reaction feature As shown in following formula:
<mrow> <mi>ln</mi> <mi>t</mi> <mo>=</mo> <mfrac> <mi>E</mi> <mrow> <mi>R</mi> <mi>T</mi> </mrow> </mfrac> <mo>-</mo> <mi>A</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mi>T</mi> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mi>E</mi> <mi>R</mi> </mfrac> <mo>-</mo> <mi>A</mi> </mrow>
Wherein, E represents pyrolysis activation energy, and T represents temperature, and t represents the time, and R is gas constant, and A is constant to be determined.
CN201610811514.4A 2016-09-08 2016-09-08 Method and system for establishing oil shale in-situ conversion temperature, time and conversion rate chart Active CN107808019B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610811514.4A CN107808019B (en) 2016-09-08 2016-09-08 Method and system for establishing oil shale in-situ conversion temperature, time and conversion rate chart

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610811514.4A CN107808019B (en) 2016-09-08 2016-09-08 Method and system for establishing oil shale in-situ conversion temperature, time and conversion rate chart

Publications (2)

Publication Number Publication Date
CN107808019A true CN107808019A (en) 2018-03-16
CN107808019B CN107808019B (en) 2021-06-15

Family

ID=61576240

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610811514.4A Active CN107808019B (en) 2016-09-08 2016-09-08 Method and system for establishing oil shale in-situ conversion temperature, time and conversion rate chart

Country Status (1)

Country Link
CN (1) CN107808019B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106940331A (en) * 2017-03-21 2017-07-11 山东京博石油化工有限公司 Intermittent liquid phase bulk polymerization prepares the measuring method of polyacrylic polymerisation conversion
CN111610119A (en) * 2019-02-25 2020-09-01 中国石油天然气股份有限公司 Crude oil activation energy determination method
CN112632738A (en) * 2019-10-09 2021-04-09 中国石油天然气股份有限公司 Method and device for determining hydrocarbon generation activation energy in shale sample in-situ conversion
CN113756779A (en) * 2021-08-30 2021-12-07 陕西省煤田地质集团有限公司 Method for monitoring and evaluating degree of in-situ pyrolysis reaction of oil-rich coal
CN115221735A (en) * 2022-09-20 2022-10-21 中国石油大学(华东) Scale increasing method for numerical simulation of in-situ conversion of electrically heated oil shale

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101045667A (en) * 2006-03-31 2007-10-03 中国石油化工股份有限公司 Combined catforming of high-production low carbon alkene
US20100252483A1 (en) * 2006-11-06 2010-10-07 Omer Refa Koseoglu Process for removal of nitrogen and poly-nuclear aromatics from fcc feedstocks
CN104007083A (en) * 2014-06-11 2014-08-27 中国石油大学(北京) Pyrolysis measurement method of oil shale based on light tracking
CN104156593A (en) * 2014-08-08 2014-11-19 中国石油大学(华东) Shale oil and gas yield evaluation model building and parameter calibrating method applied under closed system
CN104849172A (en) * 2014-02-18 2015-08-19 中国石油化工股份有限公司 Oil shale oil content test experiment apparatus and test method thereof

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101045667A (en) * 2006-03-31 2007-10-03 中国石油化工股份有限公司 Combined catforming of high-production low carbon alkene
US20100252483A1 (en) * 2006-11-06 2010-10-07 Omer Refa Koseoglu Process for removal of nitrogen and poly-nuclear aromatics from fcc feedstocks
CN104849172A (en) * 2014-02-18 2015-08-19 中国石油化工股份有限公司 Oil shale oil content test experiment apparatus and test method thereof
CN104007083A (en) * 2014-06-11 2014-08-27 中国石油大学(北京) Pyrolysis measurement method of oil shale based on light tracking
CN104156593A (en) * 2014-08-08 2014-11-19 中国石油大学(华东) Shale oil and gas yield evaluation model building and parameter calibrating method applied under closed system

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
FENGTIAN BAI 等: "Thermal and kinetic characteristics of pyrolysis and combustion of three oil shales", 《ENERGY CONVERSION AND MANAGEMENT》 *
关德范 等: "烃源岩生排烃理论研究与泥页岩油气", 《中外能源》 *
徐峰: "油页岩热解及燃烧反应的化学动力学研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 *
杨帆: "油页岩的热转化利用", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106940331A (en) * 2017-03-21 2017-07-11 山东京博石油化工有限公司 Intermittent liquid phase bulk polymerization prepares the measuring method of polyacrylic polymerisation conversion
CN106940331B (en) * 2017-03-21 2019-09-06 山东京博石油化工有限公司 Intermittent liquid phase bulk polymerization prepares the measurement method of polyacrylic polymerisation conversion
CN111610119A (en) * 2019-02-25 2020-09-01 中国石油天然气股份有限公司 Crude oil activation energy determination method
CN112632738A (en) * 2019-10-09 2021-04-09 中国石油天然气股份有限公司 Method and device for determining hydrocarbon generation activation energy in shale sample in-situ conversion
CN113756779A (en) * 2021-08-30 2021-12-07 陕西省煤田地质集团有限公司 Method for monitoring and evaluating degree of in-situ pyrolysis reaction of oil-rich coal
CN115221735A (en) * 2022-09-20 2022-10-21 中国石油大学(华东) Scale increasing method for numerical simulation of in-situ conversion of electrically heated oil shale

Also Published As

Publication number Publication date
CN107808019B (en) 2021-06-15

Similar Documents

Publication Publication Date Title
CN107808019A (en) Oil shale in-situ conversion temperature, time, conversion ratio plate method for building up and system
Daggupati et al. Laboratory studies on cavity growth and product gas composition in the context of underground coal gasification
Sajdak et al. Use of plastic waste as a fuel in the co-pyrolysis of biomass: part III: optimisation of the co-pyrolysis process
Li et al. Compaction property prediction of mixed gangue backfill materials using hybrid intelligence models: A new approach
CN1575374A (en) Seismic monitoring of in situ conversion in a hydrocarbon containing formation
Zhang et al. Kinetic parameters estimation and reaction model modification for thermal degradation of Beizao oil shale based on thermogravimetric analysis coupled with deconvolution procedure
Meijssen et al. In-situ oil shale development in jordan through icp technology
CN107688669A (en) A kind of oil shale in-situ exploits converting rate Forecasting Methodology
Han et al. Heat transfer characteristic of oil shale particle during the retorting
Samdani et al. Kinetics of heterogeneous reactions with coal in context of underground coal gasification
Khalidy et al. Assessment of geochemical modeling applications and research hot spots—a year in review
Guo et al. Mechanism and reservoir simulation study of the autothermic pyrolysis in-situ conversion process for oil shale recovery
Gao et al. Numerical and experimental investigation of production performance of in-situ conversion of shale oil by air injection
Qiao et al. Research on cobalt-rich crust cutting modes by deep-sea mining vehicle and cutting performance evaluation
Zhang et al. An efficient modeling approach to simulate heat transfer rate between fracture and matrix regions for oil shale retorting
Lakshmikanthan et al. Studies on characterization of mechanically biologically treated waste from Bangalore city
Knyazeva et al. A two-phase model of shale pyrolysis
Bandyopadhyay et al. Thermal deactivation studies of coconut shell pyrolysis
Gregg et al. Laboratory and modelling investigation of a Colorado oil-shale block heated to 900 C
Sennoune et al. Toward the control of the smoldering front in the reaction-trailing mode in oil shale semicoke porous media
Wang Multicriteria sustainability ranking of biohydrogen systems
Wang et al. Spontaneous combustion of coal
Qu et al. Simulation Study of Hydrocarbon Generation Characteristics in Siliceous Rock
Kačur et al. Investigation of underground coal gasification in laboratory conditions: A review of recent research
Selormey et al. Development of mathematical model for predicting methane-to-carbon dioxide proportion in anaerobic biodegradability of cattle blood and rumen content

Legal Events

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