山西省2017年农田灌溉水有效利用系数测算分析研究

以山西省172个样点灌区为研究区域,按照《全国农田灌溉水有效利用系数测算分析技术指导细则》测算方法及要求,测算样点灌区农田灌溉水有效利用系数,在此基础上,以全省灌区统计资料为权重,测算分析得出山西省2017年大、中、小和纯井灌区农田灌溉水利用系数分别为0.484 6、0.501 0、0.475 8和0.639 8,全省系数为0.538。通过对测算成果分析,认为测算结果真实反映了当年山西省灌区农田灌溉用水效率,相关成果对山西省灌区工程建设、用水管理和规划编制等有一定指导意义。     

不同定义的湿位涡分析及在台风中的诊断

本文对不同定义的湿位涡做了理论分析,并利用1522号台风“彩虹”的数值模拟结果对各种湿位涡进行了诊断。主要结论有:经典湿位涡、广义湿位涡和改进湿位涡的差异主要是由不同定义的位温造成的,相当位温、广义位温和修改位温的构成均是在位温基础上添加一显含水汽的附加量;经典湿位涡、广义湿位涡和改进湿位涡的构成均能分为干、湿分量两部分,其干分量表达式相同,都与Ertel干位涡的定义一样,水物质相变潜热的影响隐含在位温中;不同定义湿位涡的本质差异表现在不同的湿分量上,湿分量的表达式中显含了水物质的作用。对台风的诊断分析发现,改进湿位涡分布与Ertel干位涡非常相似,呈现中空分布的位涡塔结构,大值区对应眼墙内侧,改进湿位涡湿分量与经典湿位涡的湿分量分布相似,只是湿分量的绝对值更小,这反映了改进湿位涡既能保持干位涡的分布特征,其分布和演变可反映台风的结构和演变,又能合理地体现水汽分布的影响,所以在台风诊断中有更广泛的应用前景。经典湿位涡在低层表现为负值,这与水汽梯度的分布关系很大,但与垂直速度、潜热加热大值区等都没有很好的匹配关系,用其分析台风结构和演变具有一定局限性;广义湿位涡其形式较复杂,仅在近饱和区域才能发挥其诊断优势。     

高低频系统相互作用对Rossby波破碎过程的影响

利用美国国家环境预报中心和能源部（NCEP/DOE）的逐日再分析资料（NCEP-DOE AMIP-Ⅱ）,对2010年12月20日发生在北太平洋一次典型的反气旋式波破碎（AWB）事件进行研究,分析了波破碎过程中等熵位涡场的演变特征,揭示了波破碎过程中高频扰动以及低频信号的逐日演变特征,并对2010年冬季350 K等熵面上逐日高频位涡（PV）扰动和低频变化做经验正交函数（EOF）分析,得到了其主要模态,并从等熵位涡方程出发研究了波破碎过程中位涡高、低频变化的原因。研究表明,波破碎过程中高频低PV空气从北太平洋西部日本附近沿东北方向向对流层上层侵入,而来自阿拉斯加湾附近的高频高PV空气向对流层下层侵入。高频位涡场EOF分解得到的前两个模态共同描述了北太平洋中纬度地区自西向东移动的天气尺度波列;低频位涡场EOF分解的第一模态在北太平洋呈弧形波列结构。天气尺度波在传播过程中受到低频场的平流作用逐渐偏离其传播主要模态的位置,并发生破碎,同时高频流场对高频位涡的平流可以产生低频变化,使得低频变化的空间形态向其冬季主要模态转变。     

球面距离计算方法及精度比较

球面距离(角间距)计算是天文或地理学中极常用的计算之一,也是目标查找、锥形检索、交叉证认等方法的基础。数学上,通过球面几何可以直接计算出两点的距离,前人已经推导出了多个复杂程度不一的计算方法。但是由于计算机的精度有限,在进行数值计算时有舍入误差,导致公式计算结果出现偏差。对几个常用的球面距离计算公式进行了考察,测试并对比它们在不同计算环境下的精度与优缺点。此外还展示并比较了几种常用天文软件包、数据库的球面距离计算方法,以期有助于天文工作者选择适合自己当前需要的计算方法。     

An efficient hybrid integral-equation method for point-absorber wave energy converters with a vertical axis of symmetry

To assist in the prototyping and controller design of point-absorber wave energy converters (WECs), an easy-to-implement hybrid integral-equation method is presented for computing the frequency-domain hydrodynamic properties of bodies with a vertical axis of symmetry in waves. The current hybrid method decomposes the flow domain into two parts: an inner domain containing the body and an outer domain extending to infinity. The solution in the inner domain is computed using the boundary-element method, and the outer-domain solution is expressed using eigenfunctions. Proper matching at the domain boundary is achieved by enforcing continuity of velocity potential and its normal derivative. Body symmetry allows efficient computation using ring sources in the inner domain. The current method is successfully applied to three different body geometries including a vertical truncated floating cylinder, the McIver toroid, and the coaxial-cylinder WEC being developed in the authors’ laboratory. In particular, the current results indicate that, by replacing the flat bottom of the coaxial-cylinder WEC with the Berkeley-Wedge (BW) shape, viscous effect can be significantly reduced with only minor negative impact on wave-exciting force, thus increasing WEC efficiency. Finally, by comparing to experimental measurements, the current method is demonstrated to accurately predict the heave added mass and wave-exciting force on the coaxial-cylinder WEC with BW geometry. If a viscous damping correction factor is used, the heave motion amplitude can also be accurately computed.     

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.     

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.     

胶州湾增养殖海域营养状况与赤潮形成的初步研究

根据对胶州湾女姑山增养殖海域1998年5月～9月的连续监测资料,参照潜在性富营养化的概念,应用NQI指数对该海域的营养状况进行分析。认为该海域水质富营养化是7月3～8日Skeletonemacostatum和Biddulphiaaurita混合型赤潮形成的基础,磷、硅营养盐的消耗是赤潮消亡的主要原因;赤潮消亡之后浮游植物群落发生演替,水体叶绿素a仍保持较高含量,NQI指数也相应较高,水质表现为磷限制潜在性富营养化,由于磷酸盐的限制没有发展为赤潮。     

逐步逼近曲化平方法研究

位场资料的常规处理解释方法是建立在平面数据理论之上的,然而实际的位场数据大多为曲面数据,若把曲面数据当成平面数据进行处理,必然导致很大的误差,因此有必要进行“曲化平”处理. 曲化平目前存在的主要问题是计算精度低、计算量大,因此,研究快速、精度高且适合大数据量处理的曲化平方法具有重要的价值. 本文在已有的泰勒级数曲化平方法基础上提出了逐步逼近技术和平均平面技术,使得曲面位场资料处理的精度得到了很大提高. 最后通过理论模型和实际资料的处理验证了该方法的应用效果.     

位场曲化平的插值-迭代法

将起伏曲面B上的位场向下延拓至曲面最低点的平面A的插值-迭代法步骤是：1）将曲面B上的场值放置在水平面A上具有相同水平坐标的点上,作为A上的初值;2）用若干水平面切割B,从A的初值,用快速傅里叶变换法(FFT)向上延拓出这些平面的场值,用插值的方法从这些平面的场值计算曲面B的场值;3）根据B上的实测值与计算值的差值,对A上的值进行加权改正;4）重复步骤2）和3）,直到B上的差值小到可以忽略.这种插值-迭代法具有高的计算速度,比通常的FFT法延拓得更深,可以超过10倍点距.文中给出计算实例.