丹麦法稳健估计在隧道椭圆拟合中的应用

隧道断面在设计时为一标准圆形,受到施工及外部环境的影响变形为偏心率很小的椭圆。减小粗差影响、精确确定椭圆的几何参数是其研究的重点。通过对上海某一隧道断面的测量试验,得出用丹麦法稳健估计可以有效减小粗差的影响,而且拟合值可以很好地逼近无粗差情况下的最小二乘拟合值。采用Matlab 语言实现算法,并详细介绍了用二项式变换求解椭圆参数的方法。     

Statistical evaluation of tectonomagmatic discrimination diagrams for granitic rocks and proposal of new discriminant-function-based multi-dimensional diagrams for acid rocks

The four tectonic discrimination diagrams of Pearce et al. [Journal of Petrology, v. 25, p. 956–983] for granitic rocks were first evaluated using the literature cited by these authors as well as from our new database. The first diagram (Y?Nb) cannot discriminate volcanic-arc and collision settings. Both Y?Nb and Yb?Ta diagrams have an overlapping field for within-plate and ocean-ridge granitoids. The remaining two diagrams (Y?+?Nb?Rb and Yb?+?Ta?Rb) use a mobile element (Rb) in their y-axis. Although these diagrams successfully discriminate volcanic-arc and within-plate granites, they perform less well for collision tectonics. Besides, felsic or acid rocks are scarce in ocean-ridge settings, which limits the usefulness of these diagrams for this geological environment. Therefore, using an extensive database, we proposed a set of five new discriminant-function-based multi-dimensional diagrams for acid magmas from four tectonic settings (island arc, continental arc, continental rift, and collision). The very similar tectonic settings of island and continental arcs are discriminated for the first time. These diagrams are based on correct statistical treatment of compositional data, because they use natural logarithm transformation of major-element ratios and linear discriminant analysis (LDA). The use of discordant outlier-free samples prior to LDA improved the success rates by about 3–5%. Success rates of these diagrams as inferred from a testing set were between 76% and 88% for island arc, 60% and 92% for continental arc, and 72% and 84% for both continental rift and collision settings. Finally, application of these new diagrams to case studies not compiled in our initial database used for constructing these diagrams provided the following results: a collision setting for the Himalayas at about 30 Ma; an island arc setting for Quaternary acid rocks from geothermal boreholes in El Salvador; an island- or continental-arc setting for northern Italy at 35–52 Ma; a continental-arc setting for the Italy–Austria border at about 30 Ma; either a rift or a collision setting for northern Nigeria at about 164 Ma; a collision setting for central Nigeria at about 144 Ma and for the Cretaceous Masirah ophiolites of Oman; and an island arc setting for the Cretaceous Semail ophiolites of Oman. In spite of the relative mobility of major elements, these applications suggest utility of the new discrimination diagrams for all four tectonic settings.     

Identification of Archaean plate tectonic processes from multidimensional discrimination diagrams and probability calculations

We applied our group's previously published multidimensional diagrams in 2006–2012 and corresponding probability estimates in 2011–2012 to geochemical data for Archaean rocks compiled from cratons in Australia, South Africa, Brazil, Canada, and India. Tectonic processes similar to present-day plate tectonics evidently were active at least since the Palaeoarchaean (?3570 Ma). This seems to be true in spite of a presumably hotter Earth at that time. For the eastern part of the Pilbara craton (Australia), a Palaeoarchaean (3570–3450 Ma) and Mesoarchaean (2900 Ma) continental arc setting apparently evolved to a collision (Col) setting during the Neoarchaean (2600 Ma). We infer an island arc (IA) environment for Kambalda (Australia) during the Neoarchaean (2700 Ma). For the Barberton belt (South Africa), a transition from a mid-ocean ridge regime during the older part of the Palaeoarchaean (3470 Ma) to an IA setting during the younger part (3300–3260 Ma) is likely. We inferred an arc environment for the São Francisco craton (Brazil) and the Rio Maria terrane (Brazil) during the Mesoarchaean (3085–2983 Ma and 2870 Ma, respectively), whereas a within-plate setting is clearly indicated for the Carajás metallogenic province (Brazil) during the Neoarchaean (2740–2700 Ma). We also recognize an IA regime for the Mesoarchaean (3000 Ma) North Caribou and Neoarchaean (2700 Ma) Abitibi greenstone belts (Canada), and for the Gadwal greenstone belt (India) during the Neoarchaean (2700–2500 Ma). A Col setting was inferred for the Archaean sanukitoid suite (Canada) and the Kolar suture zone (India) during the Neoarchaean (2700–2660 Ma and 2630–2520 Ma, respectively).     

裸眼系列完井方式下叠式分支水平井产能预测

