电磁脉冲数值异常初探

用ＤＷＭＪ－８９Ａ型电磁波记录仪记录的异常与地震３要素的关系，讨论了判别电磁波脉冲数值异常的方法与地震预报的问题。

基于线损数值异常波动的反窃电方法研究

线损数据计算和监测能够为反窃电工作提供重要支持,文章立足电网线损数据分析进行反窃电研究,意在明确线损的影响因素,并以此为基础提出反窃电的科学方法。窃电的主要方法包括欠/施压法、断/施压法和移项法三种。应用线损理论数值计算方法,可根据线损理论计算最大值和实时线损率波动异常解决窃电问题。技术人员需要运用科学的计算方法计算线损数值,精准判断窃电情况,并提出科学的反窃电方法。

激波诱导燃烧模拟中轴对称边界数值异常研究

对钝体激波诱导燃烧模拟中引起对称轴附近出现数值异常现象的几个影响因素进行了细致讨论.在飞行速度低于爆轰速度情形的模拟中,在一定条件下靠近对称轴的位置会出现明显的数值异常——燃烧阵面扭曲.分别从边界处理方法、计算格式、化学反应机制和计算网格等4个方面对该异常产生的原因进行了详细地讨论.研究表明4个因素都会导致数值异常产生,深入的分析表明计算格式是其中最根本的因素,而这种燃烧阵面扭曲的数值异常应该与粉刺问题(Carbuncle)类似.针对如何克服轴对称边界数值异常提出了一些合理的措施,并给出了验证算例的对比与分析.

空间波数混合域磁异常场积分解三维数值模拟

本文提出一种空间波数混合域磁异常场三维数值模拟方法.该方法利用磁位三维空间域积分为卷积的特点,沿水平方向进行二维傅里叶变换,把空间域磁位满足的三维积分问题转化为不同波数之间相互独立的垂向一维积分问题.保留垂向为空间域,优势之一在于便于浅层单元剖分可适当加密,随着深度增加,单元剖分适当稀疏,可以准确模拟任意复杂地形和磁性体的磁异常,兼顾了计算精度与计算效率;优势之二在于一维积分垂向可离散为多个单元积分之和,每个单元采用二次形函数表征磁化强度,可得出单元积分的解析表达式,计算精度高、效率高.该方法充分利用一维形函数积分的高效和高精度、快速傅里叶变换的高效性及算法高度并行性,实现了磁异常场高效、高精度的数值模拟.设计棱柱体模型,将模型解析解与空间波数混合域法的数值解对比,结果表明该方法计算精度高、效率高.设计了组合棱柱体复杂模型,对比分析了标准FFT扩边法与Gauss-FFT法的计算精度与计算效率,总结了标准FFT的扩边系数选取策略.针对任意复杂地形条件下的磁异常模拟问题,本文提出一种适用于起伏地形条件下的磁异常场快速计算方法,并对其有效性进行了验证.

借助EGM2008的海域重力异常反距离加权插值算法

构建高精度高分辨率的海域重力异常模型,是将海洋重力测量数据应用于全球高程基准统一和水下重力匹配导航等技术领域的关键步骤之一。针对反距离加权插值算法仅利用了海域重力观测值的空间信息,没有考虑重力观测值物理特性的缺陷,提出一种借助EGM2008重力场模型构建海域重力异常模型的反距离插值算法。以DTU10重力异常数值模型为基础数据,设计了海域重力测量数据格网化实验,结果表明本文算法能有效改进海域重力异常格网模型的计算精度。

应用地震数据体结构特征法预测油层分布规律

塔河油田八区的主产层为碳酸盐岩裂缝-缝洞型储层,具有很强的非均质性,给储层的油气预测带来了一定的困难。应用传统的储层预测方法很难准确预测出油区内油藏的油水界面,从而很难弄清储层的分布、发育规律。针对塔河油田八区特殊的地质条件,采用了地震数据体结构特征法对该区奥陶系储层进行了油气预测,并对油区内、外的区域进行了有利区块的划分。在有利区块中设计的18口钻探井位,经钻井证实均获得了较高的油气产量,取得了很好的经济效益。

电力大数据支持下VOCs污染综合防治方法研究

