基于法方程病态的谱修正迭代算法的探讨

本文通过分析法方程病态问题产生的原因,结合谱修正迭代算法原理,探讨法方程病态对参数估值的影响,对谱修正迭代的改进算法的适用范围进行扩展,从理论上证明修正因子r取大于零的实数时,迭代改进算法的逆矩阵二范数值随迭代次数的增加而趋近于零,从而得到平差参数接近真值的估值。经过分析发现,谱修正迭代改进算法的迭代速度主要取决于修正因子r和迭代初值的取值。本文从理论和实例证明了迭代速度与修正因子取值保持线性的变化规律,同时用实例证明了不同初值对迭代速度的影响。     

CBERS-02B卫星数据在涪陵区土地利用分类中的应用

为推进国产卫星数据的应用,本文以CBERS-02B卫星的CCD数据和HR数据为数据源,以重庆市涪陵区为实验区,将DEM、NDVI数据作为补充,采用分层分类的方法,对实验区土地利用/土地覆盖现状信息进行提取,分析CBERS-02B数据的特点及在应用中存在的问题。研究结果表明,CBERS-02B数据可满足复杂地区土地利用宏观监测的需要。     

应用稀疏非负矩阵分解聚类实现高光谱影像波段的优化选择

针对高光谱影像数据高维性、高度相关性和冗余性等特点,提出应用稀疏非负矩阵分解聚类实现高光谱影像波段的优化选择。通过稀疏非负矩阵分解方法对高光谱影像进行稀疏化表示,同时顾及其可聚类的特性,在保留所选波段物理意义的基础上,得到波段选择后的高光谱影像降维数据。通过该方法对PHI-3高光谱影像进行波段选择的试验分析,应用聚类特征有效性分析波段聚类结果,并采用波段子集的信息量、相关性和可分性3类评价指标来验证方法的效果。最终,从运行效率和分类精度两方面证明了基于无监督聚类的稀疏非负矩阵分解对高光谱影像的波段选择的实用性。     

耦合遥感观测和元胞自动机的城市扩张模拟

在传统元胞自动机(CA)模型中,静态的模型参数和模型误差不能释放是影响城市扩张模拟效果的两个重要原因。文中引入集合卡尔曼滤波方法到CA模型中,提出了基于联合状态矩阵的地理元胞自动机。该模型在模拟过程中可以通过同化遥感观测数据,动态地调整模型参数和纠正模拟结果,使模型参数能够反映转换规则的时空变化,同时也能较好地释放积累的模型误差。将模型应用于东莞市的城市扩张模拟中,实验结果表明,模型能够准确地调整模型参数使之符合城市发展模式,同时也能有效地控制模型误差,其模拟的空间格局与真实情况吻合。     

稀疏存储的显式有限元三角网格地震波数值模拟及其PML吸收边界条件

有限元法是复杂介质地震模拟的有力工具,它能比较客观地反映地震波的传播,比较细致地再现地震图像.但是,为了获得较精确的结果,有限元法模拟地震波的传播需要的网格点数多,具有计算量大和消耗内存多的缺点.针对上述缺点,本文对刚度矩阵采用压缩存储行(CSR)格式,以减少计算量并节省内存;采用集中质量矩阵得到对角的质量矩阵以提高有限元法(显式有限元)的计算效率;时间离散采用保能量的Newmark算法以提高有限元法的计算精度;采用变分形式(弱形式)的PML吸收边界条件对人工截断边界进行处理.通过与高精度的数值方法--谱元法的数值试验的对比表明,上述方法的引入可使有限元法在计算精度和计算效率方面均可取得比较显著的改进.为了获得相当的计算精度,相比于7阶谱元法,显式有限元法需要更精细的网格.然而,显式有限元法的计算速度比前者快近2倍,而内存需求仅为谱元法的1/4~1/6.     

Effective Correction of Mass Bias for Rhenium Measurements by MC-ICP-MS

