CN109579685A - A kind of frost penetration calculation method and device - Google Patents
A kind of frost penetration calculation method and device Download PDFInfo
- Publication number
- CN109579685A CN109579685A CN201811544248.9A CN201811544248A CN109579685A CN 109579685 A CN109579685 A CN 109579685A CN 201811544248 A CN201811544248 A CN 201811544248A CN 109579685 A CN109579685 A CN 109579685A
- Authority
- CN
- China
- Prior art keywords
- soil
- layer
- frequency
- current time
- degree
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/02—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring length, width or thickness
- G01B7/06—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring length, width or thickness for measuring thickness
- G01B7/08—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring length, width or thickness for measuring thickness using capacitive means
- G01B7/082—Height gauges
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
The present invention relates to frozen soil field of measuring technique, especially a kind of frost penetration calculation method and device.The information that the sensor for detecting frost penetration detects is obtained, the frequency of every layer of soil of current time is obtained；The frequency of every layer of soil of current time is made comparisons with setting reference frequency, obtain every layer of soil of current time freezes degree；Freeze degree according to every layer of soil of current time, and freeze the actual depth of the corresponding soil of degree, obtains the depth of freezing of current time soil.The corresponding frequency of soil is measured by sensor, the degree of freezing of soil is obtained compared with setting reference frequency, and carries out the calculating of depth, this method accuracy is higher, can accurately analyze frost penetration, and be simply easily achieved.
Description
Technical field
The present invention relates to frozen soil field of measuring technique, especially a kind of frost penetration calculation method and device.
Background technique
Frozen soil refers to zero degrees celsius hereinafter, and various soil and rock containing ice.Frozen soil is boundless in China's distribution,
The total score cloth area of ever-frozen ground and seasonal frozen ground almost accounts for 2/3rds of entire national territorial area.The variation of temperature decides
Many processes such as formation, development and the degeneration and ablation of frozen soil more influence variation and the earth frozen soil circle of thickness of frozen earth layer
Develop, as global warming temperature increases, permafrost degradation will become inevitable, in Frozen Ground Area the stability of engineering project,
Weather conditions variation just becomes the hot topic of research.Precise measurement frost penetration, studies it with the development law of temperature change,
The recession situation of movement of tundra is grasped, can be studied for engineering in permafrost regions construction, agricultural structure adjustment, climate change etc.
Basic data is provided and supports that there is great scientific meaning.
In China, frozen soil observation in meteorology automatic Observation field is not able to achieve automation so far.Currently, the jelly used both at home and abroad
Native observation method mainly has three classes, first is that freezing meter observation method, surveys the rubber equipped with distilled water for being placed on underground by manually touching
Sebific duct freezes situation and carrys out measuring frozen, since water and specific heat capacity difference are larger, adds the soil texture, the content of organic matter and acid
The influence of the factors such as basicity, soil freezing temperature and pure water be not identical, in addition, carrying out freezing meter observation taking out rubber tube
When, it will receive outside air temperature influence, destroy thermal balance in rubber tube, therefore, freezing meter observation frozen soil is easy to appear deviation,
It is different surely to reflect practical frost penetration；Second is that the direct method of measurement (such as hand digging or drilling), has measurement directly, intuitively
Advantage, but pure manual operations is wanted, large labor intensity is observed, data observation density is inadequate, is unable to real-time monitoring soil freezing-thawing depth
And its development and change, and observation person's subjective factor is affected；Third is that remote sensing method, it is realized to the region even big ruler in the whole world
Frost penetration observation is spent, but since frozen soil is influenced by ambient temperature, is often existed with multilayer form, remote sensing method is difficult to observe
Thawing situation inside frozen soil layer, hardly possible have met the frozen soil observation requirements of fining.
It is existing in order to know the frozen state of soil on the basis of summarizing various frozen soil detection method advantage and disadvantage at present
Occurs the freeze thawing state that soil is detected using frozen soil sensor in technology.Frozen soil sensor is by multi-layer capacity sensor group
At, multiple groups plane capacitance is arranged in the same plane according to a certain distance, using its edge effect, in addition measuring circuit,
Soil is divided into dry ground, wet soil according to soil moisture content, if the soil moisture reduces to a certain extent, might have frozen soil, soil
When earth freezes or thaws, internal moisture can be undergone phase transition, and soil overall dielectric constant can be caused to change, when sensor external is situated between
When matter changes, the state of each layer sensor surrounding soil just can determine that, can carry out the differentiation of frozen soil accordingly.But only root
The state of soil can not accurately be gone out by sensor Direct Recognition according to the variation of dielectric constant, and then can not accurately be identified
The depth of frozen soil is unfavorable for the use of the frozen soil sensor.
Summary of the invention
The object of the present invention is to provide a kind of frost penetration calculation method and device, the frozen soil to solve the prior art is deep
Spend the problem of calculation method inaccuracy.
In order to realize that frost penetration calculates, the frost penetration calculation method inaccuracy of the prior art is solved the problems, such as.This hair
It is bright that a kind of frost penetration calculation method is provided, include the following steps:
1) information that the sensor for detecting frost penetration detects is obtained, the frequency of every layer of soil of current time is obtained；
2) frequency of every layer of soil of current time is made comparisons with setting reference frequency, obtains every layer of soil of current time
Freeze degree；
3) degree is freezed according to every layer of soil of current time, and freezes the actual depth of the corresponding soil of degree, obtained
To the depth of freezing of current time soil.
Beneficial effect is to measure the corresponding frequency of soil by sensor, obtain soil compared with setting reference frequency
Freeze degree, and carry out the calculating of depth, this method accuracy is higher, can accurately analyze frost penetration, and simple
It is easily achieved.
Further, in order to enable the setting reference frequency arrived is more reasonable, and meet the natural law, described in step 2)
The acquisition process for setting reference frequency is as follows:
Obtain for detect frost penetration sensor detect information, obtain different moments different layers soil temperature and
Frequency；
Look for and meet frequency corresponding to the soil horizon of following temperature condition: the temperature of soil is in setting solidification point range
It is interior, and the temperature at continuous N number of moment is downward trend and the temperature at this N number of moment is greater than the first set temperature value；
Frequency corresponding to all soil horizons for meeting the temperature condition is corrected, and then obtains the reference frequency.
Further, described to correct all soil for meeting the temperature condition in order to more accurately obtain reference frequency
Layer corresponding to frequency, and then obtain the reference frequency include: will be corresponding to all soil horizons for meeting the temperature condition
Frequency be fitted, be fitted to frequency corresponding to critical-temperature, seek being averaged for frequency corresponding to all critical-temperatures
Value, obtains the reference frequency.
Further, for better, simpler differentiation frozen soil state, the process of step 2) includes:
If the frequency of current time layer soil is more than or equal to the first setting value of reference frequency, the current time layer
The degree of freezing of soil is fully charge；
If the frequency of current time layer soil is less than or equal to the second setting value of reference frequency, the current time layer
The degree of freezing of soil is not freeze；
Current time remainder layer soil freezes between degree do not freeze between fully charge and.
It further, further include according to difference after step 2) in order to more accurately distinguish the degree of freezing of every layer of soil
The temperature of every layer of soil of moment freezes the step of degree is modified to every layer of soil of obtained current time: if current time
The temperature of certain layer of soil is greater than the second set temperature value, then the degree of freezing for correcting this layer of soil of current time is not freeze；If
The temperature of current time layer soil is less than third set temperature value, then the degree of freezing for correcting this layer of soil of current time has been
Freeze entirely.
It further, further include root after step 2) in order to which the degree of freezing more accurately to every layer of soil distinguishes
The step of degree is modified is freezed to every layer of soil of obtained current time according to center is freezed:
Freeze degree according to every layer of soil, obtain freezing the corresponding number of plies in center and freezes total depth；
The frozen state for setting the soil horizon in freezing range centered on freezing center is corrected as fully charge；Its
In, the upper limit for setting freezing range, which is subtracted to freeze the corresponding number of plies in center, freezes the setting multiple of total depth, described to set
The lower limit for determining freezing range is to freeze the corresponding number of plies in center to add the setting multiple for freezing total depth.
It further, further include according to upper a period of time after step 2) in order to more accurately correct the frozen state of every layer of soil
The frozen state of the upper layer of soil of certain layer of frozen state and this layer of soil is carved to freeze to every layer of soil of obtained current time
The step of degree is modified: if the frozen state of last moment layer soil be do not freeze and this layer of soil of current time it is upper
The frozen state of layer soil is not freeze, and the frozen state of amendment this layer of soil of current time is not freeze.
It further, further include according to every layer of current time in step 3) in order to realize the differentiation to soil freezing region
Soil freezes degree, judges the frozen region of the soil for only some frozen region or has two parts frozen region
Step: if occurring melting region behind inline freezing region, melting region, occurring frozen region again determines that soil has later
Two parts frozen region；Wherein, in frozen region every layer of soil freeze degree be between fully charge and between not freezing or
Fully charge, the degree of freezing for melting every layer of soil in region is not freeze.
Further, abnormal frequency signal occurs in order to prevent, improves and calculates accuracy, further includes the every layer of soil that will be obtained
The step of frequency and temperature of earth are modified: if the frequency and temperature of certain layer of soil are in setting normal frequency range and setting temperature
It spends except range, using the frequency of this layer of soil of last moment and temperature as the frequency and temperature of this layer of soil of current time.
The present invention provides a kind of frost penetration computing device, including memory and processor, and the processor is for executing
Instruction stored in memory is to perform the steps of
1) information that the sensor for detecting frost penetration detects is obtained, the frequency of every layer of soil of current time is obtained；
2) frequency of every layer of soil of current time is made comparisons with setting reference frequency, obtains every layer of soil of current time
Freeze degree；
3) degree is freezed according to every layer of soil of current time, and freezes the actual depth of the corresponding soil of degree, obtained
To the depth of freezing of current time soil, the corresponding frequency of soil is measured by sensor, is obtained compared with setting reference frequency
Soil freezes degree, and carries out the calculating of depth, and this method accuracy is higher, can accurately analyze frost penetration, and
And it is simply easily achieved.
Further, in order to enable the setting reference frequency arrived is more reasonable, and meet the natural law, the step in the device
It is rapid 2) described in setting reference frequency acquisition process it is as follows:
Obtain for detect frost penetration sensor detect information, obtain different moments different layers soil temperature and
Frequency；
Look for and meet frequency corresponding to the soil horizon of following temperature condition: the temperature of soil is in setting solidification point range
It is interior, and the temperature at continuous N number of moment is downward trend and the temperature at this N number of moment is greater than the first set temperature value；
Frequency corresponding to all soil horizons for meeting the temperature condition is corrected, and then obtains the reference frequency.
Further, in order to more be accurately obtained reference frequency, the amendment in the device is all to meet the temperature
Frequency corresponding to the soil horizon of degree condition, and then obtaining the reference frequency includes: to meet the temperature condition for all
Frequency corresponding to soil horizon is fitted, and is fitted to frequency corresponding to critical-temperature, is sought corresponding to all critical-temperatures
Frequency average value, obtain the reference frequency.
Further, for better, simpler differentiation frozen soil state, the process of the step 2) in the device includes:
If the frequency of current time layer soil is more than or equal to the first setting value of reference frequency, the current time layer
The degree of freezing of soil is fully charge；
If the frequency of current time layer soil is less than or equal to the second setting value of reference frequency, the current time layer
The degree of freezing of soil is not freeze；
Current time remainder layer soil freezes between degree do not freeze between fully charge and.
Further, it in order to which more accurately distinguish every layer of soil freezes degree, after the step 2) in the device, also wraps
It includes and the step of degree is modified is freezed to every layer of soil of obtained current time according to the temperature of every layer of soil of different moments:
If the temperature of current time layer soil is greater than the second set temperature value, the degree of freezing for correcting this layer of soil of current time is
Do not freeze；If the temperature of current time layer soil is less than third set temperature value, the jelly of this layer of soil of current time is corrected
Knot degree is fully charge.
Further, in order to which the degree of freezing more accurately to every layer of soil distinguishes, the step 2) in the device
Afterwards, further include the steps that basis is freezed center and is modified to the degree of freezing of every layer of soil of obtained current time:
Freeze degree according to every layer of soil, obtain freezing the corresponding number of plies in center and freezes total depth；
The frozen state for setting the soil horizon in freezing range centered on freezing center is corrected as fully charge；Its
In, the upper limit for setting freezing range, which is subtracted to freeze the corresponding number of plies in center, freezes the setting multiple of total depth, described to set
The lower limit for determining freezing range is to freeze the corresponding number of plies in center to add the setting multiple for freezing total depth.
Further, in order to more accurately correct the frozen state of every layer of soil, after the step 2) in the device, further include
According to last moment layer frozen state and the frozen state of the upper layer of soil of this layer of soil come to every layer obtained of current time
Soil freezes the step of degree is modified: if the frozen state of last moment layer soil be do not freeze and current time should
The frozen state of the upper layer of soil of layer soil is not freeze, and the frozen state of amendment this layer of soil of current time is not freeze.
It further, further include that basis is worked as in the step 3) in the device in order to realize the differentiation to soil freezing region
Every layer of soil of preceding moment freezes degree, judges the frozen region of the soil for only some frozen region or has two parts
The step of frozen region.
Further, further include in the step of abnormal frequency signal occurs in order to prevent, improves and calculates accuracy, the device
The step of frequency of every layer of obtained soil and temperature are modified: if the frequency and temperature of certain layer of soil are in the normal frequency of setting
Except rate range and set temperature range, using the frequency of this layer of soil of last moment and temperature as this layer of soil of current time
Frequency and temperature.
Detailed description of the invention
Fig. 1 is a kind of main flow chart of frost penetration calculation method of the invention；
Fig. 2 is that a kind of reference frequency of frost penetration calculation method of the invention updates flow chart；
Fig. 3 is a kind of degree calculation flow chart of frost penetration calculation method of the invention；
Fig. 4 is a kind of temperature adjustmemt flow chart of frost penetration calculation method of the invention；
Fig. 5 is to freeze degree schematic diagram in the somewhere period of the invention；
Fig. 6 is to freeze degree schematic diagram after abnormality processing of the invention；
Fig. 7 is that a kind of exception of occurring of sensor of the invention detection freezes degree schematic diagram；
Fig. 8 is the result figure before a kind of smoothing processing of frost penetration calculation method of the invention；
Fig. 9 is the result figure after a kind of smoothing processing of frost penetration calculation method of the invention.
Specific embodiment
The present invention will be further described in detail with reference to the accompanying drawing.
The present invention provides a kind of frost penetration calculation method, as shown in Figure 1, including the following steps:
1) obtain for detect frost penetration sensor detect information, obtain every layer of soil of current time frequency and
Temperature.Wherein, if the frequency of certain layer of soil is being set except normal frequency range, using the frequency of this layer of soil of last moment as
The frequency of this layer of soil of current time.Herein, normal frequency range may be configured as [40,60].Certainly, it is also desirable to acquisition
Temperature is modified: if the temperature of certain layer of soil is being set except normal temperature range, by the temperature of this layer of soil of last moment
Temperature as this layer of soil of current time.Herein, setting normal temperature range may be configured as [- 40,50].
The present invention use for detecting the sensor of frost penetration for multiple groups plane capacitance is arranged according to a certain distance
The multi-layer capacity sensor constituted in the same plane is arranged, in order to avoid the interference of multiple sensors from each other, sensor work
Time sharing power supply mode is used when making, synchronization only has a working sensor.When sensor external medium changes, electricity
Capacitance can change, and can cause the variation of frequency on oscillating circuit, just can determine that the state of each layer sensor surrounding soil.
The data that sensor obtains contain the temperature and frequency information of every layer of soil.This part sensor is existing technology, no
Do excessive introduction.
2) frequency of every layer of soil of current time is made comparisons with setting reference frequency, obtains every layer of soil of current time
Freeze degree.
Reference frequency updates:
Ideal reference frequency is the frequency of soil freezing eve, which can not directly acquire, therefore it is desirable that
It obtains close to the soil frequency for freezing not freeze also but.
The acquisition process of setting reference frequency (i.e. reference frequency) is as follows in the step 2), as shown in Fig. 2, passing through the benchmark
Frequency process carries out reference frequency update, but is not limited to the process:
1, the information that the sensor for detecting frost penetration detects is obtained, the temperature of different moments different layers soil is obtained
And frequency.
2, look for and meet frequency corresponding to the soil horizon of following temperature condition: the temperature of soil is in setting solidification point model
In enclosing, and the temperature at continuous N number of moment is downward trend and the temperature at this N number of moment is greater than the first set temperature value.
Soil freezing temperature is less than 0 degree, but due to sensor error, 0 degree may have already appeared and freeze, by field trial number
According to analysis, temperature can guarantee that soil does not freeze at 1 degree, therefore select 1 degree of frequency as reference frequency.Due to each
Data of hour acquisition, the variation of temperature are discrete especially in high level, may span across 1 degree, therefore temperature range is arranged
[0.7,1.5], when temperature reaches this range, acquisition current frequency.It should be noted that existing when temperature is 1 degree by subzero jelly
Knot state is warming up to 1 degree and soil still freezes this case at this time, it is therefore desirable to be screened.For example, screening foundation are as follows: close
8 groups of data temperatures are downward trend, and nearly 8 groups of data temperature mean values are greater than 1.4.
3, frequency corresponding to all soil horizons for meeting temperature condition is corrected, and then obtains reference frequency.
Wherein, frequency corresponding to all soil horizons for meeting temperature condition is corrected, and then obtains reference frequency.It can incite somebody to action
Frequency corresponding to all soil horizons for meeting temperature condition is fitted, and is fitted to frequency corresponding to critical-temperature, is sought
The average value of frequency corresponding to all critical-temperatures, obtains reference frequency.
The amendment is directed to the frequency for meeting temperature condition obtained after above-mentioned screening.Since taken frequency is corresponding
Temperature may any temperature between [0.7,1.5], there is very strong otherness, therefore we are modified frequency, lead to
Over-fitting function is 0.7 degree of corresponding frequency frequency translation, to improve accuracy.Since soil frequency is by precipitation, evaporation, water
The factors such as migration are divided to influence, identical temperature lower frequency there may also be larger difference, therefore to the frequency for meeting condition occurred
Rate is averaged as reference frequency.
Frequency is usually increased with the reduction of temperature, we are 0.7 degree corresponding frequency translation by fitting function
Frequency, fitting function in Fig. 2 are as follows:
F′add=Fadd+(Tt,i-0.7)*0.08
Wherein, FaddFor the frequency for meeting temperature condition, F 'addFor revised frequency, for updating reference frequency, Tt,i
For Current Temperatures, analyze at 1 degree or so by observation data, it is about 0.08 that frequency, which changes corresponding temperature variation ratio, current
Temperature frequency substitute into formula can be by the frequency that is fitted when obtaining 0.7 degree.
Freeze the calculating of degree:
Firstly, further include the division to degree is freezed during above-mentioned step 2), as shown in figure 3, dividing as follows:
If the frequency of current time layer soil is more than or equal to the first setting value of reference frequency, the current time layer
The degree of freezing of soil is fully charge；If the frequency of current time layer soil is set less than or equal to reference frequency second
Value, then the degree of freezing of this layer of soil of current time is not freeze；Current time remainder layer soil freezes degree between complete
Between freezing and not freezing.
For example, second setting value is 0.45, F as shown in figure 3, the first setting value is 0.75dFor benchmark frequency, f is rememberedminFor
Fd+ 0.45, remember fmaxFor Fd+ 0.75, remember Ft,iIndicate the frequency of certain i layers of soil at current time, St,iIndicate frozen state.If working as
Preceding frequency then thinks that reaching freeze frequency soil starts to freeze beyond reference frequency value 0.45, even Ft,i≤fmin, St,i=0；If
Fully charge, even F are thought when beyond freeze frequency 0.75t,i≥fmin, St,i=1；Frequency values between [0.45,075] are corresponding
Freeze degree [0,1], specially St,i=(Ft,i-fmin)/(fmax-fmin)。
The above process is by the way that the frequency of every layer of soil of current time to be compared with reference frequency, to the jelly of every layer of soil
Knot state has carried out preliminary judgement.But due to sensor is unstable, reference frequency data are influenced by the external world, water translocation etc. because
Element influence, the frozen state that previous step may be made to obtain is not accurate enough, can follow the steps below come to obtain it is current when
The frozen state for carving certain layer of soil is modified.
Temperature adjustmemt:
In order to which more accurately distinguish every layer of soil freezes degree, as shown in figure 4, further including after above-mentioned step 2)
The step of degree is modified is freezed to every layer of soil of obtained current time according to the temperature of every layer of soil of different moments: if
The temperature of current time layer soil is greater than the second set temperature value, then the degree of freezing for correcting this layer of soil of current time is not
Freeze；If the temperature of current time layer soil is less than third set temperature value, freezing for this layer of soil of current time is corrected
Degree is fully charge.
For example, being that (horizontal axis is time, the longitudinal axis for the depth of freezing schematic diagram of somewhere period as shown in Figure 5 and Figure 6
For depth), there is the apparent layer (i.e. a plurality of black horizontal line in Fig. 5) not freezed but in subzero in Fig. 5.To abnormal number
According to being handled, the second set temperature value of setting is 3, and third set temperature value is -3, it may be assumed that if the Current Temperatures of certain layer of soil are big
Think that this layer of soil must not freeze at this time in 3；This layer of soil at this time is thought if the Current Temperatures of certain layer of soil are less than -3
It must freeze, freeze schematic diagram as shown in fig. 6, it will be evident that the exception in Fig. 5 is freezed from the Fig. 6 after abnormality processing
Layer has been corrected.
In addition, needing to carry out deep anomalies processing since soil can only be freezed from upper layer to lower layer: if last moment layer
The frozen state of soil is not freeze and the frozen state of the upper layer of soil of this layer of soil of current time is not freeze, and amendment is current
The frozen state of this layer of soil of moment is not freeze, as shown in fig. 7, the white point in circle is a kind of abnormal conditions.
In addition, every annual precipitation, the Weather Elements such as temperature are variant, and the same season soil moisture of the same area is in not the same year
Part is likely to occur larger difference, therefore for better measurement result, further includes what frequency reference was reset in the process of this method
Process: when continuously melting number of days up to 30 days is 720 small, then it is assumed that enter thawing season, reference frequency is zeroed.
3) degree is freezed according to every layer of soil of current time, and freezes the actual depth of the corresponding soil of degree, obtained
To the depth of freezing of current time soil.
Depth of freezing calculating is to freeze degree according to every layer, calculates the bound of the frozen soil depth of freezing this moment.Soil exists
When actually freezing, the cold and hot alternating of weather leads to the alternating of freeze thawing, in fact it could happen that the phenomenon that multilayer is freezed.But for many years to observation station
Observation is divided into part frozen region (second part freezing zone data analysis found that two parts region at most occurs in frozen region
Domain) and lower part frozen region (first part frozen region), and upper partial region thickness is generally relatively thin, 9 observation stations top
The depth of freezing of point frozen region is up to 5cm, and frequency of occurrence is lower.
It mainly includes freezing the contents such as center acquisition, bound calculating, depth smooth that the depth of freezing, which calculates,.
Below with reference to specific example, each step being related to be described in further detail.
For example, the original degree of freezing obtained through the above steps is as shown in the A column of table 1, total number of plies is 35 layers.
Table 1
A) it obtains and freezes center.Degree is freezed according to every layer of soil first, obtains freezing total depth, freezes total depth
28, the number of plies where freezing center, 28/2=14 are then found, then freezing layer corresponding to center is 21 (i.e. central cores)；
Wherein, freeze degree number of plies where freezing center to sum upwards to be 14, summation downwards is also 14.
B) freeze immediate vicinity amendment.After center is freezed in acquisition, closer to freeze center occur melt layer probability it is smaller,
Therefore distance is freezed center and is set as fully charge less than the layer for freezing 40% specific gravity of total depth by we.For example, freezing center
The number of plies at place is 21 (i.e. central cores), then extent of amendment is (10,32) (21-28*40%, 21+28*40%) i.e., C in table 1
It is classified as and corrected freezes degree.
C) judge that still there is two parts frozen region in the frozen region of the soil for only some frozen region, to carry out
The depth of freezing calculates.
It calculates first part frozen region and freezes upper and lower limit, when first part frozen region is when to freeze the upper limit be 0,
Then some frozen region, when not being 0, then there may be two parts frozen regions.Then traversal is current from top to bottom
Moment each layer soil freezes degree, occurs melting region after encountering inline freezing region, melts region and occurs again later
Frozen region then thinks that soil has two parts frozen region, and the frozen region initially encountered is as second part frozen region, meter
It calculates the depth of freezing of the second layer region and updates the depth of freezing of first part frozen region.
Firstly, being superimposed degree is freezed to central core, as shown in D column in table 1.Then it negates, as shown in the E column of table 1.So
Bound is settled accounts by following equation afterwards.It should be noted that freezing the reality that degree is soil corresponding to 1 in the embodiment
Depth is 2.5.
Mid=21, layer=35 are substituted into:
So the upper limit F_Start for obtaining first part frozen region is respectively 12.5, the lower limit F_ that first layer freezes
Stop is 87.5, and the two is all greater than 0, therefore there may be second part frozen regions.
Then, it is traversed from top to bottom by C column, meets layer of the degree greater than 0 for the first time as start, by start downward time
It goes through and encounters 0 for stop, start is 1, stop 6 in figure.Then the upper limit S_Start of second part is (start-1) * 2.5=
0, second part lower limit S_Stop are (stop-1) * 2.5=12.5.
It is 87.5 that the lower limit F_Stop of first part frozen region, which is remained unchanged, the upper limit F_ of first part frozen region
Start should be moved down:
F_Start=12.5+ (1+0.8+0.8+0.2+0.2-1) * 2.5=17.5
Finally, ultimate depth are as follows:
F_Start=17.5, F_Stop=87.5
S_Start=0, S_Stop=12.5
It should be noted that first part frozen region here is the lower region of two layers of frozen region, second part freezes
Tie region is the top area of two layers of frozen region.
Using the above method, the depth of different moments soil can be obtained, as a result as shown in Figure 8.Due to sensor is unstable,
The factors such as water translocation influence, and gained depth of freezing data jump is obvious, do not conform to the actual conditions, therefore carry out to the gained depth of freezing
Smoothing processing.The smoothing processing that the present invention uses is by asking the mean value of preceding 3 groups of depth of freezing as the current depth of freezing, smoothly
Treated, and the depth of freezing is as shown in Figure 9.Above-mentioned method can be realized by the program in processor run memory,
Constitute a kind of device for realizing the above method.
Specific embodiment of the present invention is presented above, but the present invention is not limited to described embodiment.
The technological means in above-described embodiment is converted by the way of being readily apparent that those skilled in the art, is replaced,
Modification, and play the role of with the present invention in relevant art means it is essentially identical, realization goal of the invention it is also essentially identical,
The technical solution formed in this way is to be finely adjusted to be formed to above-described embodiment, and this technical solution still falls within protection of the invention
In range.
Claims (10)
1. a kind of frost penetration calculation method, which comprises the steps of:
1) information that the sensor for detecting frost penetration detects is obtained, the frequency of every layer of soil of current time is obtained；
2) frequency of every layer of soil of current time is made comparisons with setting reference frequency, obtains freezing for every layer of soil of current time
Degree；
3) degree is freezed according to every layer of soil of current time, and freezes the actual depth of the corresponding soil of degree, worked as
The depth of freezing of preceding moment soil.
2. frost penetration calculation method according to claim 1, which is characterized in that set reference frequency described in step 2)
Acquisition process it is as follows:
The information that the sensor for detecting frost penetration detects is obtained, the temperature and frequency of different moments different layers soil are obtained
Rate；
Look for and meet frequency corresponding to the soil horizon of following temperature condition: the temperature of soil is being set within the scope of solidification point,
And the temperature at continuous N number of moment is downward trend and the temperature at this N number of moment is greater than the first set temperature value；
Frequency corresponding to all soil horizons for meeting the temperature condition is corrected, and then obtains the reference frequency.
3. frost penetration calculation method according to claim 2, which is characterized in that the amendment is all to meet the temperature
Frequency corresponding to the soil horizon of condition, and then obtaining the reference frequency includes: by all soil for meeting the temperature condition
Frequency corresponding to earth layer is fitted, and is fitted to frequency corresponding to critical-temperature, is sought corresponding to all critical-temperatures
The average value of frequency obtains the reference frequency.
4. frost penetration calculation method according to claim 1, which is characterized in that the process of step 2) includes:
If the frequency of current time layer soil is more than or equal to the first setting value of reference frequency, this layer of soil of current time
Freeze degree be fully charge；
If the frequency of current time layer soil is less than or equal to the second setting value of reference frequency, this layer of soil of current time
Freeze degree be do not freeze；
Current time remainder layer soil freezes between degree do not freeze between fully charge and.
5. frost penetration calculation method according to claim 4, which is characterized in that further include according to difference after step 2)
The temperature of every layer of soil of moment freezes the step of degree is modified to every layer of soil of obtained current time: if current time
The temperature of certain layer of soil is greater than the second set temperature value, then the degree of freezing for correcting this layer of soil of current time is not freeze；If
The temperature of current time layer soil is less than third set temperature value, then the degree of freezing for correcting this layer of soil of current time has been
Freeze entirely.
6. frost penetration calculation method according to claim 4, which is characterized in that further include that basis is freezed after step 2)
Freeze the step of degree is modified to every layer of soil of obtained current time in center:
Freeze degree according to every layer of soil, obtain freezing the corresponding number of plies in center and freezes total depth；
The frozen state for setting the soil horizon in freezing range centered on freezing center is corrected as fully charge；Wherein, institute
It states to set the upper limit of freezing range and subtract to freeze the corresponding number of plies in center and freezes the setting multiple of total depth, the setting is freezed
The lower limit of range is to freeze the corresponding number of plies in center to add the setting multiple for freezing total depth.
7. frost penetration calculation method according to claim 4, which is characterized in that further include according to upper one after step 2)
The frozen state of the upper layer of soil of moment layer frozen state and this layer of soil is come the jelly to every layer of soil of obtained current time
The step of knot degree is modified: if the frozen state of last moment layer soil is not freeze and this layer of soil of current time
The frozen state of upper layer of soil is not freeze, and the frozen state of amendment this layer of soil of current time is not freeze.
8. frost penetration calculation method according to claim 4, which is characterized in that further include according to current in step 3)
Every layer of soil of moment freezes degree, judges the frozen region of the soil for only some frozen region or has two parts jelly
The step of tie region: if occurring melting region behind inline freezing region, occur frozen region again after melting region and sentence
Determining soil has two parts frozen region；Wherein, the degree of freezing of every layer of soil is not freeze between fully charge and in frozen region
Between knot or fully charge, the degree of freezing for melting every layer of soil in region is not freeze.
9. frost penetration calculation method according to claim 1, which is characterized in that further include the every layer of soil that will be obtained
The step of frequency and temperature are modified: if the frequency and temperature of certain layer of soil are in setting normal frequency range and set temperature model
Except enclosing, using the frequency of this layer of soil of last moment and temperature as the frequency and temperature of this layer of soil of current time.
10. a kind of frost penetration computing device, which is characterized in that including memory and processor, the processor is for executing
Instruction stored in memory is to realize the step as described in any one of claim 1-9 in frost penetration calculation method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811544248.9A CN109579685B (en) | 2018-12-17 | 2018-12-17 | Frozen soil depth calculation method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811544248.9A CN109579685B (en) | 2018-12-17 | 2018-12-17 | Frozen soil depth calculation method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109579685A true CN109579685A (en) | 2019-04-05 |
| CN109579685B CN109579685B (en) | 2020-11-13 |
Family
ID=65929814
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811544248.9A Active CN109579685B (en) | 2018-12-17 | 2018-12-17 | Frozen soil depth calculation method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109579685B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113588920A (en) * | 2021-07-22 | 2021-11-02 | 北京安赛博技术有限公司 | Method, system, equipment and medium for identifying and monitoring soil freezing and thawing process |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH093940A (en) * | 1995-06-20 | 1997-01-07 | Tokyo Gas Co Ltd | Constructing method for underground tank |
| CN102607498A (en) * | 2012-03-02 | 2012-07-25 | 河南省气象科学研究所 | Frozen soil and dry soil layer measuring transducer |
| CN204536237U (en) * | 2015-04-10 | 2015-08-05 | 长安大学 | A kind of Freezing Soils hydro-thermal power comprehensive test system |
| CN107272073A (en) * | 2017-04-26 | 2017-10-20 | 北京交通大学 | A kind of method that application GPR calculates frozen soil relative water content |
-
2018
- 2018-12-17 CN CN201811544248.9A patent/CN109579685B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH093940A (en) * | 1995-06-20 | 1997-01-07 | Tokyo Gas Co Ltd | Constructing method for underground tank |
| CN102607498A (en) * | 2012-03-02 | 2012-07-25 | 河南省气象科学研究所 | Frozen soil and dry soil layer measuring transducer |
| CN204536237U (en) * | 2015-04-10 | 2015-08-05 | 长安大学 | A kind of Freezing Soils hydro-thermal power comprehensive test system |
| CN107272073A (en) * | 2017-04-26 | 2017-10-20 | 北京交通大学 | A kind of method that application GPR calculates frozen soil relative water content |
Non-Patent Citations (1)
| Title |
|---|
| 李鹏等: "电容式冻土测量传感器设计", 《传感器与微系统》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113588920A (en) * | 2021-07-22 | 2021-11-02 | 北京安赛博技术有限公司 | Method, system, equipment and medium for identifying and monitoring soil freezing and thawing process |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109579685B (en) | 2020-11-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Nandintsetseg et al. | Seasonal change of soil moisture in Mongolia: its climatology and modelling | |
| CN106770418B (en) | The device and method that Rock And Soil internal moisture migrates in monitoring frozen-thaw process in real time | |
| Zhou et al. | Monitoring and modeling the influence of snow pack and organic soil on a permafrost active layer, Qinghai–Tibetan Plateau of China | |
| CN110031608A (en) | Ever-frozen ground frozen swell and melt settlement characteristic test apparatus and its test method | |
| Shang et al. | Energy exchange of an alpine grassland on the eastern Qinghai-Tibetan Plateau | |
| CN101836562B (en) | Index method for diagnosing degree of water shortage of crop | |
| CN103743771A (en) | Freeze-thawing circulation detection apparatus for nature saline soil | |
| CN103760623B (en) | Full-automatic water surface evaporation capacity monitoring system | |
| Iwata et al. | The role of permafrost in water exchange of a black spruce forest in Interior Alaska | |
| CN203587589U (en) | Freezing-melting circulation detecting device suitable for natural saline soils | |
| KR20210045253A (en) | Irrigation system and method using water shortage diagnosis of crops | |
| Kutikoff et al. | Heat storage and its effect on the surface energy balance closure under advective conditions | |
| CN113642269B (en) | Precise irrigation method and irrigation system | |
| CN117147007B (en) | Facility agriculture low-temperature cold damage forecasting system | |
| CN106323374A (en) | Method for detecting glacial ablation based on tracer technique | |
| Hou et al. | Assessing crop coefficients for Zea mays in the semi-arid Hailiutu River catchment, northwest China | |
| CN112699621A (en) | Potential evapotranspiration estimation method considering change of underlying surface condition | |
| CN109579685A (en) | A kind of frost penetration calculation method and device | |
| CN110657834A (en) | Method and device for judging icing state of power grid transmission line in micro-terrain area | |
| CN206573506U (en) | The device that Rock And Soil internal moisture is migrated in monitoring frozen-thaw process in real time | |
| CN107748187A (en) | Traditional building wall rate-of-loss of coolant method of testing | |
| CN113030431B (en) | Method and system for measuring soil field water retention | |
| CN116659800A (en) | Automatic determination method for mine ventilation resistance | |
| CN200996055Y (en) | Electrolyte overheat online inspecter for aluminum electrolyzer | |
| Li et al. | Crop model data assimilation with particle filter for yield prediction using leaf area index of different temporal scales |
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 |