不同种类的挥发性有机化合物(volatile organic compounds, VOCs)污染在性质、挥发度和毒性上存在差异,确定适当的检测量限制数值具有较大难度。为此,引入电力大数据,研究新的VOCs污染综合防治方法。判断污染数据的检测阈值,记录设备中污染物的基础排放量,通过预警段区间信息判断污染数据检测量限制数值,实现VOCs污染数据检测,通过计算光源强度数值和阳光辐射强度参数建立污染传播模型,实现污染数据定量分析,通过电力大数据技术检测不同参数,实现污染综合防治。实验结果表明,研究方法应用后,半年内的VOCs污染月排放量在60万吨以下,防治优化效果明显。

SDZ-7130故障原因分析解决

本文对SDZ-7130双感应-八侧向测井仪的故障进行原因分析,从八侧向电极总成故障、八侧向电路板电容断裂、感应数值异常三个方面进行论述,针对各故障现象,提出解决方案,确保仪器性能稳定,提高测井质量和测井时效。1引言SDZ-7130双感应-八侧向测井仪应用频繁,针对小队在测井过程中就反映测井异常现象,对该仪器半年来常见故障进行分析。

Efficient and Effective 4D Trajectory Data Cleansing

As the rapid development of aviation industry and newly emerging crowd-sourcing projects such as Flightradar24 and FlightAware,large amount of air traffic data,particularly four-dimension(4D)trajectory data,have become available for the public.In order to guarantee the accuracy and reliability of results,data cleansing is the first step in analyzing 4D trajectory data,including error identification and mitigation.Data cleansing techniques for the 4D trajectory data are investigated.Back propagation(BP)neural network algorithm is applied to repair errors.Newton interpolation method is used to obtain even-spaced trajectory samples over a uniform distribution of each flight’s 4D trajectory data.Furthermore,a new method is proposed to compress data while maintaining the intrinsic characteristics of the trajectories.Density-based spatial clustering of applications with noise(DBSCAN)is applied to identify remaining outliers of sample points.Experiments are performed on a data set of one-day 4D trajectory data over Europe.The results show that the proposed method can achieve more efficient and effective results than the existing approaches.The work contributes to the first step of data preprocessing and lays foundation for further downstream 4D trajectory analysis.

A new algorithm for estimating gillnet selectivity

The estimation of gear selectivity is a critical issue in fishery stock assessment and management.Several methods have been developed for estimating gillnet selectivity,but they all have their limitations,such as inappropriate objective function in data fitting,lack of unique estimates due to the difficulty in finding global minima in minimization,biased estimates due to outliers,and estimations of selectivity being influenced by the predetermined selectivity functions.In this study,we develop a new algorithm that can overcome the above-mentioned problems in estimating the gillnet selectivity.The proposed algorithms include minimizing the sum of squared vertical distances between two adjacent points and minimizing the weighted sum of squared vertical distances between two adjacent points in the presence of outliers.According to the estimated gillnet selectivity curve,the selectivity function can also be determined.This study suggests that the proposed algorithm is not sensitive to outliers in selectivity data and improves on the previous methods in estimating gillnet selectivity and relative population density of fish when a gillnet is used as a sampling tool.We suggest the proposed approach be used in estimating gillnet selectivity.

Numerical Investigation of Transmission of Low Frequency Sound Through a Smooth Air-water Interface

It is the traditional belief that sound transmission from water to the air is very weak due to a large contrast between air and water impedances. Recently, the enhanced sound transmission and anomalous transparency of air-water interface have been introduced. Anomalous transparency of air-water interface states that the sound generated by a submerged shallow depth monopole point source localized at depths less than 1/10 sound wavelength, can be transmitted into the air with omni-directional pattern. The generated sound has 35 times higher power compared to the classical ray theory prediction. In this paper, sound transmission through air-water interface for a localized underwater shallow depth source is examined. To accomplish this, two-phase coupled Helmholtz wave equations in two-phase media of air-water are solved by the commercial finite element based COMSOL Multiphysics software. Ratios of pressure amplitudes of different sound sources in two different underwater and air coordinates are computed and analyzed against non-dimensional ratio of the source depth （D） to the sound wavelength （λ）. The obtained results are compared with the experimental data and good agreement is displayed.

NUMERICAL SIMULATION OF SSTA IMPACTS OVER THE GLOBAL OCEAN ON THE ANOMALOUS CIRCULATION OVER EURASIA IN JANUARY 2008

