改进的最小二乘法在水文分析计算中的应用

在水分析计算中，经常涉及到变量之间的线性或非线性拟合，而在拟合各种特性曲线时，通常应用以实测资料与拟合曲线间的误差平方和最小作为目标函数的方法——最小二乘法，但这种方法忽视了所有实测点应与拟合曲线间的相对误差尽量不超过某一百分比的原则，为了达到上述要求，提出了非线性的加权最小二乘法及线性相关方程的最小距离平方和法，探讨改进了传统的最小二乘法达到优化的效果。最小距离平方和法与常用的图解法相比，本法所得成果较为客观；与传统的单方向(x或y方向)最小二乘回归法相比，所求线性方程不会因坐标系的选取而改变。最后应用算例进行了初步讨论。     

油气化探综合异常模糊划分动态聚类分级评价

通过选取描述油气化探综合异常的各种特征参量，借助模糊数学工具，利用最优聚类中心所具有的代表特征，建立分级指标，对油气化探综合进行分级评价。鄂尔多斯盆地油气化探综合异常的评价实例结果表明，利用油气化探综合异常各种特征参量进行动态聚类分级评价是可行的。     

任意面积储量计算方法研究

这里介绍了一种任意面积储量计算方法，该方法的核心技术是充分利用现有的石油勘探、开发数据库资料，由计算机来自动确定任意含油面积，进而通过容积法来计算石油储量。该方法解决了用求积仪来求取面积时的大量手工操作，在储量计算过程中自动化程度高，具有方便、灵活、实用的特点。     

浅水方程数值计算方法的研究

求解以水位为变量的连续方程,并根据Navier-Stokes方程压力修正算法的基本思想,建立了浅水方程的水位修正算法,放宽了对离散时间步长的限制.通过对离散方程系数矩阵的重新构造,建立了高分辨率有限元格式,该格式既具有较高的离散精度又避免了数值解的伪振荡.对动量方程的阻力项做负坡线性化处理,提高了露滩计算的稳定性.数值模拟结果与解析解吻合良好,表明所建立的数值计算方法是正确的和可靠的.     

地震参数计算与数据传输程序设计

设计了一种地震参数计算与数据传输程序，能够通过交互方式自动完成mb、mB、MS7、MS、ML五个震级及其他相关地震参数的计算，自动生成大震速报文件，直接连接PCU速报。不仅适用于数字地震仪同时也可以用于模拟仪器(如DD—1，DK—1)解决了SSB大震速报软件数据不易改动、容易出错、速度慢等问题。同时也完善了SSDP软件参数计算功能，提高了数字测震台的效能和速报水平，充分发挥了数字地震仪的优越性。     

Compression of the North Hemisphere derived from space geodesy

Introduction With the development of global tectonics and overall detections for global tectonics with multi-geophysical methods, ones can roundly study on the geological tectonics of sampling and magnetic stripe image, so as to summarize and interpret the geometrical and kinematical charac-teristics for the distribution of the ocean and the land, and spreading state of the global tectonics in a global scale. From a comprehensive view, the South and North hemispheres are clearly unsym-metrical…     

The inversion of anelastic coefficient,source parameters and site respond using genetic algorithm

It gradually becomes a common work using large seismic wave data to obtain source parameters, such as seismic moment, break radius, stress drop, with completingof digital seismic network in China (Hough, et al, 1999; Bindi, et al, 2001). These parameters are useful on earthquake prediction and seismic hazard analysis.Although the computation methods of source parameters are simple in principle and the many research works have been done, it is not easy to obtain the parameters accurately. There are two factors affecting the stability of computation results. The first one is the effect of spread path and site respond on signal. According to the research results, there are different geometrical spreading coefficients on different epicenter distance. The better method is to introduce trilinear geometrical spreading model (Atkinson, Mereu, 1992; Atkinson, Boore, 1995; WONG, et al, 2002). In addition, traditional site respond is estimated by comparing with rock station, such as linear inversion method (Andrews, 1982), but the comparative estimation will introduce some errors when selecting different stations. Some recent research results show that site respond is not flat for rock station (Moya, et al, 2000; ZHANG,. et al, 2001; JIN, et al, 2000; Dutta, et al, 2001). The second factor is to obtain low-frequency level and corner frequency fromdisplacement spectrum. Because the source spectrum model is nonlinear function,these values are obtained by eye. The subjectivity is strong. The small change of corner frequency will affect significantly the result of stress drop.     

观测有效性的度量方法

基于Hilbert空间理论 ,通过对观测方程的系数阵实施格拉姆_施密特 (Gram_Schmidt,G_S)正交化过程 ,分析观测的两个特性 ,从一个新的角度论证了参数估计的必要观测问题 ,以及多余观测与必要观测的关系 ;用欧氏范数‖·‖2 和Frobenius‖·‖ F 范数作为基本度量工具对观测提供参数估计的信息量进行度量分析 ,研究观测的有效性和参数估计的可靠性。采用观测空间分析方法研究了测绘学科中必要观测、多余观测等基本概念 ;提出了观测信息量、观测有效性的概念和相应的度量方法 ;从观测的空间结构出发分析参数估计系统的可靠性 ,从一个侧面扩展和丰富了测绘学中的可靠性理论 ;丰富了研究测绘学中一些基本问题的分析方法 ,也为信息学中的数据挖掘问题提供了一种分析思路     

DLT程序设计与应用

在近景摄影测量直接线性变换原理的基础上，加入附加条件建立附加有平行条件的直接线性变换函数模型，探讨了在解算过程中用数据探测法剔除粗差，用验后补偿的方法补偿系统误差的实现。鉴于非量测像机在近景摄影洲量中应用导致直接线性变换解法的重要性和实用性及数据量大、解算过程过于复杂的情况。用VC 编制了直接线性变换的解算程序，对其关键的方法和用法进行了说明。     

A Linear Programming Approach to Determine the Normative Composition of Sedimentary Rocks

The analysis of modal and normative composition of sedimentary rocks is widely used for studying their sources and tectonic settings. The normative calculation of the mineral composition of rocks in this study is formulated as a linear programming problem and is solved by means of the simplex method. This enables both simultaneous and successive subtraction of a set of basic minerals from a rock sample represented by its chemical composition {SiO₂...LOI}. Such an approach provides a more exact calculation of the contribution of basic minerals in the rock. This mathematical approach is used to study two representative sets of sandstones and fine-grained rocks from a Meso- to Neoproterozoic marginal basin of southeastern Siberia (Uchur–Maya region, Yakutia) and a Pennsylvanian-Lower Permian uplifted continental block in Colorado, USA. The calculated normative mineral compositions of the Siberian sandstones are consistent with the observed modal compositions. These sandstones vary from K- Feldspar rich arkoses at the base of the sequence (the Uchur Group, lower Riphean) to quartz arenites or lithic sandstones and wacke in transgressive successions of the middle-upper Riphean. Arkoses and quartz arenites are dominant in Meso- to Neoproterozoic Siberia. These samples represent craton interior uplifted basement and quartzose, recycled orogen provenance of a stable craton in Rodinia. There are higher but consistent discrepancies between the calculated and observed compositions for the Pennsylvanian to Lower-Permian arkoses and quartz arenites (Sangre de Cristo, Belden, and Maroon Formations). The differences between the predicted and observed mineralogy may be due to uncertainties in the modes in the matrix and cement of the sandstones. This normative program should supplement modal calculations and provide better genetic constructions, especially in case of matrix-rich sandstones.