CN105547957B - Soil bacterial diversity wetted front suction and rainfall spurt volume computational methods - Google Patents
Soil bacterial diversity wetted front suction and rainfall spurt volume computational methods Download PDFInfo
- Publication number
- CN105547957B CN105547957B CN201510917405.6A CN201510917405A CN105547957B CN 105547957 B CN105547957 B CN 105547957B CN 201510917405 A CN201510917405 A CN 201510917405A CN 105547957 B CN105547957 B CN 105547957B
- Authority
- CN
- China
- Prior art keywords
- mrow
- infiltration
- soil
- msub
- mfrac
- 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
- 239000002689 soil Substances 0.000 title claims abstract description 45
- 230000001580 bacterial effect Effects 0.000 title 1
- 238000000205 computational method Methods 0.000 title 1
- 230000008595 infiltration Effects 0.000 claims abstract description 58
- 238000001764 infiltration Methods 0.000 claims abstract description 58
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 33
- 230000001186 cumulative effect Effects 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 238000002474 experimental method Methods 0.000 claims abstract description 8
- 229920006395 saturated elastomer Polymers 0.000 claims abstract description 7
- 230000008859 change Effects 0.000 claims abstract description 5
- 238000000034 method Methods 0.000 claims description 15
- 241000428199 Mustelinae Species 0.000 claims description 5
- 238000005070 sampling Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- 238000009736 wetting Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000006424 Flood reaction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000012212 insulator Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
- G01N15/0893—Investigating volume, surface area, size or distribution of pores; Porosimetry by measuring weight or volume of sorbed fluid, e.g. B.E.T. method
Landscapes
- Chemical & Material Sciences (AREA)
- Dispersion Chemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
Abstract
一种土壤下渗湿润峰吸力和降雨初损量计算方法,包括以下步骤:通过双环入渗实验,测量并记录供水水量随入渗时间的变化,计算实时累积入渗量;基于入渗时间和实时累积入渗量,拟合必要参数;基于拟合参数和固定入渗水头高度,计算湿润峰吸力;基于拟合参数、土壤初期含水量和土壤孔隙度,计算土壤饱和导水率;基于湿润峰吸力、土壤饱和导水率、降雨强度、土壤初期含水量和土壤孔隙度,计算降雨初损量。
A calculation method for soil infiltration peak suction and rainfall initial loss, comprising the following steps: through a double-loop infiltration experiment, measuring and recording the change of water supply volume with infiltration time, and calculating the real-time cumulative infiltration amount; based on the infiltration time and Real-time cumulative infiltration, fitting necessary parameters; based on fitting parameters and fixed infiltration head height, calculating wet peak suction; based on fitting parameters, initial soil water content and soil porosity, calculating soil saturated hydraulic conductivity; based on wet Peak suction, soil saturated hydraulic conductivity, rainfall intensity, initial soil water content and soil porosity, to calculate initial rainfall loss.
Description
技术领域technical field
本公开涉及水文学降雨产流模型的土壤参数确定,特别涉及一种基于双环入渗的土壤下渗湿润峰吸力和降雨初损量计算方法。The disclosure relates to the determination of soil parameters of a hydrological rainfall runoff model, in particular to a calculation method for soil infiltration peak suction and initial rainfall loss based on double-ring infiltration.
背景技术Background technique
我国是暴雨洪水最频发的国家之一,对洪水灾害进行及时准确的预报预警对于降低灾害损失十分重要。洪水灾害预报预警的关键问题之一是准确模拟场次降雨产流过程。Green-Ampt模型是物理机制清晰的产流模型,然而其湿润峰固定吸力参数难以确定,严重限制了该模型的应用。传统的双环入渗实验仅仅能够测试土壤饱和导水率,也不能确定Green-Ampt模型中的湿润峰固定吸力参数,难以满足Green-Ampt模型应用于流域水文学对土壤参数的需求。因此,亟需一种可靠实用的方法,能够快速计算Green-Ampt模型中的湿润峰固定吸力参数和降雨初损量。my country is one of the countries where rainstorms and floods occur most frequently. Timely and accurate forecasting and early warning of flood disasters is very important for reducing disaster losses. One of the key issues in flood disaster forecasting and early warning is to accurately simulate the process of rainfall and runoff generation. The Green-Ampt model is a runoff model with a clear physical mechanism, but it is difficult to determine the fixed suction parameters of the wet peak, which severely limits the application of the model. The traditional double-ring infiltration experiment can only test the soil saturated hydraulic conductivity, and cannot determine the fixed suction parameters of the wet peak in the Green-Ampt model, which is difficult to meet the soil parameter requirements of the Green-Ampt model applied to watershed hydrology. Therefore, there is an urgent need for a reliable and practical method that can quickly calculate the fixed suction parameters of the wet peak and the initial rainfall loss in the Green-Ampt model.
发明内容Contents of the invention
为了弥补现有技术的缺失,本公开提供一种土壤下渗湿润峰吸力和降雨初损量计算方法。In order to make up for the deficiencies in the prior art, the present disclosure provides a calculation method for soil infiltration peak suction and initial rainfall loss.
本公开采用以下解决方案:The present disclosure employs the following solutions:
一种土壤下渗湿润峰吸力和降雨初损量计算方法,包括以下步骤:A calculation method for soil infiltration peak suction and rainfall initial loss, comprising the following steps:
步骤1:通过双环入渗实验,测量并记录供水水量Q随入渗时间t的变化,根据公式(1)计算实时累积入渗量h:Step 1: Through the double-loop infiltration experiment, measure and record the change of the water supply Q with the infiltration time t, and calculate the real-time cumulative infiltration h according to the formula (1):
其中,S表示采用的双环入渗仪的内环面积;Among them, S represents the area of the inner ring of the double-ring infiltration instrument used;
步骤2:基于所述入渗时间t和实时累积入渗量h,根据公式(2)拟合参数A和参数B:Step 2: Based on the infiltration time t and the real-time cumulative infiltration amount h, fitting parameter A and parameter B according to formula (2):
步骤3:根据公式(3)计算湿润峰吸力hf:Step 3: Calculate the wet peak suction h f according to formula (3):
hf=B-h0 (3)h f =Bh 0 (3)
其中,h0表示固定入渗水头高度;Among them, h 0 represents the fixed infiltration head height;
步骤4:根据公式(4)计算土壤饱和导水率K:Step 4: Calculate soil saturated hydraulic conductivity K according to formula (4):
其中,θ0表示土壤初期含水量,θs表示土壤孔隙度;Among them, θ 0 represents the initial water content of the soil, and θ s represents the soil porosity;
步骤5:根据公式(5)计算降雨初损量Ia:Step 5: Calculate the initial rainfall loss I a according to the formula (5):
其中,P表示降雨强度。Among them, P represents the rainfall intensity.
优选地,通过环刀取样测量所述土壤孔隙度θs。Preferably, said soil porosity θ s is measured by ring knife sampling.
优选地,通过土壤含水量测试测定所述土壤初期含水量θ0。Preferably, the initial soil water content θ 0 is determined by a soil water content test.
优选地,通过最小二乘法拟合所述参数A和参数B。Preferably, said parameter A and parameter B are fitted by a least square method.
优选地,所述供水水量Q是供水马氏瓶的水量。Preferably, the water supply quantity Q is the water quantity of the water supply Martens flask.
本公开的有益效果是基于双环入渗计算土壤下渗湿润峰吸力和降雨初损量,能够快速准确地确定Green-Ampt模型所涉及的土壤参数,进而进行洪水灾害预报预警。The beneficial effect of the present disclosure is that the calculation of soil infiltration peak suction and initial rainfall loss based on the double-ring infiltration can quickly and accurately determine the soil parameters involved in the Green-Ampt model, and then perform flood disaster forecasting and early warning.
附图说明Description of drawings
通过结合附图对本公开示例性实施例进行更详细的描述,本公开的上述以及其它目的、特征和优势将变得更加明显,其中,在本公开示例性实施例方式中,相同的参考标号通常代表相同部件。The above and other objects, features and advantages of the present disclosure will become more apparent by describing the exemplary embodiments of the present disclosure in more detail with reference to the accompanying drawings, wherein, in the exemplary embodiments of the present disclosure, the same reference numerals are generally represent the same part.
图1示出了根据示例性实施例的土壤下渗湿润峰吸力和降雨初损量计算方法的流程图;Fig. 1 shows the flow chart of the method for calculating the peak suction of soil infiltration and initial rainfall loss according to an exemplary embodiment;
图2示出了根据示例性实施例的入渗量与入渗时间的曲线图。FIG. 2 shows a graph of infiltration amount versus infiltration time, according to an exemplary embodiment.
具体实施方式Detailed ways
下面将参照附图更详细地描述本公开的优选实施例。虽然附图中显示了本公开的优选实施例,然而应该理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本公开更加透彻和完整,并且能够将本公开的范围完整地传达给本领域的技术人员。Preferred embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. Although preferred embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure can be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
图1示出了根据示例性实施例的土壤下渗湿润峰吸力和降雨初损量计算方法的流程图,该方法包括以下步骤:Fig. 1 shows the flow chart of the calculation method of soil infiltration peak suction and rainfall initial loss according to an exemplary embodiment, the method includes the following steps:
步骤1,通过双环入渗实验,测量并记录供水水量Q随入渗时间t的变化,根据以下公式(1)计算实时累积入渗量h:Step 1. Through the double-loop infiltration experiment, measure and record the change of the water supply Q with the infiltration time t, and calculate the real-time cumulative infiltration h according to the following formula (1):
其中,S表示采用的双环入渗仪的内环面积。Among them, S represents the area of the inner ring of the double-ring infiltration instrument used.
在实际操作中,优选采用双环入渗仪进行双环入渗实验。双环入渗仪是本领域中常用的测量水渗入土壤的渗透速度的设备。标准的双环入渗仪包括直径不同的外环和内环。垂直渗透水流向边缘时,外环可以起到隔离的作用。由于内环中的水是垂直流动的,因而测量仅限于在内环中进行。In actual operation, it is preferable to use a double-ring infiltration instrument to carry out double-ring infiltration experiments. A double-loop infiltration meter is a device commonly used in the art to measure the rate of infiltration of water into soil. A standard dual-ring infiltrator includes outer and inner rings of different diameters. The outer ring acts as an insulator when the vertically percolating water flows towards the edge. Since the water in the inner ring flows vertically, measurements are limited to the inner ring.
作为一种优选方案,采用马氏瓶作为供水设备,因此供水水量Q是指供水马氏瓶的水量,当采用其他供水设备时,供水水量Q是指所采用的供水设备的水量。As a preferred solution, the Martens bottle is used as the water supply equipment, so the water supply quantity Q refers to the water quantity of the water supply Martens flask, and when other water supply equipment is used, the water supply quantity Q refers to the water quantity of the water supply equipment adopted.
步骤2,基于步骤1中的入渗时间t和实时累积入渗量h,根据以下公式(2),拟合参数A和参数B:Step 2, based on the infiltration time t and the real-time cumulative infiltration h in step 1, according to the following formula (2), fitting parameter A and parameter B:
在实际实施过程中,可以通过最小二乘法编制计算机程序或者利用通用数据处理软件(例如TableCurve 2D),拟合参数A和参数B。In the actual implementation process, the computer program can be compiled by the least square method or the general data processing software (such as TableCurve 2D) can be used to fit the parameters A and B.
步骤3,根据以下公式(3)计算湿润峰吸力hf:Step 3, calculate the wet peak suction hf according to the following formula (3):
hf=B-h0 (3)h f =Bh 0 (3)
其中,h0表示固定入渗水头高度。Among them, h 0 represents the fixed infiltration head height.
步骤4,根据以下公式(4)计算土壤饱和导水率K:Step 4, calculate the soil saturated hydraulic conductivity K according to the following formula (4):
其中,θ0表示土壤初期含水量,可通过土壤含水量测试来测定,θs表示土壤孔隙度,可通过环刀取样来测定。Among them, θ 0 represents the initial water content of soil, which can be determined by soil moisture test, and θ s represents the soil porosity, which can be measured by ring knife sampling.
步骤5,根据以下公式(5)计算降雨初损量Ia:Step 5, calculate the initial rainfall loss I a according to the following formula (5):
其中,P表示降雨强度。根据实际应用,P可以通过雨量计测量,也可以是设计的降雨强度。Among them, P represents the rainfall intensity. Depending on the practical application, P can be measured by a rain gauge, or it can be the designed rainfall intensity.
应用示例Application example
本示例提供一种基于野外双环入渗实验的土壤下渗湿润峰吸力和降雨初损量计算方法,包括以下步骤:This example provides a calculation method for soil infiltration peak suction and initial rainfall loss based on field double-ring infiltration experiments, including the following steps:
步骤1:进行双环入渗实验,以1秒为采样周期,利用挂式电子秤连接电脑串口测量并记录供水马氏瓶的水量随入渗时间t的变化,根据公式(1)计算实时累积入渗量h,得到的实时累积入渗量h和入渗时间t的曲线图如图2所示。Step 1: Carry out a double-loop infiltration experiment, with a sampling period of 1 second, use a hanging electronic scale to connect to the serial port of the computer to measure and record the change of the water volume of the water supply Martens flask with the infiltration time t, and calculate the real-time cumulative infiltration according to formula (1). Infiltration amount h, the obtained real-time cumulative infiltration amount h and infiltration time t curves are shown in Figure 2.
步骤2:基于步骤1中的实时累积入渗量h和入渗时间t的数据,根据公式(2),通过最小二乘法拟合公式(2)中的参数A和参数B,得到A=32.45s/cm,B=25.91cm,R2=0.9997,其中R2表示拟合确定系数。Step 2: Based on the data of real-time cumulative infiltration h and infiltration time t in step 1, according to formula (2), fit parameter A and parameter B in formula (2) by least square method, and obtain A=32.45 s/cm, B=25.91 cm, R 2 =0.9997, where R 2 represents the fitting determination coefficient.
步骤3:测量得到固定入渗水头高度h0=10cm,根据公式(3)计算湿润峰吸力hf=15.91cm。Step 3: Measure the fixed infiltration head height h 0 =10 cm, and calculate the wetting peak suction h f =15.91 cm according to the formula (3).
步骤4:通过土壤含水量测试测量得到土壤初期含水量θ0=0,通过环刀取样测得的土壤孔隙度θs=36%,根据公式(4)计算土壤饱和导水率K=399.22mm/hr。Step 4: Obtain the initial soil water content θ 0 =0 through the soil moisture test measurement, the soil porosity θ s =36% measured by the ring knife sampling, and calculate the soil saturated hydraulic conductivity K=399.22mm according to the formula (4) /hr.
步骤5:根据公式(5)计算降雨初损量:Step 5: Calculate the initial rainfall loss according to formula (5):
上述技术方案只是本发明的一种实施例,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施例所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。The above-mentioned technical solution is only an embodiment of the present invention. For those skilled in the art, on the basis of the application methods and principles disclosed in the present invention, it is easy to make various types of improvements or deformations, and is not limited to The methods described in the above specific embodiments of the present invention, therefore, the above-described methods are only preferred and not limiting.
Claims (5)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510917405.6A CN105547957B (en) | 2015-12-10 | 2015-12-10 | Soil bacterial diversity wetted front suction and rainfall spurt volume computational methods |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510917405.6A CN105547957B (en) | 2015-12-10 | 2015-12-10 | Soil bacterial diversity wetted front suction and rainfall spurt volume computational methods |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105547957A CN105547957A (en) | 2016-05-04 |
CN105547957B true CN105547957B (en) | 2018-05-04 |
Family
ID=55827299
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510917405.6A Expired - Fee Related CN105547957B (en) | 2015-12-10 | 2015-12-10 | Soil bacterial diversity wetted front suction and rainfall spurt volume computational methods |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105547957B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108426803A (en) * | 2018-03-02 | 2018-08-21 | 江西理工大学 | A kind of assay method of ion type rareearth One-dimensional Vertical Infiltration rule |
CN109444016B (en) * | 2018-11-05 | 2024-05-31 | 北京市水科学技术研究院 | Simple soil infiltration rate measuring device and measuring method |
CN109887241A (en) * | 2019-04-08 | 2019-06-14 | 河北省水利水电勘测设计研究院 | A kind of mountain flood weather warning calculation method and system |
CN110029647B (en) * | 2019-05-10 | 2020-10-09 | 长沙理工大学 | A method for calculating the ultimate bearing capacity of homogeneous foundation under pre-peak rainwater infiltration |
CN110400014A (en) * | 2019-07-23 | 2019-11-01 | 华东师范大学 | A Numerical Simulation Method of Multi-source Floods in Coastal Cities Based on GIS Raster Operation |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101738356A (en) * | 2010-01-15 | 2010-06-16 | 鲁东大学 | Method for measuring and calculating soil wetting front suction |
CN101762445A (en) * | 2010-01-15 | 2010-06-30 | 鲁东大学 | Soil saturation hydraulic conductivity measuring and calculating method based on infiltration time characteristic parameters |
CN103063820A (en) * | 2013-01-04 | 2013-04-24 | 张振华 | Method and special device for measuring soil hydrodynamic parameters in situ in field |
CN203365275U (en) * | 2013-04-12 | 2013-12-25 | 新疆维吾尔自治区水利水电科技信息中心 | Point source infiltration wetting front testing device based on Markov bottle |
-
2015
- 2015-12-10 CN CN201510917405.6A patent/CN105547957B/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101738356A (en) * | 2010-01-15 | 2010-06-16 | 鲁东大学 | Method for measuring and calculating soil wetting front suction |
CN101762445A (en) * | 2010-01-15 | 2010-06-30 | 鲁东大学 | Soil saturation hydraulic conductivity measuring and calculating method based on infiltration time characteristic parameters |
CN103063820A (en) * | 2013-01-04 | 2013-04-24 | 张振华 | Method and special device for measuring soil hydrodynamic parameters in situ in field |
CN203365275U (en) * | 2013-04-12 | 2013-12-25 | 新疆维吾尔自治区水利水电科技信息中心 | Point source infiltration wetting front testing device based on Markov bottle |
Non-Patent Citations (3)
Title |
---|
A combined rainfall infiltration model based on Green-Ampt and SCS-curve number;李军 等;《Hydrological Processes》;20141212;第29卷;第2628-2634页 * |
微咸水入渗条件下Philip模型与Green-Ampt模型参数的对比分析;史晓楠;《土壤学报》;20070331;第44卷(第2期);第360-363页 * |
水循环综合模拟系统的降雨产流模型研究;刘昌明 等;《河海大学学报(自然科学版)》;20150930;第43卷(第5期);第378-382页第2-4节 * |
Also Published As
Publication number | Publication date |
---|---|
CN105547957A (en) | 2016-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105547957B (en) | Soil bacterial diversity wetted front suction and rainfall spurt volume computational methods | |
CN104568694B (en) | Method for testing gas-water relative permeability of dense core | |
CN203929586U (en) | A kind of soil infiltration automatic water supplement and measurement mechanism | |
CN101762445B (en) | Soil saturation hydraulic conductivity measuring and calculating method based on infiltration time characteristic parameters | |
CN107340219B (en) | Oil reservoir dynamic capillary effect evaluation experiment data processing method | |
CN103091226A (en) | Device and method for detecting porosity of saturated soil | |
CN206430670U (en) | A kind of river bottom mud measurer for thickness | |
CN103885078B (en) | Sampled measurements topsoil precipitation rate of radon obtains the method and device of soil potentiality radon consistence | |
CN104569350B (en) | Test method for original water saturation of irregular full-diameter rock samples from sealed coring | |
CN105841847B (en) | A kind of method for estimating Surface latent heat fluxes | |
CN101526515B (en) | Method and device for measuring disseminated parameter of solute migration | |
CN106769771A (en) | A kind of measuring method of the unsaturated soil infiltration coefficient based on low-field nuclear magnetic resonance technology | |
CN101738356B (en) | Method for measuring and calculating soil wetting front suction | |
CN106525753A (en) | Convenient and simple remote-sensing soil moisture monitoring method | |
CN102323178B (en) | Method and device for measuring physical property indexes of soil body | |
CN204649586U (en) | A kind of Oil in Super-low Permeability rock core leads to nitrogen water two-phase relative permeability determinator | |
CN106709168B (en) | River-based flow prediction method | |
CN106202002A (en) | A kind of for detecting the method whether series of hydrological parameter makes a variation | |
CN204855299U (en) | Density detector | |
CN206479579U (en) | A kind of axle sample resistivity experimental provision of temperature control three | |
CN109960861B (en) | Earth surface and underground water source segmentation method based on full-accumulation flow-producing mode | |
CN208109744U (en) | A kind of Soil salinity detection recorder | |
CN105095669B (en) | The evaluation method of slope ground surface roughness and infiltration parameter under the conditions of current scour | |
CN116932996A (en) | Estimation method of soil organic matter content based on improved Hapke model | |
CN205787135U (en) | GPS static observation quality of data real-time detecting system |
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 |
Granted publication date: 20180504 Termination date: 20201210 |
|
CF01 | Termination of patent right due to non-payment of annual fee |