一种模拟非均匀介质中弹性波传播新的k-space方法          下载免费PDF全文

黄继伟  刘洪 《地球物理学报》2020,63(8):3091-3104

传统的伪谱(PS)方法,采用傅里叶变换(FT)计算空间导数具有很高的精度,每个波长仅需要两个采样点,而时间导数采用有限差分(FD)近似因而精度较低.当采用大时间步长时,由于时空精度不平衡,PS法存在不稳定性问题.原始的k-space方法可以有效地克服这些问题但是却无法适用于非均匀介质.为了提高原始k-space方法模拟非均匀介质波动方程的精度,我们提出了一种新的k-space算子族.它是用非均匀介质的变速度代替原k-space算子中的常数补偿速度构造得到,引入低秩近似可以高效求解.我们将构造的新的k-space算子应用于耦合的二阶位移波动方程,而不是交错网格一阶速度应力波动方程,使模拟弹性波的计算存储量减少.我们从数学上证明了基于二阶波动方程的k-space方法与基于一阶波动方程的k-space方法是等价的.数值模拟实验表明,与传统的PS、交错网格PS和原始的k-space方法相比,我们的新方法可以在时间和空间步长较大的均匀和非均匀介质中,为弹性波的传播提供更精确的数值解.在保持稳定性和精度的同时,采用较大的时空采样间隔,可以大大降低数值模拟的计算成本.  相似文献   

A nearly analytic symplectic partitioned Runge–Kutta method based on a locally one-dimensional technique for solving two-dimensional acoustic wave equations     

Chol Sim  Chunyou Sun  Nam Yun 《Geophysical Prospecting》2020,68(4):1253-1269

In this paper, we develop a new nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique for numerically solving two-dimensional acoustic wave equations. We first split two-dimensional acoustic wave equation into the local one-dimensional equations and transform each of the split equations into a Hamiltonian system. Then, we use both a nearly analytic discrete operator and a central difference operator to approximate the high-order spatial differential operators, which implies the symmetry of the discretized spatial differential operators, and we employ the partitioned second-order symplectic Runge–Kutta method to numerically solve the resulted semi-discrete Hamiltonian ordinary differential equations, which results in fully discretized scheme is symplectic unlike conventional nearly analytic symplectic partitioned Runge–Kutta methods. Theoretical analyses show that the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique exhibits great higher stability limits and less numerical dispersion than the nearly analytic symplectic partitioned Runge–Kutta method. Numerical experiments are conducted to verify advantages of the nearly analytic symplectic partitioned Runge–Kutta method based on locally one-dimensional technique, such as their computational efficiency, stability, numerical dispersion and long-term calculation capability.  相似文献   

Structure-from-motion photogrammetry analysis of historical aerial photography: Determining beach volumetric change over decadal scales     

Rafael C. Carvalho  David M. Kennedy  Yakufu Niyazi  Chloe Leach  Teresa M. Konlechner  Daniel Ierodiaconou 《地球表面变化过程与地形》2020,45(11):2540-2555

Historical aerial photographs are an invaluable tool in shoreline mapping and change detection in coastal landscapes. We evaluate the extent to which structure-from-motion (SfM) photogrammetric methods can be applied to quantify volumetric changes along sandy beaches, using archival imagery. We demonstrate the application of SfM-derived digital surface models (DSMs) at East Beach and Lady Bay in southwest Victoria, Australia, using photographic datasets taken in 1969, 1977 and 1986, and compare them to LiDAR-derived DSMs acquired at both sites in 2007. The SfM approaches resulted in two entire and two partial suitable DSMs out of six datasets. Good-quality DSMs were spatially continuous with a good spread of ground control points (GCPs) near the beach at Lady Bay, whereas unsuitable DSMs were mostly restricted by poor distribution and number of GCPs in spatially segmented areas of East Beach, due to limited overlapping of images, possible poor quality of GCPs and also the propagation of errors in the derived point clouds. A volume of approximately 223 000 ± 72 000 m³ was deposited at Lady Bay between 1969 and 2007, despite minimal erosion observed near the breakwater. The partially suitable dataset of East Beach indicated that beach erosion of at least 39 m³ m⁻¹ occurred immediately to the east of the seawall after 1977. We also discuss the drawbacks and strengths of SfM approaches as a benchmark of historical erosion assessments along sandy beaches. © 2020 John Wiley & Sons, Ltd.  相似文献   

利用CEEMD分析中国沿海GNSS站高程时序变化     

刘豪  周东旭  王盼龙  张化疑  王永康  姜传苓 《测绘通报》2020,(3):39-43,77