The geochemistry of Re-Os and the recent use of Re as a non-traditional stable isotope both need accurate and precise quantification of ¹⁸⁷Re/¹⁸⁵Re ratios. This paper reports rhenium isotopic data obtained from the analysis of a standard solution and geological samples by MC-ICP-MS. We show that measured isotopic ratios are modified by matrix effects that cannot be accounted for by the standard solution bracketing technique. The bias resulting from measurements on a spiked (¹⁸⁵Re-enriched) sample is shown to alter the apparent Re concentration by several percent. When spiking samples and calibrators with tungsten, simultaneous measurement of tungsten and rhenium isotopes compensates for the matrix-induced modification of mass bias. Rhenium and tungsten are shown to have different fractionation factors. This may be due to the fact that the two elements fractionate in a different but systematic way, or that the reference isotopic ratios used for elemental Re and W are incoherent with one another. The consistency of fractionation through time can be used to obtain an empirical relationship between W and Re measured ratios from a standard solution to obtain a sample's fractionation-corrected ¹⁸⁷Re/¹⁸⁵Re spiked ratio on samples containing pg g⁻¹ levels of Re, even if some matrix capable of affecting mass bias remains in the final solution.     

An Improved Atmospheric Vector Radiative Transfer Model Incorporating Rough Ocean Boundaries

The radiative transfer model (RT3), a vector radiative transfer (VRT) scheme in a plane-parallel atmosphere, was bounded by a rough ocean surface in this study. The boundary problem was solved using a Fourier series decomposition of the radiation field as a function of the azimuth. For the case of a rough ocean surface, the decomposition was obtained by developing both the Fresnel reflection matrix and the probability distribution of the water facet orientation as Fourier series. The effect of shadowing by ocean surface waves was also considered in the boundary condition. The VRT model can compute the intensity and degree of polarization of the light at the top of the atmosphere (TOA), the ocean surface, and any level of the atmosphere in the ocean-atmosphere system. The results obtained by our model are in good agreement with those computed by Ahmad’s model. The simulated results showed that the shadow effects of wave facets on the intensity and the degree of polarization are negligible except at the ocean surface near the grazing angle, possibly because we did not consider the effect of white caps.     

矩阵方程A^TXA=D扰动分析

利用矩阵分块与矩阵范数的性质,研究矩阵方程A^TXA=D,获得该方程的扰动界,这些结果可用于模型修正中的数值计算。     

热压高磷铁基胎体性能及其金刚石钻头研究

针对在硬而致密岩石中钻进时效低和钻头使用寿命短以及传统金刚石钻头胎体材料价格高等众多难题,采用混料回归试验设计方法,进行了替代碳化钨基金刚石钻头胎体材料的新型铁基胎体配方的试验研究;分析了热压钻头的金刚石出刃与岩石研磨性等岩性之间的内在联系,认为钻头的胎体成分及其性能是关键。因此,从热压金刚石钻头的胎体成分与性能研究入手,经过反复的试验研究与优化组合,研究出热压高磷铁基孕镶金刚石钻头的胎体材料,并试制出金刚石钻头。在室内对可钻性Ⅸ级岩石进行钻进的试验表明：钻进时效平均达到1．91m,钻头平均工作寿命达到60．26m。试验结果表明,高磷铁基眙体是一种性能良好的金刚石钻头胎体材料,且该类型胎体的热压金刚石钻头是一种成本明显降低、具有广谱性能的金刚石钻头。     

坚硬致密岩层用绳索取心钻头的研制与应用

为了克服绳索取心钻头钻进坚硬致密岩层出现的进尺缓慢、钻头寿命短的难题,着重从钻头胎体配方、金刚石参数、钻头水口等方面进行优化设计、试验并生产应用,使研制的绳索取心钻头(S77)具备良好的适应性,钻进坚硬～硬、中～弱研磨性岩层克服了“打滑”现象并延长了钻头的寿命,取得了一定的经济效益。