面向矿产资源信息的空间关联性分析

关联性分析是综合多源、多类数据,发现和挖掘数据中潜在的相关关系,提取和挖掘数据之间的关联性。该方法对于充分利用地质大数据、发现地质要素之间的共生组合规律具有重要的意义。本文首先基于数据的空间位置,将不同类型的地质空间数据建立联系,形成空间属性数据库;其次应用统计方法,对不同来源数据中的属性特征进行分析,发现热液型金矿的形成与火山作用存在明显相关性;最后基于Apriori算法提取热液型金矿的伴生矿与侵入岩的频繁项集,发现伴生矿与侵入岩酸碱性二者密切相关。因此,在今后的工作中,有望应用空间关联性的方法对地球物理、地球化学、遥感等多源地质数据开展深入的研究分析。     

Brittleness index predictions from Lower Barnett Shale well-log data applying an optimized data matching algorithm at various sampling densities

The capability of accurately predicting mineralogical brittleness index(BI)from basic suites of well logs is desir-able as it provides a useful indicator of the fracability of tight formations.Measuring mineralogical components in rocks is expensive and time consuming.However,the basic well log curves are not well correlated with BI so correlation-based,machine-learning methods are not able to derive highly accurate BI predictions using such data.A correlation-free,optimized data-matching algorithm is configured to predict BI on a supervised basis from well log and core data available from two published wells in the Lower Barnett Shale Formation(Texas).This transparent open box(TOB)algorithm matches data records by calculating the sum of squared errors be-tween their variables and selecting the best matches as those with the minimum squared errors.It then applies optimizers to adjust weights applied to individual variable errors to minimize the root mean square error(RMSE)between calculated and predicted(BI).The prediction accuracy achieved by TOB using just five well logs(Gr,pb,Ns,Rs,Dt)to predict BI is dependent on the density of data records sampled.At a sampling density of about one sample per 0.5 ft BI is predicted with RMSE～0.056 and R2～0.790.At a sampling density of about one sample per 0.1 ft BI is predicted with RMSE～0.008 and R2～0.995.Adding a stratigraphic height index as an additional(sixth)input variable method improves BI prediction accuracy to RMSE～0.003 and R2～0.999 for the two wells with only 1 record in 10,000 yielding a BI prediction error of＞±0.1.The model has the potential to be applied in an unsupervised basis to predict BI from basic well log data in surrounding wells lacking mineralogical measure-ments but with similar lithofacies and burial histories.The method could also be extended to predict elastic rock properties in and seismic attributes from wells and seismic data to improve the precision of brittleness index and fracability mapping spatially.     

Plio-Pleistocene sedimentation in coastal grabens,north-east Rhodes,Greece

The stratigraphy of marine Plio-Pleistocene sediments from northeast Rhodes (Dodecanese islands, Greece) is revised in the light of facies mapping and the recognition of three major transgression–regression cycles. Before late Pliocene submergence, metamorphosed Mesozoic limestones formed a high-relief karstic landscape. During subsequent transgressions, subsiding basins with spectacular cliffed and bioeroded margins were infilled by a mosaic of carbonate-dominated sediments, sourced by high rates of carbonate productivity on narrow shelves rimming the basins and their drowning islands. Periods of relative sea-level fall superimposed a complex series of ‘fossil’ coastal geomorphological features such as cliffs, abrasion platforms, surf caves, notches, boulder beaches and palaeokarsts. These created the rugged present day topography of the island's northeast coast. Three main lithostratigraphic units, the Kritika, Rhodes and Lindos Acropolis Formations, are defined. The Rhodes and Lindos Acropolis Formations are subdivided into mappable facies groups rather than members, as the lithofacies are strongly diachronous. The Kritika and Rhodes Formations were each subaerially exposed and partly eroded before re-transgression. Although there is a general pattern of lithological succession, details vary across the study area and some facies groups are restricted in development, indicating differences in tectonic behaviour and palaeo-geomorphology between neighbouring basins. A preliminary process-response model is presented for the sedimentation history of the Rhodes and Lindos Acropolis Formations, using sedimentological, palaeontological, palaeoecological and ichnological data.     

使用因子分析方法研究九江5.7级地震前的地震活动异常

使用2005年11月26日江西九江5.7级地震前震中附近地区的地震资料, 选择反映地震活动时、 空、 强特征的地震频次N(M_L≥3.0)、 b值、 η值、 A(b)值、 Mf值、 Ac值、 C值和D值参量进行因子分析。 这些参量之间具有一定的相关性, 各参量在不同时段的变化各有所异, 预报效果并不理想。 但是根据因子分析可以得到反映地震活动时、 空、 强特征的综合指标W_fa值, 该指标在九江5.7级地震前出现明显的异常变化。 表明综合指标W_fa值可以较好地反映地震活动的异常特征。 文中还对与因子分析结果的有关问题进行了讨论。     

干旱生态系统净初级生产力估算及变化探测

