Intel MKL在Bernese GNSS数据处理中的应用

为提高Bernese GNSS software数据处理效率，将英特尔数学核心函数库(math kernel library，MKL)应用于Bernese 精密定轨数据处理，对比分析多个MKL矩阵求逆函数与Bernese 原有程序的计算效率。使用2019-03全球200个测站北斗/GNSS数据进行实验分析，结果表明，采用参数预消除策略时，参数预消除步骤消耗时间明显大于矩阵求逆，使用MKL处理数据效率提升不明显；而未采用参数预消除策略时，使用MKL矩阵求逆函数可显著提高矩阵求逆效率，其中dpotri函数矩阵求逆计算效率最高，消耗时间平均值为133 s，相比Bernese 原有程序计算速度可提高13倍。     

用传输矩阵方法研究多层介质增透膜的透射特性

利用麦克斯韦方程组和边界条件,推导出由高低折射率介质组成周期性的多层增透膜系的传输矩阵,用传输性矩阵求出增透膜系的反射率和透射率表示式,采用Maple语言和数值计算方法得到4层和8层高反射膜系的反射率随入射光波长变化曲线,分析和讨论入射光的入射角对反射率和透射率的影响.得到的结果有利于对多层介质增透膜系结构的理解,为光学仪器镀膜的设计提供了理论依据.     

基于弱势人群个体可达性评价的社区服务设施配置研究

当前我国人口形势面临老龄化与少子化风险,抓好“一老一小”、确保老有所养和幼有所育成为我国城乡发展中特别关注的方面。在这一背景下,如何营造有利于养老和生育的社会环境,如何针对老年人、育龄妇女和儿童等弱势人群提供更为合理的设施配置成为民生工程的重要内容。社区服务设施作为城市公共服务的空间载体,其配置模式的完善与提升是落实以人为本、集中体现社会公平的重要路径。当前,从网络地图规划路径API获取的出行时间成本矩阵能够为可达性研究提供接近真实的出行时间数据,此外个体化的人口数据如实有人口数据能够以其丰富的属性信息为识别各类弱势人群提供依据,同时精确刻画人口的空间分布。本文基于弱势人群的需求特征与相关规范,对老年人、育龄妇女和儿童三类弱势人群提出了明确定义,并构建了面向各类弱势人群个体使用需求的社区服务设施供给标准。在此基础上,利用实有人口数据、网络地图出行时间成本矩阵和POI数据,综合考虑社区服务设施的服务容量,提出了基于弱势人群个体可达性评价的社区服务设施供给评价与布局优化方法,并以广州市人民南社区为例进行了应用实践。研究表明,人民南社区的弱势人群服务设施供给存在不足,优化配置后设施的供需情况更加合理,服务人口覆盖率亦明显提升。本文提出的方法能够从独立个体尺度精确评价社区弱势人群服务设施的布局合理性,在拓展数据类型及应用方法方面为当前人口形势下精细化研究服务设施的空间分布和供需问题提供新的思路。     

SCATTERING PHASE MATRICES OF ICE CRYSTALS WITH HEXAGONAL PRISM AND TRIANGULAR PYRAMID FORM—A VECTOR RAY TRACING METHOD

The work presented previously by the authors(Cai and Liou,1982)has been extended in this paper.By making use of our improved model the calculations on scattering phase matrices of hexagonal prismice crystals(HPIC)have been conducted for monodisperse and polydisperse systems.Compared with themodel of Cai and Liou,the required computational quantity is decreased by about two orders of magni-tude and the errors of results are less for the new model.Meanwhile,the scattering phase matrices of triangularpyramid ice crystals(TPIC)are also computed in the paper,and the comparison between the scatterings of thetwo forms of ice crystals is performed.     

扩展SSH Hubbard模型的晶格振动矩阵

在计及最近邻电子相互作用的扩展ＳＳＨＨｕｂｂａｒｄ模型时，得到了晶格振动矩阵，它对电子的正负自旋态对称和非对称的各类非线性元激发情形均适用。     

湖相泥页岩基质储层直井压裂选层方法研究——以松辽盆地三肇凹陷松页油3井为例

湖相泥页岩基质储层是中国页岩油赋存的主要载体,相对于夹层型和裂缝型储层更难被压裂改造和实现工业开发,因而压裂层位的选择至关重要. 本文以松辽盆地三肇凹陷松页油3井为例,建立了一套湖相泥页岩基质储层直井压裂选层方法. 通过录井、测井解释优选出青山口组一段中的I类页岩油层,对I类页岩油层的电性、物性、含油气性、烃源岩性、脆性、可压性6项指标进行综合分析,基于上述6项指标,利用层次分析法对所有I类页岩油层进行综合排序,确定了27、25和23号层为最有利页岩油层. 松页油3井针对上述页岩油层进行了压裂试油,实现了三肇凹陷页岩油工业油流突破,证实选层方法可行有效,可推广应用于湖相基质型页岩油的甜点层优选和水平井部署.     

煤样渗透率有效应力与基质收缩双重耦合效应及其与煤阶关系