In this paper, we discussed the features of atmospheric circulations over Eurasia as a response to sea surface temperature anomalies （SSTAs） over the tropical Indian Ocean, the equatorial Pacific, Kuroshio and the North Atlantic. Our results are shown as follows： （1） CAM3.0, driven by the combined SSTAs over the four oceanic regions, can simulate well the features of anomalous atmospheric circulations over Eurasia in January 2008, indicating that the effects of the SSTAs over these four regions were one of the key causes of the anomalous systems over Eurasia. （2） The SSTAs over each key region contributed to the intensification of blocking over the Urals Mountains and a main East Asian trough. However, the influence of the SSTAs over individual oceanic regions differed from one another in other aspects. The SSTAs over the North Atlantic had an impact on the 500-hPa anomalous height （Z500A） over the middle-high latitudes and had a somewhat smaller effect over the low latitudes. For the warm SSTAs over Kuroshio, the subtropical high was much stronger, spread farther north than usual, and had an anomalous easterly that dominated the northwest Pacific Ocean. The warm SSTAs over the tropical Indian Ocean could have caused a negative Z500A from West Asia to Middle Asia, a remarkably anomalous southwesterly from the Indian Ocean to the south of China and an anomalous anticyclone circulation over the South China Sea-Philippine Sea region. Because of the La Nifia event, the winter monsoon was stronger than normal, with an anomalously cooler northerly over the southeastern coastal areas of China. （3） The combined effects of the SSTAs over the four key regions were likely more important to the atmospheric circulation anomalies of January 2008 over Eurasia than the effects of individual or partly combined SSTAS. This unique SSTA distribution possibly led to the circulation anomalies over Eurasia in January 2008, especially the atmospheric circulation anomalies over the subtropics, which were more similar to those of the winter E1 Nifio events than to the circulation anomalies following La Nifia.

NUMERICAL MODELING STUDY OF EFFECTS OF EASTERN PACIFIC WARM POOL ON ENSO

In this study, sensitivity experiments were conducted with the Zebiak-Cane ocean-atmosphere coupled model forced by the wind stress anomaly from the U.S. National Centers for Environmental Prediction/National Center for Atmospheric Research （NCEP/NCAR） reanalysis data to study the impacts of eastern Pacific warm pool on the formation and development of ENSO events. The effects of climatological mean sea surface temperature of the warm pool on forecast skill during the ENSO events of 1982-1999 are more considerable that those of climatological mean meridional winds and ocean currents. The forecast skill for the 1997/1998 E1 Nifio event is characterized by sensitivity to climatological mean sea surface temperature and anomalies of northerly winds and currents. The forecast skill is found insensitive to climatological mean northerly meridional winds and currents.

中小城市连续运行参考站脱离水准测量应用研究

文章结合泉州市连续运行参考站的建立实例,通过加入精化水准面模型和未加精化水准面模型两种实验数据比对,分析在建设CORS过程中,用精化水准面模型为基准才能得到更加精确的高程数据,中小城市在没有精化水准面模型前提下用坐标联测并进行高程数值拟合的方法求解出近似似大地水准面而得出的高程在测图方面也是可行的。

放射工作人员个人剂量监测异常增高调查分析

目的分析佛山市顺德区某医疗机构放射工作人员2014年第2季度个人剂量监测数值异常增高的原因。方法收集该医疗机构2014年度放射工作人员的个人剂量监测数据,以出现个人剂量监测数值异常增高的第2季度数据作为调查对象,采用现场职业卫生调查法,结合个人剂量计收发和日常管理情况进行分析。结果 2014年第2季度共检测16名放射工作人员的个人剂量计,人均有效剂量为38.51 m Sv,集体有效剂量为616.09人·m Sv,其中6名放射工作人员的有效剂量〉20.00 m Sv,年剂量当量为98.53~109.20 m Sv/a。现场职业卫生调查结果显示,各放射机房工作场所泄漏的辐射水平,以及个人剂量计收发期间存放的工作场所的辐射水平均接近本底值,不可能导致个人剂量监测结果异常。该医疗机构个人剂量计的日常存放存在漏洞,结合本次超标的6个个人剂量检测结果相近,推测本次个人剂量计异常增高的结果可能为人为统一进行照射所致。结论应加强放射防护知识的宣传教育和放射工作人员管理,制定科学合理的管理制度。