针对中国西部干旱区特有的气候-植被特征,利用卫星遥感混合像元分解技术对干旱区植被进行了光合作用植被和非光合作用植被区分和组分解析,尝试通过对干旱覆被植被灌层结构进行解析而解决因植被区系差异和环境背景干扰的问题,并参考国际上遥感-生态模型GLO-PEM和CASA,借助遥感生态反演的物理分析,初步构建起基于遥感与生态过程的干旱区适用的光能利用率模型NPP-PEM,并以中国西部干旱区喀什地区叶尔羌-喀什噶尔河流域山地-绿洲-荒漠生态系统为案例,利用AVHRR/NOAA气象卫星遥感数据和气候资料估算了1992年和1998年中国西部喀什地区叶尔羌-喀什噶尔河流域山地-绿洲-荒漠生态系统1 km分辨率年净第一性生产力,并进行了变化探测分析。模拟检验结果精度较好,生态系统碳吸收的空间异质性特征明显。结果表明,考虑了干旱植被生理特征和灌层结构的光能利用模型,模拟结果较为合理,也为引入其他生态模型应用到干旱区生态系统研究提供了借鉴,从而为干旱区陆地生态系统碳循环研究开辟了途径。     

三参数小波及其在地震资料分析中的应用

提出了一类含有三个可调参数的分析小波（称为三参数小波）,通过选择这三个参数,可使其最适合于所给定的问题．Morlet小波及改进的Morlet小波是其特例．通过对三参数小波中的参数加以约束,可以获得一类新的近似解析小波．论证了通过恰当的选择参数,三参数小波不仅适合于分析包含慢变频率和振幅分量的信号,而且也适合于包含快变分量的信号．以三参数小波为分析小波,提出了一种用于薄互层地震资料分析的方法．该方法能够刻画薄互层的沉积旋回,也能研究薄互层内部的局部结构,并采用模型资料验证了方法的有效性．     

立式储液罐地震破坏快速评估方法

本文采用大型有限元程序ANSYS对立式储液罐的地震反应进行了非线性数值分析,考虑了罐体和内部液体的相互作用,得到了3个不同容积的储罐,在不同峰值的地震动时程输入下的反应结果。结合地震灾害中立式储罐的地震破坏现象和特点,以及现有的立式储罐地震破坏预测方法,提出了立式储罐的地震破坏快速评估方法。     

利用气象要素计算五道梁地区土壤热流量的试验

根据五道梁地区1994—1997年的热流板测量土壤热流值,初步分析了该地区土壤热流随时间的变化特征。利用1994—1996年的月平均资料分别分析了总辐射、地表温度、地表温度和80 cm深度温度差和2 cm、10 cm深度的土壤热流的相关性,并建立了相应的回归方程;利用得到的回归方程对1997年的土壤热流作计算检验,就月平均资料来说,用总辐射和地温差来计算浅层土壤热流结果是比较可靠的。     

利用奇异谱分析提高极移中长期预报精度

由于空间大地观测数据传输耗时及处理过程复杂, 导致极移测量值的获取存在时延, 无法满足对高精度的极移预报值有重大需求的应用领域. 针对极移复杂的时变特性, 提出一种基于奇异谱分析(singular spectrum analysis, SSA)的预报方法. 首先用SSA分离提取极移时序中的高频组分与低频组分; 其次建立最小二乘(least square, LS)外推与自回归(AutoreGressive, AR)模型对极移高频和低频组分进行组合预报. 结果表明, SSA方法能够准确地分离和提取极移低频和高频组分, 对低频和高频组分组合预报可以显著改善极移的中长期(30--365d)预报精度, 与国际地球自转和参考系服务局(International Earth Rotation and Reference Systems Service, IERS)提供的A公报中的极移预报值相比, SSA方法对极移PMX分量(本初子午线方向)和PMY分量(西90$^\circ$子午线方向)的中长期预报精度改进最高分别可达45.97%和62.44%. 研究结果验证了SSA方法对极移中长期预报改进的有效性.     

大型射电望远镜主反射面调整方法研究

随着毫米波天文学和空间通信的重要性日益提高, 对天线性能提出了越来越高的要求, 而天线性能往往受到其反射器表面精度的限制. 微波全息技术是一种快速有效的检测反射面天线表面轮廓的测量技术. 通过微波全息测量得到天线口径场, 计算天马65m射电望远镜反射面与理想抛物面的偏差. 天马65m射电望远镜的主反射面板是放射状的, 有14圈. 面板的每个角都固定在面板下方促动器的螺栓上进行上下移动, 且相邻面板交点处的拐角共用一个促动器. 采用平面拟合的方法可以计算各块面板拐角处的调整值, 但是同一个促动器会得到4个不同的调整量. 通过平面拟合, 同时以天线照明函数为权重的平差计算方法得到相邻面板拐角的一个平差值, 即天马65m射电望远镜1104个促动器的最佳调整值. 通过多次调整和新算法的应用, 天马65m射电望远镜反射面的面形精度逐渐提高到了0.24mm.