CN108520101A - Geothermal well well casing scaling prediction method - Google Patents

Geothermal well well casing scaling prediction method Download PDF

Info

Publication number
CN108520101A
CN108520101A CN201810204981.XA CN201810204981A CN108520101A CN 108520101 A CN108520101 A CN 108520101A CN 201810204981 A CN201810204981 A CN 201810204981A CN 108520101 A CN108520101 A CN 108520101A
Authority
CN
China
Prior art keywords
state parameter
group
well casing
deviation
geothermal
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
CN201810204981.XA
Other languages
Chinese (zh)
Other versions
CN108520101B (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.)
Guangzhou Institute of Energy Conversion of CAS
Original Assignee
Guangzhou Institute of Energy Conversion of CAS
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 Guangzhou Institute of Energy Conversion of CAS filed Critical Guangzhou Institute of Energy Conversion of CAS
Priority to CN201810204981.XA priority Critical patent/CN108520101B/en
Publication of CN108520101A publication Critical patent/CN108520101A/en
Application granted granted Critical
Publication of CN108520101B publication Critical patent/CN108520101B/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
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/40Searching chemical structures or physicochemical data

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Evolutionary Computation (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Computing Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a kind of methods of geothermal well well casing scaling prediction, including:Step 101:Obtain the original state parameter of stratum depths and well casing;Step 102:Original state parameter input multiphase flow coupling geochemical model is simulated, different moments side-wall begriming thickness is obtained.The geothermal well well casing scaling prediction method of the present invention, the state parameter that geochemical model simulates geothermal well different moments is coupled by multiphase flow, namely simulate the system status of different moments, can obtain when fouling and the fouling thickness of different moments, to for rationally and effectively reduction fouling and the implementation of scale removal scheme guidance is provided.

Description

Geothermal well well casing scaling prediction method
Technical field
The present invention relates to geothermal well technical fields, and in particular to a kind of geothermal well well casing scaling prediction method.
Background technology
Geothermal energy resources are the clean energy resourcies of earth interior generally existing, with Development of China's Urbanization and job development.Underground heat provides Source becomes very popular green energy resource.It is contemplated that coming years China increases the explosion type for welcoming geothermal exploitation, application It is long.
However, there is also many restraining factors for the utilization of geothermal energy, wherein geothermal well and pipe fitting scale problems is each The primary problem often faced in class geothermal applications.No matter shallow-layer or deep geothermal resources, either generate electricity, plant, support It grows, almost all of geothermal applications mode can all encounter pipeline and equipment scaling problem.Once fouling, system operation cost meeting It improves, resource utilization reduces.When serious, system breakdown can not be even continuing with.
Currently, for geothermal well well casing scale problems research still in the primary stage, for may fouling time, The researchs such as position, rate are still weak, lack and suitable estimate means and method.Therefore, fouling processing is still with removing after fouling Based on dirt, complicated for operation and waste of resource.
Invention content
In view of the above-mentioned drawbacks in the prior art or insufficient, the present invention provides a kind of geothermal well well casing scaling prediction method, Can Accurate Prediction well casing fouling condition, to take measures to prevent fouling in advance.
To achieve the above object, the technical solution adopted in the present invention is:
A kind of method of geothermal well well casing scaling prediction, including:
Step 101:Obtain the original state parameter of stratum depths and well casing;
Step 102:Original state parameter input multiphase flow coupling geochemical model is simulated, is obtained not Side-wall begriming thickness in the same time.
The step 102 includes as an improvement of the present invention:
Step 201:Original state parameter input geochemistry data library is simulated, obtains and exports first Group state parameter;
Step 202:First group of state parameter input multi- scenarios method model is simulated, obtains and exports second Group state parameter;
Step 203:First group of state parameter is compared with second group of state parameter, obtains deviation;
Step 204:Judge whether the deviation is less than or equal to required precision;
If so, being defined as second group of state parameter to meet the state parameter of required precision, it is transferred to step 205;
If it is not, second group of state parameter is then inputted the geochemistry data library again, repeat step 201 to 204 deviation corrects subcycle;
Step 205:The state parameter for meeting required precision is inputted the geochemistry data library to simulate, Obtain the side-wall begriming situation at this moment;
Step 206:The state parameter for meeting required precision is inputted the multi- scenarios method model to simulate, is obtained To subsequent time state parameter, step 201 is returned to, is recycled into next round.
Further, the step 204 further includes:
If after having carried out the deviation correction subcycle of n wheels, the deviation is still greater than required precision, then according to institute It states deviation and corresponds to the adjustment original state parameter, wherein n >=1.
Further, the multi- scenarios method model includes:
Diabatic process governing equation
Wherein, u is interior energy (kJ kg-1), ρ is density (kg m-3), v is fluid flow rate (scalar-unit m s-1), and λ is heat Conductance (J (smK)-1), T is temperature (K), and t is the time (s);
The mass transfer equation of diffusion and transmittance process control
Wherein, c indicates molar density (mol m-3), D is diffusion coefficient (m2 s-1), R indicates component geochemical reaction Rate (mol m-3 s-1), i indicates component.
Further, the original state parameter includes temperature, pressure, phase, component and concentration of component.
Compared with prior art, the present invention advantage is:
Geochemical model is coupled by multiphase flow and simulates the state parameter of geothermal well different moments, that is, is simulated not System status in the same time, can obtain when fouling and the fouling thickness of different moments, to rationally and effectively to reduce Fouling and the implementation of scale removal scheme provide guidance.
Description of the drawings
Fig. 1 is the flow chart of one embodiment of geothermal well well casing scaling prediction method of the present invention;
Fig. 2 is the flow chart of another embodiment of geothermal well well casing scaling prediction method of the present invention.
Specific implementation mode
The application is described in further detail with reference to the accompanying drawings and examples.It is understood that this place is retouched The specific embodiment stated is used only for explaining related invention, rather than the restriction to the invention.It also should be noted that being Convenient for description, is illustrated only in attached drawing and invent relevant part.
Fig. 1 shows a kind of geothermal well well casing scaling prediction method according to the embodiment of the present application comprising:
Step 101:Obtain the original state parameter of stratum depths and well casing.
In the present embodiment, state parameter can be temperature, pressure, phase, component, concentration of component, flow velocity, solubility, ratio Hold the physical quantity with description state behaviors such as chemical reaction equilibrium parameters.Temperature can be the temperature of stratum different depth, also may be used The temperature for thinking liquid in well casing can also be the temperature etc. of gas in well casing.Similarly, pressure can be liquid in well casing Pressure can also be the partial pressure of gas with various in well casing.Phase can be any one in gas, liquid, solid three-phase or multiphase.Component Can be in well casing different phase at being grouped as situation, such as gaseous component may include air, carbon dioxide, liquid phase group It may include dissolved substance etc. in water and water to divide, and solid phase components may include the mineral of Precipitation, such as calcium carbonate. Concentration of component can be the ion that liquid is included in well casing, the concentration of compound, as calcium ion concentration, carbanion are dense Degree.For normal conditions in order to ensure that can obtain subsequently judging state parameter needed for system fouling, original state parameter can be down to Include temperature, pressure, phase, component and concentration of component less.
Original state parameter can be the state parameter of well casing different depth position, i.e. original state parameter can be more The state parameter of group different depth position, depth location can be specifically configured according to forecast demand, e.g., original state ginseng Number can be well casing at interval of the state parameter at 100 meters, or can also be every at interval of 10 meters of state parameter The state parameter of a point.Original state parameter can be obtained by field survey, such as the temperature of well head, pressure, flow velocity; It can be obtained by local existing data information, such as the temperature of the well casing of different depth, pressure, phase, dissolved constituent composition And concentration etc.;It can also be obtained by the data operation that existing data and/or field survey obtain, such as the temperature of depth of stratum Degree, pressure etc.;It can also estimate to obtain according to experimental simulation or other modes.That is, can be obtained according to state parameter Complexity and specific forecast demand obtain original state parameter in different ways.
Step 102:Original state parameter input multiphase flow coupling geochemical model is simulated, when obtaining different Carve side-wall begriming thickness.
Multiphase flow coupling geochemical model can be different gas-liquid-solid phase change multi-phase migration reaction geochemistry mould Type such as can be water-gas-salt-ore body system multiphase flow coupling geochemical model, can also be water-gas-salt system multiphase flow Geochemical model is coupled, can also be water-salt-ore body system multiphase flow coupling geochemical model, can also be that other can be with The multiphase flow for simulating geothermal well system couples geochemical model.Original state parameter input multiphase flow is coupled into geochemistry Model can simulate heterogeneous system in a certain moment well casing by the model, and obtain the moment corresponding state ginseng Number, so as to according to the chemical reaction equilibrium condition of the moment system, judge whether there is Precipitation and precipitation is precipitated Quality or volume;It is based on the moment state parameter simultaneously, coupling geochemical model by multiphase flow is simulated, and can be obtained To the state parameter of subsequent time system, so as to continue to judge Precipitation and the quality or volume of precipitation is precipitated.It is logical Multiphase flow coupling geochemical model is crossed persistently to simulate each moment at well casing different depth position in geothermal well system, it can To obtain whether each moment at well casing different depth position has the quality or volume of Precipitation and precipitation precipitation, to Can obtain the borehole wall at well casing different location when fouling and fouling thickness.
In above-described embodiment, the state parameter that geochemical model simulates geothermal well different moments is coupled by multiphase flow, Namely simulate different moments system status, can obtain when fouling and the fouling thickness of different moments, to for close Reason is effectively reduced fouling and the implementation of scale removal scheme provides guidance.
Fig. 2 shows the geothermal well well casing scaling prediction methods according to another embodiment of the application comprising:
Step 200:Obtain the original state parameter of stratum depths and well casing.
Step 201:Original state parameter input geochemistry data library is simulated, obtains and exports first group of shape State parameter.
In the present embodiment, geochemistry data library belongs to a part for multiphase flow coupling geochemical model, is to integrate The applied database set up on the basis of scientific research data and data for many years can be comprehensive geochemistry data Library, can also be with the relevant thematic data base of geothermal well, different geochemistry numbers can be specifically selected according to specific requirements According to library.After original state parameter inputs geochemistry data library, it can be obtained with the relevant data of original state parameter by inquiring To first group of state parameter, first group of state parameter can also be obtained by calculation, can also be obtained by inquiry and calculations incorporated To first group of state parameter.That is, input state parameter is to geochemistry data library, which can be with According to the rule and algorithm of setting, obtains and export corresponding state parameter.State parameter that is acquired and exporting, i.e., first Group state parameter, may include the original state parameter after the correction simulated by geochemistry data library, can also Include the other new state parameters simulated by geochemistry data library, which meets geochemical Act on.For example, inputting at certain moment well casing different depth position temperature, pressure and component volume of data to geochemistry Database is simulated, and temperature, pressure and a series of numbers of component after being corrected at the moment different depth position can be obtained According to, or the concentration of component at the moment well casing different depth position and temperature, pressure and component after correction can be obtained Volume of data.
In the present embodiment, by the way that state parameter is inputted the simulation of geochemistry data library, it not only may be implemented to join state Several supplements and extension, can also realize the correction to state parameter, and it is bright in simulation geochemistry process to solve single model There is a problem of that deviation is larger when aobvious complex system and can not be applicable in.
Step 202:First group of state parameter is inputted multi- scenarios method model to simulate, obtains and exports second group of shape State parameter.
In the present embodiment, multi- scenarios method model also belongs to a part for multiphase flow coupling geochemical model, Ke Yiwei Gas-liquid-solid phase change multi-phase migration reaction model can also be fluid stream such as the coupling of fluid flowing-heat transfer-mass transfer three Dynamic-coupling of mass transfer two, can also be that fluid flowing-heat transfer-mass transfer-power four couples, can also be other geothermal well systems Relevant multiple physical field coupling model.It is appreciated that multi- scenarios method model can be from physical fields such as heat, quality and momentum The mathematical model for the description fluid state real-time change that angle is established.Second group of state parameter can be according to multi- scenarios method First group of state parameter after the correction that modeling obtains, second group of state parameter meet heat conservation, the conservation of mass and The laws such as the conservation of momentum.For example, the first of the well casing different depth position that input is obtained by upper step geochemistry data library The at most field coupling models of group state parameter T1, P1, c1, i1 and ρ 2, can obtain after being corrected at well casing different depth position Second group of state parameter T2, P2, c2, i2 and ρ 2, wherein T are temperature, and P is pressure, and c is molar density, and i is component, and ρ is close Degree.
Preferably, multi- scenarios method model includes:The mass transfer side of diabatic process governing equation and diffusion and transmittance process control Journey.
In the present embodiment, according in water-gas-salt-ore body system in geothermal well fluid flowing, heat transfer and mass transfer process, The diabatic process governing equation of heat, momentum and quality change situation and diffusion during describing fluid heat transferring can be established With the mass transfer equation of transmittance process control.Due to the mass transfer equation of diabatic process governing equation and diffusion and transmittance process control Involved in fluid velocity correlated condition parameter, thus may not need be added flow process governing equation, so that it may to obtain fluid Second group of state parameter.
For example, the diabatic process governing equation of description well casing thermal change situation can be as shown in formula (1):
Wherein, u is interior energy (kJ kg-1), ρ is density (kg m-3), v is fluid flow rate (scalar-unit m s-1), and λ is heat Conductance (J (smK)-1), T is temperature (K), and t is the time (s).
The diffusion and transmittance process governing equation for describing well tube fluid (dirt) can be as shown in formulas (2):
Wherein, c indicates molar density (mol m-3), D is diffusion coefficient (m2s-1), R indicates component geochemical reaction speed Rate (mol m-3 s-1), i indicates component.
A series of T1, P1, c1, i1 and ρ 1 of the well casing different depth positions obtained by geochemistry data library is defeated Enter above formula (1)-(2), and by simultaneous above formula (1)-(2), can obtain a series of t2, p2 after the correction of multi- scenarios method model, C2, i2 and ρ 2.
In more embodiments, diabatic process governing equation and the mass transfer equation that is controlled with transmittance process of diffusion can be adopted With the equation of other forms, as long as can be used for describing the process of thermal change, momentum change and mass change in geothermal well i.e. The technique effect of the state parameter after being corrected can be achieved.
Step 203:First group of state parameter is compared with second group of state parameter, obtains deviation.
In the present embodiment, first group of state of the well casing different depth position obtained by geochemistry data library is joined It counts and is compared by second group of state parameter that multi- scenarios method is simulated, the deviation of the two can be obtained.The deviation It can be indicated with the absolute value of the difference of the two, can also be obtained with the absolute value of the difference of the two and geochemistry data library the The ratio that one group of state parameter is divided by indicates, can also can be used to indicate that deviation feelings between the two with other forms The numerical value of condition indicates.The state parameter being compared can be a certain item state parameter, and such as concentration can also be a few shapes State parameter, such as temperature, pressure and concentration.When being that a few state parameters are compared, one of which state can first be selected to join Number is compared, and after this state parameter meets the requirements, then other is selected to be compared.
Step 204:Whether judgment bias value is less than or equal to required precision;
If so, being defined as simulated by multi- scenarios method second group of state parameter to meet the state of required precision Parameter is transferred to step 205;
If it is not, second group of state parameter is then inputted geochemistry data library again, step 201 is repeated to step 204, Subcycle is corrected into next round deviation.
In the present embodiment, required precision can be preset according to requirements such as prediction accuracies, can also be according to other need It asks and is set.Deviation can be compared with the value set by required precision, whether judgment bias value is less than or equal to Value set by required precision.If the deviation is less than or equal to the value set by the required precision, illustrate multi- scenarios method mould Quasi- obtained state parameter meets geochemistry data library and requiring simultaneously, i.e., the state parameter not only met each physical field conservation but also Meet geochemistry process, then the state parameter can be made to the state parameter of this moment system.If the deviation is more than precision It is required that then illustrating that the state parameter that geochemistry database simulation obtains is differed with the state parameter that multi- scenarios method is simulated It is larger, which can again be inputted to geochemistry data library and simulated to obtain third group state parameter, and by Three groups of state parameters input multi- scenarios method model and are simulated to obtain the 4th group of state parameter again, by third group state parameter It is compared with the 4th group of state parameter, judges whether its deviation is less than or equal to required precision, if so, the 4th group of shape State parameter is to meet the state parameter of required precision, is followed if otherwise continuing return to step 201 and carrying out next round deviation syndrome Ring, until obtaining the status parameter values for meeting required precision." first ", " second " and similar word in the embodiment of the present application Language is not offered as any sequence, quantity or importance, and is used only to distinguish different groups of state parameter.
Preferably, if after having carried out n wheel deviation correction subcycles, deviation is still greater than required precision, then according to deviation Corresponding adjustment original state parameter, wherein n >=1.
In the present embodiment, since original state parameter can be according to existing data or experimental simulation or other modes Estimation obtains, thus there may be certain deviations.Suitable cycle-index n can be set according to specific requirements.Preferably, n More than 1000.If the two groups of state parameters obtained after being taken turns by geochemistry data library and multi- scenarios method mold cycle operation n Still greater than required precision, i.e., the state parameter obtained after deviation aligning step cycle n wheels can not meet geochemical deviation simultaneously When learning database and multi- scenarios method model, then correction original state parameter can be corresponded to according to the deviation, to ensure in system The accuracy of state parameter.Original state parameter to be adjusted can be obtained by field survey, can also pass through experimental simulation It obtains, can also be obtained by other suitable modes.
Subcycle is corrected by deviation to be corrected original state parameter, it is possible to prevente effectively from because of original state parameter The inaccurate situation of deviation is larger and occurs simulation, improves the accuracy of modeling, at the same can also effectively avoid be System enters endless loop state because being unable to get satisfactory state parameter.
Step 205:The state parameter input geochemistry data library of required precision will be met, obtains side-wall begriming situation.
The state parameter of required precision will be met, i.e., meet the shape of geochemistry database and multi- scenarios method model simultaneously State parameter, input geochemistry data are simulated in library, can obtain well casing different depth position and be in the state parameter item Under part when system chemical reaction equilibrium scale-forming ion concentration, by the concentration of scale-forming ion and chemical reaction in the state parameter Concentration when balance compares, judge the borehole wall under the conditions of the state parameter at system different location whether fouling (it is heavy to be precipitated Form sediment), and the quality or volume of fouling (precipitation is precipitated), that is, side-wall begriming situation are calculated, obtain side-wall begriming thickness.
Step 206:The state parameter for meeting required precision input multi- scenarios method model is simulated, according to flow velocity or Other state parameters with time correlation can obtain the system subsequent time state parameter, i.e. next round recurrent state is joined Number arrives step 201, is recycled into next round.
Since the state parameter in geothermal well system is in the process of real-time change, next round can be recycled shape State parameter inputs geochemistry data library and obtains third group state parameter, and the shape that the geochemistry data library is obtained again State parameter input multi- scenarios method model obtains the 4th group of state parameter, and two groups of state parameters are compared, judge its deviation With the size of required precision, the state parameter for meeting required precision is obtained again.Such persistent loop geochemistry data library and The simulation process of multi- scenarios method model, the case where each moment side-wall begriming at well casing different depth position can be obtained.
In the present embodiment, by using the mode in geochemistry data library and multi- scenarios method models coupling, to state parameter Simulation in real time and correction, may be implemented the real-time simulation to geothermal well multiphase flow complex system, and it is multiple to solve single modeling Deviation is larger when miscellaneous system and the scope of application has the problem of limitation;And by the modeling, obtain different moments knot Dirty situation provides reasonable guidance for design, operating and the scale removal maintenance etc. of geothermal device.
Above description is only the preferred embodiment of the application and the explanation to institute's application technology principle.Art technology Personnel should be appreciated that invention scope involved in the application, however it is not limited to skill made of the specific combination of above-mentioned technical characteristic Art scheme, while should also cover in the case where not departing from the inventive concept, by above-mentioned technical characteristic or its equivalent feature into Other technical solutions of the arbitrary combination of row and formation.

Claims (5)

1. a kind of geothermal well well casing scaling prediction method, which is characterized in that the method includes:
Step 101:Obtain the original state parameter of stratum depths and well casing;
Step 102:Original state parameter input multiphase flow coupling geochemical model is simulated, different moments are obtained Side-wall begriming thickness.
2. geothermal well well casing scaling prediction method according to claim 1, which is characterized in that the step 102 includes:
Step 201:Original state parameter input geochemistry data library is simulated, obtains and exports first group of state Parameter;
Step 202:First group of state parameter input multi- scenarios method model is simulated, obtains and exports second group of state Parameter;
Step 203:First group of state parameter is compared with second group of state parameter, obtains deviation;
Step 204:Judge whether the deviation is less than or equal to required precision;
If so, being defined as second group of state parameter to meet the state parameter of required precision, it is transferred to step 205;
If it is not, second group of state parameter is then inputted the geochemistry data library again, step 201 is repeated to 204 Deviation corrects subcycle;
Step 205:The state parameter for meeting required precision is inputted the geochemistry data library to simulate, obtains this The side-wall begriming situation at moment;
Step 206:The state parameter for meeting required precision is inputted the multi- scenarios method model to simulate, is obtained next Moment state parameter, returns to step 201, is recycled into next round.
3. geothermal well well casing scaling prediction method according to claim 2, which is characterized in that the step 204 further includes:
If after having carried out the deviation correction subcycle of n wheels, the deviation is still greater than required precision, then according to the deviation Value is corresponding to adjust the original state parameter, wherein n >=1.
4. geothermal well well casing scaling prediction method according to claim 2, which is characterized in that the multi- scenarios method model packet It includes:
Diabatic process governing equation
Wherein, u is interior energy (kJ kg-1), ρ is density (kg m-3), v is fluid flow rate (scalar-unit m s-1), λ is thermal conductivity (J(smK)-1), T is temperature (K), and t is the time (s);
The mass transfer equation of diffusion and transmittance process control
Wherein, c indicates molar density (mol m-3), D is diffusion coefficient (m2s-1), R indicates component geochemical reaction rate (mol m-3s-1), i indicates component.
5. according to any geothermal well well casing scaling prediction methods of claim 1-4, which is characterized in that the original state Parameter includes temperature, pressure, phase, component and concentration of component.
CN201810204981.XA 2018-03-13 2018-03-13 Method for predicting scaling of geothermal well pipe Active CN108520101B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810204981.XA CN108520101B (en) 2018-03-13 2018-03-13 Method for predicting scaling of geothermal well pipe

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810204981.XA CN108520101B (en) 2018-03-13 2018-03-13 Method for predicting scaling of geothermal well pipe

Publications (2)

Publication Number Publication Date
CN108520101A true CN108520101A (en) 2018-09-11
CN108520101B CN108520101B (en) 2021-11-16

Family

ID=63433639

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810204981.XA Active CN108520101B (en) 2018-03-13 2018-03-13 Method for predicting scaling of geothermal well pipe

Country Status (1)

Country Link
CN (1) CN108520101B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109740203A (en) * 2018-12-18 2019-05-10 新疆贝肯能源工程股份有限公司 Orientation trajectory design method for underground heat well development
CN110017129A (en) * 2019-05-21 2019-07-16 水利部交通运输部国家能源局南京水利科学研究院 A kind of karst GEOTHERMAL WATER Scaling Tendency Prediction method considering sour gas degassing
CN112100850A (en) * 2020-09-16 2020-12-18 西南交通大学 High-concentration organic wastewater gathering pipeline scaling prediction method based on system dynamics
CN115060870A (en) * 2022-08-11 2022-09-16 中国长江三峡集团有限公司 Geothermic fluid scaling prediction method and device and laboratory reaction equipment

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090119076A1 (en) * 2004-06-03 2009-05-07 Madatov Arkady G Method for Generating a 3D Earth Model
US20120022692A1 (en) * 2010-07-23 2012-01-26 Chevron Phillips Chemical Company Lp Prediction and control solution for polymerization reactor operation
CN102402184A (en) * 2011-10-28 2012-04-04 中国石油集团川庆钻探工程有限公司 Shaft pressure model prediction system controlling method
US20140097338A1 (en) * 2012-10-10 2014-04-10 California Institute Of Technology Mass spectrometer, system comprising the same, and methods for determining isotopic anatomy of compounds
CN104361228A (en) * 2014-11-06 2015-02-18 中国科学院广州能源研究所 New method for calculating terrestrial heat resource quantity of dry heat rock
CN204989151U (en) * 2015-10-10 2016-01-20 马致远 Experimental device for simulation geothermol power tail water recharges chemistry and blocks up
CN106481339A (en) * 2016-09-14 2017-03-08 中国石油大学(华东) Gas hydrates based on electromagnetism thermoacoustic effect are surveyed with probing and analogy method
US20170074093A1 (en) * 2015-09-16 2017-03-16 King Fahd University Of Petroleum And Minerals Well-bore and reservoir monitoring process by logging temperature and resistivity

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090119076A1 (en) * 2004-06-03 2009-05-07 Madatov Arkady G Method for Generating a 3D Earth Model
US20120022692A1 (en) * 2010-07-23 2012-01-26 Chevron Phillips Chemical Company Lp Prediction and control solution for polymerization reactor operation
CN102402184A (en) * 2011-10-28 2012-04-04 中国石油集团川庆钻探工程有限公司 Shaft pressure model prediction system controlling method
US20140097338A1 (en) * 2012-10-10 2014-04-10 California Institute Of Technology Mass spectrometer, system comprising the same, and methods for determining isotopic anatomy of compounds
CN104361228A (en) * 2014-11-06 2015-02-18 中国科学院广州能源研究所 New method for calculating terrestrial heat resource quantity of dry heat rock
US20170074093A1 (en) * 2015-09-16 2017-03-16 King Fahd University Of Petroleum And Minerals Well-bore and reservoir monitoring process by logging temperature and resistivity
CN204989151U (en) * 2015-10-10 2016-01-20 马致远 Experimental device for simulation geothermol power tail water recharges chemistry and blocks up
CN106481339A (en) * 2016-09-14 2017-03-08 中国石油大学(华东) Gas hydrates based on electromagnetism thermoacoustic effect are surveyed with probing and analogy method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
E.STAMATAKIS等: ""Scale prediction in liquid flow through porous media: A geochemical model for the simulation of CaCO3 deposition at the near-well region"", 《JOURNAL OF GEOCHEMICAL EXPLORATION》 *
张恒等: ""水文地球化学模拟技术在康定某高温地热井结垢研究中的应用"", 《新能源进展》 *
王晓星等: ""增强型地热系统开发过程中的多场耦合问题"", 《水文地质工程地质》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109740203A (en) * 2018-12-18 2019-05-10 新疆贝肯能源工程股份有限公司 Orientation trajectory design method for underground heat well development
CN110017129A (en) * 2019-05-21 2019-07-16 水利部交通运输部国家能源局南京水利科学研究院 A kind of karst GEOTHERMAL WATER Scaling Tendency Prediction method considering sour gas degassing
CN110017129B (en) * 2019-05-21 2023-07-21 江苏省环境科学研究院 Karst geothermal water scaling trend prediction method considering acid gas degassing
CN112100850A (en) * 2020-09-16 2020-12-18 西南交通大学 High-concentration organic wastewater gathering pipeline scaling prediction method based on system dynamics
CN115060870A (en) * 2022-08-11 2022-09-16 中国长江三峡集团有限公司 Geothermic fluid scaling prediction method and device and laboratory reaction equipment
CN115060870B (en) * 2022-08-11 2022-11-29 中国长江三峡集团有限公司 Geofluorine fluid scaling prediction method and device and laboratory reaction equipment

Also Published As

Publication number Publication date
CN108520101B (en) 2021-11-16

Similar Documents

Publication Publication Date Title
CN108520101A (en) Geothermal well well casing scaling prediction method
US20130304679A1 (en) System and method for using an artificial neural network to simulate pipe hydraulics in a reservoir simulator
Tzatchkov et al. Advection-dispersion-reaction modeling in water distribution networks
CN108241777B (en) Method for calculating seepage velocity field in hydrate sediment based on unstructured grid finite element method
McBride et al. Accurate modeling of river pollutant transport
AU2011357784A1 (en) System and method for using an artificial neural network to simulate pipe hydraulics in a reservoir simulator
WO2012082273A1 (en) Method and system for coupling reservoir and surface facility simulations
CN111859748B (en) Ocean internal wave simulation method based on vertical mixed coordinates
CN104598701B (en) A kind of river pollution emulation deduction method of three-dimensional scenic segment iteration
Zhou et al. Route optimization of pipeline in gas-liquid two-phase flow based on genetic algorithm
Waldrop et al. Three‐dimensional computation of buoyant plumes
CN109827079B (en) Oil spill source detection method based on submarine petroleum pipeline information physical system
CN103617563A (en) Water supply pipeline non-monitoring node pressure determining method based on ground statistics space analysis theory
CN109210268B (en) Big data processing method based on ultralow-power electromagnetic valve
Demin et al. A numerical method of calculation of currents and sea surface topography in multiply connected domains of the ocean
CN106292289A (en) Process industry controls the mixing elite Stochastic search optimization method of loop object
Trenty et al. A sequential splitting strategy for CO2 storage modelling
CN115293037A (en) Hydrodynamic modeling method of river and lake composite system based on machine learning
Liu et al. A comparative study on hydrodynamics and hydrochemistry coupled simulations of drainage pipe crystallization blockage in karst tunnels
Ravindran et al. Multivariate optimization of production systems: The time dimension
LU500800B1 (en) A method for predicting scaling of geothermal well pipes
Przulj Computational modelling of vortex shedding flows
CN1936907A (en) Method for dynamically setting up mould of parallely operated machine group steam mather tube
CN110717295B (en) Method for tracking streamline distribution of tight sandstone reservoir by using finite element method
Tran et al. New features to Barbaro’s heat exchanger network algorithm: heat exchanger technologies and waste heat flow representation

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