综合运用保角变换、镜像反映愿理,势的叠加原理以及等值渗流阻力等数学和渗流力学分析方法,推导出裸眼系列完井方式下叠武分支水平井的产能解析公式,当分支数目为1时,所推导的解析公式退化为著名的Borisov公式.分析了影响其产能的几个重要因素,从而为该类型分支井的完井方案设计提供参考依据.经研究发现,完井筛管挡砂精度越高,相比水平井,叠式分支水平井产能损失相对越少;叠武分支水平井的产能受其在上下封闭边界油藏中垂向位置的影响较小;相比微元段模型,解析模型预测的产能偏高,引入校正系数后可减小偏差.     

酒泉盆地营尔凹陷异常高压形成机理

酒泉盆地营尔凹陷下白垩统砂泥岩普遍发育超压。通过对地层压实特征、黏土矿物转化及有机质演化特征的分析,揭示了凹陷超压形成的机理。结果表明:①营尔凹陷纵向上发育常压带、浅层超压带、过渡带和深层超压带4个压力带,超压带内具有明显的泥岩密度降低、电阻率减小和声波时差增大的欠压实特征,表明不均衡压实作用是形成该区超压的主要原因;②凹陷内发育3个蒙脱石向伊利石转化带,第一迅速转化带对应于浅层超压带,第二迅速转化带对应于深层超压带,反映黏土矿物转化与凹陷超压有一定的联系;③深层有机质演化程度较高(镜质体反射率Ro>0.6%),与浅层超压带相比,生烃作用对深层超压带的超压贡献更大。由此认为,营尔凹陷的下白垩统地层超压是不均衡压实作用、黏土矿物脱水作用和生烃作用共同影响的结果,并且隆起带和洼陷带的成因机理有所差异。     

基于移不变全方向角提升的遥感图像降噪

在分析遥感图像结构特征及其与噪声之间主要区别的基础上,利用图像信号的方向信息,提出基于移不变全方向角提升小波（TI—OL）抑制遥感图像噪声的方法。该方法在方向提升小波变换的基础上并利用循环平移,Gabor小波滤波器和图像旋转技术改进了方向提升小波在图像去噪过程中存在的三个弊端：缺乏移不变性质,图像局部方向信息判方法断缺乏噪声鲁棒性和变换方向分布有限。消除去噪结果中的吉布斯效应,提高图像方向信息判断的准确性并保证图像纹理方向始终落在方向提升能最优表示的方向区间内。试验结果证明所提方法在处理遥感图像的过程中能在去噪的同时保留图像的细节和边缘信息,对遥感图像中的边缘信息如道路和桥梁有较好的刻画性能,较传统方法去噪性能（PSNR）和主观视觉效果（SSIM）均有较大提高。     

基于小波-HIS的遥感影像融合方法技术研究

小波是建立在泛函分析、调和分析、傅立叶分析基础上的时频原子,在时域、频域同时具有良好的局部化特性和多分辨率特性。遥感影像融合是遥感影像处理中的重要环节,是人们处理多源遥感影像数据的主要手段之一,而融合方法的选择对融合结果来说至关重要。将传统的融合方法与小波结合已经成为研究的热点,本文将小波变换与HIS融合方法相结合,并...     

高斯曲线优化能见度与气溶胶光学厚度转换模型

大气气溶胶是影响对地观测定量精度的最主要不确定性因素。随着定量遥感的发展,对气溶胶光学厚度数据的精度提出了更高要求。在广泛应用的基于辐射传输模型大气校正研究中,需要输入气溶胶光学厚度等关键参数,但与对地观测影像数据同时相的气溶胶光学厚度获取较难,而水平能见度作为表征气溶胶光学特性的间接参数可通过广泛分布的气象台站获得,可将能见度转换得到的气溶胶光学厚度数据作为同时相数据输入传输模型进行大气校正计算。本文以实测的能见度和气溶胶光学厚度数据为基础,通过拟合气溶胶标高其随时间的变化对Peterson模型进行了修正。对修正后的模型进行精度验证得到RMSE为0.254,结果表明优化的模型对精度有较大提升。     

加权总体最小二乘在三维基准转换中的应用

对比研究加权总体最小二乘(weighted total least-squares,WTLS)方法和混合最小二乘(LS-TLS)方法、最小二乘(least-squares,LS)方法在三维空间小角度直角坐标转换中的适用性。在两套坐标系下坐标测量值均存在误差时,用WTLS方法不但可以对观测向量y和系数矩阵A同时修改、将坐标先验精度引入平差计算,而且引入的权阵PA对系数阵A起到固定常数元素而只修改必要数据元素的作用,以得到更合适的参数解。     

基于改进SQPM算法的三维直角坐标转换

给出了七参数三维直角坐标转换的严密表达式,将改进SQPM算法用于大旋转角的三维直角坐标转换。实验表明该方法能在大范围内实现收敛且结果稳定、精度较高,适合于大旋转角的三维直角坐标转换。