本发明专利技术公开了一种喀斯特非连续性土壤与岩石分布的地质雷达探测方法，其步骤：A、在待测区域内选择长方形样方，移除样方中土壤表面的植物；B、在样方中按等间距设置测线；C、向土壤中插入铁钎初步探测样方内土壤深度范围；D、探地雷达型号为MALA—ProEx；E、选择测线进行测量；F、得到样方内所有测线上的GPR数据；G、在样方内按S型布点，采用比重计法测定土壤质地；H、测定的土壤质地类别和含水量；J、根据关系式计算出电磁波速；K、将样方内所有测线的GPR数据和样方内电磁波波速输入3D软件，得到样方内土壤与岩石分布的三维图像。该方法探测土壤与岩石分布，其精度高、反映直观、快速无损且适用范围较广。
【技术实现步骤摘要】
本专利技术属农业生态与环境地质
，更具体涉及一种，它适用于农业生态与环境、科研、资源勘探、地质勘查等部门对喀斯特非连续性土壤与岩石分布进行探测。
技术介绍
世界喀斯特面积约2200万km2，约占全球陆地总面积的15%。我国西南地区岩溶总面积约55万km2。喀斯特地区由于其成土过程缓慢，土层浅薄，一方面岩溶作用控制的土壤沿溶槽、溶沟垂向或侧向迁移导致土壤剖面三维空间不连续、土壤与岩石交错分布、岩石裸露率高；另一方面，由于岩溶作用，土壤二元结构发育，石-土界面明显，土壤与岩石电导率差异显著，为地质雷达的探测石-土界面提够了可行性。探明喀斯特地区土壤与岩石的分布情况，是该区水土流失防控、生态恢复、碳酸盐岩资源利用的最基本前提。以往对喀斯特地区土壤与岩石分布的探测都是采用剖面开挖的方法，该方法对喀斯特生态环境破坏大、估算不准确、费时费力。地质雷达(Ground penetrating radar,简称GPR)是基于地下介质的电导率、介电常数等电性参数的差异，利用高频电磁脉冲波的反射探测地下介质分布的一种物探手段。GPR已经被广泛应用于地质、工程、环境和考古学调查。目前，在喀斯特地区，地质雷达主要被应用到洞穴的考古学调查中，尚缺乏针对喀斯特这种特殊地质环境下非连续性土壤与岩石分布的探测方法。因此，开发一种高精度、高效、环境友好且适合喀斯特非连续性土壤与岩石分布探测的方法，是喀斯特地区土壤生态相关研究和碳酸盐岩资源调 查的迫切需要。
技术实现思路
本专利技术的目的是在于提供了一种，方法易行，操作简便，误差小(在岩土介质中500 MHz和800 MHz中心频率天线的垂直分辨率分别为:3~8 cm,2^5 cm)，能够得到喀斯特复杂地形下土壤与岩石分布的三维图。为了实现上述的目的，本专利技术采用以下技术方案: 一种，其步骤是: A.在待测区域内选择典型样方(长方形或正方形地块)，移除样方中土壤表面的植物及其凋落物； B.在样方中按等间距(O~0.8 m)设置测线，每条测线的走向和长度须一致，且起点在同一直线上； C.向土壤中插入铁钎初步探测样方内土壤深度范围(一般，坡上〈1.4m，洼地〉1.4m),便于选择不同频率的GPR天线； D.探地雷达型号:瑞典MALA公司生产ProEx系统主机与800和500MHz屏蔽天线；当土壤深度为〈1.4 m时，可选用800 MHz或500 MHz中心频率天线，土壤深度为1.4~2.2m时，采用500 MHz中心频率天线；在使用500 MHz和800 MHz中心频率天线时，其探测参数分别设置为:时窗62.09心和38.68 ns，采样频率7730.26 MHz和12409.1 MHz，天线偶极子间距分别设为0.18 m和0.14 m ;道间距(一条测线上紧邻两个测点之间的距离)按测线长短设置(测线由短到长可设置为0.01~0.1 m:测线长I~5米时，按每条测线上有100个测点为宜，测线长5~10米时道间距可设置为0.05 m，测线长>10 m时，道间距可设置为0.10 m),电磁波原始波速设置为0.1 m/ns ； E.选择第一条测线进行测量:测线较短(〈5m)且测线上每个点的位置要求精确时，可选用点测法(ProEx主机自带，按道间距测量各个测点)按设置的道间距长度移动天线；测线较长(> 5 m)时，采用轮测法(Ρι?Εχ主机自带，连续采集测线上数据)进行测量；测量测线上有岩石出露的地方在采集软件上进行标记(数据显示在命名为mrk的文件中)，得到这条测线上的GPR数据包括:命名为rd3 (数据文件，记录采集的电磁波信息)、rad (道头文件，记录采集过程参数设置信息)、mrk (标记文件，记录各测线上标记的物体)的三个文件； F.按步骤E测定样方内其它测线，得到样方内所有测线上的GPR数据； G.在样方内按“S”型布点，采用土钻按0-10cmU0-20 cm,20-30 cm、30_50 cm,50-70cm、70-100 cmUOOcm-基岩层，钻取各层次土壤样品，带回实验室测定其含水量，并计算样方内土壤平均含水量；同时，按多点采样法在样方内采集一个表层土壤的混合样品带回实验室，采用比重计法测定土壤质地(GB 7845-87)； H.根据步骤G中测定的土壤质地类别和含水量，选择不同的波速和质量含水量关系式，算出样方内电磁波的平均波速，鉴于喀斯特土壤质地粘重，重点考虑粘质土壤的电磁波波速；当测定的土壤为砂质粘壤土时，电磁波波速(V，单位m/ns)和质量含水量(β )的经验关系式为 V = 5.9184 0J - 2.1403 Θm 2 - 0.1841 Θm + 0.1546 (决定系数R2=0.9889)，当土壤为粘壤土时，关系式为 K= 9.7339 Θm 3 - 4.9935 ~ 2 + 0.4912 +0.102 (R2=0.9933)，当土壤为粉(砂)质粘土时，关系式为 κ= 12.745 Θm 3 - 9.5112 Θm 2+ 2.0415 Θm - 0.028 (R2=0.9914)； 1.分别将各条测线上采集的GPR数据(rd3文件)和计算得到的电磁波速输入Reflexw软件(探地雷达系统自配软件)，经过以下处理:一维滤波/去直流漂移(ID-Filter/Subtract-Dc-Shift)、静校正 / 移动开始时间(Static Correction/Move start time)、增益 / 能量衰减(Gain/Energy decay)、二维滤波 / 抽取平均值(2D_Filter/Subtractingaverage)、一维滤波 / 巴特沃斯带通滤波(lD-Filter/bandpassbutterworth)、二维滤波 /滑动平均(2D-Filter/Running average),得到各条测线下方岩土分布的二维图像，图像可以反映各条测线上土壤深度情况，并可自动生成所有测点的土壤深度； J.将样方内所有测线的GPR数据(rd3文件)和样方内电磁波波速输入Easy 3D软件(探地雷达系统自配软件)，经过以下处理:DC校正(DC adjustment)、自动增益控制(AGC)、删除平均道(Delete mean trace)、时域的带通滤波(FIR)、在二维点阵取平均值(Average),得到样方内土壤与岩石分布的三维图像。上述的技术方案中，所述步骤G中，100 cm深度以下不宜利用土钻采集土壤样品用于测定含水量，但该深度以下含水量变化极小，可以能钻取到的最大深度处土壤含水量代替100 Cm深度以下土壤含水量。上述的技术方案中，所述步骤H中，列举了喀斯特地区最典型三种质地土壤适用的GPR电磁波速与土壤含水量的经验关系式，这些拟合的关系式经过了多次室内模拟和原位探测试验验证。该方法探测土壤与岩石分布，其精度高、反映直观、快速无损且适用范围较广。本专利技术与现有技术相比，具有以下优点和效果: a)误差小、精度高:使用500 MHz和800 MHz中心频率天线在岩土工程探测中的垂直分辨率(误差)分别为3?8 cm、2?5 cm，可满足对一般土壤生态系统探测的需要。b)直观:该探测方法最终以三维图的形式直观反映喀本文档来自技高网...
【技术保护点】
一种喀斯特非连续性土壤与岩石分布的地质雷达探测方法，其步骤是： A、在待测区域内选择长方形或正方形地块，移除样方中土壤表面的植物及其凋落物；B、在样方中按等间距：0～0.8 m设置测线，每条测线的走向和长度一致，起点在同一直线上；C、向土壤中插入铁钎初步探测样方内土壤深度范围，坡上< 1.4 m，洼地> 1.4 m，选择不同频率的GPR天线；D、探地雷达型号：瑞典MALA公司生产ProEx系统主机与800和500 MHz屏蔽天线；土壤深度为< 1.4 m，选用800 MHz或500 MHz中心频率天线，土壤深度为1.4～2.2 m，采用500 MHz中心频率天线；使用500 MHz和800 MHz中心频率天线时，其探测参数分别设置为：时窗62.09 ns和38.68 ns，采样频率7730.26 MHz和12409.1 MHz，天线偶极子间距分别设为0.18 m和0.14 m；道间距：一条测线上紧邻两个测点之间的距离，按测线长短设置：测线由短到长设置为0.01～0. 1 m：测线长1～5米时，按每条测线上有100个测点为宜，测线长5～10米时道间距设置为0.05 m，测线长>10 m，道间距设置为0.10 m，电磁波原始波速设置为0.1 m/ns；E、选择第一条测线进行测量：测线较短< 5 m测线上每个点的位置，选用点测法按设置的道间距长度移动天线；测线长> 5 m时，采用轮测法进行测量；测量测线上有岩石出露的地方在采集软件上进行标记，得到测线上的GPR数据包括：命名为rd3、rad、mrk的三个文件；F、按步骤E测定样方内其它测线，得到样方内所有测线上的GPR数据；G、在样方内按S型布点，采用土钻按0‑10 cm、10‑20 cm、20‑30 cm、30‑50 cm、50‑70 cm、70‑100 cm、100cm‑基岩层，钻取各层次土壤样品，带回实验室测定其含水量，并计算样方内土壤平均含水量；按采样法在样方内采集一个表层土壤的混合样品带回实验室，采用比重计法测定土壤质地； H、根据步骤G中测定的土壤质地类别和含水量，选择不同的波速和质量含水量关系式，算出样方内电磁波的平均波速，测定的土壤为砂质粘壤土时，电磁波波速和质量含水量的关系式为ν = 5.9184 θm3 – 2.1403θm 2 ‑ 0.1841θm + 0.1546 ，土壤为粘壤土时，关系式为ν= 9.7339θm 3 – 4.9935θm 2 + 0.4912θm + 0.102，土壤为粉质粘土时，关系式为ν= 12.745θm 3 – 9.5112θm 2 + 2.0415θm + 0.028；I、分别将各条测线上采集的GPR数据和计算得到的电磁波速输入Reflexw软件，经过以下处理：一维滤波/去直流漂移、静校正/移动开始时间、增益/能量衰减、二维滤波/抽取平均值、一维滤波/巴特沃斯带通滤波、二维滤波/滑动平均，得到各条测线下方岩土分布的二维图像，自动生成所有测点的土壤深度； J、将样方内所有测线的GPR数据和样方内电磁波波速输入Easy 3D软件，经过以下处理：DC校正、自动增益控制、删除平均道、时域的带通滤波、在二维点阵取平均值，得到样方内土壤与岩石分布的三维图像。...
【技术特征摘要】
1.一种喀斯特非连续性土壤与岩石分布的地质雷达探测方法，其步骤是: A、在待测区域内选择长方形或正方形地块，移除样方中土壤表面的植物及其凋落物； B、在样方中按等间距:0~0.8 m设置测线，每条测线的走向和长度一致，起点在同一直线上； C、向土壤中插入铁钎初步探测样方内土壤深度范围，坡上〈1.4m，洼地〉1.4 m，选择不同频率的GPR天线； D、探地雷达型号:瑞典MALA公司生产ProEx系统主机与800和500MHz屏蔽天线；土壤深度为〈1.4 m，选用800 MHz或500 MHz中心频率天线，土壤深度为1.4~2.2 m，采用500 MHz中心频率天线；使用500 MHz和800 MHz中心频率天线时，其探测参数分别设置为:时窗62.09 ns和38.68 ns，采样频率7730.26 MHz和12409.1 MHz，天线偶极子间距分别设为0.18 m和0.14 m ;道间距:一条测线上紧邻两个测点之间的距离，按测线长短设置:测线由短到长设置为0.01~0.1 m:测线长I~5米时，按每条测线上有100个测点为宜，测线长5~10米时道间距设置为0.05 m，测线长>10 m，道间距设置为0.10 m,电磁波原始波速设置为0.1 m/ns ； E、选择第一条测线进行测量:测线较短〈5m测线上每个点的位置，选用点测法按设置的道间距长度移动天线；测线长〉5 m时，采用轮测法进行测量；测量测线上有岩石出露的地方在采集软件上进行标记，得到测线上的GPR数据包括:命名为rd3、rad、mrk的三个文件； F、按步骤E测定...
【专利技术属性】
技术研发人员：苏以荣，陈香碧，夏银行，黎蕾，
申请(专利权)人：中国科学院亚热带农业生态研究所，
类型：发明
国别省市：湖南;43
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。