最小二乘定权方法对流动重力计算结果的影响

在处理流动重力观测数据时,采用最小二乘方法解算由权系数和观测值组成的超定方程组,得到最优权值,并利用最优权计算重力值。结果表明,等权方法计算得到的结果点值精度不高,但包含了较全的重力异常地震分析预报信息;不等权方法计算得到的结果点值精度较高,但削弱了有用的重力异常地震分析预报信息;最小二乘方法计算得到的结果不但有较高的点值精度,对M_S4.0及以上地震的分析预报效果也最优。     

扩展共线方程结合车载组合全景相机影像区域网平差

利用扩展共线方程实现组合相机单片影像的空中三角测量,通过该方法来提高全景影像的定位精度。由于组合全景相机中各子相机刚性固定,相机间的相对位姿参数可以精确标定,因此利用扩展共线方程模型,可以将组合全景相机的构像单元从单个子相机影像扩展为整个组合相机获取的多张影像。与已有的球模型方法,包括球理想模型和球严密模型方法相比,本文所用的扩展共线方程方法使用各子相机影像像点作为观测值,避免了球模型方法中由各子相机摄影中心不重合导致的球模型投影误差和摄站内重叠区域融合误差。试验结果表明:与球模型方法相比,本文方法无需模型投影和摄站内重叠区域融合等繁琐的处理过程,避免了由其引起的精度损失,保证了空三网稳定性的同时提高了定位精度。     

电力系统元件非线性微分-代数子系统模型的逆系统控制:一种新算法

针对电力系统元件非线性微分-代数子系统模型,本文提出一种新算法研究其逆系统控制问题.所提出的新算法不需要对控制输出及其高阶导数做复杂的变换,具有更好的应用性.本文的逆系统控制方法主要分为两步：第1步,利用所提出的新算法来判断被控元件的可逆性,若可逆,则基于状态反馈与动态补偿,构造出元件的α阶积分右逆系统,实现复合系统的线性化和解耦;第2步,利用线性控制的理论和方法设计闭环控制器,使得元件被控对象满足期望的性能指标.最后按照本文所提出的方法,研究了多机电力系统的分散非线性汽门控制问题.仿真结果验证了本文所提方法的有效性.     

斜压海洋动力学的一种三维数值模式──Ⅱ.温度、盐度和垂直涡动粘性系数的计算

根据海水温度和盐度平流扩散方程给出一种数值计算方案并采用混合长度理论给出垂直涡动粘性系数的计算方法。对于温度和盐度方程,其平流过程采用了Lax－Wendroff格式,水平扩散采用显格式,垂直扩散采用隐格式。时间步长主要受平流过程的Courant-Friedrichs-Lewy条件限制。垂直涡动粘性系数计算依据Prandtl混合长度理论,并考虑了海水层化的抑制作用,因而其数值与流场及密度场结构有关。温度、盐度及垂直涡动粘性系数的计算与动力方程中内模态的计算同步进行。应用本模式模拟渤、黄、东海由潮流、密度流和风海流迭加而成的综合海流,得出了良好的结果。     

官厅水库水质模型研究

根据官厅水库狭长、水浅和水文、水力、水质迁移扩散等特点与河流相似的事实,采用一维St.Venant方程组模拟其水流状况;采用完全混合反应器概念模拟其水质变化,两者的藕合就构成了官厅水库的水质模型。该模型在模拟官厅水库NH₃-N浓度时空变化时,精度令人满意,并已用于官厅水库的水质预测和水质管理。     

Refraction Effect in Geothermal Heat Flow Due to a 3-D Prism in a Two-Layered Earth

A mathematical model of the disturbance of the stationary geothermal field due to a three-dimensional perturbing body embedded near the surface of a two layered earth. The theoretical analysis is based on the generalized theory of the double-layer potential, similar to the boundary integral method likely to direct current geoelectricity problems. Special attention is paid to cases in which a pan of the body's surface is tangential to the planar boundaries of the substratum or of the surface of the earth The numerical calculations are performed for a three-dimensional prism using methods similar to the boundary element method (BEM). Numerous graphs are shown for the disturbance of the temperature field and of the heat flow on the surface of the earth and along vertical profiles (boreholes).     

Some aspects of head‐variance evaluation

We compare two methods of evaluating head covariance for two‐dimensional steady‐state flow in mildly heterogeneous bounded rectangular aquifers. The quasi‐analytical approach, widely used in stochastic subsurface hydrology, is based on the Green's function representation, and involves numerical four‐fold integration. We compare this approach with a numerical solution of the two‐dimensional boundary‐value problem for head covariance. We show that the finite differences integration of this problem is computationally less expensive than numerical four‐fold integration of slowly‐convergent infinite series. This revised version was published online in July 2006 with corrections to the Cover Date.     

