CN105426666B - Gas hydrate dissociation rate of gas release computational methods and its device - Google Patents
Gas hydrate dissociation rate of gas release computational methods and its device Download PDFInfo
- Publication number
- CN105426666B CN105426666B CN201510744846.0A CN201510744846A CN105426666B CN 105426666 B CN105426666 B CN 105426666B CN 201510744846 A CN201510744846 A CN 201510744846A CN 105426666 B CN105426666 B CN 105426666B
- Authority
- CN
- China
- Prior art keywords
- hydrate
- gas
- decomposition
- model
- mol
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Organic Low-Molecular-Weight Compounds And Preparation Thereof (AREA)
Abstract
The present invention provides a kind of gas hydrate dissociation rate of gas release computational methods and device, the method to include：According to decomposition of hydrate theory basis, establish hydrate slurry decomposition kinetics model and be fitted relevant parameter numerical value in decomposition of hydrate model；Phase flash distillation is carried out according to gas component, calculates and obtains natural gas relevant parameter；Real time temperature pressure data after the amount of the substance of gas hydrates and decomposition reaction start when acquisition gas hydrate dissociation is initial；According to the real time temperature pressure data of hydrate slurry decomposition kinetics model, decomposition of hydrate model parameter numerical value, natural gas related property parameter, decomposition of hydrate after the amount of the substance of gas hydrates and decomposition reaction start when initial, gas hydrate dissociation rate of gas release is obtained.Present invention is mainly applied to the directions such as Gas Hydrate exploitation, the removal of flow assurance gas hydrates frozen block, gas hydrates storaging and transport technology and hydrate separation.
Description
Technical field
The present invention relates to fluid calculation field, espespecially a kind of calculating for calculating gas hydrate dissociation rate of gas release
Method and its device are mainly used in Gas Hydrate, it can be achieved that the prediction of gas hydrates rate of gas release
The directions such as exploitation and flow assurance, the removal of gas hydrates frozen block, gas hydrates storaging and transport technology and hydrate separation.
Background technology
With the development of petroleum industry, the fields such as drilling well, oil recovery and oil-gas transportation extend to deep-sea.In offshore oil work
In industry, oil gas is post-processed after offshore mining comes out by well head through pipeline to FPSO, offshore platform or land.Sea
Bottom pipeline especially seabed gathering line, the medium of conveying is mostly without separating treatment, usually oil gas water multiphase flow pipeline.
Submarine pipeline operating pressure is higher, and running temperature is relatively low, is easier to form hydrate in pipeline.Hydrate is in deep-sea pit shaft and pipeline
Middle generation, decomposition, so as to the trend for hindering flowing.In the case where hydrate generation is more serious, flow can be significantly
It reduces or even stops production.In addition, the formation and decomposition of hydrate, can influence to transport positioned at the valve in seabed and the safety of relevant device
Row.
During drilling well, oil recovery and oil-gas transportation, the formation of gas hydrates and decompose to pit shaft, submarine pipeline and
It is influenced caused by undersea device etc. and its requirement to technology and cost is increasingly taken seriously.Due to hydrate traditional at present
Suppression technology has a cost expenses of apparent technical limitation and great number in ocean petroleum developing, hydrate slurry conveying technology by
It is gradually of interest by petroleum industry and domestic and international researcher.In the oil and gas multiphase flow system of multiphase conveying, hydrate can be with pipe
The variation of pressure and temperature in road and formed or decomposed.It accurately to predict the hydrate amount of having in pipeline, not only need to count
Calculate the synthesis speed of hydrate, it is also necessary to grasp the Forecasting Methodology of decomposition of hydrate rate.At present, it is mixed both at home and abroad for multiphase
The research of decomposition of hydrate mechanism is less in defeated system, some mostly existing decomposition researchs are for porous media or reaction kettle
The decomposition of middle hydrate, and it is fewer and fewer for the decomposition of hydrate research in multi-phase mixed delivering system under flow regime.
Invention content
It is a kind of for calculating the skill of rate of gas release during gas hydrate dissociation present invention aims at providing
Art, i.e., under to multi-phase mixed delivering system in High-Voltage Experimentation circuit on the basis of the decomposition mechanism research of hydrate, Binding experiment number
According to gas hydrate dissociation kinetic model is established, so as to calculate gas hydrates rate of gas release.
In order to achieve the above object, the present invention specifically provides a kind of gas hydrate dissociation rate of gas release computational methods,
The method includes：According to the theory basis of decomposition of hydrate, establish hydrate slurry decomposition kinetics model and be fitted
Relevant parameter numerical value in decomposition of hydrate model；Phase flash distillation is carried out according to gas component, calculates and obtains natural gas correlation ginseng
Number；Real time temperature after the amount of the substance of gas hydrates and decomposition reaction start when acquisition gas hydrate dissociation is initial
Pressure data；According to relevant parameter numerical value, the natural gas phase in the hydrate slurry decomposition kinetics model, the model
Real-time temperature after the amount of the substance of gas hydrates and the decomposition reaction start when related parameter, the decomposition of hydrate are initial
Pressure data is spent, obtains gas hydrate dissociation rate of gas release.
In above-mentioned gas hydrate dissociation rate of gas release computational methods, it is preferred that described according to hydrate point
The theory basis of solution is established hydrate slurry decomposition kinetics model and is included：According to the intrinsic kinetics of decomposition of hydrate
And the influence of heat transfer factor, mass transfer factor to decomposition of hydrate rate, it establishes comprising intrinsic, heat transfer, mass transfer and coalescence characteristic
Hydrate slurry decomposition kinetics model.
In above-mentioned gas hydrate dissociation rate of gas release computational methods, it is preferred that described according to hydrate point
The theory basis of solution is fitted relevant parameter numerical value in decomposition of hydrate model and includes：Intended using high pressure circuit experimental data
Close the relevant parameters numerical value such as reaction constant, carry-over factor and the activation energy of eigen decomposition reaction.In above-mentioned gas hydrates
It decomposes in rate of gas release computational methods, it is preferred that described to be included according to gas component progress phase flash distillation：According to natural
The component of natural gas of gas, and phase flash distillation is carried out based on thermodynamics Phase Equilibrium theory.
In above-mentioned gas hydrate dissociation rate of gas release computational methods, it is preferred that the natural gas correlation ginseng
Number is including at least under the physical parameter table of natural gas, hydrate formation curve, three-phase equilibrium fugacity and different temperatures pressure state
Fugacity.
In above-mentioned gas hydrate dissociation rate of gas release computational methods, it is preferred that the hydrate slurry point
It solves kinetic model and includes the following formula：
In above formula, nd(t)For the gas burst size of t moment, mol；nd(end)Gas consumption at the end of for hydrate formation
Amount, mol；T be the decomposition of hydrate time, s；Kd0For the rate constant of eigen decomposition reaction, mol/ (MPasm2)；T is day
Right gas hydrate decomposing solution real time temperature, K；R is mol gas constant, 8.314J/ (molK)；AsFor decomposition of hydrate
Total surface area, m2；fo,iFor gas phase fugacity of the component a certain under experiment condition at hydrate shell outer boundary, MPa；fe,iIt is a certain
The hydrate phase balance fugacity of component, MPa；α is carry-over factor；Δ E be activation energy, J/mol.
The present invention also provides a kind of gas hydrate dissociation rate of gas release computing device, described device includes：It builds
Form unit for the theory basis according to decomposition of hydrate, establishes hydrate slurry decomposition kinetics model and is fitted mould
Relevant parameter numerical value in type；Flash evaporation unit for carrying out phase flash distillation according to gas component, calculates and obtains natural gas correlation ginseng
Number；Collecting unit, for obtain gas hydrate dissociation it is initial when gas hydrates substance amount and decomposition reaction open
Real time temperature pressure data after beginning；Data processing unit, for according to the hydrate slurry decomposition kinetics model, described
Natural gas is hydrated when decomposition of hydrate Parameters in Mathematical Model numerical value, the natural gas relevant parameter, the decomposition of hydrate are initial
The amount of the substance of object and the decomposition reaction start after real time temperature pressure data, obtain gas hydrate dissociation gas release
Put rate.
In above-mentioned gas hydrate dissociation rate of gas release computing device, it is preferred that the modeling unit is also wrapped
Containing fitting module, the fitting module is used for according to relevant parameter number in high pressure circuit fitting experimental data decomposition of hydrate model
Value.
In above-mentioned gas hydrate dissociation rate of gas release computing device, it is preferred that the modeling unit includes：
According to the theory basis of decomposition of hydrate, following hydrate slurry decomposition kinetics model is established：
In above formula, nd(t)For the gas burst size of t moment, mol；nd(end)Gas consumption at the end of for hydrate formation
Amount, mol；T be the decomposition of hydrate time, s；Kd0For the rate constant of eigen decomposition reaction, mol/ (MPasm2)；T is day
Right gas hydrate decomposing solution real time temperature, K；R is mol gas constant, 8.314J/ (molK)；AsFor decomposition of hydrate
Total surface area, m2；fo,iFor gas phase fugacity of the component a certain under experiment condition at hydrate shell outer boundary, MPa；fe,iIt is a certain
The hydrate phase balance fugacity of component, MPa；α is carry-over factor；Δ E be activation energy, J/mol.
The advantageous effects of the present invention are：It can be realized to gas during gas hydrate dissociation using this method
The calculating and prediction of rate of release, in this way based on the software write can contribute to Monitoring Gas Pipeline safety and steady
Operation can also be used as and carry out the exploitation of hydrate ore deposit in a deep going way, using hydrate as medium storing with transporting gas, purifying with detaching
The basis of the technical research such as gas and application hydrate slurry conveying gas.
Description of the drawings
Attached drawing described herein is used to provide further understanding of the present invention, and forms the part of the application, not
Form limitation of the invention.In the accompanying drawings：
Fig. 1 is gas hydrate dissociation rate of gas release computational methods flow chart provided by the present invention；
Fig. 2 is gas hydrate dissociation rate of gas release computing device schematic diagram provided by the present invention；
Fig. 3 is hydrate formation curve of the obtained natural gas of phase flash calculation under three kinds of different phases；
Fig. 4 A be inhibitor concentration 3% moisture content, 15% pressure 6.0MPa when, it is different in flow rate under decomposition of hydrate gas
Burst size analogue data and experimental data comparison schematic diagram；
When Fig. 4 B are inhibitor concentration 3% moisture content, 15% flow velocity 0.91m/s, the decomposition of hydrate under different initial pressures
Gas burst size analogue data and experimental data comparison schematic diagram；
When Fig. 4 C are 3% flow velocity 0.91m/s pressure 5.0MPa of inhibitor concentration, the decomposition of hydrate gas under different water cut
Body burst size analogue data and experimental data comparison schematic diagram；
When Fig. 4 D are 20% flow velocity 0.91m/s pressure 5.0MPa of moisture content, the decomposition of hydrate under different inhibitor concentrations
Gas burst size analogue data and experimental data comparison schematic diagram.
Specific embodiment
Purpose, technical scheme and advantage to make the embodiment of the present invention are more clearly understood, with reference to embodiment and attached
Figure, is described in further details the present invention.Here, the illustrative embodiments of the present invention and their descriptions are used to explain the present invention,
It is but not as a limitation of the invention.
Shown in please referring to Fig.1, the present invention specifically provides a kind of gas hydrate dissociation rate of gas release computational methods,
The method includes：S101 establishes hydrate slurry decomposition kinetics model simultaneously according to the theory basis of decomposition of hydrate
It is fitted relevant parameter numerical value in decomposition of hydrate model；S102 carries out phase flash distillation according to gas component, calculates and obtains naturally
Gas phase related parameter；The amount of the substance of gas hydrates and decomposition reaction start when S103 acquisition gas hydrate dissociations are initial
Real time temperature pressure data afterwards；S104 is according to the hydrate slurry decomposition kinetics model, the decomposition of hydrate model
The amount of the substance of gas hydrates and described when parameter values, the natural gas relevant parameter, the decomposition of hydrate are initial
Decomposition reaction start after real time temperature pressure data, obtain gas hydrate dissociation rate of gas release.
In above-mentioned steps S101, the theory basis according to decomposition of hydrate establishes hydrate slurry decomposition
Kinetic model specifically includes：According to the intrinsic kinetics of decomposition of hydrate and heat transfer factor, mass transfer factor to decomposition of hydrate
The hydrate slurry decomposition kinetics model for including intrinsic, heat transfer, mass transfer and coalescence characteristic is established in the influence of rate.
In above-mentioned steps S101, the theory basis according to decomposition of hydrate is fitted decomposition of hydrate model
Middle relevant parameter numerical value includes：The rate constant K reacted according to high pressure circuit fitting experimental data eigen decompositiond0, carry-over factor
The parameter values such as α, activation energy Δ E.
It is described to be included according to gas component progress phase flash distillation in above-mentioned steps S102：According to the makings of natural gas
Component, and phase flash distillation is carried out based on thermodynamics Phase Equilibrium theory.
In the above-described embodiments, the natural gas relevant parameter includes at least the physical parameter table of natural gas, hydrate life
Fugacity under curve, three-phase equilibrium fugacity and different temperatures pressure state.Wherein, the hydrate slurry decomposition kinetics mould
Type includes the following formula：
In above formula, nd(t)For the gas burst size of t moment, mol；nd(end)Gas consumption at the end of for hydrate formation
Amount, mol；T be the decomposition of hydrate time, s；Kd0For the rate constant of eigen decomposition reaction, mol/ (MPasm2)；T is day
Right gas hydrate decomposing solution real time temperature, K；R is mol gas constant, 8.314J/ (molK)；AsFor decomposition of hydrate
Total surface area, m2；fo,iFor gas phase fugacity of the component a certain under experiment condition at hydrate shell outer boundary, MPa；fe,iIt is a certain
The hydrate phase balance fugacity of component, MPa；α is carry-over factor；Δ E be activation energy, J/mol.
Shown in please referring to Fig.2 again, the present invention also provides a kind of gas hydrate dissociation rate of gas release computing device,
Described device includes：Modeling unit 201 for the theory basis according to decomposition of hydrate, establishes hydrate slurry decomposition
Relevant parameter numerical value in kinetic model and model of fit；Flash evaporation unit 202, for carrying out phase sudden strain of a muscle according to gas component
It steams, calculates and obtain natural gas relevant parameter；Collecting unit 203, for obtain gas hydrate dissociation it is initial when gas water
Close object substance amount and decomposition reaction start after real time temperature pressure data；Data processing unit 204, for according to
Relevant parameter numerical value, the natural gas relevant parameter, the hydrate in hydrate slurry decomposition kinetics model, the model
Real time temperature pressure data after the amount of the substance of gas hydrates and the decomposition reaction start when decomposing initial obtains day
Right gas hydrate decomposes rate of gas release.
In a preferred embodiment of the invention, the modeling unit 201 also comprising fitting module, use by the fitting module
According to high pressure circuit fitting experimental data relevant parameter numerical value.
In the above-described embodiments, the modeling unit includes：According to the theory basis of decomposition of hydrate, establish following
Hydrate slurry decomposition kinetics model：
Wherein, correlation function is as defined above, and this will not be detailed here.
Below by way of specific example, above-mentioned gas hydrate dissociation rate of gas release computational methods are made furtherly
Bright, relevant technical staff in the field is when it is found that it is only to help to understand that technical scheme of the present invention does not protect model to the present invention
It encloses and limits.
First, the decomposition of hydrate kinetic model for considering intrinsic, heat transfer, mass transfer and coalescence characteristic is established；
The present invention on the basis of the microcosmic resolution characteristic of the defeated system hydrate slurry of pipe and decomposition kinetics experiment law,
The defeated system hydrate slurry decomposition kinetics mould of pipe is proposed by the influence for considering intrinsic kinetics, heat and mass transfer
Type.The model is on the basis of Kim decomposes intrinsic kinetics model, using the promotion of fugacity difference characterization decomposition of hydrate reaction
Power, while influence the effects that consider heat transfer, mass transfer, particle coalescence to decomposition of hydrate process, and will heat transfer, effect of mass transmitting coupling
Close apparent decomposition rate constant KdIn, the influence of particle/drop coalescence is introduced into decomposable process in phase interfacial parameter, into
And obtain the decomposition of hydrate kinetic model for considering above-mentioned influence factor.
Specific modeling process is as follows, is first depending on the semibatch that Canadian scientist Kim is directed to the stirring of constant temperature and pressure band
The intrinsic power of methane hydrate decomposition reaction for not considering heat and mass transport and being influenced on decomposition of hydrate that reaction kettle system proposes
Learn model.The motive force of model selection fugacity difference characterization decomposition of hydrate reaction, the fugacity in decomposable process in reaction kettle
Consistently lower than balance fugacity.Its expression formula is：
In above formula：nHFor the amount of methane hydrate substance, mol；T be the decomposition of hydrate time, s；KdIt is anti-for eigen decomposition
The rate constant answered, mol/ (MPasm2)；AsFor the total surface area of decomposition of hydrate, m2；F is that methane exists under experiment condition
Gas phase fugacity at hydrate shell outer boundary, MPa；feFor methane gas fugacity, MPa under condition of three-phase balance.
In view of influence and practical pipe defeated system of the factors such as heat transfer, mass transfer to decomposition of hydrate rate in the case of practical pipe is defeated
Interaction violent between particle/drop during the complexity and decomposition of hydrate of flow velocity rule.Therefore, the present invention will pass
The defeated system decomposition reaction rate constant K of practical pipe is coupled in heat, influence of the mass transfer factor to decomposition of hydrate ratedIn, table
It is up to formula：
This experiment loop using heating by the way of carry out decomposition of hydrate, in decomposable process in circuit natural gas fugacity
It is consistently higher than it and balances fugacity, with reference to above formula, obtain the natural gas suitable for High-Voltage Experimentation circuit under multi-phase mixed delivering system and be hydrated
The decomposition model of object：
In above formula：nd(t)For the gas burst size of t moment, mol；nd(end)Consumption at the end of for hydrate formation
Tolerance, mol；T be the decomposition of hydrate time, s；Kd0For the rate constant of eigen decomposition reaction, mol/ (MPasm2)；T is
Gas hydrate dissociation system real time temperature, K；R is mol gas constant, 8.314J/ (molK)；AsFor decomposition of hydrate
Total surface area, m2；fo,iFor gas phase fugacity of the component a certain under experiment condition at hydrate shell outer boundary, MPa；fe,iFor certain
The hydrate phase balance fugacity of one component, MPa；α is carry-over factor；Δ E be activation energy, J/mol.
In the parameter needed for above calculate, the air consumption n at the end of hydrate formationd(end), decomposition of hydrate when
Between t, gas hydrate dissociation system real time temperature T, decomposition of hydrate total surface area AsIt can be measured by experimental data
And it is simply calculated；The rate constant K of eigen decomposition reactiond0, carry-over factor α, activation energy Δ E can be tested by high pressure circuit
Data are fitted to obtain.
Then, phase flash calculation day is carried out based on ripe thermodynamics Phase Equilibrium theory according to the component of natural gas of natural gas
So fugacity under the physical parameter table of gas, hydrate formation curve, three-phase equilibrium fugacity and different temperatures pressure state；
Secondly, the amount of the substance of gas hydrates and decomposition are anti-in system when acquisition gas hydrate dissociation is initial
Real time temperature pressure data after should starting；
Finally, the decomposition of hydrate intrinsic-kinetics equation established before being combined using each parameter determined above is carried out
It calculates, obtains gas hydrate dissociation rate of gas release.
Above-mentioned gas hydrate dissociation rate of gas release computational methods are verified with an example below；Please refer to table
Shown in 1, table 1 is decomposition of hydrate model parameter under the conditions of the differential responses obtained according to high pressure circuit fitting experimental data；
Table 1
According to upper table, below using deionized water, -20# diesel oil and Shan capital line natural gas as test medium, in height
Compacting, which is tested in circuit, carries out the heating decomposition of hydrate experiment of oil gas water three phase constant volume.Pass through data collecting system the real time measure, note
Record temperature, the pressure and other parameters in circuit.By handling experimental data, obtain the gas during decomposition of hydrate and release
High-volume ndWith dimensionless decomposition rate rd。
(1) gas burst size nd, can be obtained by the following formula：
Wherein, nd(t2)For t2The gas burst size at moment, mol；P(t1) it is t1Pressure in moment circuit, MPa；P(t2) be
t2Pressure in moment circuit, MPa；T(t1) it is t1Temperature in moment circuit, K；T(t2) it is t2Temperature in moment circuit, K；R is gas
Body constant, 8.314J/ (molK)；V be circuit volume, m3；Z is compressibility factor.Compressibility factor Z under different pressures, according to
BWRS equations acquire, and concrete numerical value is please referred to shown in table 2, and table 2 is the compressibility factor Z under different pressures.
Table 2
(2) dimensionless decomposition rate rd, can be obtained by the following formula：
Wherein, dimensionless decomposition rate rd, 1/h；nd(t)For the gas burst size of t moment, mol；nd(end)It is given birth to for hydrate
Into the air consumption at the end of process, mol；T is the time, h.
Above-mentioned experiment acquires the dry gas gas sample before experiment, the gas sample of equilibrium state and hydration after hydrate generation respectively
Gas sample after the completion of object decomposition, is analyzed using gas chromatograph, determines its makings composition as shown in table 3 to table 5.Then,
The physical property that phase flash calculation natural gas is carried out based on ripe thermodynamics Phase Equilibrium theory according to the component of natural gas of natural gas is joined
Fugacity under number table, hydrate formation curve, three-phase equilibrium fugacity and different temperatures pressure state.As table 3 to table 5 is shown
Carried out in High-Voltage Experimentation circuit decomposition of hydrate experiment in, decomposition experiment start before, hydrate generation after equilibrium state with
And the makings composition of the natural gas after the completion of decomposition of hydrate.The obtained natural gas of phase flash calculation is illustrated in figure 3 to exist
Hydrate formation curve under three kinds of different phases.
Table 3
Table 4
Table 5
Above-mentioned experiment carries out decomposition of hydrate experiment using the method for constant volume heating, according to experimental data, utilizes the present invention
Middle established decomposition of hydrate intrinsic-kinetics equation, calculates under following different experimental conditions, decomposition of hydrate reaction process
The gas burst size of each time phase, as a result as shown in Fig. 4 A to Fig. 4 D.Referring again to table 6, data shown in table 6 are under different condition
Decomposition of hydrate Model Error Analysis is as a result, according to the absolute deviation of data display simulation result of calculation and experimental data in table 6
Between 7%~27%, but it was found from the data analysis in Fig. 4 A to Fig. 4 D：Model provided by the invention can describe well
The variation tendency of gas burst size during decomposition of hydrate.Experimental data is not filtered, is that simulation is caused to calculate
As a result with the relative deviation of experimental data it is bigger than normal the main reason for.In general, this modeling result of calculation and experimental data
The goodness of fit is in tolerance interval.
Table 6
The results show that utilize gas hydrate dissociation rate of gas release computational methods provided by the present invention and device
Calculating and prediction to rate of gas release during gas hydrate dissociation can be achieved, in this way based on write it is soft
Part can contribute to the operation of Monitoring Gas Pipeline safety and steady, can also be used as and carry out the exploitation of hydrate ore deposit in a deep going way, using hydration
Object is as medium storing with transporting gas, purifying with detaching the technical research such as gas and application hydrate slurry conveying gas
Basis.
Particular embodiments described above has carried out the purpose of the present invention, technical solution and advantageous effect further in detail
Describe in detail it is bright, it should be understood that the above is only a specific embodiment of the present invention, the guarantor being not intended to limit the present invention
Range is protected, all within the spirits and principles of the present invention, any modification, equivalent substitution, improvement and etc. done should be included in this
Within the protection domain of invention.
Claims (7)
1. a kind of gas hydrate dissociation rate of gas release computational methods, which is characterized in that the method includes：
According to the theory basis of decomposition of hydrate, establish related in hydrate slurry decomposition kinetics model and model of fit
Parameter values；
Phase flash distillation is carried out according to gas component, calculates and obtains natural gas relevant parameter；
Real-time temperature after the amount of the substance of gas hydrates and decomposition reaction start when acquisition gas hydrate dissociation is initial
Spend pressure data；
According to relevant parameter numerical value, the natural gas correlative in the hydrate slurry decomposition kinetics model, the model
Real-time temperature after the amount of the substance of gas hydrates and the decomposition reaction start when property parameter, the decomposition of hydrate are initial
Pressure data is spent, obtains gas hydrate dissociation rate of gas release；
The hydrate slurry decomposition kinetics model includes the following formula：
In above formula, nd(t)For the gas burst size of t moment, mol；nd(end)Air consumption at the end of for hydrate formation,
mol；T be the decomposition of hydrate time, s；Kd0For the rate constant of eigen decomposition reaction, mol/ (MPasm2)；T is natural gas
Decomposition of hydrate system real time temperature, K；R is mol gas constant, 8.314J/ (molK)；AsSummary table for decomposition of hydrate
Area, m2；fo,iFor gas phase fugacity of the component a certain under experiment condition at hydrate shell outer boundary, MPa；fe,iFor a certain component
Hydrate phase balance fugacity, MPa；α is carry-over factor；△ E be activation energy, J/mol.
2. gas hydrate dissociation rate of gas release computational methods according to claim 1, which is characterized in that described
According to the theory basis of decomposition of hydrate, establish hydrate slurry decomposition kinetics model and include：According to decomposition of hydrate
The influence to decomposition of hydrate rate of intrinsic kinetics and heat transfer factor, mass transfer factor, establish comprising intrinsic, heat transfer, mass transfer
And the hydrate slurry decomposition kinetics model of coalescence characteristic.
3. gas hydrate dissociation rate of gas release computational methods according to claim 1, which is characterized in that described
According to the theory basis of decomposition of hydrate, it is fitted relevant parameter numerical value in decomposition of hydrate model and includes：According to high compression ring
Road gas hydrate dissociation fitting experimental data relevant parameter numerical value.
4. gas hydrate dissociation rate of gas release computational methods according to claim 1, which is characterized in that described
Phase flash distillation is carried out according to gas component to include：According to the component of natural gas of natural gas, and based on thermodynamics Phase Equilibrium theory into
Row phase flashes.
5. gas hydrate dissociation rate of gas release computational methods according to claim 1, which is characterized in that described
Natural gas relevant parameter includes at least physical parameter table, hydrate formation curve, three-phase equilibrium fugacity and the non-equality of temperature of natural gas
Spend the fugacity under pressure state.
6. a kind of gas hydrate dissociation rate of gas release computing device, which is characterized in that described device includes：
Modeling unit for the theory basis according to decomposition of hydrate, establishes hydrate slurry decomposition kinetics model simultaneously
Relevant parameter numerical value in model of fit；
Flash evaporation unit for carrying out phase flash distillation according to gas component, calculates and obtains natural gas relevant parameter；
Collecting unit, for obtain gas hydrate dissociation it is initial when gas hydrates substance amount and decomposition reaction open
Real time temperature pressure data after beginning；
Data processing unit, for according to relevant parameter numerical value in the hydrate slurry decomposition kinetics model, the model,
The amount of the substance of gas hydrates and the decomposition reaction are opened when the natural gas relevant parameter, the decomposition of hydrate are initial
Real time temperature pressure data after beginning obtains gas hydrate dissociation rate of gas release；
The modeling unit includes：According to the theory basis of decomposition of hydrate, establish following hydrate slurry and decompose power
Learn model：
In above formula, nd(t)For the gas burst size of t moment, mol；nd(end)Air consumption at the end of for hydrate formation,
mol；T be the decomposition of hydrate time, s；Kd0For the rate constant of eigen decomposition reaction, mol/ (MPasm2)；T is natural gas
Decomposition of hydrate system real time temperature, K；R is mol gas constant, 8.314J/ (molK)；AsSummary table for decomposition of hydrate
Area, m2；fo,iFor gas phase fugacity of the component a certain under experiment condition at hydrate shell outer boundary, MPa；fe,iFor a certain component
Hydrate phase balance fugacity, MPa；α is carry-over factor；△ E be activation energy, J/mol.
7. gas hydrate dissociation rate of gas release computing device according to claim 6, which is characterized in that described
Also comprising fitting module, the fitting module is used for the theory basis according to decomposition of hydrate, utilizes high pressure modeling unit
Circuit fitting experimental data relevant parameter numerical value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510744846.0A CN105426666B (en) | 2015-11-05 | 2015-11-05 | Gas hydrate dissociation rate of gas release computational methods and its device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510744846.0A CN105426666B (en) | 2015-11-05 | 2015-11-05 | Gas hydrate dissociation rate of gas release computational methods and its device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105426666A CN105426666A (en) | 2016-03-23 |
| CN105426666B true CN105426666B (en) | 2018-06-08 |
Family
ID=55504874
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510744846.0A Active CN105426666B (en) | 2015-11-05 | 2015-11-05 | Gas hydrate dissociation rate of gas release computational methods and its device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105426666B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107620590B (en) * | 2017-08-08 | 2018-06-22 | 广州海洋地质调查局 | A kind of sea bottom hydrate recovery process balances each other dynamic method for visualizing and device |
| CN109243540B (en) * | 2018-08-13 | 2022-02-08 | 中国石油大学(北京) | Method and device for calculating release rate of hydrate decomposition gas |
| CN109344446B (en) * | 2018-09-06 | 2021-05-25 | 中国石油大学(北京) | Method and device for calculating viscosity of hydrate slurry |
| CN109723970B (en) * | 2019-01-04 | 2019-11-08 | 中山大学 | The Monitoring and management system and method that gas hydrate slurry body transports |
| CN115408889B (en) * | 2022-11-01 | 2023-03-24 | 中国石油大学(华东) | Method for exploiting sea natural gas hydrate by combining heat injection, fracturing and depressurization |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1529238A (en) * | 2003-10-17 | 2004-09-15 | 大庆油田有限责任公司 | Three-dimensinal composite drive computer simulating method for low-concentration surface active agent and phase state joint |
| CN101477093A (en) * | 2008-12-29 | 2009-07-08 | 中国科学院广州能源研究所 | Gas hydrate kinetic analysis apparatus |
| CN103257203A (en) * | 2012-02-16 | 2013-08-21 | 中国科学院物理研究所 | High-pressure integrated experiment apparatus for novel energy hydrate |
| CN103676652A (en) * | 2013-12-06 | 2014-03-26 | 杭州电子科技大学 | Dynamic modeling method of flash evaporator |
| CN104462684A (en) * | 2014-12-02 | 2015-03-25 | 辽宁石油化工大学 | Method for predicting generation of hydrate in nature gas pipeline |
-
2015
- 2015-11-05 CN CN201510744846.0A patent/CN105426666B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1529238A (en) * | 2003-10-17 | 2004-09-15 | 大庆油田有限责任公司 | Three-dimensinal composite drive computer simulating method for low-concentration surface active agent and phase state joint |
| CN101477093A (en) * | 2008-12-29 | 2009-07-08 | 中国科学院广州能源研究所 | Gas hydrate kinetic analysis apparatus |
| CN103257203A (en) * | 2012-02-16 | 2013-08-21 | 中国科学院物理研究所 | High-pressure integrated experiment apparatus for novel energy hydrate |
| CN103676652A (en) * | 2013-12-06 | 2014-03-26 | 杭州电子科技大学 | Dynamic modeling method of flash evaporator |
| CN104462684A (en) * | 2014-12-02 | 2015-03-25 | 辽宁石油化工大学 | Method for predicting generation of hydrate in nature gas pipeline |
Non-Patent Citations (4)
| Title |
|---|
| 天然气水合物分解动力学模型研究;雷昊等;《重庆科技学院学报（自然科学版）》;20150831;第17卷(第4期);第85-89页 * |
| 天然气水合物形成条件预测及防止技术(续);李长俊等;《管道技术与设备》;20020430(第2期);正文第10页第1栏，表6，图6 * |
| 流动体系水合物分解研究进展;史博会等;《油气储运》;20140725;第33卷(第7期);摘要，正文第2大节 * |
| 管道流动体系天然气水合物分解动力学模型研究进展;王淼等;《天然气化工（C1化学与化工）》;20150425;第40卷(第2期);第82-87页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105426666A (en) | 2016-03-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105426666B (en) | Gas hydrate dissociation rate of gas release computational methods and its device | |
| Zhan et al. | Experimental study on characteristics of methane hydrate formation and dissociation in porous medium with different particle sizes using depressurization | |
| WO2021227384A1 (en) | Supergravity hydrate research test system and method | |
| US9970267B2 (en) | Experimental device for simulating exploitation of natural gas hydrate in permeable boundary layer | |
| US10095819B2 (en) | Simulation experiment system and simulation method of entire natural gas hydrate exploitation process | |
| WO2018112902A1 (en) | Experimental device and method for researching relationship between sand production behavior and porous medium radial deformation in natural gas hydrate mining process | |
| CN105301200B (en) | Testing apparatus for characteristics of sand production during mining of natural gas hydrate | |
| Li et al. | Experimental study on gas production from methane hydrate in porous media by huff and puff method in pilot-scale hydrate simulator | |
| Li et al. | Experimental investigation into gas production from methane hydrate in sediment by depressurization in a novel pilot-scale hydrate simulator | |
| WO2017088226A1 (en) | Experimental apparatus and method for simulating stratum deformation in natural gas hydrate exploitation process | |
| Zhou et al. | Investigation on the kinetics of carbon dioxide hydrate formation using flow loop testing | |
| WO2017008354A1 (en) | Experimental device and experimental method for studying porous medium skeleton change in natural gas hydrate decomposition process | |
| CN102865066B (en) | Experiment device and method for deepwater shaft multiphase flow containing natural gas hydrate phase changes | |
| CN106544070B (en) | Method and device for generating and replacing natural gas hydrate | |
| WO2018112899A1 (en) | Experimental device and method for conducting multiphase separation on natural gas hydrate well drilling liquid | |
| CN103233704A (en) | Method and device for simulating experiment of permafrost region natural gas hydrate mining through replacement of CO2/N2 | |
| WO2014176794A1 (en) | Three-dimensional strata stability simulation device for natural gas hydrate extraction | |
| Li et al. | Kinetic behaviors of methane hydrate formation in porous media in different hydrate deposits | |
| CN107503715A (en) | A kind of device for simulating parallel water horizontal well seawater injection exploitation of gas hydrate | |
| CN102539296A (en) | Method and special device for hydrate deposit seepage flow test | |
| Almenningen et al. | Experimental investigation of critical parameters controlling CH4–CO2 exchange in sedimentary CH4 hydrates | |
| CN108051354A (en) | A kind of hypotonic hydrate sediment permeability survey method and apparatus based on impulse attenuation analysis | |
| Zhang et al. | Replacement of CH4 in hydrate in porous sediments with liquid CO2 injection | |
| CN109243540A (en) | The calculation method and device of decomposition of hydrate rate of gas release | |
| Sayani et al. | Comprehensive review on various gas hydrate modelling techniques: Prospects and challenges |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |