MATLAB语言在非线性最小二乘估计中的应用

伴随着测量事业的发展,要编制测量参数估计的矩阵运算程序越来越复杂,难度也越来越大。MATLAB与其他编程语言相比,具有简单直观、速度快等特点。本文将MATLAB语言运用于非线性模型参数估计,并给出了单纯形法进行非线性参数估计的函数程序。     

1:100万中华人民共和国数字地质图空间数据库

中华人民共和国1∶100万数字地质图数据库是根据统一标准和要求编制的专题数字地质图,以中华人民共和国1∶50万数字地质图数据库为数据源,在综合研究的基础上,充分吸收近年来地质调查1∶25万、1∶20万、1∶5万等区调工作新成果、新资料,应用地质编图新技术、新理论和新方法编制而成。以年代地层单位为主、辅以岩石地层的表示方法,侵入岩按“岩性加年代”表示方法,进一步提高了基础地质的研究程度。地质图数据库内容丰富,信息量大,数据量约为1.2 GB,标示了岩石地层单位5347个图例,侵入体“岩性加时代”单位1780个图例,跨省区重要断裂93条,各省（市、自治区）内重要断层558条,同位素年龄数据1545个（组）,有代表性的钻孔382个。所有地质体的面元及线元、同位素年龄和钻孔都建立了相应的属性,相邻图幅之间进行了接图处理,编写了编图说明书和元数据。数据库采用统一的线型库、符号库、色标库等,成图过程全部采用MapGIS6.5平台计算机辅助成图,成图精度高,质量好,符合设计要求,全国64幅图采用分4个片区的工作方法,最后统一编制而成。该数字地质图是目前中国资料最全、内容最新的1∶100万地质图,是中国第一份应用GIS技术的1∶100万数字地质图的最新成果,充分反映了中国地质构造特点和当前地质研究的新水平。     

均匀充电矩形薄板电场复形调优反演计算

针对充电电场反演计算参数较多以及复形调优算法处理约束条件的直接有效性,作者就均匀充电矩形薄板反演问题编制程序,计算地质体理深、倾角及大小。尽管迭代次数较高,计算结果颇为理想,取得了较好的效果。     

Study of red tide prediction model for the Changjiang Estuary

ＩＮＴＲＯＤＵＣＴＩＯＮＡｑｕｉｃｕｌｔｕｒｅａｒｅａｒｅｄｔｉｄｅｓｉｇｎａｌｓｏｃｅａｎｐｏｌｌｕｔｉｏｎａｎｄｄｉｓａｓｔｅｒ.Ａｌｔｈｏｕｇｈｅｎｖｉｒｏｎｍｅｎｔｓｃｉｅｎｔｉｓｔｓｈａｖｅｃｏｎｄｕｃｔｅｄｍａｎｙｉｎ ｄｅｐｔｈｒｅｓｅａｒｃｈｅｓｏｎｔｈｅｃａｕｓｅａｎｄｆｏｒｍａｔｉｏｎｍｅｃｈａｎｉｓｍｏｆｒｅｄｔｉｄｅ ,ｍａｎｙｏｆｉｔｓａｓｐｅｃｔｓａｒｅｓｔｉｌｌｕｎｋｎｏｗｎｂｅｃａｕｓｅｒｅｄｔｉｄｅｉｓａｃｏｍｐｌｅｘｐｈｅｎｏｍｅｎｏｎａｎｄｔｈｅｅｃｏｌｏｇｉｃ…     

1∶50万地质图数据库的研建

全国1∶50万地质图数据库以国产软件MAPGIS作为基础平台，以各省（区）区域地质志附图和香港特别行政区1∶20万地质图为基础资料，采用现代地质学、地层学、岩石学等的新理论和新的表示方法，按岩石地层单位、花岗岩谱系单位，侵入岩按时代加岩性编制而成。该库由数字地质图库和地理底图库构成，可按任意空间范围、不同比例尺、不同投影方式和不同地质、地理内容等单项或多项任意检索；通过检索能自动形成具有图例、图外整饰的完整的地质图。     

胶西北新城金矿床硫同位素地球化学

新城金矿床是胶西北金矿集区中典型的破碎带蚀变岩型金矿床,其热液成矿作用可划分为四个阶段:黄铁矿-石英-绢云母阶段(I)、石英-黄铁矿阶段(II)、石英-多金属硫化物阶段(III)和石英-方解石阶段(IV),其中金主要赋存于II和III阶段的黄铁矿内。该矿床赋矿围岩为郭家岭岩体,岩性为石英二长岩和二长花岗岩,主要为胶东群变质基底经部分熔融形成。胶东群变粒岩硫同位素较为均一(δ34S值介于6.9‰~9.4‰,均值为8.0‰);郭家岭岩体的δ34S值介于6.0‰~16.0‰,均值为8.6‰,反映了其硫同位素组成总体上继承了变粒岩的硫同位素特征;长英质脉岩的δ34S值变化范围为0.8‰~8.5‰(均值为6.7‰),其中4件样品的δ34S值变化范围为7.4‰~8.5‰,反映了其硫主要源自于郭家岭岩体,变粒岩可能提供了部分硫源;而其中1件样品的δ34S值仅为0.8‰,符合岩浆硫来源特征,表明深部岩浆可能也提供了部分硫源。新城金矿床矿石中硫化物δ34S值变化范围较大(4.3‰~10.6‰,均值为8.3‰),表明矿石硫可能源于郭家岭岩体、变粒岩和长英质脉岩,最终主要来源于胶东群变质基底。I阶段黄铁矿颗粒较小(5~600μm),晶形主要为立方体,反映其处于温度较高(300~350℃)、成矿流体的过饱和度较低、低氧逸度和硫逸度、冷却快速、物质供应不足的成矿环境;黄铁矿δ34S值变化范围为8.4‰~10.6‰,均值为9.7‰,反映了矿石硫可能源自于δ34S值较高的郭家岭岩体和变粒岩。II和III阶段黄铁矿粒径变化较大(3μm~2.5mm),晶形主要为五角十二面体,反映其处于中-低温度(200~300℃)、成矿流体过饱和度高、高氧逸度和硫逸度、缓慢冷却同时物质供应充分的成矿环境。其中,II阶段黄铁矿δ34S值变化范围为7.7‰~9.7‰,均值为8.7‰,表明矿石硫源除郭家岭岩体和变粒岩外,δ34S值较低的长英质脉岩可能也提供了部分硫源。III阶段硫化物δ34S值变化范围较大(4.3‰~9.4‰,均值为7.1‰),闪锌矿-方铅矿硫同位素热力学平衡温度范围为180~282℃,氧逸度约为10-37.3~10-36.8,反映了矿石硫源自于郭家岭岩体、变粒岩和长英质脉岩,硫化物δ34S值变化范围较大可能是硫同位素分馏达到平衡的结果。IV阶段黄铁矿粒径最小(1~5μm),为晶形完好、表面光滑的立方体晶形,表明其处于较低温度(200℃),成矿流体的过饱和度较低、氧逸度和硫逸度较低、物质供应不足的成矿环境。     

安庆大别山区矿泉水化学特征及成因模式

安庆大别山区位于大别山造山带东侧,调查研究了该区饮用天然矿泉水6处(锶偏硅酸型2处,锶型1处,偏硅酸型3处),理疗天然矿泉水硅酸水3处.在分析矿泉水的水化学特征和环境同位素特征的基础上,进一步探讨其成因模式.冷泉(井)水(AH1、AH3、AH4、AH5、AH8和AH9)的主要离子为Ca2+、Mg2+、Na+和HCO3-...     

全国1∶100万地质图空间数据库建设进展

地质图空间数据库的建立对于基础地质图件的快速更新、信息资源共享以及基础地质科学研究具有重要意义。全国1∶100万地质图空间数据库(2005年版本)是以1∶50万地质图数据库为基础,充分吸收近年来1∶25万、1∶5万区域地质调查的新资料新成果,应用当代最先进的地质科学理论、现代信息技术和先进制图学理论,使用中国自主开发的MapGIS平台,结合ArcGIS、Access等软件进行研究建立的。数据库包括地理底图库、地质图库、元数据库、外挂岩石地层库,数据量约1.6GB。数据库管理系统具有较强的图幅检索、查询、自动成图、自动生成图例等功能。     

FY-2E晴空风矢同化对台风分析和预报的影响研究

为了研究“二阶差分法”反演的晴空区风矢同化在台风分析和预报中的作用,以1509号台风“灿鸿”和1211号台风“海葵”为例,首先利用WRF-3DVAR系统对晴空风矢进行同化,探讨了晴空风矢的引入对模式初始场的影响。然后利用WRF模式对两个个例分别进行48 h的预报试验。通过对比控制试验和同化试验,结果表明,同化晴空风矢资料能够对初始风场和位势高度场进行合理的调整,在台风周围引导气流的作用下,台风路径与实况更靠近,从而提高了台风路径的预报效果。除此之外,同化晴空风矢对台风强度以及风场预报也有一定的改善作用,还可更准确地预报出降水的落区及雨强,提高降水预报质量。因此,晴空风矢的引入,有利于改善模式的初始场,从而提高WRF模式对台风的预报能力。      

基于FY-2E分裂窗差值法的台风外围晴空区风场反演研究

利用MODTRAN4大气辐射传输模式,结合FY-2E光谱响应函数,分析热带大气廓线下FY-2E分裂窗亮温差(BTD=IR1-IR2)对下垫面温度(Ts)、大气水汽含量(WV)及温度廓线(TP)的敏感性。分析表明,在水汽充沛的热带洋面,水汽是引起BTD变化的主要因素,Ts和TP对分裂窗亮温差的影响不大。其中,大气水汽含量WV在2.0~5.0 g/cm2内变化0.42 g/cm2时,能够引起FY-2E分裂窗亮温差0.2 K的变化,达到卫星NEdT值。以2012年11号台风“海葵”、2013年7号台风“苏力”和12号台风“潭美”为例,利用分裂窗差值法反演了台风外围晴空区风场,结果表明,分裂窗差值法能够反演晴空区风场,并弥补云导风产品的不足,其反演结果与NCEP再分析资料低层(800~900 hPa)风场具有很好的一致性,且精度与云导风的风矢测量精度相当。