天津地面相对湿度资料的非均一性检验及订正

研究中利用NCEP/DOEAMIP-ⅡReanalysis(R-2)比湿再分析资料和基础观测资料两种参考序列,对天津历史月平均相对湿度数据中存在的非均一性因素进行了检验和订正。结果显示,有9个地面站的相对湿度序列存在显著断点,占台站总数的69%。查阅台站元数据得到导致其突变的主要原因是自动站业务化,其次为迁站和仪器变更。从订正量来看,负订正量比例约占96. 3%,其中分布范围在-5.0%～-1.5%,约占总订正量的80%以上。对比分析订正前后月平均相对湿度序列的方差和趋势变化,均一性订正基本修正了非均一性因素造成的突变影响,减缓相对湿度序列长期出现异常偏干的趋势变化。同时,与朱亚妮等(2015)研制的同类数据产品误差分析发现,两类数据的误差平均值较小,但由于一些主观因素影响使得个别台站MAE、SE误差值达到3. 0%以上,造成MAE、SE误差值范围在0～2%的台站比例分别仅在84. 6%和76. 9%及两者以上。从而,一定程度上能够反映出本研究所采用的数据处理技术手段相对以往具有优越性。     

Homogenization of the Daily Land Surface Temperature over the Mainland of China from 1960 through 2017

Land surface temperature(LST) is one of the most important factors in the land-atmosphere interaction process. Raw measured LSTs may contain biases due to instrument replacement, changes in recording procedures, and other non-climatic factors. This study attempts to reduce the above biases in raw daily measurements and achieves a homogenized daily LST dataset over China using 2360 stations from 1960 through 2017. The high-quality land surface air temperature(LSAT)dataset is used to correct the LST warming biases especially evident during cold months in regions north of 40°N due to the replacement of observation instruments around 2004. Subsequently, the Multiple Analysis of Series for Homogenization(MASH) method is adopted to detect and then adjust the daily observed LST records. In total, 3.68 × 103 effective breakpoints in 1.65 × 106 monthly records(about 20%) are detected. A large number of these effective breakpoints are located over large parts of the Sichuan Basin and southern China. After the MASH procedure, LSTs at more than 80% of the breakpoints are adjusted within +/– 0.5°C, and of the remaining breakpoints, only 10% are adjusted over 1.5°C.Compared to the raw LST dataset over the whole domain, the homogenization significantly reduces the mean LST magnitude and its interannual variability as well as its linear trend at most stations. Finally, we perform preliminary analysis upon the homogenized LST and find that the annual mean LST averaged across China shows a significant warming trend [0.22°C(10 yr)–1]. The homogenized LST dataset can be further adapted for a variety of applications(e.g.,model evaluation and extreme event characterization).     

Application of a semi-spectral cloud water parameterization to cooling tower plumes simulations

In order to simulate the plume produced by large natural draft cooling towers, a semi-spectral warm cloud parameterization has been implemented in an anelastic and non-hydrostatic 3D micro-scale meteorological code. The model results are compared to observations from a detailed field experiment carried out in 1980 at Bugey (location of an electrical nuclear power plant in the Rhône valley in East Central France) including airborne dynamical and microphysical measurements. Although we observe a slight overestimation of the liquid-water content, the results are satisfactory for all the 15 different cases simulated, which include different meteorological conditions ranging from low wind speed and convective conditions in clear sky to high wind and very cloudy. Such parameterization, which includes semi-spectral determination for droplet spectra, seems to be promising to describe plume interaction with atmosphere especially for aerosols and cloud droplets.     

Influence of Thermal Effects on the Wind Field Within the Urban Environment

Micrometeorological conditions in the vicinity of urban buildings strongly influence the requirements that are imposed on building heating and cooling. The goal of the present study, carried out within the Advance Tools for Rational Energy Use towards Sustainability (ATREUS) European research network, is the evaluation of the wind field around buildings with walls heated by solar radiation. Two computational fluid dynamics (CFD) codes were validated against extensive wind-tunnel observations to assess the influence of thermal effects on model performance. The code selected from this validation was used to simulate the wind and temperature fields for a summer day in a specific region of the city of Lisbon. For this study, the meteorological data produced by a non-hydrostatic mesoscale atmospheric model (MM5) were used as boundary conditions for a CFD code, which was further applied to analyze the effects of local roughness elements and thermodynamic conditions on the air flow around buildings. The CFD modelling can also provide the inflow parameters for a Heating, Ventilation and Air Conditioning (HVAC) system, used to evaluate the building energy budgets and to predict performance of the air-conditioning system. The main finding of the present three-dimensional analyses is that thermal forcing associated with the heating of buildings can significantly modify local properties of the air flow.     

Adaptation of Pressure Based CFD Solvers for Mesoscale Atmospheric Problems

General purpose Computational Fluid Dynamics (CFD) solvers are frequently used in small-scale urban pollution dispersion simulations without a large extent of ver- tical flow. Vertical flow, however, plays an important role in the formation of local breezes, such as urban heat island induced breezes that have great significance in the ventilation of large cities. The effects of atmospheric stratification, anelasticity and Coriolis force must be taken into account in such simulations. We introduce a general method for adapting pressure based CFD solvers to atmospheric flow simulations in order to take advantage of their high flexibility in geometrical modelling and meshing. Compressibility and thermal stratification effects are taken into account by utilizing a novel system of transformations of the field variables and by adding consequential source terms to the model equations of incompressible flow. Phenomena involving mesoscale to microscale coupled effects can be analyzed without model nesting, applying only local grid refinement of an arbitrary level. Elements of the method are validated against an analytical solution, results of a reference calculation, and a laboratory scale urban heat island circulation experiment. The new approach can be applied with benefits to several areas of application. Inclusion of the moisture transport phenomena and the surface energy balance are important further steps towards the practical application of the method.     

