项目技术总结报告
深圳桑达百利电器有限公司编号：QF7.3-048-PM
SHENZHENSANGDABAILIELECTRONICCO.,LTD.
项目技术总结报告
项目名称：部
门：起止年月：项目经理：
大功率LED智能控制恒流电源
桑达百利电源研发部
年2月10日起---201*年12月30止吴云波
桑达百利研发部制20深圳桑达百利电器有限公司编号：QF7.3-048-PM
SHENZHENSANGDABAILIELECTRONICCO.,LTD.
一．背景二．系统慨述三.技术分析与总结深圳桑达百利电器有限公司编号：QF7.3-048-PM
SHENZHENSANGDABAILIELECTRONICCO.,LTD.
一.背景随着LED技术的不断发展，推动了白光LED的问世，产业开始了绿色照明时代。由于LED能耗少、热辐射低、发光效率高，是一种节能、环保、经济、安全的新型照明器件，符合”十二五”节能减排的规划要求,大功率LED是低电压、大电流的驱动器件,目前户外照明灯具广泛使用大功率LED,因此，加快LED户外灯具配套大功率LED智能控制恒流电源技术研究成为当今首要问题。大功率LED要成为照明业的主体，其中安全、高效智能控制恒流电源技术驱动研究是推广应用大功率LED的关键,而且本项目符合国家的产业政策，是国家鼓励发展的项目。产品市场前景广阔，经济效益和社会效益显著，符合国家质量标准，所以开发本项目。二.系统慨述LED恒流电源驱动是LED电源的一种,是采用开关电源变换器，做成隔离型的恒流电源，其输出电流恒定且智能可调，设计时还要注意输入功率因数要高。主要原因是：1.避免驱动电流超出最大额定值，影响其可靠性。2.获得预期的亮度要求，并保证LED灯具各个LED亮度、色度的一致性3.智能接收外界各种控制器信号,使LED灯具开、关、亮度灰度可控.三.技术分析和总结1.恒流精度较难控制我们采用了精确的恒流控制IC,恒流取样电阻采用了1%的精密电阻而且温度系统变化小的电阻,DC-DC输出储能电感的电感量控制在5%.2.效率提升到90%以上较难整体效率由AD-DC和DC-DC恒流效率组成,我们AC-DC的效率调试到大于95%,DC-DC恒流效率大于96%,所有整体效率做到了90%以上.深圳桑达百利电器有限公司编号：QF7.3-048-PM
SHENZHENSANGDABAILIELECTRONICCO.,LTD.
3.在环境温度较高时的稳定性不理想设计时我们选用工业品元器件,特别在定制元件变压器设计时选用CLASSF等级,变压器温度等级可以达到150度.4.环境温度低温-40℃难以启动.设计时使用日本进口电容,能够在低温-40℃长期工作,特别是VCC电容的选择.5.高功率小体积的电源较难设计提高开关电源的频率,从而减少变压器的体积,PCB上尽量使用贴片元器件,PCBLAYOUT优化设计,PCBA底部增加绝缘片可以使PCB和外壳之间的距离减少从而减少了外壳的体积.6.电源的使用寿命难以达到5万小时.设计时选用了长寿命的南高温电解电容,因为电解电容在电源器件中寿命是最短的.项目负责人（签字）：吴云波开发部经理（签字）:范传玺年月日
扩展阅读：科技发展项目技术总结报告
优秀青年科技人才计划项目总结报告
项目名称：委托单位：承担单位：受资助人：起止时间：
介质含水率与探地雷达信号关系研究国土资源部科技与国际合作司中国地质科学院物化探研究所方慧
201*年10月～201*年11月
国土资源部
二○○九年十一月二十日目录
一、研究领域及资助研究项目概况............................................................................31、研究领域..............................................................................................................32、资助研究项目概况..............................................................................................3二、研究领域国内外发展趋势和前沿......................................................................41、探地雷达技术应用现状......................................................................................42、测量物质含水率的主要方法及存在问题..........................................................53、探地雷达探测物质含水率研究领域现状及发展趋势......................................6三、研究工作总结........................................................................................................91、研究项目实施情况..............................................................................................92、研究工作取得的主要成果和创新点................................................................10四、经费使用情况......................................................................................................26五、所在单位审核意见..............................................................................................261、单位对受资助者给予的支持............................................................................262、单位学术委员会对受资助者研究成果的评价................................................263、对受资助者资助期间总体工作情况的评价....................................................27主要参考文献..............................................................................................................29
一、研究领域及资助研究项目概况
1、研究领域
⑴专业领域：勘探地球物理
⑵主要研究方向：探地雷达在土壤、建筑材料含水率监测领域的应用技术研究
2、资助研究项目概况
⑴项目名称：介质含水率与探地雷达信号关系研究⑵起止时间：201*年12月201*年11月
⑶目标任务：应用数值模拟技术和物理实验技术开展介质含水率与探地雷达信号关系研究，建立探地雷达测量信号与介质含水率之间准确合理的数学关系模型，推动探地雷达方法技术进步。
⑷主要研究内容：探地雷达三维正演技术研究；介质含水率与探地雷达信号关系物理实验研究；介质含水率与探地雷达信号关系数值模拟研究。
⑸工作成果：两年来，针对上述研究内容开展了较系统的研究工作，基本实现了设计的任务目标，取得如下主要成果和认识。
①利用时间域有限差分方法实现了探地雷达三维正演计算，编制三维正演软件。解决了目前探地雷达常见软件无法模拟介质孔隙度及含水率变化的问题，为应用数值模拟技术开展介质含水率与探地雷达信号关系问题研究提供了有力工具。
②首次应用数值模拟方法研究了在介质孔隙度及含水率发生变化时，对雷达信号传播特征的影响规律。并与物理模型实验结果进行了对比分析，证明了数值模拟方法的有效性。
③根据数值模拟结果，分析了当介质的孔隙度或含水率发生改变时，介质的等效介电常数及雷达波幅值、传播速度等参数随孔隙度及含水率的变化规律，并对这些参数相对介质孔隙度或含水率变化的灵敏程度进行了对比分析。④在介质含水率数值模拟中，不仅考虑了介质电导率的影响，也考虑了高频电磁场条件下，由于极化滞后效应造成的介电损耗，并利用物理模型实验结果与数值模拟结果的对比分析，证明了这种思路的合理性。
⑤通过物理模型实验，分析了石英砂、沥青等材料的等效介电常数、雷达波频谱及传播速度等参数随介质含水率的变化规律。在此基础上，提出了介质含水率与等效介电常数之间的数学关系式。这些数学关系模型，经进一步实验验证后，可作为探地雷达检测公路材料含水性的基础。
二、研究领域国内外发展趋势和前沿
1、探地雷达技术应用现状
探地雷达是一种高频电磁法。与探空雷达相似，探地雷达利用发射天线以宽频短脉冲形式向地下发射高频电磁波，电磁波在介质电磁性质不同的界面处会产生反射，并被接收天线所接收，通过分析电磁波在时间、空间的传播特性实现探测地下目标体的空间位置、规模和物理性质等目的。探地雷达技术具有分辨率高、无损、高效等特点。
探地雷达技术的应用最早可追溯到上世纪初。早在1910年，德国的G.Leimback和Lwyc曾以专利形式阐明了这一现象。第一次正式应用是在1929年用以确定冰河的深度（Stern,1929,1930），之后这种技术几乎消失。直到1950年因有飞机失事掉进格陵兰岛的冰缝中，才再次采用探地雷达技术。受仪器性能和理论研究等因素的限制，探地雷达初期的应用仅限于波吸收很弱的冰层、盐岩矿等介质中(Cook,1964;Barringer,1965;Lundien,1966)。随着仪器信噪比的大大提高和数据处理技术的进步，七十年代以后，探地雷达的实际应用范围迅速扩大。1972年更被阿波罗号宇宙飞船带上了月球(Simmonsetal.,1972)。目前，探地雷达已广泛应用于工程勘察、考古、环境、军事等领域。在仪器制造方面，国际几大著名厂商相继推出适于不同应用的多种仪器系统。在理论研究方面，主要集中在信号处理和正反演研究等方面(Olhoeft,201*)。每两年召开一次的国际探地雷达会议基本反映出探地雷达技术的研究及应用现状。
我国探地雷达研究始于七十年代初期，原地质矿产部物探研究所、煤炭部煤炭科学院等科研单位开展过探地雷达仪器研制和野外实验工作。目前，我国已有几百家单位拥有探地雷达设备和有关技术人员，行业覆盖地质、冶金、煤炭、水利、交通、建筑、考古、环境及军事等。主要使用进口仪器，也有少量国产仪器在销售和使用。在研究方面主要是针对信号处理技术，少数大学等科研机构开展了正反演方法研究。
2、测量物质含水率的主要方法及存在问题
实际生活中，常常需要研究或了解天然物质和人工材料的孔隙度和含水率等参数随空间或时间的变化，如了解土壤、岩石的孔隙度及含水率的大小在地质灾害预防，海侵程度监测，冻土层调查，赋水层位的确定，公路、机场跑道危险隐患调查，建筑地基状况的评估及种植业管理等方面都是十分重要的指标。了解建筑材料孔隙度及含水率情况同样是评价材料质量的重要指标之一。如沥青是一种广泛用于铺设公路路面的材料，影响沥青公路质量的一个重要指标是沥青材料的孔隙度大小。其原因是由于孔隙中可能充满空气、水、冰或者它们的混合物，它们的存在会严重影响沥青材料的整体物理性质，进而造成材料质量发生变化。目前，检测物质含水量的常见方法主要有烘干法、电阻法、中子仪法、γ射线(透射)法、时域反射仪法(TDR)法等。这些方法原理不同、各具特色，有些方法简便、经济，有些方法测量精度很高。但这些方法普遍存在如下缺点：一是只能采用定点测量方式，无法实现空间上的连续测量，若开展大面积测量，成本较高；二是测量结果受采集样品或测量探头附近物质的含水状态影响较大，其测量结果有时不能准确代表物质整体含水情况；三是有些方法要求测量探头埋设在测量物质中，因此只能适用于土壤等非固结物质，无法对沥青、混凝土等建筑材料的含水情况进行检测；四是需要采样测试的方法会对检测物质造成一定程度的破坏。因此，研究精确、高效、无损的探测技术正在成为上述领域的需要。3、探地雷达探测物质含水率研究领域现状及发展趋势
目前，探地雷达的应用还主要集中在探测目标体的空间位置、几何形态等方面。在资料处理和解释中通常假设目标体及其周围介质是均匀的，然而事实上，大自然中常见物质，如土壤、岩石及人工合成建筑材料等，都是由多种成分组成的，因此这些物质的电磁性质既取决于物质组成成分的物理性质，也受物质的结构及其孔隙度、含水率、温度等多种因素的影响。在多数情况下，这些因素在一定尺度范围内并非处处相同，因此天然物质存在着不均匀性，雷达波在其中的传播特性会因此发生一定程度的改变。特别是，由于水具有较高的介容率，又是有极分子，不仅会改变物质的电导率，更会改变物质整体介电常数。因此，物质的含水状况对雷达波的传播速度和能量损耗都会产生很大影响，使雷达波在介质中的传播特性对介质含水率的变化十分敏感，因此应用探地雷达技术探测物质含水情况具有良好的地球物理前提，可以应用探地雷达技术探测介质含水率的变化情况。重要的是，相对上述几种检测物质含水率的常见方法，应用探地雷达技术测量介质含水率变化情况具有如下主要优点：
①与测量介质非接触。因此既可以应用于土壤等非固结状物质，也可以应用于岩石、沥青、混凝土等固结状物质，且不会对介质造成任何破坏，是真正的无损检测；
②现代雷达提供了高密度采样测量方式，因此可以对测量介质实现空间或时间上连续测量，相对传统的定点测量方式，更能准确地对介质整体含水状况进行评估。
③探地雷达不仅具有较高的横向分辨能力，也具有较高的纵向分辨能力和一定的勘探深度，因此可以同时对地下不同物质层（如路基不同基层）同时进行检测，且可以同时对不同层位的含水情况分别进行评估；
④测量效率高。如采用车载雷达，可以50公里/小时的速度进行测量，因此更适合开展大范围的监测。正是因为探地雷达技术具有上述优势，应用探地雷达检测介质含水率成为近年来探地雷达技术新的研究方向。然而，介质含水率的变化与雷达波的传播特性之间的关系十分复杂，需要通过开展理论和实验研究建立起两者之间的定量关系，才能使探地雷达技术真正应用于探测物质含水率这一领域。在这一方面已有一些学者开展了相关研究。
Hasted（1973）通过实验获得了25℃条件下水的介电常数随电磁波频率变化情况（图2.1），可见在高频电磁场作用下水的极化特性表现出较强的频散特征。
Topp（1980）通过实验给出的土壤介电常数与土壤含水率之间的近似关系式：
图2.1水的介电常数随频率变化曲线
5.31022.92102b5.5104b24.3106b3
（0.022m3m-3，Jacobsen和Schjonning，1994）。
（2.1）
实践证明上述实验公式可以在不同类型、成分的土壤条件下取得较高的精度
更多的学者（如Dobson1985，Roth1990，Friedman1998，Jones和Friedman201*）则是建立土壤不同组成成分的介电常数和含量多少（包括含水率的多少）与土壤整体介电常数之间的关系模型。在这些模型中，土壤整体介电常数与土壤颗粒、土壤中含水率以及土壤中的空气含量之间的关系可以用如下CRIM（ComplexRefractionIndexModel）模型来描述
1b(w(1n)s(n)a)
（2.2）其中，b为土壤整体介电常数，为土壤含水率，n为土壤的孔隙度，s、w及a分别为土壤颗粒、水及空气的介电常数。系数α与电场方向和土壤构造的相对关系有关。
上述表明前人的研究主要集中在土壤含水率与土壤整体介电常数之间的关系，这些关系式在测量土壤含水率方面取得较好效果，但不能直接用于描述其它介质的含水率与介质介电常数之间的关系。
在其它介质研究方面，S.Laurens等研究了混凝土材料中含水率变化与雷达信号之间的关系，讨论了雷达信号的速度、介电常数、幅值及相位与含水率之间的关系。LanboLiu等研究了沥青材料中孔隙度、含水率与整体介电常数之间的关系，研究表明：在干燥条件下，介质孔隙度的变化对整体介电常数影响不大；随着含水率的增加介质整体介电常数明显增大。但研究尚不够深入。
在应用方面，美国在应用探地雷达技术监测葡萄园土壤水分状况方面取得了很好的应用效果；一些国家在应用探地雷达技术监测高速公路路基含水情况方面也开展了实验和实际应用工作。
从国内刊物发表的文章看，国内只有少数科技人员开展了部分研究工作，如杨厚荣等开发了WPRT-1型原油持(含)水率雷达测井仪，巧妙地利用了雷达探测技术,可将原油持水率的测量范围扩大到0～100%,测量精度达1%,有效地解决了油田高含水率生产的测量问题。
在其它方面，如能否应用数值模拟技术研究介质含水率变化与雷达信号之间的关系，以及其它因素如环境温度、介质孔隙的大小等因素对测量结果的影响等方面的研究工作开展得还很少。
研究介质含水率与探地雷达信号之间的关系可以利用数值模拟和物理模拟两种方法。数值模拟方法方便、灵活，但由于数值模拟过程中进行了一定程度的近似处理，因此数值模拟结果与物理模拟结果会有一定差别，可以用于研究一般性规律。物理模拟更能准确地反映特定环境下物质含水率变化对雷达波的影响规律，对探地雷达资料精确解释是必不可少的，因此后者应用较为普遍。虽然，探地雷达数值模拟技术最近十年得到了较大发展，国际上一些大学如荷兰Delft理工大学、美国俄亥俄州立大学、科罗拉多矿业学院及一些商业公司开展了大量研究工作，已有商业或免费二、三维正演软件出现。国内中国矿业大学等科研机构也有学者开展了数值模拟技术研究。但是，在数值模拟研究领域，除算法研究外，应用研究主要集中在研究均匀介质中局部不均匀体的响应，雷达天线极化特性，介质频散特性等方面。由于这些研究中都假设介质是均匀的，因此目前大多数软件不能用于模拟介质不均匀变化，也就难以直接用来模拟介质含水情况。也很少见到利用数值模拟技术研究介质含水率变化对雷达波的传播特性影响规律等方面的报道。
总之，应用探地雷达技术可以对大规模的测量介质无损、高效、低成本地实现横向与纵向空间连续观测，是近年来探地雷达技术研究的新方向。但是，在建立探地雷达信号与介质含水率之间的定量关系方面工作开展得还很不够，需要开展更系统的物理实验；在数值模拟方面还需要研制开发更适合的计算软件。
三、研究工作总结
1、研究项目实施情况
根据合同要求，按计划全面开展了研究工作，主要包括探地雷达三维正演技术研究；介质含水率与探地雷达信号关系物理实验研究；介质含水率与探地雷达信号关系数值模拟研究等研究内容，全面完成了合同规定的研究任务。完成的主要工作量包括：
①研制一套探地雷达三维正演软件；
②完成了石英、沥青砂等不同材料的探地雷达检测含水率的物理模型试验；③完成了部分理论模型的数值模拟研究；
④完成了物理实验和数值计算数据的整理和分析；
⑤初步建立了适于沥青、石英砂等材料的含水率-介电常数数学关系式。2、研究工作取得的主要成果和创新点
（1）受资助期间取得的主要研究成果
成果一利用时间域有限差分法实现了探地雷达三维正演计算，研制了一套探地雷达三维正演软件，解决了目前探地雷达常见软件无法模拟介质孔隙度、含水率变化的问题，为应用数值模拟方法开展介质含水率问题研究提供了有力工具。
自然条件下，水赋存于介质内部的孔隙中，而且在非饱和状态下，一部分孔隙含水，另一部分孔隙则充填为空气。因此，开展介质含水率数值模拟时，要求正演程序要能实现均匀介质中随机分布有不同比例的自由水或空气等“杂质”，而目前常见的正演软件均假设介质是均匀的，无法模拟介质含水状态，要应用数值模拟技术研究介质含水率与探地雷达信号关系，首先要研制一套合适的探地雷达三维正演软件。
1.1基本理论
实验表明，所有的电磁现象都服从麦克斯韦方程，在时间域中，麦克斯韦方程有如下形式：
EHBt
(3.1)
(3.2)
DJtDB0
(3.3)
(3.4)
在导电介质中，对于单色电磁波，利用DE，BH，jE的关系及场矢量EE(r)eit、BB(r)e，很容易推导出频率域波动方程的解为
EE0esrei(krt)(3.5)sri(krt)BB0eeit其中k和s与介质物性参数和电磁波频率有关。值得指出的是，水分子在交变电场情况下容易被极化，偶极距随交变电场不断改变方向。受分子的惯性影响，偶极子的取向需要一定时间（驰豫时间），出现极化滞后现象，即电场和感应偶极矩之间出现了相位差，这时水的介电常数实际为复数，即：Re()iIm()，此时，3.5式中k和s分别为：
1Im()k"1()21Re()2(3.6)1Im()212s"1()12Re()因子S决定振幅随传播距离的减速度，故称S为衰减系数。电磁波在介质中的
12传播速度v由因子K决定
vk112(3.7)
Im()21"1()1Re()2可以看出，介电常数的实部为介质的固有极化，虚部为介电损耗，并与导电
率合并成为介质的有效导电率。实验表明介电常数的实部和虚部还随频率变化，变化规律可以用Debye公式来描述（Debye，1929）
(f)s1(iffrel)(3.8)
其中，s、为直流和极高频状态下介质的介电常数，frel为弛豫频率。对于25℃条件下的自然水，其Debye参数为εs=80.1，ε∞=4.2，frel=17.1GHz（Hasted,1973）。根据Debye模型可以计算出当雷达波频率为1GHz时，由于水分子的迟豫作用产生的附加电导率约为0.265s/m，这与自然状态下常见物质的电导率在同一量级或更高一些，因此在介质含水率模型正演计算中，必须考虑水分子所产生的附加电导率的影响。
1.2数值模拟的实现
探地雷达三维正演计算最常采用的方法是时间域有限差分法。时间域有限差分法最早由Yee于1966年提出，是一种对麦克斯韦方程进行离散化的简单实用技术。Yee巧妙地在剖分单元内使电场分量与磁场分量在时间和空间上相互分离，利用中心差商代替微商把连续变量离散化，使含时间变量的麦克斯韦旋度方程离散以后构成显式差分方程，从而可以在时间上迭代求解，而不需要进行矩阵求逆运算。由给定相应电磁问题的初始值，应用时间域有限差分法就可以逐步推进地求得以后各个时刻空间电磁场的分布。
具体计算流程如下：
否根据发射脉冲信号设置下一时间步长发射天线的电场值是是否完成设计时间步长？应用吸收边界条件计算边界上的电场值根据各单元磁场值计算各单元下一时间步长的电场值设定发射天线位置、长度及极化方向设定发射脉冲信号设计时间步长设计正演模型：计算空间尺度；剖分单元尺度；各剖分单元电/磁性参数赋值。为计算空间各单元电场分量赋初始值根据各单元电场值计算各单元磁场值
1.3数值计算结果检验
根据电磁场理论可以推导出球坐标系下描述自由空间中电偶极子辐射场时空分布的解析公式如下:
r2cc2cc2Er,t22cose22sinpt(3.9)er4rrtrtrtrc结束输出数据文件式中p(t)为偶极子的电偶极矩，c为光速。
为方便求解，发射信号选为高斯脉冲，其函数的时域形式为：
4tt02(3.10)Eitexp2其中τ为常数，决定了高斯脉冲的宽度，脉冲峰值出现在tt0时刻，如图3.1所示，实际计算中τ=2ns。图3.2、3.3分别为电偶极子下方10和25厘米处平行偶极子方向电场的解析解和时间域有限差分法三维正演对比结果。
对比结果显示，在电偶极子下方10厘米处，解析解和数值解十分接近；25厘米处的结果显示在信号晚时两者存在少许差别，其原因可能是边界效应引起的。数值解与解析解的一致性证明研制的正演软件是可靠的。
0.050.040.030.020.010-0.01-0.0210.8E(v)0.60.40.201*2345t(ns)图3.1高斯脉冲信号
解析解0.5解析解数值解
0.4幅值(v/m)0.30.20.100-0.1
24681012幅值(v/m)数值解时间（ns）
时间（ns）
图3.2天线下方10厘米处数值解和
解析解计算结果图3.3天线下方25厘米处数值解和
解析解计算结果
成果二首次应用数值模拟方法进行了介质孔隙度及含水率变化对雷达信号传播特征影响规律的研究。在介质含水率数值模拟中，不仅考虑了介质电导率的影响，也考虑了高频电磁场条件下，由于极化滞后效应造成的介电损耗，并利用物理模型实验结果与数值模拟结果的对比分析，证明了这种思路的合理性。研究结果表明应用数值模拟方法研究介质孔隙度及含水率的变化对探地雷达信号传播特性的影响规律是可行的。
2.1介质孔隙度变化与雷达信号传播特性关系2.1.1数值模拟模型设计：
数值模拟模型由相对介电常数（r）、电导率（/sm-１）横向均匀的三层介质组成，其中第一层介质为空气；模型及其坐标系统见图3.4，即：各层介质物性参数：r1＝1.0、1＝0.0s/m，r2＝5.4、2＝0.0s/m、r3＝15.0、3→∞；
模型整体尺寸：0.50×0.50×0.50m，第一界面位于XOY平面上0.25m，第二界面在0.45m处；
雷达波场源为电偶极子天线，位于第一界面中心上方0.01m处；电偶极子极化方向平行X轴，主频1GHz；
剖分单元尺寸：0.005×0.005×0.005m。为研究介质孔隙度变化对雷达信号传播的影响：①在第二层介质中加入随机分布，相对介电常数r＝1.0、电导率=0.0s/m的异常单元，借以模拟均匀介质中分布有充满空气的孔隙；②异常单元所占模型剖分单元的体积百分比分别为：2％，4％，，10％。
图3.4模型示意图
模型参数：第一层：εr1＝1.0，σ1=0.0(s/m)，厚度h1=25cm
第二层：εr2＝5.4，σ2=0.0(s/m)，厚度h2=20cm第三层：εr3＝15.0，σ3→∞
空气天线xy=5.40s/mz
2.1.2模型正演结果分析
介质的孔隙可看作均匀介质中加入具有另一种介电常数的“杂质”。雷达波的传播特性会随着“杂质”含量的多少发生改变。
图3.5为介质孔隙度分别为2%、4%、、10%时，通过天线中心，平行电偶极子极化方向剖面上t=2.0ns时刻电场分量Ex空间分布状况，可见随着孔隙度增大，电场的背景“噪声”明显增强。
计算结果表明：随着介质孔隙度增大，由于孔隙内充填的空气使介质的等效介电常数减小（图3.6），雷达波传播速度增大（图3.7），反射信号的走时明显减小，信号的幅值也有小幅增大（图3.8）。雷达波传播速度由孔隙度为零时的12.91cm/ns近似线性地增加到孔隙度为10%时的13.61cm/ns，即孔隙度每增加
图3.5t=2.0ns时孔隙度不同的介质中垂直电偶极子方向电场Ex空间分布图
1：孔隙度＝0%；2：孔隙度＝2%；3：孔隙度＝4%；
6：孔隙度＝10%4：孔隙度＝6%；5：孔隙度＝8%；
1%，反射波传播速度平均增大0.54%。介质等效介电常数由孔隙度为零时的5.29减小到孔隙度为10%时的4.86，即孔隙度每增加1%，等效介电常数平均减小0.8%。
不同的观测参数相对介质孔隙度变化的灵敏程度也不一样。图3.9为不同观测参数的灵敏度曲线。可以看出介质的等效介电常数对于介质孔隙度变化最灵敏，反射波速度次之，这是因为介电常数与电磁波传播速度呈平方反比关系。信号的幅值相对灵敏度较低，这是因为信号幅值的改变主要是因为介质孔隙度的变化改变了介质的等效介电常数，从而改变了介质与相邻介质间的反射和折射系数，使电磁波的传播发生改变，但这种改变是有限的。另一方面，介质孔隙度的存在又有可能使雷达波产生散射从而降低幅值的信噪比，因此实际工作中若根据反射信号幅值反演介质孔隙度可能较其它两个参数反演精度低。
5.413.8速度（cm/ns）等效介电常数13.65.213.413.25134.80246810
12.80246810
0.0278孔隙度(%)孔隙度（%）
图3.6介质等效介电常数随孔隙度变化曲线
图3.7雷达波传播速度随孔隙度变化曲线
1.2反射波幅值（v/m）等效介电常数雷达波速度灵敏度（dB）0.02760.80.02740.40.02720反射信号幅值
0.0270246810
-0.40246810孔隙度（%）孔隙度（%）
图3.8反射波幅值随孔隙度变化曲线图3.9不同观测参数随孔隙度变化灵敏度曲线2.2介质含水率变化与雷达信号传播特性关系
介质含水率变化对雷达信号的影响有时被认为是干扰，有时又可被利用。第一章中已经提到由于水的特殊性质，即具有高介电常数和相对较大的附加电导率，使电磁波在含水介质中的传播变得更加复杂，因此，介质含水率模型可以看作是均匀介质中加入同时具有高介电常数和高电导率“杂质”的模型。事实上，介质中含水率大小与介质孔隙度是密切相关的。介质在饱和状态下，孔隙度越大，介质含水率越高。当介质处于不饱和状态时，介质中既含有水也含有空气。这里将计算两类模型，模型一中只考虑介质含水率变化，模型二中同时考虑含水率和孔隙度的变化。
2.2.1数值模拟模型一设计
数值模型仍由相对介电常数（r）、电导率（/sm-１）横向均匀的三层介质组成，其中第一层介质为空气；模型及其坐标系统见图3.10。
空气天线x=5.40s/myz
图3.10模型示意图
模型参数：第一层：εr1＝1.0，σ1=0.0(s/m)，厚度h1=35cm
第二层：εr2＝5.4，σ2=0.0(s/m)，厚度h2=10cm第三层：εr3＝15.0，σ3→∞
介质中存在的自由水可以看作是均匀介质中加入一种同时具有高介电常数和高电导率的“杂质”。雷达波的传播特性会随着“杂质”含量的多少发生改变。为研究介质含水率变化对雷达信号传播的影响，在第二层均匀介质中加入一定体积百分比，呈随机分布的高介电常数、高电导率的异常单元。异常单元：①所占模型中剖分单元的体积百分比分别为：2％，4％，，10％；②由于水为有极分子，在高频电磁场作用下，水分子的驰豫作用会产生附加电导率；因此在数值模拟计算时，必须考虑其附加电导率对雷达波的影响。依据第一章提到的Debey模型，在天线（雷达波场源）主频为1GHz的条件下，异常单元的物性参数取：
相对介电常数r＝80.1；
水自身电导率为＝0.1/sm-1；驰豫作用产生的附加电导率Im()＝0.265/sm-1，并假设附加电导率在雷达波主频附近的小范围内不随频率而改变。
2.2.2模型正演结果分析
图3.11为介质含水率分别为2%、4%、、10%时，通过天线中心，分别平行电偶极子极化方向的剖面上t=2.0ns时刻电场分量Ex空间分布状况。可见：
①随着含水率增大，雷达波场产生越来越强的散射现象，电场的背景“噪声”明显增强；
②介质等效介电常数由不含水时的5.4，增大到含水率为10%时的9.8（图3.12），即含水率每增加1%，等效介电常数平均增大7.7%。随着含水率增大，雷达波的传播速度明显减小（图3.13）。由不含水时的12.9cm/ns近似线性地减小到含水率为10%时的9.57cm/ns，即含水率每增加1%反射波传播速度平均减小2.9%。信号的绝对幅值则由不含水时的114.4mv/m减小到含水率为10%时的84.9mv/m（图3.14），即含水率每增加1%，信号绝对幅值平均减小1.9%，从而造成探地雷达的有效勘探深度也随之减小。
总之，介质含水率的改变会使雷达波的空间传播特性发生强烈改变。作为干扰源，介质含水率的变化，会降低探地雷达的有效勘探深度和对弱不均匀体的空间分辩能力；对目标层（体）埋藏深度的估计可能出现偏差。
正因为电磁波对介质含水率变化如此灵敏，所以探地雷达技术更适合用于探测含水率分布不均匀介质情况。根据模型正演得到的反射波传播速度、幅值及介质的等效介电常数与介质含水率之间的关系式可用来指导探地雷达的资料解释。
图3.11t=2.0ns时含水率不同的介质中垂直电偶极子方向剖面电场
Ex分量空间分布图
1：含水率＝0%：2：含水率＝2%；含水率＝4%；4,含水率＝6%5：含水率＝8%；6：含水率＝10%图3.15为不同观测参数的灵敏度曲线。可以看出介质的等效介电常数对于介质含水率变化最为灵敏；反射波速度与信号幅值的灵敏度相当。
此外，由于水在高频电磁场作用下能产生较强的附加电导率，致使信号的幅值随介质含水率变化发生明显改变；水又具有较高的介电常数使得电磁波的传播速度随介质含水率的变化十分敏感。但自然界中的常见物质大都不具有这样双重性质，这也就为应用探地雷达方法区分水与非水提供了十分有利的前提。
等效介电常数1013速度（cm/ns）012345678910
9128711610590246810
0.12含水率(%)
含水率(%)
图3.12等效介电常数随含水率变化曲线图3.13雷达波速度随含水率变化曲线
信号幅值（v/m）灵敏度（dB）6543210等效介电常数传播速度信号幅值
0.11
0.10.090.0802468100246810
含水率（%）含水率（%）
图3.14反射波幅值随含水率变化曲线图3.15不同测量参数灵敏度曲线
2.2.3数值模拟模型二设计
模型二是为了与物理实验结果进行对比，以检验数值模拟技术的有效性。物理实验结果取自法国作者S.Lauren发表在“Non-destructiveTestinginCivilEngineering201*”会议论文集中的文章中。作者为研究混凝土含水率变化与介质等效介电常数及雷达波传播特性的关系，进行了物理实验。实验模型为立方体，长和宽均为25cm，高为7cm。雷达天线放在模型的上方，天线主频为1.5GHz。模型的底面放在一金属板上，以使雷达波产生全反射（图3.16）。实验结果得到了介质等效介电常数、传播速度随介质含水率变化的关系曲线（图3.17）。根据实验数据进行了CRIM模型分析，在CRIM模型中，介质固体材料的介电常数取为4.0，孔隙度取14%。
图3.16混凝土实验模型（左）及测量装置（右）
图3.17实验结果
图3.18CRIM模型与实测结果对比
为了模拟上述物理实验结果，数值模型仍由相对介电常数（r）、电导率（/sm-１
）横向均匀的三层介质组成，其中第一层介质为空气，第二层介质的相对介电常
数和厚度与物理模型相同，分别取r2＝4.2，h2=7cm，雷达天线主频设计为1.5GHz，第三层取电导率为无穷大，以模拟金属板。模型及其坐标系统见图3.19。空气天线=4.20s/mxyz
图3.19模型示意图
模型参数：第一层：εr1＝1.0，σ1=0.0(s/m)，厚度h1=25cm
第二层：εr2＝4.2，σ2=0.0(s/m)，厚度h2=7cm第三层：εr3＝15.0，σ3→∞
为研究介质孔隙度及含水率变化对雷达信号传播的影响，孔隙度及含水率变化异常场的设计，是在第二层均匀介质中加入呈随机分布的异常单元：
①异常单元所占模型剖分单元的体积百分比为14%；
②异常单元中一部分充填水，另一部分充填空气，二者在空间上均为随机分布。填水单元所占模型剖分单元的体积百分比分别为0％，2％，4％，，14％，同时，充填空气的异常单元所占模型剖分单元的体积百分比分别为14％，12％，，2％，0％。
③依据Debey模型，在天线（雷达波场源）主频为1.5GHz的条件下，异常单元的物性参数取：
相对介电常数r＝80.1；
水自身电导率为＝0.1/sm-1；弛豫作用产生的附加电导率Im()＝
0.6/sm-1，并假设附加电导率在雷达波主频附近的小范围内不随频率而改变。
图3.20为数值模型计算结果与物理实验结果对比情况。可以看出，无论是反射信号幅值，还是根据速度换算出的介质等效介电常数，其数值模拟结果与物理实验结果都吻合较好，特别是当介质含水率较低时（含水率小于6%），两种结果更为接近。随着含水率增大，介质的不均匀程度增强，干扰程度随之加大，计算结果相对较分散。但数值模拟技术能较好地计算出反射信号幅值的相对变化情况。
数值模拟结果与物理模型实验结果对比分析证明了数值模拟技术的可行性和有效性，也说明了在数值模拟计算中由于引入了附加电导率，使其结果更为合理。因此，数值模拟计算技术可以作为研究介质不均匀性对电磁波传播特性影响规律的有效手段。
反射信号幅值相对改变量（dB）
介质等效介电常数含水率（%）
物理实验数据
含水率（%）
物理实验数据
数值模拟计算曲线
数值模拟计算曲线
图3.20物理实验与数值模拟结果对比
成果三开展了石英砂、公路沥青等材料的探地雷达检测含水率物理模型实验。通过根据物理实验结果，分析了这些材料的等效介电常数、幅值、传播速度、频谱特征等参数随介质含水率的变化规律；用数值模拟方法和前人的试验数据对实验结果进行了检验，证明了实验结果的可靠性，试验结果表明应用探地雷达技术检测介质含水率是可行的，其无损、非接触、高效的特性表现出其它检测介质含水率的传统技术所不具备的优势。
（2）资助期间发表论文、论著、报告及专利情况资助期间发表论文4篇，完成研究报告7部。
⑴ExperimentaldeterminationofbulkdielectricpropertiesandporosityofporousasphaltandsoilsusingGPRandacyclicmoisturevariationtechnique，Geophysics，Vol.71,No.4，第三作者；⑵StructuralfeaturesoftheCoqenbasinincentralTibetbymagnetotelluricsounding，JournalofChinaUniversityofGeosciences，Vol.18,SpecialIssue,第一作者；
⑶青藏高原措勤盆地大地电磁测量初步结果，物化探计算技术，201*年增刊，第一作者；
⑷介质含水率与探地雷达信号关系数值模拟研究，物探与化探，Vol.33,No.5,201*第一作者；
⑸羌塘盆地石油地质走廊大剖面综合地质调查（玉盘湖双湖段）大地电磁测深调查报告，第一作者；
⑹措勤盆地石油地质走廊大剖面综合地质调查（洞错东措勤段）大地电磁测深调查报告，第二作者；
⑺羌塘盆地龙尾湖区块连续电磁剖面（CEMP）测量成果报告，第一作者；⑻羌塘盆地关键地段连续电磁剖面测量成果报告，第一作者；⑼青藏高原石油地质调查大地电磁测量综合研究报告，第一作者；⑽青藏高原非地震油气勘探方法技术综合研究成果报告，第三作者；⑾青藏高原油气资源战略选区调查与评价成果报告，主笔一章。（3）负责和参加相关科研项目情况
资助期间，除完成资助课题外，还承担了如下科研工作。
⑴参加国家油气专项“青藏高原油气资源战略选区调查与评价（XQ201*06）”过程中，负责完成了如下专题：
①洞错东措勤石油地质走廊大剖面大地电磁测量（科油[201*]10号），201*年，经费130万，第一负责人；
②羌塘盆地龙尾错区块连续电磁剖面（CEMP）测量（科油[201*]15号），201*年，经费50万元，第一负责人；③羌塘盆地重点地区电磁阵列剖面测量（科油[201*]07号），201*年，经费256万元，第一负责人；
④青藏高原油气勘探方法技术综合研究（科油[201*]13号、科油[201*]09号、科油[201*]06号、科油[201*]07号），201*-201*年，经费366万元，第二负责人。
⑵国土资源部公益性行业科研专项“重要成矿区带地球物理深部探测与研究示范（201*11039）”，201*-201*，总经费304万，第一负责人。
⑶国土资源大调查项目“龙门山及邻近构造带综合地球物理调查（121201*914049）”，201*-201*年，201*年度经费200万元，第一负责人。
⑷国土资源大调查与物化探所基本科研业务费联合资助项目“永久冻土区天然气水合物电磁异常特征研究”，201*年，总经费80万元，第一负责人。
⑸国土资源大调查项目（121201*980046）“松辽外围中新生代盆地群油气地质综合调查”，201*-201*年，201*年度经费180万元，第三负责人。
四、经费使用情况
见决算报告。
五、所在单位审核意见
1、单位对受资助者给予的支持
为完成本项目，中国地质科学院物化探研究所在人力、财力和物力方面给与了受资助人大力支持。不仅提供了用于开展物理模拟实验所需的实验室、探地雷达及计算机等设备，还免收了项目参加人员的工资福利等费用。
2、单位学术委员会对受资助者研究成果的评价
探地雷达是一种快速高效地对浅层目标体进行精细探测的手段，其应用领域十分宽广。受资助者根据实际应用中发现的问题，以开发新的检测介质含水率的方法技术，拓宽探地雷达应用领域为目的，以介质含水率对雷达波传播特性影响规律研究为切入点，选定“介质含水率与探地雷达信号关系研究”为题，既具有理论意义，又具有实用价值。
研究工作在以下几方面具有创新性：
⑴利用时间域有限差分方法实现了探地雷达三维正演计算。解决了目前探地雷达常见软件无法对模型介质加入随机扰动参数问题，为应用数值模拟技术开展介质不均匀性问题研究提供了有力工具。
⑵采用数值模拟方法研究了介质含水率变化对雷达波传播特性的影响规律，并通过与物理实验结果的对比分析，检验了所实现的数值模拟方法的有效性。
⑶在研究介质含水率影响的数值模型中，不仅考虑了介质导电率的影响，也考虑了高频电磁场条件下由于极化滞后效应造成的介电损耗，并通过物理模型实验数值模拟结果的合理性。
⑷通过探地雷达与传统介质含水率检测方法（TDR法）物理实验结果的对比研究，证明了探地雷达是检测介质含水率的一种有效方法，而其本身所具有的无损、非接触、高效的特性使探地雷达可能成为一种实用的检测介质含水率的方法。
⑸通过物理模型实验，分析研究了石英砂、沥青公路材料的等效介电常数随介质含水率的变化规律，在此基础上，提出了介质含水率与等效介电常数之间的数学关系式。
总之，受资助者完成的研究工作很有意义，并取得了一些具有创新性的进展和认识，较圆满地完成了设计的工作任务。
中国地质科学院地球物理地球化学勘查研究所学术委员会主任：胡平
3、对受资助者资助期间总体工作情况的评价
方慧同志任我所油气与深部地球物理研究室主任，是物化探所国家科技创新基地地球物理研究领域的科研骨干。受资助期间，不仅较好地完成了受资助的科研工作，同时，作为负责人还先后承担了国家专项“青藏高原油气资源战略选区调查与评价”项目中的4个专题，国土资源部公益性行业科研专项1项，地质大调查工作项目2项，其它项目1项。带领科研团队多次进入青藏高原藏北无人区、青海祁连山地区和川西龙门山地区开展工作，在探测羌塘盆地等青藏重点含油气盆地构造格架方面，在研究青藏高原油气资源战略调查与评价有效物化探方法技术组合方面，在探索冻土区天然气水合物调查物化探有效技术方面，在探测龙门山构造带深部结构及地震活动性等方面取得了不少新的认识和成果，成绩显著。
该同志基础理论扎实，学术思想敏锐，有强烈的创新意识和很好的职业道德以及团队合作精神，具有良好的科研素质。
中国地质科学院地球物理地球化学勘查研究所所长：韩子夜主要参考文献
1.Annan,A.P..Radarinterferometrydepthsounding:partI-theoreticaldiscussion.Geophysics,
1973,38(3):557~580.
2.B.EngquistandA.Majda,“Absorbingboundaryconditionsforthenumericalsimulationof
waves,”Math.Comput.,1977,31:629~651
3.D.J.Daniels,SurfacePenetratingRadar,ShortRunPressLtd，1988
4.Dobson,M.C.,F.T.Ulaby,M.T.Hallikainen,andM.A.El-Rayes,Microwavedielectric
behaviourofwetsoil-partII:dielectricmixingmodels.IEEETrans.onGeoscienceandRemoteSensing,1985,GE-23:35~46.
5.FawwazT.,ThomasH.etal,Microwavedielectricpropertyofdryrocks,IEEE.Trans.on
GeoscienceandRemoteSensing,1990,28(3):325~336.
6.Friedman,S.P.Asaturationdegree-dependentcompositespheresmodelfordescribingthe
effectivedielectricconstantofunsaturatedporousmedia.Wat.Resour.Res.1998,34:2949~2961.
7.GaryR.Olhoeft,Maximizingtheinformationreturnfromgroundpenetratingradar,Journalof
appliedgeophysics,201*,43:175~187
8.G.Mur,Absorbingboundaryconditionsforthefinite-differenceapproximationofthe
time-domainelectromagnetic-fieldequations,IEEETrans.Electromagn.Compat.,1981,EMC-23:377~382.
9.Grote,K.,S.S.Hubbard,andY.Rubin.GPRmonitoringofvolumetricwatercontentinsoils
appliedtohighwayconstructionandmaintenance.LeadingEdgeofExpl.201*.21:482~485.10.Hubbard,S.S.,K.Grote,andY.Rubin,Mappingthevolumetricsoilwatercontentofa
Californiavineyardusinghigh-frequencyGPRgroundwavedata,.LeadingEdgeofExpl.201*,21:552~559.
11.Huisman,J.A.,C.Sperl,W.Bouten,andJ.M.Verstraten.Soilwatercontentmeasurementsat
differentscales:accuracyoftimedomainreflectometryandgroundpenetratingradar.J.ofHydrol.201*,245:48~58.
12.J.A.Huismanetal.Measuringsoilwatercontentwithgroundpenetratingradar:areview,
VadoseZoneJournal,201*,2:476~491.
13.J.R.Wang,Thedielectricpropertyofsoil-watermixturesatmicrowavefrequencies,Radio
Science,1980,GE-18:288~214.J.R.Wang,ThomasJ.Schmugge,Anempiricalmodelforthecomplexdielectricpermittivity
ofsoilsasafunctionofwatercontent,IEEETrans.Electromagn.Compat.,1981,EMC-23:377~382.
15.J.L.Davis,A.P.Annan,Ground-penetratingradarforhigh-resolutionmappingofsoiland
rockstratigraphy,Geophysicalprospecting,1989,37:531~551
16.Jones,S.B.,andS.P.Friedman,Particleshapeeffectontheeffectivepermittivityof
anisotropicorisotropicmediaconsistingofalignedorrandomlyorientedellipsoidalparticles.Wat.Resour.Res.201*.36:2821~2833.
17.J.Q.Shangetal.Dielectricconstantandrelaxationtimeofasphaltpavementmaterials,
Journalofinfrastructuresystem,1999,5(4):135~142
18.K.S.Yee,“NumericalsolutionofinitialboundaryvalueproblemsinvolvingMaxwell’s
equationinisotropicmedia,”IEEETrans.AntennasPropagat.,1966,14:302~30719.K.Shlageretal.,“ASurveyoftheFinite-DifferenceTime-DomainLiterature”,
20.LanboLiu,Ground-Penetratingradar:whatcanittellaboutthemoisturecontentofthehotmix
asphaltpavement?，TRB201*annualmeeting，201*
21.Ledieu,J.,P.DeRidder,P.DeClerck,andS.Dautrebande.Amethodofmeasuringsoil
moisturebytimedomainreflectometry.J.ofHydrol.1986.88:319~328.
22.MarttiT.H.,FawwazT.U.etal.,MicrowavedielectricbehaviorofwetsoilPartI:
Empiricalmodelsandexperimentalobservations,IEEE.Trans.onGeoscienceandRemoteSensing,1985,GE-23:25~34
23.S.Laurens,J.P.Balayssac,etal.NondestructiveevaluationofconcretemoisturebyGPR
technique:experimentalstudyanddirectmodeling,Internationalsymposium:Non-destructivetestingincivilengineering201*
24.Topp,G.C.,J.L.Davis,andA.P.Annan,Electromagneticdeterminationofsoilwater
content:measurementsincoaxialtransmissionlines,WaterResourcesResearch,1980,16(3):574-~582,
25.王为民，电磁场论，长春地质学院，1982
26.葛德标等，电磁波时域有限差分方法，西安电子科技大学出版社，201*.
27.米萨克N纳比吉安主编，赵经祥等译。勘查地球物理-电磁法，第一卷。北京：地质
出版社，1992。
28.何兵寿，等．FDTD法在地质雷达正演模拟中的应用，中国煤田地质，1999(4):81~8329.谭捍东，大地电磁三维正反演问题研究，博士学位论文，201*
30.黄南晖，有耗媒质中电磁波的传播特性，地球科学，1993,3(18):257~265..
友情提示：本文中关于《项目技术总结报告》给出的范例仅供您参考拓展思维使用，项目技术总结报告：该篇文章建议您自主创作。
来源：网络整理 免责声明：本文仅限学习分享，如产生版权问题，请联系我们及时删除。
公文素材库友情提示：本网站所有内容为共享上传提供，不涉及任何商业利益，本站对此不承担任何保证责任！
Copyright © www.bsmz.net Corporation, All Rights Reserved 共享时代 共享你我他 版权所有