利用完备经验模态分解方法（CEEMD）对我国沿海地区6个GNSS基准站（2010—2018）的高程时序数据进行了处理分析。结果表明：CEEMD在高程时间序列分析中具有一定的优越性，可准确分解出各GNSS站高程时序中存在的周、月、季节、年等变化周期项，其中周年运动是主要贡献项，各站高程时间序列的短周期变化与潮汐变化周期具有密切关联性；沿海GNSS站的地面沉降既具有区域的一致性，又存在区域间差异性，其中D区DBJO、DZJJ站呈现先下降后上升的趋势，N区NZUH、NWZU站呈下降趋势，B区的BZMW呈上升趋势，而同海区的BLHT站则呈显著的下降趋势。  相似文献   

基于时相变化的晴空条件下Himawari-8 AHI火点检测     

鄢俊洁  瞿建华  冉茂农  张芳芳 《遥感学报》2020,24(5):571-577

Himawari-8静止气象卫星具有高空间分辨率、高观测频次和高时效特点,对于火点检测具有很强优势。对Himawari-8卫星的3.9μm和11.2μm两通道亮温进行了连续时相变化研究,得出两通道的亮温在时间上的变化差值稳定且规律明显。根据两通道的亮温时相特征,考虑白天可见光对3.9μm通道的影响,并结合火点产生时引起的亮温变化特征,提出了适用于晴空条件下改进的火点检测算法。在多处进行了此算法的实验,例如2018-11-27 T 16:40(UTC时)河北张家口市桥东区一化工厂附近发生的严重爆炸起火事件以及2019-02-28澳大利亚西南部发生的火灾事件,均快速有效的检测到了火点。实验表明,改进的火点检测算法能很好的进行火点检测,并能解决晨昏交界、冰雪下垫面、常规火源点、太阳耀光等火点检测的难题。  相似文献   

城市“腾笼换鸟”产业转型升级的困境与出路——基于北京市的调查研究     

王颖  段霞  吴康 《地理科学》2020,40(5):786-792

剖析北京“腾笼换鸟”产业转型升级的现状,在针对“新鸟”进笼“老鸟”去哪问题,“老笼”空间结构合理优化问题,“老鸟”涅槃“新鸟”培育问题,“老鸟”和“新笼”承接问题的分析基础上,提出应完善京津冀产业链协作,“腾笼换鸟”拓展区域联动发展空间;东西城合并成首都特区,“腾笼换鸟”保护北京历史文化名城;切实发挥政府的职能作用,“腾笼换鸟”提升企业内在动力机制;规划建设“微中心”小城镇,“腾笼换鸟”促进北京人口有效疏解的建议。  相似文献   

1966—2018年河西走廊春季寒潮频次及影响因子分析          下载免费PDF全文

吕青松  辛存林  王瑞  张勃  张华 《干旱区地理》2020,43(4):946-954

基于甘肃河西走廊15个国家基准基本站点提供的1966年1月1日~2018年12月31日春季（3~5月）逐日最低气温值数据,利用线性回归分析方法、Spline空间插值法和Mann-Kendall趋势和突变检验法,探讨了河西走廊地区1966—2018年春季（3~5月）寒潮频次变化及其影响因子。研究表明：（1） 1966—2018年河西走廊地区春季（3~5月）单站寒潮频次总体呈现下降趋势〔–0.098次·（10 a）^–1〕,其中1980—2010年寒潮频次呈显著下降趋势,2010年之后下降趋势变缓,未通过显著性检验;区域寒潮53 a来频次总体呈缓慢下降趋势〔–0.015次·（10 a）^–1〕。（2） 近53 a河西走廊地区春季三个月中,单站寒潮总量4月>3月>5月,其中4月、5月寒潮频次下降不明显,3月频次下降显著。（3） 空间上,大致以北大河和黑河干流为界,两河中间区域春季寒潮频次低,而北大河以北和黑河干流以东区域则是寒潮高值区,走廊外围地区寒潮频次较高,且大多呈显著下降趋势,寒潮频次与气温距平存在明显负相关关系,内部地区变化趋势不明显。（4） 河西地区春季寒潮频次受气候变暖、地形和大气环流的影响,寒潮频次变化趋势存在地区差异。研究可提高对甘肃河西走廊寒潮演化过程的认知,为河西走廊气候变化的进一步研究奠定基础。  相似文献   

毛乌素沙地东南缘全新世气候不稳定性     

李想  苏志珠  马义娟  张彩霞  柳苗苗 《中国沙漠》2020,40(2):109-117