Numerical modeling of the anomalous post-rift subsidence in the Baiyun Sag,Pearl River Mouth Basin

The Baiyun Sag is the deepest sag in the Pearl River Mouth Basin in northern continental margin of South China Sea,with the maximum sediment thickness over 12.5 km above the basement including >6.5 km sediments above the 30 Ma breakup unconformity.According to the theoretical models for the rifted basins,the post-rift subsidence is driven solely by the thermal contraction and can be calculated as the function of the lithospheric stretching factor.A method combining the forward modeling and reverse backstripping was designed to estimate lithospheric stretching factor.Using the 2D forward modeling based on the flexural cantilever model,we simulated the multi-rifting process of the Baiyun Sag with constrain of the backstripped profiles.By doing this the lithospheric stretching factor was obtained,and then the theoretical post-rift subsidence was calculated.The calculated theoretical subsidence was much smaller than the observed subsidence given by backstripping.Along the 1530 line in the Baiyun Sag,the anomalous post-rift subsidence is over 2 km in the sag center,and varies slightly to the north and south edges of the sag.This suggests that the anomalous post-rift subsidence continues beyond the sag both in the continental shelf to the north and in the continental slope to the south.The sensitivity tests in the forward modeling process indicate that only the use of low-angle faults (≤13°) can we simulate the shape of the backstripped profile.

Numerical simulations on origin of Galilean moons' magnetic anomalies

Galileo mission detected the magnetic anomalies originated from Galilean moons.These anomalies are likely generated in the moons' interiors,under the influence of a strong ambient Jovian field.Among various possible generation mechanisms of the anomalies,we focus on magneto-convection and dynamos in the interiors via numerical simulation.To mimic the electromagnetic environment of the moons,we introduce in our numerical model an external uniform magnetic field B0 with a fixed orientation but varying field strength.Our results show that a finite B0 can substantially alter the dynamo processes inside the core.When the ambient field strength B0 increases to approximately 40% of the field generated by the pure dynamo action,the convective state in the core changes significantly:the convective flow decreases by 80% in magnitude,but the differential rotation becomes stronger in much of the fluid layer,leading to a stronger field generated in the core.The field morphologies inside the core tend to align with the ambient field,while the flow patterns show the symmetry-breaking effect under the influence of B0.Furthermore,the generated field tends to be temporally more stable.

Dynamics of saltwater intrusion in the Modaomen Waterway of the Pearl River Estuary

Observations indicate an abnormal characteristic of saltwater intrusion in the upper Modaomen Waterway of the Pearl River Estuary,i.e.,the maximum salinity occurs during the neap tide or the coming moderate tide instead of during the spring tide.To explore the associated dynamic mechanisms,a high resolution three-dimensional numerical model was set up based on the Finite Volume Coastal Ocean Model(FVCOM),which covered the entire river network,the Pearl River Estuary,and the adjacent sea.Numerical experiments illustrated that the upper Modaomen Waterway is significantly influenced by the saltwater intrusion from the Hongwan Waterway,a narrow and shallow channel connecting the Modaomen Waterway to the sea.Specific topography,spring-neap tidal variation,local wind stress,and their interaction drive an up-estuary residual current in the Hongwan Waterway,which is much stronger during the neap tide than during the spring tide.As a result,more saltwater in the Hongwan Waterway is spilled over into the Modaomen Waterway during the neap tide or the coming moderate tide.This is the inherent dynamic mechanism why the saltwater intrusion in the upper Modaomen Waterway reaches its maximum during the neap tide or the coming moderate tide.Besides,we also found that the winter prevailing wind can pronouncedly enhance the saltwater intrusion in the Modaomen Waterway.

Robust structured total least squares algorithm for passive location

A new approach called the robust structured total least squares(RSTLS) algorithm is described for solving location inaccuracy caused by outliers in the single-observer passive location. It is built within the weighted structured total least squares(WSTLS)framework and improved based on the robust estimation theory.Moreover, the improved Danish weight function is proposed according to the robust extremal function of the WSTLS, so that the new algorithm can detect outliers based on residuals and reduce the weights of outliers automatically. Finally, the inverse iteration method is discussed to deal with the RSTLS problem. Simulations show that when outliers appear, the result of the proposed algorithm is still accurate and robust, whereas that of the conventional algorithms is distorted seriously.