一种近岸区波浪破碎模型

从波浪破碎的能量关系入手,以紊流能量方程为基础,考虑破碎区内单个波在不同破碎阶段所提供的紊动能量强度的变化过程,提出了一种波浪破碎模式.通过将这一模型引入Boussinesq方程中,初步建立了一种近岸区波浪变形数学模型,并用波浪水槽实验资料对模型模拟波高和平均水位的情况进行了初步验证,得到了良好的结果.     

港湾中海啸流的模拟及其时空演变特征对输入条件不确定性的响应分析

一直以来,海啸波特征作为表征海啸潜在破坏性的参数指标得到了广泛应用,特别是针对近场极端海啸事件造成的灾害来说,这种表征具有较好的适用性.然而总结分析历史海啸事件造成的损失发现：在远场近岸及港湾系统中,海啸诱导的强流却是造成损失的主要原因.陆架或港湾振荡导致海啸波幅快速升降诱发强流,可能促使港工设施受到威胁及损害,进而对海啸预警服务及海事应急管理提出了新的挑战.因此,全面理解与评估海啸在港湾中诱发的灾害特征,探索港湾中海啸流的数值模拟方法,发展针对港湾尺度的海啸预警服务指导产品尤为迫切.受限于海啸流验证数据的缺乏及准确模拟海啸流技术方法的诸多不确定性,大部分海啸数值模拟研究工作主要是针对水位特征的研究及验证,可能导致对港湾中海啸灾害危险性认识的曲解与低估.本研究基于非线性浅水方程,针对夏威夷群岛三个典型港湾建立了精细化海啸数值模型（空间分辨率达到10 m）,并联合有限断层破裂模型计算分析了日本东北地震海啸在三个港湾及其邻近区域的海啸特征,波、流计算结果与实测结果吻合较好,精细化的海啸港湾模型模拟结果可信.模拟发现港湾中较小的波幅,同样可以产生强流.综合分析日本东北地震海啸波、流特征对输入条件不确定性的响应结果发现：港湾中海啸波-流能量的空间分布特征差异较大,这与港湾系统中海啸波的驻波特性相关;相比海啸波幅空间特征,海啸流特征具有更强的空间敏感性;海啸流时空分布特征对输入条件的不确定性响应比海啸波幅对这些不确定性的响应更强,海啸流的模拟与预报更有挑战性;不确定性对海啸流计算精度的影响会进一步传导放大港湾海啸流危险性的评估及对港工设施产生的应力作用的误差,合理的输入条件对海啸流的精确模拟至关重要.最后,希望通过本文的研究可以从海啸波-流特征角度更加全面认识近岸海啸灾害特征,拓展海啸预警服务的广度与深度,从而为灾害应急管理部门提供更加科学合理的辅助决策产品.     

An approach to analysing plot scale infiltration and runoff responses to rainfall of fluctuating intensity

Simulated rainfall of fluctuating intensity was applied to runoff plots on bare dryland soils in order to explore a new method for analysing the non‐steady‐state responses of infiltration and overland flow. The rainfall events all averaged 10 mm/h but included intensity bursts of up to 70 mm/h and lasting 5–15 min, as well as periods of low intensity and intermittency of up to 25 min. Results were compared with traditional steady‐state estimates of infiltrability made under simulated rainfall sustained at a fixed intensity of 10 mm/h. Mean event infiltration rate averaged 13.6% higher under fluctuating intensities, while runoff ratios averaged only 63% of those seen under constant intensity. In order to understand the changing soil infiltrability, up to three affine Horton infiltration equations were fitted to segments of each experiment. All equations had the same final infiltrability f_c, but adjusted values for coefficients f₀ (initial infiltrability) and K_f (exponential decay constant) were fitted for periods of rainfall that followed significant hiatuses in rainfall, during which subsurface redistribution allowed near‐surface soil suction to recover. According to the fitted Horton equations, soil infiltrability recovered by up 10–24 mm/h during intra‐event rainfall hiatuses of 15 to 20‐min duration, contributing to higher overall event infiltration rates and to reduced runoff ratios. The recovery of infiltrability also reduced the size of runoff peaks following periods of low intensity rainfall, compared with the predictions based on single Horton infiltration equations, and in some cases, no runoff at all was recorded from late intensity peaks. The principal finding of this study is that, using a set of affine equations, the intra‐event time variation of soil infiltrability can be tracked through multiple intensity bursts and hiatuses, despite the lack of steady‐state conditions. Copyright © 2016 John Wiley & Sons, Ltd.