提高任意广角波动方程计算效率及压制倏逝波干扰方法

任意广角波动方程（Arbitrarily Wide-angle Wave Equation,AWWE）是一种高精度的时空域单程波动方程,与其他时空域单程波方程相比,具有更大的成像倾角,因此是偏移成像的有力工具之一.本文改进了现有的AWWE有限差分计算方案,避免了原计算方案中对矩阵反复求逆的过程,同时减少了矩阵相乘的次数,显著地提高了AWWE数值计算的效率.数值实验证实,采用改进之后的计算方案,其计算结果与原方案完全一致.此外,本文采用了一种方法,在有倏逝波干扰的情况下,能在频率-波数（f-k）域区分出倏逝波与非倏逝波,为设计f-k滤波器提供合理的视速度门槛信息.偏移实验证实,在深度延拓时仅仅进行数次f-k滤波就能很好地抑制倏逝波干扰,提高偏移成像质量.     

基于计算岩石物理方法的页岩储层弹性参数提取

发展了应用数值计算方法获取页岩储层的速度、各向异性参数的计算岩石物理系列方法.该系列方法包括了大尺度精细地质模型数值建模、计算网格尺度的地球物理建模和地震波数值模拟提取岩石物理弹性参数.本文方法利用储层的统计数据而不是具体岩心的测量数据,可获得储层岩石物理弹性参数的变化规律.相比于基于岩心测试的岩石物理方法,本文方法可精细考虑实际储层的非均匀特征,可得到岩心测试难以求取的与尺寸效应高度相关的弹性参数,也避免了求取弹性参数变化规律时获取不同地质特征岩心的困难.本文发展了计算岩石物理方法,为计算岩石物理面临的大尺度地质建模和计算能力限制问题提供了有效的解决方案.文中以胜利罗家的页岩储层为例,求得了储层TOC含量从3%到21%变化情况下储层的P波、S波速度以及各向异性参数变化规律.

     

Methods for locating the proper position of a planned fishway entrance near a hydropower tailrace

The spatial distribution of upstream migrating Atlantic salmon (Salmo salar) spawners was studied in 2008 and 2009 in the surroundings of the tailrace from a hydropower station in the River Umeälven. This area is problematic because the fish have difficulties finding their way into the original riverbed, which prompted an investigation on the feasibility of adding a fishway in the area. Echo sounding was used in 2008 to investigate the spatial use of fish in the tailrace channel during the time of salmon migration. Presence of other fish species was so low that all echo sounding detections were assumed to be salmon. In 2009, data on wild radio-tagged salmon (n = 94) was collected in the same tailrace by an array of underwater antennas with a detection radius of approximately 10 m, to validate the results from the echo sounding. Both types of surveys showed aggregations of fish in one part of the tailrace. As a final step, Computational Fluid Dynamics (CFD) modeling was performed to analyze hydrodynamics. This CFD modeling showed a coincidence of fish detections in areas with turbulence intensities between 0.6 and 0.8, which may be an indication that the fish are holding in these areas to reduce energy expenditure during migration. A high proportion of the radio-tagged salmon were observed in the tailrace area (a median of 21 days between their first and last detection), indicating that the salmon are delayed in their upstream migration by the attraction to the tailrace in this area. The overall data on fish behavior in the tailrace led us to suggest a location for a new fishway where the fish aggregation was most pronounced. The number of detections from either technology had similar distributions over the tailrace, suggesting that the results are comparable. Thus, the split beam echo sounder can preferably be used to investigate which location is most appropriate for locations of entrances to new fishways since many more individual fishes are covered by this method compared to telemetry. Furthermore, there is no need to handle the fish in the echo sounding studies as is required in telemetry studies.     

地震波动方程的局部间断有限元方法数值模拟

近年来,随着地震波数值模拟对计算精度和效率的要求越来越高,间断有限元方法开始受到越来越多的关注.本文中,针对具有吸收边界条件的二维地震声波波动方程,作者提出了一种基于局部间断有限元方法的数值模拟算法.该算法在空间上使用局部间断有限元方法进行离散,在时间上采用了显式蛙跳格式.在这种时空离散的组合方式下,每个时间步上,此算法在空间剖分的每个单元上的求解计算是相互独立的,因而具有极高的并行性.通过数值算例,我们将该算法与连续有限元方法进行了比较.结果表明,本算法不仅具有对起伏构造的良好适应性,而且在计算效率和计算精度等方面,都具有优越性.     

SNREI地球模型负荷勒夫数数值计算的新进展

基于Ｌｏｎｇｍａｎ和Ｆａｒｒｅｌｌ的表面负荷形变理论，深入研究计算ＳＮＲＥＩ弹性地球模型表面负荷勒夫数的数值计算方法，取得新的进展：采用了一种新的代换方法─ｒ￣ｎ因子法，比传统的无量纲化法简单易行，且能达到很高的计算精度；确立了检验表面边界条件、微分方程积分线性性和计算结果的收敛性，用于评价数值计算结果精度的有效准则；给出了ＰＲＥＭ地球模型的一组高精度计算结果．