为研究煤层气在排采过程中不同煤阶煤储层渗透率动态变化规律,利用煤岩三轴应力应变（基质收缩膨胀）测试系统,对褐煤、气煤和无烟煤样开展了有效应力与基质收缩双重效应物理模拟实验。固定轴压和围压不变,改变气体平衡压力,模拟开发过程中储层压力变化特征,测试其动态渗透率。利用实验结果,分析了不同煤阶煤岩在排采过程中动态渗透率反弹特征,并对比分析煤岩动态渗透率改善效果的差异性。研究表明：气体平衡压力从5 MPa降至1 MPa过程中,在有效应力和基质收缩双重效应作用下,褐煤样的归一化渗透率依次为1. 00、0. 60、0. 57、0. 57、0. 52,气煤样依次为1.00、0. 64、0. 50、0. 54和0. 55,无烟煤样依次为1.00、0. 74、0. 58、0. 50和0. 56。随气体平衡压力下降,中阶及高阶煤样动态渗透率先下降后上升,整体呈不对称“V”型变化规律,但拐点略有不同;低阶煤样动态渗透率呈先下降后基本稳定的趋势,整体呈斜“L”型变化规律。在有效应力和基质收缩双重效应影响下,中阶及高阶煤样动态渗透率改善效果优于低阶煤样。     

常用锂同位素地质标准物质的多接收器电感耦合等离子体质谱分析研究

锂同位素研究是非传统稳定同位素地球化学研究的前沿,已广泛应用于从地表到地幔的岩石圈及流体等固体地球科学的研究领域。准确测定锂同位素比值是应用该同位素体系的前提。本文报道了国际上7种常用地质标准物质(BHVO-2、JB-2、BCR-2、AGV-2、NKT-1、L-SVEC、IRMM-016)的锂同位素组成数据。分析中采用硝酸-氢氟酸混合酸消解岩石标准样品,通过3根阳离子交换树脂(AG50W-X8,200~400目)填充的聚丙烯交换柱和石英交换柱对锂进行分离富集,利用Neptune型多接收器电感耦合等离子体质谱(MC-ICPMS)测定锂同位素比值,使用标准-样品交叉法(SSB)校正仪器的质量分馏。实验得到这7种常用地质标准物质的锂同位素组成与测试精度(2SD)分别为:δ7LiBHVO-2—L-SVEC=4.7‰±1.0‰(n=53),δ7LiJB-2—L-SVEC=4.9‰±1.0‰(n=20),δ7LiBCR-2—L-SVEC=4.4‰±0.8‰(n=8),δ7LiAGV-2—L-SVEC=6.1‰±0.4‰(n=14),δ7LiNKT-1—L-SVEC=9.8‰±0.2‰(n=3),δ7LiL-SVEC—L-SVEC=-0.3‰±0.3‰(n=10),δ7LiIRMM-016—L-SVEC=0.0‰±0.5‰(n=10),这些数据在误差范围内与国际上已发表的数据一致。Li同位素分析精度可以达到大约0.5‰,长期的分析精度即外部重现性≤±1.0‰,达到了国际同类实验室水平。7种常用地质标准物质的锂同位素组成数据的发表为锂同位素研究提供了统一的标准,使地质样品的锂同位素数据的质量监控成为可能。在基质效应的研究中,使用不同量的IRMM-016配制的标准溶液过柱,深入探讨了样品量对锂同位素测定值的影响,结果表明,在现有测试精度下,只要分析样品的锂含量达到100μg/L,且不超过树脂的承载量,样品的锂同位素组成在误差范围内与真值吻合,样品量的大小不影响锂同位素测定结果的准确性。     

泥质杂基含量对柴达木盆地昆2井孔隙演化模式影响的压实模拟实验研究

随着勘探程度的提高,油气勘探必然向深层延伸。弄清深部储层的孔隙演化规律,提高预测精度,是深入勘探的重要基础。通过对柴达木盆地昆2井岩屑和岩芯样品的系统采集,运用铸体薄片、扫描电镜和测井解释等手段,系统分析了古近系与新近系地层的岩石学、矿物学和储层物性特征,发现泥质杂基含量与砂岩的孔隙演化密切相关。为了正演孔隙的演化过程,运用分选很好的风成石英砂和粒度细小的红色粘土为材料,按照粘土杂基含量为5%、10%、15%、20%,制作了四类石英砂开展模拟实验。松散沉积物的压实模拟,从1.4MPa试验加压到210.4MPa,分别进行了21个循环实验。将实验模拟的结果,与昆2井的岩芯资料进行了对比分析。研究表明,碎屑岩泥质杂基含量越高,初始孔隙度越小,孔隙度衰减越快;泥质含量较低的净砂岩,当上覆压力达到76.5MPa时,颗粒大量破碎并重新排列,导致粒间体积减缩以及少量直线型粒内裂缝出现,储层物性明显变差。压实模拟的实验数据,能够为快速沉降区砂岩的储层质量预测提供依据。     

Generating non-stationary random fields of auto-correlated,normally distributed CPT profile by matrix decomposition method

This paper emphasises the true realisation of Cone Penetration Test (CPT) profiles considering non-stationary nature of the data. Formulation of stationary random field theory has been modified and adapted to non-stationary state in order to take into account the mean and variance variability for soil properties. Multi-variance correlation matrix along with the Cholesky decomposition technique was employed to produce realisations of non-homogenous and non-stationary random fields of CPT profiles. A piecewise and segmental data realisation according to the lithology and site class specifications acquired directly from CPT data is adopted in this study so as to render an accurate data simulation. For validation of proposed method 8 CPT test profiles collected from Urmia Lake site have been introduced and simulated by the stationary and non-stationary algorithms. The mean correlation coefficient between the actual CPT data profiles and related realisations along with some other important statistical parameters and their coefficients of variation strongly demonstrate that non-stationary random field generation technique gives quite better accuracy, by comparison to the conventional stationary random field generation scheme.