CN105842426A - 一种利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法 - Google Patents
一种利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法 Download PDFInfo
- Publication number
- CN105842426A CN105842426A CN201610391977.XA CN201610391977A CN105842426A CN 105842426 A CN105842426 A CN 105842426A CN 201610391977 A CN201610391977 A CN 201610391977A CN 105842426 A CN105842426 A CN 105842426A
- Authority
- CN
- China
- Prior art keywords
- soil
- stone
- stone column
- erosion
- width
- 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
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Food Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Medicinal Chemistry (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (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)
Abstract
本发明公开了利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法,它包括它包括:选择石柱,在所选石柱与表土层接触处做标记线a;测出石柱过标记线a水平面的宽度,以该宽度为边长在相邻土层上标记一个面积为S的正方形区域;经过一场雨的淋洗或隔年后,在所测石柱与表土层接触处做标记线a1,测量出所测石柱过标记线a1水平面的宽度,以该宽度为边长在相邻土层上标记一个面积为S’的正方形区域,测量标记线a至a1石柱的变化高度H,以该高度作为测量侵蚀或堆积土壤体积的高度;计算所测区域侵蚀或堆积的土壤体积V;测定各区域土壤的密度ρ;计算土壤侵蚀或堆积量A;解决了现有技术精度低、误差大、模型参数难以获得等技术问题。
Description
技术领域:
本发明属于喀斯特地区土壤侵蚀或堆积量计量技术,尤其涉及一种利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法。
背景技术:
现有技术常用径流小区、测针法、放射性核素示踪法、模型估算和3S结合法等来对土壤侵蚀或堆积量进行测量。径流小区:因监测手段的粗略性和薄弱性制约了地面监测水平的发展,无法提供较长时期土壤侵蚀数据,难以应用于土壤侵蚀大范围空间变化;测针法的自动化程度低,耕作活动会引起侵蚀针对周围的敏感区域不同程度的扰动,受到环境和人为等诸多弊端的限制;放射性核素示踪法:某些模型参数值难以获得,限制模型应用;某些技术应用起步较晚,有关定量模型少;某些方法仍处在探索阶段,精度有待检验;模型估算:精度易受研究区条件限制;某些模型参数不易确定,推广应用有限制性;侵蚀过程模拟有区域性,模型使用需注意适应范围和可塑性;3S结合法:信息共享困难;土壤侵蚀测报技术未标准化,测报结果带主观因素,误差较大。
发明内容:
本发明要解决的技术问题:提供一种利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法,以解决现有技术对喀斯特地区土壤侵蚀或堆积量存在的精度低、误差大、模型参数难以获得等技术问题。本发明技术方案:
一种利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法,它包括未来土壤侵蚀或堆积量的计算,它包括:
步骤1、选择地层倾角小于10°,高度居所测区域石林平均高度的石柱,所选石柱周围有边长大于石柱与土层接触面宽度的正方形土层区,在所选石柱与表土层接触处做标记线a;
步骤2、测出所选石柱过标记线a水平面的宽度,以该宽度为边长在相邻土层上标记一个面积为S的正方形区域;
步骤3、经过一场雨的淋洗或隔年后,在所测石柱与表土层接触处做标记线a1,测量出所测石柱过标记线a1水平面的宽度,以该宽度为边长在相邻土层上标记一个面积为S’的正方形区域,测量标记线a至a1石柱的变化高度H,以该高度作为测量侵蚀或堆积土壤体积的高度;
步骤4、计算所测区域侵蚀或堆积的土壤体积V;
步骤5、测定各区域土壤的密度ρ;
步骤6、计算土壤侵蚀或堆积量A。
步骤4计算所测区域侵蚀或堆积的土壤体积V的计算公式为:
步骤2所述测出所选石柱过标记线a水平面的宽度,采用经纬仪或水准仪进行测量。
步骤5所述测定各区域土壤的密度ρ的方法为采用比重瓶法。用步骤1的方法对石柱做标记,一年(△t)后用步骤2的方法测石柱变化高度为△H,用步骤5的方法测出土壤密度ρ,
则可求出在面积S1内,土壤的侵蚀模数A1,计算公式为:A1=△H/△t×S1×ρ。
本发明的有益效果:
本发明利用石林来计算喀斯特地区土壤侵蚀或堆积量,可以应用土壤侵蚀量或土壤堆积量来反映喀斯特地区石林未来演变情况,可开发石林旅游区,从而给人们带来更多的经济效益,或可为测量大量水土流失破坏生态环境提供科学依据。该方法简单易行,提供较长时期土壤侵蚀数据,能避免像测针法因耕作活动会引起侵蚀针对周围的敏感区域不同程度的扰动,受到环境和人为等诸多弊端的限制,更不需要耗费大量人力、物力和时间,测量精度比较高等;解决了现有技术对喀斯特地区土壤侵蚀或堆积量存在的精度低、误差大、模型参数难以获得等技术问题。
说明书附图:
图1为本发明测试样品标记示意图。
具体实施方式:
对未来土壤侵蚀或堆积量的测定:
(1)对区域内土壤侵蚀或堆积量测定:
步骤一,选择地层平缓的区域,地层倾角小于10°,若地层倾角大于10°,石林不发育,选取高度居该片石林平均高度的石柱,同时要满足所选石柱周围有边长大于石柱与土层接触面宽度的正方形土层区,在该石柱与表土层接触处做标记线a;
步骤二,用经纬仪或水准仪量出石柱过标记线a水平面的宽度,以该宽度为边长在相邻土层上标记一个面积为S的正方形区域;
步骤三,经过一场雨的淋洗或隔年后,在该石柱与表土层接触处做标记线a1,用步骤二的方法测量出同一石柱过标记线a1水平面的宽度,以该宽度为边长在相邻土层上标记一个面积为S’的正方形区域,同时用步骤二的方法测量标记线a-a1石柱的变化高度H,以该高度作为测量侵蚀或堆积土壤体积的高度;
步骤四,计算出各区域侵蚀或堆积的土壤体积V,由公式可得出结果;
步骤五,测定各区域土壤的密度ρ,通常采用比重瓶法,将已知质量的土样放入水中或其他液体,排尽空气,求出由土壤置换出的水的体积;再以烘干土样的质量除以土壤固相体积即得土壤密度ρ。公式:ρ=M/V1
利用环刀等工具测得的土壤的体积﹣由土壤置换出的水的体积=土壤固相体积,求得的土壤固相体积V1,
步骤六,计算土壤侵蚀或堆积量A。即是土壤体积V与土壤密度ρ的乘积,公式:A=Vρ
步骤七,用SPSS软件对该20根石柱进行聚类分析(筛选出相似性高的石柱),为分析各区域石林形态大小的相似性和差异性,反映出客观事物的内在联系,对石柱高度进行聚类分析,从而减小随机选取计算的误差,得出生长相似性更高的结果。
步骤八,选取同样石柱20根进行测量,因石林具有石柱高度、宽度、土壤密度及石柱形成形态不同等特征差异存在,需要对每根石柱测定结果求平均值,可用SPSS软件对该20根石柱进行方差分析,目的主要是检验被测量区石柱高度变化大小的差异性是由于抽样测量造成的随机误差,还是由于被测量区因地貌部位和构造岩性等差异所造成的条件误差,从而减小测量结果的误差,提高测量结果的精确度。
(2)对土壤侵蚀模数的测定:用步骤一的方法对石柱做标记,一年(△t)后用步骤二的方法测石柱变化高度为△H,用步骤五的方法测出土壤密度ρ,则可求出在面积S1内,土壤的侵蚀模数为A1。
公式:A1=△H/△t×S1×ρ。
对过去土壤侵蚀量的测定:
第一,对石柱做等距离标记如h0、h1、h2、…、hn-1、hn;
第二,用手提电动切割机或角向磨光机分别对各标记点处取样;
第三,对所取样品用核素36Cl测定土壤对石柱的侵蚀年限,其原理是假定到达该地表石柱的宇宙射线通量是个常数,即不随时间变化,地表石柱中宇宙核素累积的浓度与石柱暴露的时间相关。知道了样品中宇宙核素的浓度N和宇宙核素的产生率P,就可以求出地表石柱的暴露时间t。浓度N可通过对样品的分析测得,再测出宇宙核素的产生率P,λ为放射性宇宙核素的衰变系数,为半衰期的倒数,1/Ma根据t=-ln(1-λN/P)/λ可得出结果。若h0处所测时间为t0、h1处所测时间为t1、h2处所测时间为t2、…、hn-1处所测时间为tn-1、hn处所测时间为tn,则有h1﹣h0=△h1、h2﹣h1=△h2、…、hn﹣hn-1=△hn,且△h1=△h2=…=△hn侵蚀年限分别为△t1、△t2、…、△tn;
第四,测定土壤密度ρ;
第五,计算在面积S2内,S2为所选取的不同测量区内,不同年限石柱被侵蚀相等高度时的土壤侵蚀量A2,公式:A2=△hn/△tn×S2×ρ。
Claims (5)
1.一种利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法,它包括未来土壤侵蚀或堆积量的计算,其特征在于它包括:
步骤1、选择地层倾角小于10°,高度居所测区域石林平均高度的石柱,所选石柱周围有边长大于石柱与土层接触面宽度的正方形土层区,在所选石柱与表土层接触处做标记线a;
步骤2、测出所选石柱过标记线a水平面的宽度,以该宽度为边长在相邻土层上标记一个面积为S的正方形区域;
步骤3、经过一场雨的淋洗或隔年后,在所测石柱与表土层接触处做标记线a1,测量出所测石柱过标记线a1水平面的宽度,以该宽度为边长在相邻土层上标记一个面积为S’的正方形区域,测量标记线a至a1石柱的变化高度H,以该高度作为测量侵蚀或堆积土壤体积的高度;
步骤4、计算所测区域侵蚀或堆积的土壤体积V;
步骤5、测定各区域土壤的密度ρ;
步骤6、计算土壤侵蚀或堆积量A。
2.根据权利要求1所述的利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法,其特征在于:步骤4计算所测区域侵蚀或堆积的土壤体积V的计算公式为:
3.根据权利要求1所述的利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法,其特征在于:步骤2所述测出所选石柱过标记线a水平面的宽度,采用经纬仪或水准仪进行测量。
4.根据权利要求1所述的利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法,其特征在于:步骤5所述测定各区域土壤的密度ρ的方法为采用比重瓶法。
5.根据权利要求1所述的利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法,其特征在于:它还包括对土壤侵蚀模数的测定,其测定方法为:用步骤1的方法对石柱做标记,一年(△t)后用步骤2的方法测石柱变化高度为△H,用步骤5的方法测出土壤密度ρ,则可求出在面积S1内,土壤的侵蚀模数A1,计算公式为:A1=△H/△t×S1×ρ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610391977.XA CN105842426B (zh) | 2016-06-06 | 2016-06-06 | 一种利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610391977.XA CN105842426B (zh) | 2016-06-06 | 2016-06-06 | 一种利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105842426A true CN105842426A (zh) | 2016-08-10 |
CN105842426B CN105842426B (zh) | 2018-03-23 |
Family
ID=56576625
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610391977.XA Active CN105842426B (zh) | 2016-06-06 | 2016-06-06 | 一种利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105842426B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108491653A (zh) * | 2018-03-29 | 2018-09-04 | 中国科学院地球化学研究所 | 一种喀斯特地区降雨侵蚀力计算方法 |
CN110986767A (zh) * | 2019-12-10 | 2020-04-10 | 贵州省山地资源研究所 | 一种喀斯特洞穴石笋沉积监测方法 |
CN115166211A (zh) * | 2022-09-08 | 2022-10-11 | 北京江河惠远科技有限公司 | 输电线路全地形土壤侵蚀模数的自动计算方法及系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1664582A (zh) * | 2005-03-18 | 2005-09-07 | 广州地理研究所 | 一种坡面土壤侵蚀量的测计方法 |
CN201740777U (zh) * | 2010-04-06 | 2011-02-09 | 长江水利委员会长江科学院 | 坡面土壤侵蚀快速测定系统 |
CN102353764A (zh) * | 2011-09-16 | 2012-02-15 | 中国科学院水利部成都山地灾害与环境研究所 | 一种土壤降雨可蚀性快速测定方法 |
CN202256323U (zh) * | 2011-10-20 | 2012-05-30 | 吉林省水土保持科学研究院 | 土壤侵蚀模数快速测定小区监测装置 |
CN202305524U (zh) * | 2011-10-20 | 2012-07-04 | 辽宁省水土保持研究所 | 一种坡面水土流失监测装置 |
-
2016
- 2016-06-06 CN CN201610391977.XA patent/CN105842426B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1664582A (zh) * | 2005-03-18 | 2005-09-07 | 广州地理研究所 | 一种坡面土壤侵蚀量的测计方法 |
CN201740777U (zh) * | 2010-04-06 | 2011-02-09 | 长江水利委员会长江科学院 | 坡面土壤侵蚀快速测定系统 |
CN102353764A (zh) * | 2011-09-16 | 2012-02-15 | 中国科学院水利部成都山地灾害与环境研究所 | 一种土壤降雨可蚀性快速测定方法 |
CN202256323U (zh) * | 2011-10-20 | 2012-05-30 | 吉林省水土保持科学研究院 | 土壤侵蚀模数快速测定小区监测装置 |
CN202305524U (zh) * | 2011-10-20 | 2012-07-04 | 辽宁省水土保持研究所 | 一种坡面水土流失监测装置 |
Non-Patent Citations (2)
Title |
---|
张信宝等: "监测土壤流失的标线法及其测定犁耕运移土壤的尝试", 《中国水土保持SWCC》 * |
李月等: "径流小区法监测水土流失的百年历程(1915-2014年)", 《中国水土保持SWCC》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108491653A (zh) * | 2018-03-29 | 2018-09-04 | 中国科学院地球化学研究所 | 一种喀斯特地区降雨侵蚀力计算方法 |
CN108491653B (zh) * | 2018-03-29 | 2020-07-07 | 中国科学院地球化学研究所 | 一种喀斯特地区降雨侵蚀力计算方法 |
CN110986767A (zh) * | 2019-12-10 | 2020-04-10 | 贵州省山地资源研究所 | 一种喀斯特洞穴石笋沉积监测方法 |
CN110986767B (zh) * | 2019-12-10 | 2021-07-06 | 贵州省山地资源研究所 | 一种喀斯特洞穴石笋沉积监测方法 |
CN115166211A (zh) * | 2022-09-08 | 2022-10-11 | 北京江河惠远科技有限公司 | 输电线路全地形土壤侵蚀模数的自动计算方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN105842426B (zh) | 2018-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Schattan et al. | Continuous monitoring of snowpack dynamics in alpine terrain by aboveground neutron sensing | |
Smith et al. | The Murrumbidgee soil moisture monitoring network data set | |
San'kov et al. | Geometry and rate of faulting in the North Baikal Rift, Siberia | |
CN104035138B (zh) | 一种全球及局部海洋扰动重力的精确快速计算方法 | |
Navas et al. | Assessing soil redistribution in a complex karst catchment using fallout 137Cs and GIS | |
US10401489B2 (en) | Method for drilling and characterizing subsurface porosity | |
Yao et al. | Reevaluation of the Late Pleistocene slip rate of the Haiyuan fault near Songshan, Gansu Province, China | |
CN106846475B (zh) | 一种河流历史形态数字复原方法 | |
Hector et al. | Gravity effect of water storage changes in a weathered hard-rock aquifer in West Africa: results from joint absolute gravity, hydrological monitoring and geophysical prospection | |
Poggi et al. | Seismic site-response characterization of high-velocity sites using advanced geophysical techniques: application to the NAGRA-Net | |
Karakostas et al. | The aftershock sequence of the 2008 Achaia, Greece, earthquake: joint analysis of seismicity relocation and persistent scatterers interferometry | |
CN105842426A (zh) | 一种利用石林来计算喀斯特地区土壤侵蚀或堆积量的方法 | |
Chen et al. | Revisiting the deformed high shoreline of Lake Bonneville | |
Peterson et al. | How spatial patterns of soil moisture dynamics can explain field‐scale soil moisture variability: Observations from a sodic landscape | |
Takekawa et al. | Measuring sediment accretion in early tidal marsh restoration | |
CN102419457B (zh) | 一种利用单个竖直钻孔的电视图像确定深部岩体结构面产状的方法 | |
Nath et al. | Topographic gradient based site characterization in India complemented by strong ground-motion spectral attributes | |
Kämpf et al. | Nonvolcanic carbon dioxide emission at continental rifts: the bublak Mofette Area, Western Eger Rift, Czech Republic | |
CN101266153B (zh) | 测绘工程类陀螺全站仪精度评定方法 | |
CN107270871A (zh) | 一种大型湖泊洲滩湿地地表高程动态变化测量方法 | |
Zhang et al. | Imaging hydrological dynamics in karst unsaturated zones by time-lapse electrical resistivity tomography | |
Fores et al. | Variability of the water stock dynamics in karst: insights from surface-to-tunnel geophysics | |
Xu et al. | Vertical slip rates of normal faults constrained by both fault walls: A case study of the hetao fault system in Northern China | |
McKenzie et al. | Differential coastal uplift quantified by luminescence dating of marine terraces, central Cascadia forearc, Oregon | |
Reich et al. | Field‐scale subsurface flow processes inferred from continuous gravity monitoring during a sprinkling experiment |
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 |