GNSS全自动解算软件研发及精度评估

根据载波相位差分定位技术,利用待定点周边稳定的CORS站网形成的区域误差改正模型,对待定控制点观测数据的噪声进行自动剔除,对电离层、对流层和多路径进行差分改正,研发GNSS全自动解算软件,并利用广西卫星定位基准站与待定控制点的GNSS观测数据对软件进行测试。结果表明,GNSS全自动解算软件具有较好的内外符合精度,其解算的三维坐标精度满足D级以下等级控制点的要求,可为工程测量应用提供新的自动化解算方案。     

GNSS接收机DCB短期时变特征分析和估计方法

提出一种使用非差非组合精密单点定位(PPP)估计和分析接收机DCB短时时变特征的方法。首先利用非差非组合PPP得到包含接收机DCB的重构电离层参数估值;然后通过IGS电离层GIMs格网模型内插剥离各历元站星斜向电离层距离延迟;最后通过最小二乘约束得到各历元接收机DCB解。由于格网本身精度(2~8 TECU)和插值精度限制,解算出来的接收机DCB并不能真实反映其短期时变特征。为此,提出利用站间单差或者历元间差分的方法还原其真实的变化态势。实验结果表明,所提出的方法能够正确估计接收机DCB,并能真实还原其短期时变特征,具有良好的适用性。     

近50年石家庄地-气温差变化特征

利用石家庄17个国家气象站1972—2021年逐日地面气温、0 cm地温资料,分析了石家庄地-气温差的变化特征,结果表明：(1)石家庄地-气温差从1月开始逐渐增加,5月达到最大值5.0℃,然后开始减小,12月达到最小值-0.8℃;地-气温差在11月到次年1月为负值;春、夏、秋季均为正值,夏季最大,春季大于秋季,冬季以负值为主;(2)石家庄多年平均地-气温差在1.6～2.6℃之间,平均为2.1℃;整体上东部大于西部。(3)近50年石家庄年平均地-气温差呈显著的减小趋势,变化速率为-0.14℃/10a;夏、秋、冬三季的减小趋势均非常显著,夏季的减小趋势最强;石家庄市区和近郊站点年和四季地-气温差的减少趋势更显著。本文结论对科学认识石家庄城市生态环境的变化具有参考意义。     

三种验潮方法水位观测性能比较与统计分析

水尺、压力式验潮仪与验潮井验潮仪是海道测量水位控制和水位观测的主要手段。通过6日6个时间段的同步比对试验,计算了三种方法获取水位观测数据之间的潮时差和潮高差,分析了三种观测手段的水位观测性能以及产生差值的可能原因,提出了使用方法的建议。     

The determination of potential difference by the joint application of measured and synthetical gravity data: a case study in Hungary

In an elementary approach every geometrical height difference between the staff points of a levelling line should have a corresponding average g value for the determination of potential difference in the Earth’s gravity field. In practice this condition requires as many gravity data as the number of staff points if linear variation of g is assumed between them. Because of the expensive fieldwork, the necessary data should be supplied from different sources. This study proposes an alternative solution, which is proved at a test bed located in the Mecsek Mountains, Southwest Hungary, where a detailed gravity survey, as dense as the staff point density (~1 point/34 m), is available along a 4.3-km-long levelling line. In the first part of the paper the effect of point density of gravity data on the accuracy of potential difference is investigated. The average g value is simply derived from two neighbouring g measurements along the levelling line, which are incrementally decimated in the consecutive turns of processing. The results show that the error of the potential difference between the endpoints of the line exceeds 0.1 mm in terms of length unit if the sampling distance is greater than 2 km. Thereafter, a suitable method for the densification of the decimated g measurements is provided. It is based on forward gravity modelling utilising a high-resolution digital terrain model, the normal gravity and the complete Bouguer anomalies. The test shows that the error is only in the order of 10⁻³mm even if the sampling distance of g measurements is 4 km. As a component of the error sources of levelling, the ambiguity of the levelled height difference which is the Euclidean distance between the inclined equipotential surfaces is also investigated. Although its effect accumulated along the test line is almost zero, it reaches 0.15 mm in a 1-km-long intermediate section of the line.     

The phase difference effects on 3-D structure of wave pressure acting on a composite breakwater

In designing the coastal structures, the accurate estimation of the wave forces on them is of great importance. In this paper, the influences of the phase difference on wave pressure acting on a composite breakwater installed in the three-dimensional (3-D) wave field are studied numerically. We extend the earlier model [Hur, D.S., Mizutani, N., 2003. Coastal Engineering 47, 329–345] to simulate 3-D wave fields by introducing 3-D Navier–Stokes solver with the Smagorinsky's sub-grid scale (SGS) model. For the validation of the model, the wave field around a 3-D asymmetrical structure installed on a submerged breakwater, in which the complex wave deformations generate, is simulated, and the numerical solutions are compared to the experimental data reported by Hur, Mizutani, Kim [2004. Coastal Engineering (51, 407–420)]. The model is then adopted to investigate 3-D characteristics of wave pressure and force on a caisson of composite breakwater, and the numerical solutions were discussed with respect to the phase difference between harbor and seaward sides induced by the transmitted wave through the rubble mound or the diffraction. The numerical results reveal that wave forces acting on the composite breakwater are significantly different at each cross-section under influence of wave diffraction that is important parameter on 3-D wave interaction with coastal structures.     

Simulation of nonlinear wave run-up with a high-order Boussinesq model

This paper considers the numerical simulation of nonlinear wave run-up within a highly accurate Boussinesq-type model. Moving wet–dry boundary algorithms based on so-called extrapolating boundary techniques are utilized, and a new variant of this approach is proposed in two horizontal dimensions. As validation, computed results involving the nonlinear run-up of periodic as well as transient waves on a sloping beach are considered in a single horizontal dimension, demonstrating excellent agreement with analytical solutions for both the free surface and horizontal velocity. In two horizontal dimensions cases involving long wave resonance in a parabolic basin, solitary wave evolution in a triangular channel, and solitary wave run-up on a circular conical island are considered. In each case the computed results compare well against available analytical solutions or experimental measurements. The ability to accurately simulate a moving wet–dry boundary is of considerable practical importance within coastal engineering, and the extension described in this work significantly improves the nearshore versatility of the present high-order Boussinesq approach.     

Detection of nighttime sea fog/stratus over the Huanghai Sea using MTSAT-1R IR data

A dual channel difference (DCD) method is applied to detect nighttime sea fog/stratus over the Huanghai Sea using the infrared (IR) data of shortwave (3.5–4.0 μm) and longwave (10.3–11.3 μm) channels from the Multi-functional Transport Satellite (MTSAT)-1R, i.e., shortwave minus longwave brightness temperature difference (SLTD). Twenty-four sea fog events over the Huanghai Sea during March to July of 2006 and 2007 are chosen to determine a suitable value of SLTD for nighttime sea fog/stratus detection, and ...     

Velocity potential formulations of highly accurate Boussinesq-type models

The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water.     

关于非齐次马氏链的强极限定理

主要研究了非齐次马氏链的强极限定理.首先应用鞅差序列收敛定理给出了关于非齐次马氏链的任意k元函数一类平均值的极限定理.最后得到一系列相关状态序偶出现频率的推论.