CN109751044A - A kind of oil area heat reservori top surface temperature isogram preparation method - Google Patents
A kind of oil area heat reservori top surface temperature isogram preparation method Download PDFInfo
- Publication number
- CN109751044A CN109751044A CN201910084363.0A CN201910084363A CN109751044A CN 109751044 A CN109751044 A CN 109751044A CN 201910084363 A CN201910084363 A CN 201910084363A CN 109751044 A CN109751044 A CN 109751044A
- Authority
- CN
- China
- Prior art keywords
- isogram
- well
- temperature
- geothermal gradient
- depth
- 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.)
- Pending
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of oily area heat reservori top surface temperature isogram preparing methods comprising the steps of: A, collects well point continuous temperature measurement data, segmentation calculates geothermal gradient；B, using known well point as control condition, the corresponding geothermal gradient isogram of depth segment locating for heat reservori is obtained by Kriging regression method between well；C, structural map and geothermal gradient isogram are divided into the underground heat computing unit to fit like a glove；D, using structural map and geothermal gradient map unit as virtual log and actual well drilled data as control point, heat reservori top surface ground temperature isogram is obtained by interpolation between well.The present invention is using first segmentation establishment geothermal gradient isogram, then the method for calculating control point temperature and working out ground temperature isogram, structural relief bring deviation when avoiding direct interpolation.Cell node is evenly distributed in entire basin as virtual log, and avoiding oil gas drilling well, to be concentrated mainly on structural high part inconvenient to temperature isogram establishment bring.
Description
Technical field
The present invention relates to oily area's geothermal exploitation field, heat reservori top surface temperature isogram establishment side, specifically a kind of oil area
Method.
Background technique
Eastern China major part oil field enters mid-later development phase, and some discarded oil-water wells can be used for exploitation of geothermal.Ground
Heat is used as a kind of clean energy resource, can carry out step application according to underground water temperature difference, be respectively used to power generation, heating, cultivation
Equal fields, the water body that temperature reduces after Ground Application refill underground, may be implemented to recycle.
There is unique advantage using discarded oil-gas field development underground heat.Discarded oil/gas well can be used as hot well, drop
The low cost of geothermal exploitation, avoids the welfare that environment can be brought while the wasting of resources, and the successful Hebei of geothermal exploitation is male
The ground such as county have the condition of smokeless city substantially；Oil/gas well generally all has a well log interpretation data, corresponding aqueous layer depth,
The parameters such as temperature, reservoir thickness, porosity are substantially clear, and opening the later liquid outlet quantity in layer position can also estimate, and avoid the party
The investment of dough figurine power, economy etc., and it is convenient for the formulation of arrangement and method for construction；Having the parameters such as reservoir thickness, porosity, area
Afterwards, the technology that oilfield technology personnel calculate oil and gas reserves is very mature, and the calculating of Geothermal deposit and the technology ten tap
Closely, it is not necessary to carry out relevant technology and grope and train.
However, the result that oil and gas reserves calculates obtains the volume of oil gas, weight, and Geothermal deposit is in addition to needing in block
It provides except geothermal water volume, weight, it is also necessary to consider the temperature of underground water, which directly determines that underground water can
There is provided thermal energy number.It can thus be seen that the key foundation for calculating Geothermal deposit is the temperature point of clear geothermal water reservoir
Cloth.In oil drilling data, generally there is the temperature-measurement material of part well point, can read well point target zone temperature, however thermometric
Well point number it is limited, and be distributed in the high-order bit of construction, the temperature between well point is affected by structural belt position, therefore more
It is lower that the reservoir ground temperature distribution map confidence level obtained between temperature, well using Kriging regression is read simply by well point.For with
Upper problem, the present invention provide a kind of reservoir top surface temperature isogram preparation method, can calculate for the storage of oil field old underground heat
Indispensable basis is provided.
Summary of the invention
It is above-mentioned to solve the purpose of the present invention is to provide a kind of oily area heat reservori top surface temperature isogram preparation method
The problem of being proposed in background technique.
To achieve the above object, the invention provides the following technical scheme:
A kind of oil area heat reservori top surface temperature isogram preparation method comprising the steps of:
A, well point continuous temperature measurement data is collected, the identical position of temperature logs slope is divided into well section, segmentation calculates ground temperature
Gradient；
B, using known well point as control condition, show that depth segment locating for heat reservori is corresponding by Kriging regression method between well
Geothermal gradient isogram；
C, structure of oil field figure is collected, and is aligned in Croel DRAW software with geothermal gradient isogram, is divided complete
Identical underground heat computing unit；
D, it will be constructed using Surfer and be digitized with geothermal gradient distribution map, the corresponding plane coordinates of reading unit node,
The numerical value such as buried depth, geothermal gradient, the corresponding ground temperature of calculate node and as virtual log, while reading actual well drilled data target zone
Temperature obtains heat reservori top surface ground temperature isogram by interpolation between well using virtual log and actual well drilled data as control point.
As further technical solution of the present invention: in step A: by well point continuous temperature measurement data segment processing, same well
The section continuous temperature measurement slope of curve is the same, and corresponding geothermal gradient numerical value is constant.
As further technical solution of the present invention: the well of the geothermal gradient isogram and well point continuous temperature measurement data
Section is corresponding, and in the corresponding geothermal gradient isogram compilation process of a certain depth segment, any one control well point is in the depth segment
Geothermal gradient not with change in depth, control the geothermal gradient between well point influenced by structural relief it is less than normal, using Kriging regression
Method obtains.
As further technical solution of the present invention: geothermal gradient isogram, underground heat in the step B and step C
Computing unit divides and digitizing solution is will to calculate the corresponding geothermal gradient isogram of target zone, structural map in Croel
The alignment of DRAW software, and divide identical statistic unit；Above-mentioned two figure is digitized using Surfer, at reading unit node
Corresponding plane coordinates, depth, geothermal gradient.
As further technical solution of the present invention: the heat reservori top surface ground temperature isopleth preparation method in the step D
It is to utilize the node buried depth read in step B and step C, geothermal gradient, the corresponding target zone ground temperature of calculate node, and node is made
For virtual well point；Target zone temperature data in actual well drilled data is read, with virtual log collectively as control point, by inserting between well
Value obtains ground temperature isogram；Retain the ground temperature statistic unit in ground temperature isogram, calculates geothermal reservoirs as subsequent sub-unit
The basis of amount.
As further technical solution of the present invention: in the step A: observe the oil field continuous temperature measurement curve that is collected into
The variation tendency of depth, the identical position of the slope of curve is divided into well section, any one section of geothermal gradient G meets: G=100*
(T1-T0)/(H1-H0), wherein H1、H0For depth, T1、T0For the corresponding temperature of respective depth, by can on the continuous temperature measurement curve of well point
To read.
As further technical solution of the present invention: the step C geothermal gradient isogram, target zone top surface structural map
In division unit method, by two kinds of figures according to the same position of well point coordinate, equal proportion ruler pair in Croel DRAW mapping software
Together, figure layer grid division is separately built, the coordinate range of the same grid control fits like a glove, and the node coordinate of grid passes through Surfer
Software digital obtains, and the corresponding geothermal gradient of node can be read in geothermal gradient isogram, and depth can construct
It is read in figure, under the conditions of known to node depth, the geothermal gradient, corresponding well point target zone ground temperature T=G* (D-H)+T0Meter
It obtains, wherein G is geothermal gradient, and D is heat reservori depth, and H is constant temperature layer depth, T0For layer of constant temperature.
Compared with prior art, the beneficial effects of the present invention are: the present invention is using well point as virtual log and practical probing well
Target zone temperature is as control point in point, zone of ignorance by gram in golden difference obtain Temperature Distribution isogram.Virtual well point
Number is more, is distributed in each tectonic position in basin, avoids practical probing well and is concentrated mainly on structural high part to temperature
It is inconvenient that isogram works out bring.Meanwhile the temperature isogram worked out and structural map, geothermal gradient isogram have
There is corresponding grid cell, sub-unit calculates hot energy storage, and the cumulative hot energy storage of layer as a purpose of each unit calculated result is had
There is higher accuracy.
Detailed description of the invention
Fig. 1 is the flow chart for oil field old water-bearing layer top surface temperature isopleth preparation method of the embodiment of the present invention；
Fig. 2 is part well point (9 well points) continuous temperature measurement curve graph of the embodiment of the present invention；
Fig. 3 is basin (0~3000m) the geothermal gradient isopleth and ground temperature dividing elements figure of the embodiment of the present invention；
Fig. 4 is that the Basin Tectonic of the embodiment of the present invention and ground temperature statistic unit divide figure；
Fig. 5 is the basin heat reservori top surface temperature isogram of the embodiment of the present invention.
Specific embodiment
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation description, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, it is obtained by those of ordinary skill in the art without making creative efforts every other
Embodiment shall fall within the protection scope of the present invention.
Embodiment 1: Fig. 1-5 is please referred to, a kind of oil area heat reservori top surface temperature isogram preparation method, detailed process is such as
Under: the existing well point continuous temperature measurement data in oil field step A, is collected, observation well point temperature calculates well point with the changing rule of depth
The corresponding geothermal gradient of different depth section；Step B, it using the geothermal gradient of thermometric well point as control point, is inserted between well by Ke Lijin
Value method draws different depth section geothermal gradient distribution map；Step C, target zone top surface structural map is collected, while top surface being constructed
Figure, geothermal gradient isogram carry out dividing elements, and the corresponding cell position of two sets of figures, control range fit like a glove；It will divide
Geothermal gradient isogram, structural map after unit are digitized, and the plane coordinates at reading unit node is put into table, together
When read geothermal gradient, the depth of respective point, calculate the corresponding temperature of the corresponding points depth segment；Step D, using cell node as
The data such as virtual controlling point, well point respective layer temperature obtain target zone ground temperature etc. by Kriging regression between well as control point
It is worth line chart.
In the geothermal gradient calculation method of the well point step A, the oil field continuous temperature measurement curve that is collected into is observed with the change of depth
The identical position of the slope of curve is divided well section by change trend, and any one section of geothermal gradient G meets: G=100* (T1-T0)/(H1-
H0), wherein H1、H0For depth, T1、T0For the corresponding temperature of respective depth, by that can be read on the continuous temperature measurement curve of well point.
In the step B geothermal gradient isogram preparation method that the above embodiment of the present invention provides, control well point is selected
Geothermal gradient is the geothermal gradient being calculated after the segment processing of well point, respective depth section geothermal gradient isogram compilation process
In control well point all select corresponding well sections, it may be assumed that geothermal gradient numerical value of any one well point on corresponding interval is not by structure
Make the influence of position, Structural range, therefore arriving for Kriging regression method can be used between well, eliminate structural relief bring it is slotting
The inaccuracy of value.
It is divided in step C geothermal gradient isogram, the target zone top surface structural map that the above embodiment of the present invention provides single
In first method, two kinds of figures are aligned according to the same position of well point coordinate, equal proportion ruler in Croel DRAW mapping software, are separately built
The coordinate range of figure layer grid division, the same grid control fits like a glove, and the node coordinate of grid passes through Surfer software number
Word obtains, and the corresponding geothermal gradient of node can be read in geothermal gradient isogram, and depth can be read in structural map
It takes, under the conditions of known to node depth, the geothermal gradient, corresponding well point target zone ground temperature T=G* (D-H)+T0It is calculated,
Wherein G is geothermal gradient, D is heat reservori depth, H is constant temperature layer depth (taking 20m), T0For layer of constant temperature (taking 15 DEG C).
In the step D heat reservori top surface temperature isogram preparation method that the above embodiment of the present invention provides, grid section
The known target zone plane coordinates of point, the corresponding temperature of respective depth.Using well point as purpose in virtual log and practical probing well point
Layer temperature be used as control point, zone of ignorance by gram in gold difference obtain Temperature Distribution isogram.Virtual well point number is more,
Be distributed in each tectonic position in basin, avoid practical probing well be concentrated mainly on structural high part give temperature isogram compile
Bring processed is inconvenient.Meanwhile the temperature isogram worked out and structural map, geothermal gradient isogram have it is corresponding
Grid cell, sub-unit calculate hot energy storage, and the cumulative hot energy storage of layer as a purpose of each unit calculated result is had higher essence
Exactness.
Embodiment 2, on the basis of embodiment 1, Fig. 2 are 9 mouthfuls of well continuous temperature measurement moneys of Dong- pu Depression of the embodiment of the present invention
Material, concave interior shares 17 mouthfuls of wells, and there are continuous temperature measurement data, separately have part well point to have formation testing thermometric data, can be used as temperature
The supplement of degree control well point.Well point continuous temperature measurement curve is observed it can be found that the well point 0~3000m of depth segment temperature is with depth
Uniformly, less, in addition 3000~4300m, 4300m are opposite with deep well point geothermal gradient for the geothermal gradient variation of single well point for variation
Stablize, therefore 17 mouthfuls of wells are all divided into the ground of 0~3000m, 3000~4300m, 4300m using deep 3 well sections as well point
Warm gradient.Well section geothermal gradient G meets: G=100* (T1-T0)/(H1-H0), wherein H1、H0For depth, T1、T0For respective depth
Corresponding temperature, by that can be read on the continuous temperature measurement curve of well point.The distribution of 17 well 0~3000m well section geothermal gradients is calculated
Between 2.60-4.26 DEG C/100m, representative geothermal gradient is 3.40 DEG C/100m, and geothermal gradient is related with tectonic position, in
Entreat rumble strip Lithospheric Thickness smaller, corresponding geothermal gradient is bigger than normal.
Embodiment 3: Fig. 3 is the basin 0~3000m geothermal gradient isopleth and ground temperature dividing elements of the embodiment of the present invention
Figure.As previously mentioned, the well point 0~3000m ground temperature is uniform with change in depth, i.e., well section geothermal gradient is stablized, well interpolation not by
Structural relief influences, therefore well point geothermal gradient distribution map can be obtained by Kriging regression.Specific step is as follows:
(1) geothermal gradient figure is divided into relatively uniform grid, node 141 of grid can be used as virtual log；
(2) geothermal gradient isogram is imported into Surfer software, the optional two mouthfuls well points with coordinate, with plane coordinates
For control, figure is puted up in foundation, and corresponding well point is aligned with geothermal gradient isogram, by geothermal gradient isogram number
Change.
(3) 141 virtual log plane coordinates are read, while reading the geothermal gradient of respective point, it is standby to import Excel table
With.
Embodiment 4: Fig. 4 is that the Basin Tectonic of the embodiment of the present invention and ground temperature statistic unit divide figure.Basin Tectonic figure is
One of the indispensable graph that oil field prospecting deployment needs, can be obtained by data collection.By Basin Tectonic figure and geothermal gradient etc.
Being worth line chart importing, Croel DRAW mapping software, 2 figures adhere to different figure layers separately simultaneously, and 2 corresponding well points of figure are aligned, on
Lower layer's coordinate fits like a glove, and structural map is divided into identical unit with reference to above-mentioned geothermal gradient isopleth dividing elements method,
Cell node coordinate coincide, and reads 141 virtual well point target zone buried depths on structural map in a similar way, puts
Enter spare after the corresponding virtual log coordinate of Excel table.
Embodiment 5: Fig. 5 is the basin heat reservori top surface temperature isogram of the embodiment of the present invention.It is anti-with Dong- pu Depression T1
It penetrates for layer top surface temperature isogram establishment, it is as follows to introduce its preparation method:
(1) virtual well point geothermal gradient G, the purpose layer depth D for reading the above method substitute into formula T=G* (D-H)+T0
Target zone ground temperature T is calculated, H is constant temperature layer depth (taking 20m), T0For layer of constant temperature (taking 150C).
(2) temperature for nearby reading 141 arrived mouthful virtual log Geothermal data of calculating and 17 mouthfuls of well T1 depth is as control
It is processed, it is obtained wherein between well by Kriging regression method, T1Reflecting layer top surface ground temperature scatter chart.And be distributed ground temperature etc.
Value line chart, which is divided into, is divided into identical unit with geothermal gradient isogram, structural map, and the later period can calculate purpose with sub-unit
Layer Geothermal deposit, unit is cumulative then it can be concluded that target zone Geothermal deposit.
It is obvious to a person skilled in the art that invention is not limited to the details of the above exemplary embodiments, Er Qie
In the case where without departing substantially from spirit or essential attributes of the invention, the present invention can be realized in other specific forms.Therefore, no matter
From the point of view of which point, the present embodiments are to be considered as illustrative and not restrictive, and the scope of the present invention is by appended power
Benefit requires rather than above description limits, it is intended that all by what is fallen within the meaning and scope of the equivalent elements of the claims
Variation is included within the present invention.Any reference signs in the claims should not be construed as limiting the involved claims.
In addition, it should be understood that although this specification is described in terms of embodiments, but not each embodiment is only wrapped
Containing an independent technical solution, this description of the specification is merely for the sake of clarity, and those skilled in the art should
It considers the specification as a whole, the technical solutions in the various embodiments may also be suitably combined, forms those skilled in the art
The other embodiments being understood that.
Claims (7)
1. a kind of oil area heat reservori top surface temperature isogram preparation method, which is characterized in that comprise the steps of:
A, well point continuous temperature measurement data is collected, the identical position of temperature logs slope is divided into well section, segmentation calculates geothermal gradient；
B, using known well point as control condition, show that depth segment is correspondingly locating for heat reservori by Kriging regression method between well
Warm gradient isogram；
C, structure of oil field figure is collected, and is aligned in Croel DRAW software with geothermal gradient isogram, division fits like a glove
Underground heat computing unit；
D, using Surfer will construct with geothermal gradient distribution map digitize, the corresponding plane coordinates of reading unit node, buried depth,
Geothermal gradient numerical value, the corresponding ground temperature of calculate node and as virtual log, while actual well drilled data target zone temperature is read, it will
Virtual log and actual well drilled data obtain heat reservori top surface ground temperature isogram by interpolation between well as control point.
2. a kind of oily area heat reservori top surface temperature isogram preparation method according to claim 1, which is characterized in that step
In rapid A: by well point continuous temperature measurement data segment processing, the same well section continuous temperature measurement slope of curve is the same, corresponding geothermal gradient
Numerical value is constant.
3. a kind of oily area heat reservori top surface temperature isogram preparation method according to claim 2, which is characterized in that institute
It is corresponding with the well section of well point continuous temperature measurement data to state geothermal gradient isogram, the corresponding geothermal gradient isopleth of a certain depth segment
In figure compilation process, any one control well point geothermal gradient in the depth segment controls the ground between well point not with change in depth
Warm gradient is influenced less than normal by structural relief, is obtained using Kriging regression method.
4. a kind of oily area heat reservori top surface temperature isogram preparation method according to claim 3, which is characterized in that institute
State the geothermal gradient isogram in step B and step C, underground heat calculates dividing elements and digitizing solution is will to calculate target zone
Corresponding geothermal gradient isogram, structural map are aligned in Croel DRAW software, and divide identical statistic unit；It utilizes
Surfer digitizes above-mentioned two figure, corresponding plane coordinates, depth, geothermal gradient at reading unit node.
5. a kind of oily area heat reservori top surface temperature isogram preparation method according to claim 1, which is characterized in that institute
State the heat reservori top surface ground temperature isopleth preparation method in step D be using the node buried depth read in step B and step C,
Warm gradient, the corresponding target zone ground temperature of calculate node, and using node as virtual well point；Read target zone temperature in actual well drilled data
Degree evidence obtains ground temperature isogram by interpolation between well with virtual log collectively as control point；Retain in ground temperature isogram
Ground temperature statistic unit, as subsequent sub-unit calculate Geothermal deposit basis.
6. a kind of oily area heat reservori top surface temperature isogram preparation method according to claim 1, which is characterized in that institute
It states in step A: the variation tendency for observing the oil field continuous temperature measurement curve that is collected into depth, by the identical position of the slope of curve
Well section is divided, any one section of geothermal gradient G meets: G=100* (T1-T0)/(H1-H0), wherein H1、H0For depth, T1、T0For phase
The corresponding temperature of depth is answered, by reading on the continuous temperature measurement curve of well point.
7. a kind of oily area heat reservori top surface temperature isogram preparation method according to claim 1, which is characterized in that institute
Step C geothermal gradient isogram is stated, in the structural map division unit method of target zone top surface, in Croel DRAW mapping software
Two kinds of figures are aligned according to the same position of well point coordinate, equal proportion ruler, separately build figure layer grid division, the seat of the same grid control
Mark range fits like a glove, and the node coordinate of grid is obtained by Surfer software digital, and the corresponding geothermal gradient of node is on ground
It is read in warm gradient isogram, depth is read in structural map, under the conditions of known to node depth, the geothermal gradient, accordingly
Well point target zone ground temperature T=G* (D-H)+T0It is calculated, wherein G is geothermal gradient, and D is heat reservori depth, and H is thermostat layer
Depth, T0For layer of constant temperature.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910084363.0A CN109751044A (en) | 2019-01-29 | 2019-01-29 | A kind of oil area heat reservori top surface temperature isogram preparation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910084363.0A CN109751044A (en) | 2019-01-29 | 2019-01-29 | A kind of oil area heat reservori top surface temperature isogram preparation method |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109751044A true CN109751044A (en) | 2019-05-14 |
Family
ID=66406732
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910084363.0A Pending CN109751044A (en) | 2019-01-29 | 2019-01-29 | A kind of oil area heat reservori top surface temperature isogram preparation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109751044A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112767515A (en) * | 2021-02-19 | 2021-05-07 | 上海市政工程设计研究总院（集团）有限公司 | Method, system, equipment and medium for drawing ground temperature curve graph of test point |
| CN113835137A (en) * | 2020-06-23 | 2021-12-24 | 中国石油天然气股份有限公司 | Method and device for predicting oil-gas phase state of deep layer of basin |
| CN114016991A (en) * | 2021-09-16 | 2022-02-08 | 西安交通大学 | Method for determining well body structure of butt joint well based on geothermal field distribution characteristics |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1854762A (en) * | 2005-04-21 | 2006-11-01 | 中国石油天然气集团公司 | Rectangular net gridding method for painting contour graph containing rift geological structure |
| US20070186640A1 (en) * | 2006-02-16 | 2007-08-16 | Johnson David O | Single point and fiber optic temperature measurement for correction of a gas column weight in a well |
| CN103914618A (en) * | 2014-03-28 | 2014-07-09 | 哈尔滨工程大学 | Electronic chart based seawater depth contour line generating method |
| US20150292956A1 (en) * | 2014-04-14 | 2015-10-15 | Baker Hughes Incorporated | Distributed temperature sensor enhancement by stimulated raman suppression |
| CN106327578A (en) * | 2015-06-29 | 2017-01-11 | 北京航天试验技术研究所 | Three-dimensional temperature field interpolation method based on different physical mediums |
-
2019
- 2019-01-29 CN CN201910084363.0A patent/CN109751044A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1854762A (en) * | 2005-04-21 | 2006-11-01 | 中国石油天然气集团公司 | Rectangular net gridding method for painting contour graph containing rift geological structure |
| US20070186640A1 (en) * | 2006-02-16 | 2007-08-16 | Johnson David O | Single point and fiber optic temperature measurement for correction of a gas column weight in a well |
| CN103914618A (en) * | 2014-03-28 | 2014-07-09 | 哈尔滨工程大学 | Electronic chart based seawater depth contour line generating method |
| US20150292956A1 (en) * | 2014-04-14 | 2015-10-15 | Baker Hughes Incorporated | Distributed temperature sensor enhancement by stimulated raman suppression |
| CN106327578A (en) * | 2015-06-29 | 2017-01-11 | 北京航天试验技术研究所 | Three-dimensional temperature field interpolation method based on different physical mediums |
Non-Patent Citations (5)
| Title |
|---|
| 吕国祥等: "《计算机地质制图及应用》", 31 October 2009 * |
| 林世辉等: "冀中坳陷现今地温场分布特征", 《东华理工学院学报》 * |
| 肖卫勇等: "渤海盆地地温场研究", 《中国海上油气（地质）》 * |
| 赵苏民等: "《沉积盆地型地热田勘查开发与利用》", 31 December 2013 * |
| 龚露: "地温场二维有限元数值模拟", 《中国优秀硕士学位论文全文数据库（基础科学辑）》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113835137A (en) * | 2020-06-23 | 2021-12-24 | 中国石油天然气股份有限公司 | Method and device for predicting oil-gas phase state of deep layer of basin |
| CN113835137B (en) * | 2020-06-23 | 2023-08-22 | 中国石油天然气股份有限公司 | Method and device for predicting deep oil-gas phase state of basin |
| CN112767515A (en) * | 2021-02-19 | 2021-05-07 | 上海市政工程设计研究总院（集团）有限公司 | Method, system, equipment and medium for drawing ground temperature curve graph of test point |
| CN114016991A (en) * | 2021-09-16 | 2022-02-08 | 西安交通大学 | Method for determining well body structure of butt joint well based on geothermal field distribution characteristics |
| CN114016991B (en) * | 2021-09-16 | 2022-10-21 | 西安交通大学 | Method for determining well body structure of butt joint well based on geothermal field distribution characteristics |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105651676B (en) | A kind of reservoir heterogeneity characterizing method under horizontal well rule well pattern | |
| CN109751044A (en) | A kind of oil area heat reservori top surface temperature isogram preparation method | |
| Le Brocq et al. | A subglacial water-flow model for West Antarctica | |
| CN108846245A (en) | City-scale geothermal field multi-boreholes efficient numerical analogy method and device | |
| CN106991244B (en) | A method of Fracture Networks connectivity and seepage calculation based on graph theory | |
| CN103778298B (en) | The multi-level finite element modeling method of two dimension flow motion in the simulation porous media improved | |
| CN109614753A (en) | A kind of heat hiding engineering evaluation method for geothermal energy resources | |
| CN109003330A (en) | A kind of three dimensional contour line method based on basement rock boundary constraint | |
| CN110056346A (en) | A kind of three-dimensional original water saturation analogy method of the oil reservoir based on Long-term change trend function | |
| CN105242328A (en) | Paleo-thermal lithosphere thickness determination method and device | |
| CN102682195A (en) | Mine shaft temperature calculation method for semi-submersible type platform transient drilling well | |
| Chen et al. | A 3-D prognostic numerical model study of the Georges Bank ecosystem. Part I: physical model | |
| Bansal et al. | A strongly coupled, fully implicit, three dimensional, three phase reservoir simulator | |
| Gerdes et al. | Large‐scale circulation and water mass distribution in the Arctic Ocean from model results and observations | |
| CN106484930A (en) | For determining the method and system of heterogeneous reservoir labyrinth well production | |
| CN102288144B (en) | Quantitative method for determining average thickness of water-bearing bed | |
| Vosgerau et al. | Towards a geothermal exploration well in the Gassum Formation in Copenhagen | |
| CN105912861B (en) | The computational methods of the conceptual hydrological model of nonparametric based on GIS | |
| Pawlowicz et al. | Diagnosing vertical mixing in a two‐layer exchange flow | |
| CN105698968A (en) | Method and system for observing distribution change of water temperature after reservoir muddy water density flows into sink | |
| Qi et al. | Research on vertical distribution of water temperature in different regulation reservoirs | |
| Sahagian | Ocean temperature-induced change in lithospheric thermal structure: A mechanism for long-term eustatic sea level change | |
| Goode et al. | Kinematic implications of consequent channels on growing folds | |
| Jaquet et al. | Groundwater flow modelling under transient ice sheet conditions in Greenland | |
| Dunlap et al. | Geohydrology and model analysis of stream-aquifer system along the Arkansas River in Kearny and Finney counties, southwestern Kansas |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190514 |
|
| RJ01 | Rejection of invention patent application after publication |