CN106052891B - Tube wall dynamic temperature inner is distributed method of real-time - Google Patents
Tube wall dynamic temperature inner is distributed method of real-time Download PDFInfo
- Publication number
- CN106052891B CN106052891B CN201610614916.5A CN201610614916A CN106052891B CN 106052891 B CN106052891 B CN 106052891B CN 201610614916 A CN201610614916 A CN 201610614916A CN 106052891 B CN106052891 B CN 106052891B
- Authority
- CN
- China
- Prior art keywords
- temperature
- variogram
- tube wall
- monitoring point
- kriging
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K7/00—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
- G01K7/02—Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements using thermoelectric elements, e.g. thermocouples
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
The present invention relates to a kind of tube wall dynamic temperature inners to be distributed method of real-time, including：Acquire the temperature of each monitoring point；Experimental temperature variogram between calculating monitoring point；Spherical variogram model is selected, Experiment variogram is fitted, obtains theoretical variogram；Central point by sample plane grid division, each grid is future position；Calculate the theoretical variogram value between theoretical variogram value and each monitoring point and the future position in prediction vertex neighborhood between monitoring point；Kriging equation groups are solved, Kriging weighting coefficients are obtained；Calculate the Kriging interpolated temperatures of future position；Obtain the Temperature Distribution cloud atlas of sample plane.The present invention can realize that simple in structure, the higher tube wall interior temperature distribution of precision is acquired with relatively low cost.
Description
Technical field
The present invention relates to a kind of tube wall dynamic temperature inners to be distributed method of real-time.Present invention can apply to pipeline wall surfaces
Interior temperature distribution fields of measurement.
Background technology
Temperature is one of most important procedure parameter, and the DYNAMIC DISTRIBUTION for monitoring pipe surface and internal temperature will be helpful to depth
Enter to understand the characteristics such as the heat transfer of pipeline.Currently, existing pipe surface Temperature Distribution research is mainly realized using infrared principles
Measurement to pipe surface temperature, real-time is preferable, and technology is also more mature.And it is most at present for tube wall interior temperature distribution
By the way of being implanted into temperature sensor, this mode cost is relatively low, has a wide range of applications, but there is also some problems.
Tube wall interior temperature distribution is monitored using implantation temperature sensor, generally uses thermocouple or thermoelectricity
Resistance.Wherein thermal resistance measurement accuracy is higher, but the response time of thermal resistance is generally long compared with thermocouple, in the field for having requirement of real-time
Conjunction cannot be met the requirements；Meanwhile existing multi way temperature signal acquisition mode does not ensure that the consistency of every road signal, therefore
Obtained Temperature Distribution is not accurate enough；Meanwhile host computer can only monitor the temperature of sampled point, can not comprehensively reflect pipe in real time
The Temperature Distribution of pars intramuralis has certain limitation.
Invention content
To solve the above problems, the present invention provides a kind of tube wall dynamic temperature inner distribution method of real-time, with opposite
Lower cost realizes simple in structure, the higher tube wall interior temperature distribution acquisition of precision.In order to achieve the above objectives, the present invention adopts
With following technical scheme：
A kind of tube wall dynamic temperature inner distribution method of real-time, used array of temperature sensor includes multiple heat
Thermocouple probe is implanted into tube wall by thermocouple probe, and the tube wall plane residing for thermocouple probe is known as " sample plane ", meanwhile, dress
Set and be furnished with cold junction temperature sensor, be placed in thermocouple cold junction, for measuring thermocouple cold junction, each thermocouple probe and
The temperature signal of cold junction temperature sensor acquisition is admitted to computer, and the monitoring method includes the following steps：
(1) it is equipped with n temperature monitoring point (x altogether in a sample planei,yi), i=1,2 ..., n, the temperature of monitoring point
For T (xi,yi), if the distance between two monitoring points of i points and j points are hij, j=1,2 ..., n；
(2) experimental temperature variogram γ between calculating monitoring point*(hij)：
(3) spherical variogram model is selected, Experiment variogram is fitted, theoretical variogram γ is obtained
(hij)：
Wherein C is base station value, and a is to become journey；
(4) sample plane is divided into m grid, the central point of each grid is future position；
(5) future position (x is calculated0,y0) theoretical variogram value γ (h in neighborhood between monitoring pointij) and each monitoring
Theoretical variogram value γ (h between point and future position0j), it is abbreviated as γ respectivelyijAnd γ0j；
(6) Kriging equation groups are solved, Kriging weighting coefficients λ is obtainedi, i=1,2 ..., n；
[λ]=[M] [K]-1
Wherein
Wherein, μ is Lagrangian number multiplier, λiFor undetermined coefficient；
(7) future position (x is calculated0,y0) Kriging interpolated temperatures T*(x0,y0)：
(8) step (5) to step (7) is repeated, until traversing all future positions of sample plane.To obtain sample plane
Temperature Distribution cloud atlas.
The good effect that the present invention has on the basis of the above is：
1) tube wall dynamic temperature inner distribution method of real-time provided by the invention can accurately measure tube wall dynamic temperature
Degree has preferable consistency；Using Kriging interpolation algorithms can real-time display Temperature Distribution cloud atlas, make monitoring process more
Add it is intuitive and convenient, meet high-precision, multizone remote temperature distribution monitoring；
2) tube wall dynamic temperature inner distribution method of real-time provided by the invention realizes hardware acquisition system and meter
The height of calculation machine terminal merges, and reduces the complexity of hardware, reduces cost, reduces the system research and development period；
3) tube wall dynamic temperature inner distribution method of real-time provided by the invention is compiled using LabVIEW graphic interfaces
Journey has the functions such as data storage and temperature alarming, records data for operator and safeguards that data provide a convenient.It can be wide
It is general to be applied to pipeline wall surface and other interior of articles temperature distributing measuring fields.
Description of the drawings
Fig. 1 is that the tube wall dynamic temperature inner of the embodiment of the present invention is distributed the schematic diagram of real-time measurement apparatus.
Fig. 2 is that the tube wall dynamic temperature inner of the embodiment of the present invention is distributed the structure of real-time measurement apparatus Signal-regulated kinase
Figure.
Fig. 3 is that the tube wall dynamic temperature inner of the embodiment of the present invention is distributed the program flow diagram of method of real-time.
Fig. 4 is that the tube wall dynamic temperature inner of the embodiment of the present invention is distributed thermoelectrical potential signal processing of method of real-time
Program flow diagram.
The Kriging interpolation algorithms of the tube wall dynamic temperature inner distribution method of real-time of the positions Fig. 5 embodiment of the present invention
Flow chart.
Specific implementation mode
With reference to the accompanying drawings and examples, the invention will be further described.
The tube wall dynamic temperature inner being related to as shown in Figure 1 for the present embodiment is distributed real-time measurement apparatus schematic diagram, including
Embedded temperature detector array 1, signal selection module 2, Signal-regulated kinase 3 and multi-channel data acquisition and control module 4
And computer 5.Embedded temperature detector array 1 is for measuring tube wall temperature data, by signal selection module 2, signal
Conditioning module 3 and multi-channel data acquisition are transferred to computer 5 with control module 4, pass through the host computer LabVIEW of computer 5
Program carries out processes temperature signal, and real-time display goes out current tube wall interior temperature distribution cloud atlas, realizes the distribution of tube wall dynamic temperature
Monitoring in real time.
Embedded temperature detector array 1 is passed by the plug-in high-precision cold junction temperature of multiple thermocouple arrays 6 and one
Sensor 7 forms.Thermocouple arrays 6 are distributed according to certain distance inside tube wall, and are fixed by the bracket, and the plane of distribution claims
For sample plane.Thermoelectrical potential signal is connected by the special extended line of thermocouple with 2 input terminal of signal selection module.Thermocouple cold junction
Temperature is measured using high-precision cold end temperature sensor 7, amplifying circuit of the cold junction temperature signal in Signal-regulated kinase 3
10 are amplified processing, subsequently input computer 5, and cold junction compensation is carried out to thermoelectrical potential using software.
Multiple thermoelectrical potential signals are connected with 2 input terminal of signal selection module, and the main body of signal selection module 2 selects for data
Device 8, the address of Enable Pin input 5 Program of the computer control of data selector 8, which is switched with certain frequency, is selected
Output of the thermoelectrical potential signal of appropriate address as signal selection module 2.Signal selection module 2 makes all thermocouple signals be multiplexed
Same modulate circuit, enhances signal conformance, and distribution measuring precision is made to greatly improve.
The thermoelectrical potential signal of appropriate address port is output and then enter putting in Signal-regulated kinase 3 from signal selection module 2
Big filter circuit 9, as shown in Figure 2.Signal amplification circuit includes two-stage amplifying circuit, and level-one amplifying circuit is using differential amplification electricity
Road 12, second amplifying circuit use in-phase proportion amplifying circuit 13.Signal filter circuit 14 is filtered using Butterworth step low-pass
Wave device filters off the interference signal and power frequency component of high frequency.Meanwhile the thermocouple cold junction that temperature sensors of high precision 7 measures
In-phase proportion amplifying circuit 15 is also passed through, computer is inputted by data collecting card.
High-speed data acquisition card is connect by USB interface with computer 5, the analog quantity that can export Signal-regulated kinase 3
The digital quantity that computer 5 can identify is converted to, and the discrete control signal that computer 5 exports is output to signal behavior
Module 2 realizes software and hardware interconnection.
Computer 5 is programmed using LabVIEW, it can be achieved that signal acquisition, signal processing, real-time display and storage and alarm
Etc. functions.System main program flow chart is as shown in Figure 3, it is necessary first to cold junction temperature and environment temperature is acquired, then to thermoelectrical potential
Signal is handled, and is converted thermoelectrical potential voltage signal to temperature signal and is carried out software cold junction compensation, obtains monitoring point temperature
After data, interpolation calculation is carried out to the temperature in sample plane using Kriging interpolation algorithms, and result cloud atlas is shown in real time
Show that, in host computer interface, thermoelectrical potential signal processing subroutine flow chart is as shown in Figure 4.Wherein, Kriging interpolation algorithms flow chart
As shown in figure 5, including the following steps：
(1) it is equipped with n temperature monitoring point altogether in a sample plane, each temperature monitoring point two-dimensional coordinate is (xi,
yi), the temperature of i=1,2 ..., n, each monitoring point are T (xi,yi), then the distance between any two monitoring point hijFor：
Wherein, i=1,2 ..., n, j=1,2 ..., n.
(2) each experimental temperature variogram γ under is calculated*(hij), it is defined as the difference of the temperature value at 2 points
The half of variance：
(3) spherical variogram model is selected, Experiment variogram is fitted, theoretical variogram γ is obtained
(hij), the form of spherical variogram is：
Wherein C is base station value, and a is to become journey.Through program test, the spherical variogram model is to the tube wall in the present embodiment
The prediction of Temperature Distribution has good effect.In other embodiments, other variograms can be also selected according to actual conditions
Theoretical model, such as Gauss model, Triangle Model.
(4) sample plane is divided into m grid, the central point of each grid is future position.
(5) future position (x is calculated0,y0) theoretical variogram value γ (h in neighborhood between monitoring pointij) and each monitoring
Theoretical variogram value γ (h between point and future position0j), it is abbreviated as γ respectivelyijAnd γ0j。
(6) Kriging equation groups are solved, Kriging weighting coefficients λ is obtainedi, i=1,2 ..., n.
[λ]=[M] [K]-1
Wherein
Wherein, μ is called Lagrangian number multiplier, and numerical value can be acquired by Kriging equation groups.
(7) future position (x is calculated0,y0) Kriging interpolated temperatures T*(x0,y0)：
(8) step (5) to step (7) is repeated, until traversing all future positions of sample plane.To obtain sample plane
Temperature Distribution cloud atlas.
The present invention is exemplarily described above in conjunction with attached drawing, it is clear that the foregoing is merely the preferable implementations of the present invention
Example, specific implementation of the invention are not subject to the restrictions described above, and that is done all within the spirits and principles of the present invention is any
Modifications, equivalent substitutions and improvements etc., should all be included in the protection scope of the present invention.
Claims (1)
1. a kind of tube wall dynamic temperature inner is distributed method of real-time, used array of temperature sensor includes multiple thermoelectricity
Even probe, tube wall is implanted by thermocouple probe, and the tube wall plane residing for thermocouple probe is known as " sample plane ", meanwhile, device
Equipped with cold junction temperature sensor, it is placed in thermocouple cold junction, for measuring thermocouple cold junction, each thermocouple probe and cold
The temperature signal of end temperature sensor acquisition is admitted to computer, and the monitoring method includes the following steps：
(1) it is equipped with n temperature monitoring point (x altogether in a sample planei,yi), the temperature of i=1,2 ..., n, monitoring point are T
(xi,yi), if the distance between two monitoring points of i points and j points are hij, j=1,2 ..., n；
(2) experimental temperature variogram γ between calculating monitoring point*(hij)：
(3) spherical variogram model is selected, Experiment variogram is fitted, theoretical variogram γ (h are obtainedij)：
Wherein C is base station value, and a is to become journey；
(4) sample plane is divided into m grid, the central point of each grid is future position；
(5) future position (x is calculated0,y0) theoretical variogram value γ (h in neighborhood between monitoring pointij) and each monitoring point with it is pre-
Theoretical variogram value γ (h between measuring point0j), it is abbreviated as γ respectivelyijAnd γ0j；
(6) Kriging equation groups are solved, Kriging weighting coefficients λ is obtainedi, i=1,2 ..., n；
[λ]=[M] [K]-1
Wherein
Wherein, μ is Lagrangian number multiplier；
(7) future position (x is calculated0,y0) Kriging interpolated temperatures T*(x0,y0)：
(8) step (5) to step (7) is repeated, until traversing all future positions of sample plane, to obtain the temperature of sample plane
Cloud charts.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610614916.5A CN106052891B (en) | 2016-07-26 | 2016-07-26 | Tube wall dynamic temperature inner is distributed method of real-time |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610614916.5A CN106052891B (en) | 2016-07-26 | 2016-07-26 | Tube wall dynamic temperature inner is distributed method of real-time |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106052891A CN106052891A (en) | 2016-10-26 |
| CN106052891B true CN106052891B (en) | 2018-09-25 |
Family
ID=57197051
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610614916.5A Expired - Fee Related CN106052891B (en) | 2016-07-26 | 2016-07-26 | Tube wall dynamic temperature inner is distributed method of real-time |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106052891B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109506800A (en) * | 2018-12-18 | 2019-03-22 | 重庆邮电大学 | Dump temperature measurement system based on Thermistor Temperature Measurement |
| CN110430625A (en) * | 2019-06-10 | 2019-11-08 | 深圳桐源科技有限公司 | Aerosol generating system, heating module, temperature checking method and device |
| CN111159857B (en) * | 2019-12-13 | 2024-02-13 | 天津大学 | Two-dimensional transient temperature field reconstruction method for sonic nozzle pipe wall |
| CN112577671B (en) * | 2020-11-27 | 2022-11-01 | 武汉工程大学 | Well lid monitoring method and system by using kriging method |
| CN113344286A (en) * | 2021-06-28 | 2021-09-03 | 北京工业大学 | Method and device for predicting indoor temperature distribution |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS55149025A (en) * | 1979-05-11 | 1980-11-20 | Shisaka Kenkyusho:Kk | Internal temperature measuring method |
| JP3337306B2 (en) * | 1994-03-10 | 2002-10-21 | 三菱重工業株式会社 | Temperature measurement method |
| CN102840930B (en) * | 2012-08-21 | 2014-06-04 | 清华大学 | Pipeline internal temperature measuring device |
| CN104062034B (en) * | 2014-07-04 | 2016-05-11 | 深圳市太科检测有限公司 | A kind of based on the thermometric contactless tube side fluid thermometry of pipeline outer wall |
-
2016
- 2016-07-26 CN CN201610614916.5A patent/CN106052891B/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN106052891A (en) | 2016-10-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106052891B (en) | Tube wall dynamic temperature inner is distributed method of real-time | |
| CN108444653B (en) | A kind of face formation ultrasonic gas leak point positioning detection device and method | |
| CN201716161U (en) | System for detecting position of dam leakage channel by taking temperature as tracer | |
| CN103162901A (en) | Nonlinear calibrating method for multiple temperature points of pressure sensor | |
| CN105067125B (en) | Cannon high speed colorimetric infrared temperature measurement apparatus and method | |
| CN104697668B (en) | Temperature measuring device and method | |
| CN101858991A (en) | System and method for detecting position of dam leakage passage by using temperature as tracer | |
| CN108287030B (en) | Method for measuring surface heat flow of embedded thermocouple | |
| CN108956083A (en) | A kind of test method for fast implementing wind-tunnel and surveying pressure on a small scale | |
| CN102967624B (en) | Device for testing Seebeck coefficient | |
| CN105548917B (en) | Calibration method of magnetic sensor array in non-shielding environment | |
| CN115113126B (en) | Device and method for testing and calibrating metal Hall probe | |
| CN110274660A (en) | Liquid amount detection method, device and system and household appliance | |
| CN105548246A (en) | Heat-conductivity-coefficient measuring experiment system through steady-state method and measuring method | |
| CN110542449A (en) | calibration system and calibration method of novel temperature and humidity sensor | |
| CN208458972U (en) | A kind of manometric module for fast implementing wind-tunnel and surveying pressure on a small scale | |
| CN109633494B (en) | Spacecraft magnetic field distribution information imaging method | |
| CN111157139A (en) | Visual measurement method for temperature distribution of single-connected combustion field | |
| CN105372288B (en) | A kind of rate of heat flow measuring instrument and measuring method | |
| CN113612552A (en) | Microwave radiometer receiver port amplitude inconsistency testing system and method | |
| CN201387363Y (en) | Object internal temperature filed distribution measuring instrument | |
| Wahab et al. | Image reconstruction for solid profile measurement in ERT using non-invasive approach | |
| CN204576333U (en) | A kind of isolation of the signal for engines ground TT&C system conditioning case | |
| CN107843347A (en) | A kind of porous media three dimensional temperature distribution measurement method | |
| CN102809445B (en) | Thermocouple detection method |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180925 Termination date: 20210726 |