毛乌素沙地处于中国季风区的西北边缘,对气候变化敏感,是古气候演化和古环境变迁研究的理想场所。全新世气温回暖期,气候整体稳定但存在暖湿冷干的次级波动,不同类型沉积相组成的地层序列记录了全新世千年、百年尺度的气候不稳定性。对毛乌素沙地不同沉积序列的研究表明:(1)毛乌素沙地东南缘的沙漠/黄土边界带DLT、NB剖面粒度及Rb/Sr比值可作为气候变化的代用指标。>63 μm砂含量增加,Rb/Sr比值降低,揭示东亚冬季风增强,风成砂堆积,风沙活动加剧,气候干冷;反之,东亚夏季风增强,降水增多,风化成壤程度加大,气候暖湿。(2)毛乌素沙地东南缘沙漠/黄土边界带及低洼沟谷区年代概率密度函数的分析显示,全新世早期~8\^5 ka BP风成砂堆积,风沙活动加剧,气候表现为干冷;全新世中期8\^5 ka~3\^0 ka BP地势较高处古土壤广泛发育,地势低洼处发育湖沼相/泥炭层,但在6.0 ka BP左右存在气候转冷过程;全新世晚期3\^0 ka BP以来气候波动频繁。(3)毛乌素沙地气候变化与北大西洋冰筏冷事件、董哥洞石笋、季风边缘区湖泊、东部沙地等记录具有可对比性和相对一致性,反映出全新世千年和百年尺度上的气候变化主要受全球变化的影响,这对预测和模拟未来气候变化具有借鉴意义。  相似文献   

气候变化下基于物种与植物耦联机制的小白额雁(Anser erythropus)潜在生境及适宜性研究--以长江流域中下游为例     

向伶  高翔  彭玉辉  梁婕 《资源与生态学报(英文版)》2020,11(2):140-149

气候变化和人类活动通过改变物种生境而影响物种多样性。小白额雁是长江流域中下游的一种具有较高生态价值的食草型濒危候鸟,受气候变化和人类活动威胁。本文以小白额雁为代表性物种,定量分析了气候变化对长江流域中下游候鸟潜在生境及适宜性空间分布格局的影响。采用Maxent模型模拟了当前情景和全球环流模型(GCMs)气候场景下小白额雁潜在生境及其适宜性分布。研究结果表明,小白额雁分布特征与其栖息地周边植物分布呈显著相关关系;运用Maxent模型模拟小白额雁六种主要食源植物的分布特征,并将其结果作为环境变量,将显著改善小白额雁潜在生境及其适宜性模型的模拟性能;在两种典型浓度情景(RCP 2.6和RCP8.5)下,2070年小白额雁潜在生境适宜性面积将下降。为应对气候变化对小白额雁的影响,应采取更加合理的管理措施和保护政策,包括调整保护区的大小、形状和用途。  相似文献   

Late Quaternary sea‐level change and evolution of Belfast Lough,Northern Ireland: new offshore evidence and implications for sea‐level reconstruction     

Ruth M. K. Plets  S. Louise Callard  J. Andrew G. Cooper  Joseph T. Kelley  Daniel F. Belknap  Robin J. Edwards  Antony J. Long  Rory J. Quinn  Derek W. T. Jackson 《第四纪科学杂志》2019,34(4-5):285-298

The interplay of eustatic and isostatic factors causes complex relative sea‐level (RSL) histories, particularly in paraglacial settings. In this context the past record of RSL is important in understanding ice‐sheet history, earth rheology and resulting glacio‐isostatic adjustment. Field data to develop sea‐level reconstructions are often limited to shallow depths and uncertainty exists as to the veracity of modelled sea‐level curves. We use seismic stratigraphy, 39 vibrocores and 26 radiocarbon dates to investigate the deglacial history of Belfast Lough, Northern Ireland, and reconstruct past RSL. A typical sequence of till, glacimarine and Holocene sediments is preserved. Two sea‐level lowstands (both max. ?40 m) are recorded at c. 13.5 and 11.5k cal a bp . Each is followed by a rapid transgression and subsequent periods of RSL stability. The first transgression coincides temporally with a late stage of Meltwater Pulse 1a and the RSL stability occurred between c. 13.0 and c. 12.2k cal a bp (Younger Dryas). The second still/slowstand occurred between c. 10.3 and c. 11.5k cal a bp . Our data provide constraints on the direction and timing of RSL change during deglaciation. Application of the Depth of Closure concept adds an error term to sea‐level reconstructions based on seismic stratigraphic reconstructions.  相似文献