Reconstructing the climate and environmental characteristics of Chengdu in the Middle Tang Dynasty from Du Fu's "Thousand Autumns of Snows"
-
摘要:
解析唐诗记载的地理景观和物候信息可重建诗人所处时代的气候环境特征。本文以唐代诗人杜甫 (公元712—770年)的《绝句四首·其三》和张籍(公元767—830年)的《成都曲》为例,选取诗中“千秋雪”和“荔枝熟”作为分析要素,结合物候学及地理信息系统 (GIS)等手段定量重建中唐时期成都地区的温度、降水和雪线。基于成都处于唐朝荔枝北界这一前提,通过GIS对现代荔枝分布点进行空间分析,初步限定中唐时期成都地区的最低年均温和年降水量分别为 (18 ± 0.4) ℃和 (
1230 ± 260) mm,较现代分别高约2 ℃和250 mm;进一步依据现代雪线、温度和降水的关系,通过GIS空间分析和理论计算推知中唐时期成都地区的雪线高度为 (4600 ± 200) m,较现代低约700 m。最后,通过重建的雪线高度及GIS的可见性分析,推测杜甫诗中“千秋雪”最有可能是九顶山的主峰狮子王峰。Abstract:Through the analysis of the geographical landscapes and climatic information recorded in Tang poems, it is possible to reconstruct the characteristics of the climate and environments of the times in which the poets lived. In this paper, we take the "Four Quatrains (No. 3)" by Du Fu (712—770 AD) and "Ode to Chengdu" by Zhang Ji (767—830 AD) of the Tang Dynasty to gain insight into the climatic conditions of their time, and select "Thousand Autumns of Snows" and "Lychee Ripening" as analyzing elements, then combine them with phenology and GIS to quantitatively reconstruct the temperature, precipitation, and snowline of the Chengdu during the Middle Tang Dynasty. Based on the premise that Chengdu was located at the northern boundary of lychee cultivation during the Tang Dynasty, a spatial analysis of modern lychee distribution points was conducted through GIS. The analysis initially identified the minimum annual mean temperature and annual precipitation of Chengdu during the Middle Tang Dynasty to be (18 ± 0.4) ℃ and (
1230 ± 260) mm, which were about 2 ℃ and 250 mm higher than those in modern times. Furthermore, based on the relationship between the modern snowline, temperature and precipitation, the snowline height in Chengdu during the Middle Tang Dynasty is estimated to be (4600 ± 200) m, which is about 700 meters lower than that in modern times. Based on the reconstructed snowline height and GIS visibility analysis, it is hypothesized that the 'Thousand Autumns of Snows' in Du Fu's poem refers to the main peak of Jiu Ding Mountains, Lion King Peak.-
Keywords:
- northern boundary of Lychee /
- snowline /
- Tang poetry /
- Middle Tang Dynasty /
- Chengdu
-
天然气水合物作为未来战略能源已成为国际共识,受到世界各国尤其是发达国家和能源短缺国家的高度重视。近年来,国际天然气水合物研发态势已逐步从资源普查阶段转入勘查试采一体化阶段,迄今已在加拿大麦肯齐三角洲、美国阿拉斯加北坡、中国青海祁连山3个陆地冻土区和日本南海海槽、中国南海神狐2个海域成功实施试采[1-2]。2013年、2017年,日本南海海槽砂质储层天然气水合物试采结果表明,大量出砂是试采被迫终止的主要原因[3-4]。2017年、2020年中国在南海神狐海域分别运用垂直井和水平井成功实施了针对泥质粉砂储层的天然气水合物试采,2017年连续稳定生产60天,累计产气30.9×104 m3,2020年连续稳定生产30天,累计产气86.14×104 m3,取得了理论与技术的重大突破[5-6]。
然而,上述海域天然气水合物试采实践暴露出的部分工程、地质难题仍悬而未决,严重制约了后续试采方案的制定和产业化进程。其中,工程地质风险调控领域的重大未解难题是:日本2013年试采和2017年第一口井试采均因储层大量出砂堵塞管道而被迫终止[2],但中国南海北部神狐海域泥质粉砂天然气水合物储层在产气的同时出砂却较少[5,7],给外界造成很多困惑。初步研究认为,造成这种现象的主要原因可能是储层在开采过程中发生蠕变,封闭原有产砂通道(即出砂蚯蚓洞、裂隙在蠕变过程中逐渐闭合)[8-9]。但是,天然气水合物储层蠕变特征对出砂规律的量化控制程度如何?目前暂时悬而未决。解决该问题的前提是厘清水合物开采过程中储层蠕变主控因素及其控制机理。本文将聚焦南海北部神狐海域泥质粉砂型天然气水合物储层特征,探讨不同的工程、地质因素对储层蠕变行为的影响规律,深入分析水合物储层的蠕变特征参数,阐明水合物降压开采过程中多相多场耦合条件下储层蠕变特征的响应机理及主控因素,为解决水合物试采工程发现的工程、地质问题提供科学依据。
另外,在综述海洋天然气水合物储层破坏特征的基础上,梳理海洋天然气水合物储层蠕变特征及主控因素,厘清关键科学问题;结合最新研究成果,阐述天然气水合物储层蠕变特征多尺度表征与探测技术体系的基本内涵,简要探讨该领域的未来研究方向。
1. 水合物储层蠕变行为的主控因素与科学挑战
1.1 水合物储层蠕变行为的主控因素
天然气水合物开采过程中的储层力学动态响应特征是揭示开采诱发的工程地质风险发生临界条件、演变规律的基础。天然气水合物开采过程中储层(特别是未固结成岩储层)的破坏过程在很大程度上属于长期载荷作用下的缓慢变形和破坏,即沉积物在开采过程中主要发生流变性破坏,蠕变参数是储层长期稳定性评价的重要指标。
从微观角度分析,蠕变伴随着沉积物内部微观结构随时间的变化、调整,主要表现为微裂纹、微压实等损伤形式的不断累积、发展和扩散。因此,蠕变本身是一个动态过程,其涉及的蠕变应变量、蠕变应变速率、蠕变破坏强度等参数也均为动态指标。基于损伤带的演化特征,大部分岩土材料的蠕变过程包含衰减蠕变、稳态蠕变和加速蠕变3个基本阶段[10]。然而,初步研究结果显示,含水合物沉积物在三轴应力状态下的蠕变特征曲线并不完全遵循常规意义上的3个基本蠕变阶段[11],不同沉积物本底特征的含水合物沉积物储层表现出不同的蠕变规律。
导致上述差异的第一个重要影响因素是水合物储层的相变(天然气水合物分解)。常规岩土材料蠕变特征分析中,沉积物内部不存在动态相变过程,也不存在水气两相动态渗流过程。因此,我们将这种单纯外应力作用下沉积物内部损伤累积引起的沉积物流变行为定义为“静态”蠕变。但天然气水合物开采过程中存在相变(即天然气水合物分解相变),天然气水合物分解产生的水气处于渗流流动状态,而天然气水合物相变和多相渗流过程本身会对储层所受的应力状态产生影响,因而储层所受的应力状态处于实时动态变化中[12]。因此,天然气水合物开采过程中储层的变形行为受相变-渗流与应力状态实时耦合的影响,比单纯应力作用下的岩土蠕变特征复杂得多[13]。我们将这种应力状态与实时相变、多相渗流过程动态耦合作用下的流变特征定义为“动态”蠕变。
为克服储层相变对储层基础物性研究带来的挑战,国内外学者曾尝试应用冰的结融过程替代天然气水合物生成分解过程,以简化对天然气水合物储层基本性质的研究[14]。如现在被天然气水合物界普遍接受并认可的饱和度阿尔奇公式及其饱和度指数取值(n≈2),其实就是基于含冰沉积物模拟结果获取的[15]。由此推断,冻土融冰条件下的蠕变行为可为天然气水合物储层蠕变特性的研究提供一定的思路借鉴。在相同应力路径下,衰减蠕变、稳定蠕变及加速蠕变阶段均受融冰引起的岩石劣化程度的影响,具有典型的非线性特征,而这种非线性行为可以通过冻融损伤效应来定量评价[16]。另外,冻融循环导致冻土层从稳定蠕变到加速破坏所经历的时间更短,对应的蠕变量及蠕变速率增大,且宏观破裂形态由单一斜剪切破坏模式向共轭断面拉剪复合破坏模式演化[17-18]。
在不恰当的井底工况条件下,地层中可能存在多轮次的天然气水合物分解-生成-再分解过程,这一过程类似于冻土层中的结融冰循环过程。因此,上述冻融循环载荷对沉积物蠕变的影响对天然气水合物储层蠕变特性分析具有重要的启示意义。但冰结融过程毕竟与天然气水合物合成-分解过程存在巨大的差异:天然气水合物合成分解过程中的潜热效应(62.83 kJ/mol)远大于冰结融过程的潜热(6.0 kJ/mol),导致天然气水合物开采过程中储层温度剧烈变化。因此,天然气水合物储层蠕变规律建模不仅需要考虑相变的影响,而且要考虑相变引起的温度变化。
针对常规不含天然气水合物的黏土沉积物的蠕变试验结果表明:温度的升高会提升黏土储层的蠕变速率,延长衰减蠕变阶段的时间[19]。在降压开采过程中,由于天然气水合物分解吸热和多相渗流Joule-Thomson效应的联合作用,近井地层的温度急剧降低[20],可能导致储层结冰或天然气水合物二次生成[21-22]。因此,天然气水合物开采条件下的储层蠕变行为受相变过程与温度效应的耦合控制,天然气水合物储层蠕变特征分析必须首先考虑相变与温度效应的耦合影响[23-24]。
储层中的气水分布状态和两相渗流也可能直接影响水合物储层的蠕变行为。原位状态下海洋天然气水合物储层处于水饱和状态,水对储层产生物理、化学及力学作用,进而在微观上改变储层的强度等宏观力学性质[25-26]。天然气水合物开采过程中,储层静止的饱水状态将被打破,沉积物孔隙被处于渗流状态的气水两相饱和[27]。因此,天然气水合物开采环境下储层要承受静水压力及多相渗流产生的动水冲刷联合作用,其流变特性与静止水饱和状态的常规储层存在本质差别。
综上,海洋天然气水合物开采过程中储层的变形行为取决于水合物分解相变、储层温度扰动、储层多相渗流、井底工作制度(即外部载荷)等因素的综合作用。厘清这些主要影响因素及其耦合作用对储层蠕变行为的控制机理,是建立水合物储层蠕变本构模型,精准预测水合物储层变性行为的关键。
1.2 水合物储层蠕变本构建模思路
为了综合描述各因素对岩土蠕变效应的影响,国内外研究者从不同的角度建立了非天然气水合物储层的蠕变本构模型。常用的建模途径主要有两种:一是基于试验结果建立应力-应变-时间经验模型,二是利用基本的虎克弹性体、牛顿黏性体、圣-维南塑性体等力学模型元件,通过基本元件的串并联建立线性蠕变模型[28]。由于线性模型无法描述沉积物的加速蠕变过程,因此,前人尝试采用各种类型的非线性元件替代线性元件,进而建立了各种类型的黏塑性损伤本构模型。
然而,上述常规岩土蠕变本构模型均未考虑天然气水合物开采过程中可能存在的由相变-渗流-传热-应力综合作用引起的蠕变硬化、蠕变损伤、热损伤、动水冲刷效应,不能直接用于评价天然气水合物储层的损伤发展规律,更不能准确预测不同蠕变阶段的转化时机。因此,前人针对冻土和常规岩层建立的蠕变本构,一定程度上为天然气水合物开采储层的蠕变分析提供了借鉴,但现有模型均不能完全反映天然气水合物开采过程中储层的相变-渗流-传热-应力多场多相多组分强耦合过程,因此,不能照搬现有模型和方法对天然气水合物开采储层的蠕变行为进行预测。
1.3 关键科学挑战
天然气水合物开采过程中,储层受到应力、相变、温压场、井底工作制度变化的影响,而这些影响因素本身与时间存在复杂的对应关系。这是开采条件下天然气水合物储层蠕变行为区别于常规岩土材料的最显著特点。
与常规不含天然气水合物沉积物相比,天然气水合物的存在会对储层蠕变特性产生哪些独特影响?天然气水合物降压开采动态过程又会对储层蠕变行为产生哪些影响?这两个问题是理解天然气水合物开采过程中储层工程地质风险诱发机制的总开关,也是水合物储层蠕变行为研究的关键科学挑战。回答这两个问题,对于准确评价特定天然气水合物储层的蠕变特征、识别天然气水合物开采条件下的潜在工程地质风险并进行有效防控具有重要意义。具体阐述如下:
(1)天然气水合物的存在及其空间分布特征对含天然气水合物沉积物蠕变特征的影响可视为“静态因素”,是导致含天然气水合物储层蠕变参数演化特征区别于常规砂土沉积物蠕变特征的最根本因素。那么,含天然气水合物沉积物在特定应力加载条件下的蠕变演化特征究竟会不会像常规砂土沉积物那样遵守严格的衰减蠕变、稳态蠕变和加速蠕变模式呢?如果不是,这其中天然气水合物颗粒又扮演何种角色?具体而言,含天然气水合物沉积物是一个包含天然气水合物、沉积物两种固相和孔隙水、气等流体的多元体系。蠕变变形过程中与天然气水合物本身直接相关的因素如天然气水合物颗粒本身的碎裂、天然气水合物颗粒与沉积物颗粒之间的相对滑动、天然气水合物颗粒本身的压缩、沉积物颗粒在天然气水合物限制下发生不规则翻转错动等,这些因素都可能导致沉积物体系宏观蠕变应变、蠕变变形量及蠕变速率的改变(图1),进而影响含天然气水合物沉积物的蠕变参数演化模式。
(2)蠕变本身是一个动态过程,降压开采条件下的储层相变及伴随而来的储层多相渗流、相变潜热、应力状态改变等均为相互依存且与时间相关的必然性动态过程。因此,降压开采对含天然气水合物沉积物体系蠕变特征的影响可以视为“动态因素”。理论上而言,降压路径本身对蠕变的影响可以基于有效应力原理等效转化为储层有效应力改变的影响,而应力路径本身改变对储层蠕变行为预测造成的挑战较为有限,因此,降压路径本身带来的挑战完全可以克服。最大的挑战是降压开采过程中相变-传热-渗流与应力的多场双向耦合过程(图2)。揭示天然气水合物开采过程中多场耦合效应对蠕变行为的控制机理,是回答天然气水合物开采过程中储层潜在工程地质风险发生、发展规律的关键。
2. 蠕变行为刻画方法展望
目前,国内外获取天然气水合物储层力学参数的主要技术手段有三轴剪切[29-30]、直剪[31]、静力触探[32]等。这些技术手段主要以快速载荷作用下沉积物发生破坏作为评价含天然气水合物沉积物体系力学性能的依据,对长期开采条件下储层缓慢变形的时效性行为的分析与预测涉及较少[33-34]。由于天然气水合物对沉积物骨架颗粒的胶结、填充作用,含天然气水合物沉积物在蠕变破坏过程中的微元损伤过程可视为由线弹性应力状态向非线性应力状态的转化[35-36],具体的损伤演化形式可能包括水合物晶体结构破坏、沉积物孔隙结构变形、含水合物沉积物压实带或微裂隙等,不同的损伤演化形态所适应的模拟与探测手段不同。
总之,宏观尺度的变形破坏源自细观甚至微观尺度的损伤积累,跨尺度模拟与探测是研究天然气水合物开采储层蠕变的最佳路径,而试样内部的微元损伤实时状态及其发展演变速率则是含天然气水合物沉积物蠕变形为跨尺度表征的纽带。
从分子尺度,分子动力学模拟为水合物体系中微观结构的变化预测与行为研究提供了重要手段[37-38]。分子动力学模拟可通过研究沉积物中水合物自身微观结构(如水合物晶体颗粒的大小、颗粒之间的结构、颗粒之间是否存在气泡等,图3)和沉积物与水合物之间的微观接触行为(如氢键的相互作用、矿物表面的吸附、甲烷气泡等)来获得蠕变行为的物理化学机制。初步研究表明,天然气水合物颗粒的蠕变主要表现为晶体溶解与晶体界面滑移,主要的影响因素包括晶体颗粒大小、外界温度与加载应力路径等[39-41]。
从纳微米尺度,X-CT、LF-NMR等微观探测技术可以为天然气水合物沉积物蠕变过程中的损伤动态探测提供有效支撑。目前国内外已基于低场核磁共振(LF-NMR)开展了水合物动态聚散过程中黏土质粉砂沉积物孔隙结构变化特征的探测,获得了不同水合物饱和度条件下样品横向弛豫时间谱、孔隙尺寸信息的对应关系、横向弛豫率与水合物饱和度的关系[42]。初步结果表明,小孔隙里的水合物优先分解,孔隙尺寸变大,核磁信号量逐渐增加;在一定水合物饱和度范围内(比如小于45%的情况),随着水合物饱和度的降低,横向弛豫率小幅减小[42-43]。实际上,除了指示水合物的优先合成/分解位置,LF-NMR横向弛豫谱(T2)也可以得到沉积物样品中的孔隙、裂隙、破碎带分布[18]。为此,笔者团队提出了基于LF-NMR探测水合物沉积物中微元损伤演化的装置(图4)与方法[44],其基本的测量原理是:将低场核磁共振探测和三轴应力加载相结合,在低场核磁共振平台上同时获取水合物储层的宏观力学数据和加载过程中的微观孔隙结构演化特征,为跨尺度研究天然气水合物储层的力学性能及其变形机理提供技术保障。
另外,将计算机断层扫描技术(X-CT)与三轴剪切技术结合,能够同时获得应力-应变参数和加载过程中的微观破碎带结构演化特征[45]。笔者运用X-CT三轴测试技术,清晰地揭示了天然气水合物沉积物样品在三轴压缩过程中的微观结构演化规律,随着水合物饱和度的增大,剪切带的宽度逐渐减小,倾角逐渐增大[30],这是由于随着饱和度的增大,水合物对骨架的支撑作用增强,沉积物极限强度增大,提高了样品的抗剪切能力,剪切带宽度减小。
从介观尺度(cm),开展水合物开采过程中沉积物变形行为的模拟实验,可以分析水合物饱和度、有效围压、偏应力水平对水合物开采过程中的沉积物变形量的影响(图5)。初步结果表明,水合物分解过程中试样的蠕变大致分为3个阶段:第1阶段为变形初始阶段,沉积物试样在围压的作用下等向固结,随着时间的推移,固结变形趋于稳定;第2阶段为加剧变形阶段,回压开始降低,试样中气压降低,随着水合物的分解,土体骨架刚度变弱,在偏应力作用下,产生体变;第3阶段为变形的稳定阶段,土体发生急剧变形,随着气压的进一步降低,水合物分解完成。
除了从应变-时间关系直观了解沉积物的蠕变特征,还可以引入一些损伤探测手段与剪切破坏实验相结合,实现从介观尺度表征沉积物内部损伤的目的。声波探测在含天然气水合物沉积物模拟实验中已有非常成熟的应用[24]。蠕变变形过程中,沉积物内部会发生微观结构的变化,即微裂隙的产生及发育,会使得沉积物的声波特征发生改变,因而沉积物内部结构的变化情况可能通过某一个或者多个宏观声波特征参数反映出来。因此,通过探测不同蠕变状态下声波衰减及波形变化规律,进一步建立声波特征与沉积物损伤动参数间的关系,能够为沉积物岩心尺度的蠕变动态演化分析提供支撑。
从中试尺度(m),分布式光纤感测技术作为一种广泛应用于野外监测和物模试验的在线监测技术[46],能够直接测量模型中的温度场、应变场,间接获取应力场、渗流场、水饱和度场[47],在岩土沉降、边坡灾害分析等方面已有非常广泛的应用。理论上而言,该技术引入天然气水合物蠕变分析模型试验,达到天然气水合物饱和度场、储层应力场、温度场、渗流场、变形量场一体化的监测目的,从而为揭示天然气水合物开采储层多相多场强耦合条件下储层的蠕变机制提供依据。但目前仍需要从基础理论、硬件选型、数据解调等多方面开展大量的基础研究工作。
对于现场尺度(km),数值模拟技术可以利用大尺度模型模拟天然气水合物开采过程中的传热(T)、渗流(H)、相变(C)和应力响应(M)及其耦合特征,为矿场尺度下天然气水合物长期开采储层蠕变特征和工程地质风险评价提供重要手段。目前,国内外研发了众多的天然气水合物开采THMC数值模拟器,并对开采过程中的THMC耦合行为进行了研究[48-49]。然而,现有的数值模拟器还缺少水合物开采蠕变过程的模拟能力,本质原因在于天然气水合物储层蠕变的本构模型的不完善,且水合物开采传热、渗流、相变等过程与蠕变行为之间的耦合关系也仍未明确,导致无法建立有效的数学模型描述长期开采条件下的水合物储层蠕变行为。因此,需要开展系统的实验研究,建立水合物储层的蠕变本构模型和水合物对蠕变参数的影响模型,从而建立有效的数学模型,发展相应的数值模拟器,对水合物长期开采条件下的蠕变行为进行预测。
综上所述,在天然气水合物测试中已发挥重要作用的纳微米探测(如X-CT、LF-NMR)技术和岩心尺度的实验模拟技术为天然气水合物储层蠕变时效性参数评价提供了重要的支撑,光纤感测技术也有望在中试尺度为揭示天然气水合物多场耦合条件下的储层蠕变机理提供一定的支撑,分子模拟和数值模拟技术则能分别从更微观(分子尺度)和更宏观(矿场尺度)的尺度弥补物理探测与分析技术的不足,为天然气水合物储层蠕变行为研究提供重要的研究手段,但目前基于上述先进技术进行天然气水合物储层蠕变特征评价处于起步阶段,亟需形成分子-纳微米-厘米-米-千米级跨尺度的天然气水合物储层蠕变评价技术体系。
3. 结论与展望
(1)本文首次提出储层蠕变现象在水合物开采过程中不仅是客观存在的,而且对工程地质灾害的发生发展都有非常重要的作用,摸清海洋天然气水合物储层力学性质及其随水合物开采过程的动态变化规律,是准确认识天然气水合物开采过程中工程地质风险发生、发展及其潜在危害的前提。海洋天然气水合物储层的蠕变行为因储层地质背景、开采模式而异,厘清这些因素是评价特定储层蠕变特征、识别水合物开采的潜在工程地质风险的前提。特别是对于南海北部神狐海域泥质粉砂型储层而言,其蠕变过程与砂质储层存在本质区别,是导致两者试采过程中工程地质风险演化差异的根本原因。
(2)天然气水合物的存在及其空间分布特征对含天然气水合物沉积物蠕变特征的影响可视为“静态因素”,是导致含天然气水合物储层蠕变参数演化特征区别于常规砂土沉积物蠕变特征的最根本因素。蠕变本身是一个动态过程,降压开采条件下的相变(天然气水合物分解)及伴随而来的储层多相渗流、相变潜热、应力状态改变等均为相互依存且与时间相关的必然性动态过程。因此,降压开采对含天然气水合物沉积物体系蠕变特征的影响可以视为“动态因素”。
(3)针对上述问题,我们需要从天然气水合物储层的力学性能弱化特征及蠕变各阶段的时效参数两方面入手,以我国南海北部天然气水合物储层样品为重点研究对象,将LF-NMR、X-CT、超声探测、光纤多场监测技术与三轴蠕变控制技术有机结合,并利用分子模拟和数值模拟手段向更微观和更宏观两个尺度拓展,建立天然气水合物储层蠕变行为的跨尺度研究方法体系,探索实际天然气水合物开采储层蠕变行为的主控因素;最终建立适合南海北部储层的、考虑天然气水合物开采过程中相变-渗流-潜热-应力综合作用的蠕变本构,从而为天然气水合物开采条件下储层流变学理论框架的建立奠定基础,并为天然气水合物开采工程安全设计提供理论支撑。
1 在成都遥望雪山. 由两幅照片看杜甫“窗含西岭千秋雪”的指向——成都的自然人文之旅 [缩略版]. https://mp.weixin.qq.com/s/QKaOFQYBwHyfb1196ghfEQ, 2022-05-02. -
图 2 现代荔枝北界及中唐成都气候相似区
a:中国现代荔枝北界及中唐成都气候相似区,b:西藏自治区山南市的中唐成都气候相似区。
Figure 2. Modern northern boundary of lychee growing and similar climatic regions in ancient Chengdu area during the Middle Tang Dynasty
a: Modern northern Boundary of lychee in China and the similar climatic regions in ancient Chengdu area during the Middle Tang Dynasty, b: The similar climatic regions of ancient Chengdu in the Middle Tang Dynasty in Shannan City, Xizang Autonomous Region.
图 3 视域分析栅格图
a:基于杜甫草堂观测点进行视域分析后的结果图,视域分析结果为“可见”与“不可见”,另外图中给出了海拔高于4600 m的区域范围、各个山峰的海拔高度及距杜甫草堂的距离;b:九顶山狮子王峰区域的视域分析结果图,图中给出了海拔高于4600 m且可见区域的可见次数,九顶山狮子王峰最多被看见15次;c:杜甫草堂西北角观测位置示意图(底图来源于百度地图)。
Figure 3. Raster map of viewshed analysis
a: The result of viewshed analysis based on the observation point of Du Fu Thatched Cottage, which can be classified into "visible" and "invisible" regions. In addition, the figure also shows the regions with an altitude above 4600 m, the altitude of each mountain and the distance from Du Fu Thatched Cottage, b: The viewshed analysis results of the Lion King Peak area of Jiu Ding Mountains. The figure shows the number of visible areas with an altitude above 4600 m, and the Lion King Peak of Jiu Ding Mountains is seen for a maximum of 15 times, c: Observation position map of the northwest corner of Du Fu Thatched Cottage ( base map is from Baidu map).
表 1 研究数据及其来源汇总
Table 1 Data for this research and their sources
数据名称 数据来源 荔枝分布点数据 《中国果树志 荔枝卷》[15],四川省植物资源信息网(http://www.scpri.ac.cn),《四川荔枝种植公布的历史考证》[30] 现代荔枝北界所处地2016—2020年年均温、年降水量数据 四川统计年鉴,来源于四川省统计局 (http://tjj.sc.gov.cn)[31-35] 中国2022年地市行政区划边界数据 中国多年度地市行政区划边界数据,数据来源于资源环境科学数据注册与出版系统 (https://www.resdc.cn)[36] 中国2022年省级行政区划边界数据 中国多年度省级行政区划边界数据,数据来源于资源环境科学数据注册与出版系统 (https://www.resdc.cn)[37] 2016—2020年中国1km分辨率年降水量数据 中国1km分辨率年降水量数据 (1982—2022年),数据来源于国家科技资源共享服务平台—国家地球系统科学数据中心 (http://www.geodata.cn) 2016—2020年中国1km分辨率年平均气温数据 中国1km分辨率年平均气温数据 (1982—2022年),数据来源于国家科技资源共享服务平台—国家地球系统科学数据中心 (http://www.geodata.cn) 现代雪线高程数据 《中国山河全书》[38],《青藏高原地区雪线时空变化动态研究》[39] GDEMV3 30M分辨率数字高程数据 中国科学院计算机网络信息中心地理空间数据云平台(https://www.gscloud.cn) 成都西面群山山峰点数据 基于谷歌地球卫星影像对山峰拾取坐标点 注:本研究所有数据采用统一的GCS_WGS_1984坐标系。 表 2 《全唐诗》中四川荔枝诗
Table 2 Lychee of Sichuan mentioned in the poems of Tang Dynasty
地点 时间 作者 诗名 提及“荔枝”的诗句 文献来源 成都 中唐 张籍 成都曲 “锦江近西烟水绿,新雨山头荔枝熟” 中唐 李端 送何兆下第还蜀 “袅猿枫子落,过雨荔枝香。劝尔成都住,文翁有草堂” [52] 779年 卢纶 送从舅成都县丞广归蜀 “晚程椒瘴热,野饭荔枝阴” 乐山 中唐 卢纶 送张郎中还蜀歌 “邛竹笋长椒瘴起,荔枝花发杜鹃鸣” 807年 薛涛 忆荔枝 “近有青衣连楚水,素浆还得类琼浆” [52] 866年 薛能 荔枝诗 “颗如松子色如樱,未识蹉跎欲半生” 宜宾 765年 杜甫 宴戎州杨使君东楼 “重碧拈春酒,轻红擘荔枝” [52] 泸州 767年 杜甫 解闷十二首·其十 “忆过泸戎摘荔枝,青峰隐映石逶迤” 891年 郑谷 荔枝树 “肠断渝泸霜霰薄,不教叶似灞陵红” [52] 891年 郑谷 将之泸郡旅次遂州遇裴晤员外谪居于此话旧凄凉因寄二首·其二 “我拜师门更南去,荔枝春熟向渝泸” 万县
(唐朝忠州)819年 白居易 重寄荔枝与杨使君,时闻杨使君欲种植故有落句之戏 “香连翠叶真堪画,红透青笼实可怜” 819年 白居易 种荔枝 “十年结子知谁在,自向庭中种荔枝” [52] 820年 白居易 荔枝楼对酒 “荔枝新熟鸡冠色,烧酒初开琥珀香” 表 3 现代荔枝分布点信息
Table 3 Information of modern lychee distribution sites
序号 地名 地区 位置 年均温/℃ 年降水量/mm 1 洪雅县 四川省眉山市 29.90°N、103.37°E 16.6 1435.5 2 峨眉山市 四川省乐山市 29.60°N、103.48°E 17.2 1555.3 3 犍为县 四川省乐山市 29.21°N、103.94°E 17.5 1141.3 4 市中区 四川省乐山市 29.55°N、103.76°E / 1082.0 5 沙湾区 四川省乐山市 29.41°N、103.54°E / 1530.3 6 夹江县 四川省乐山市 29.74°N、103.57°E 17.1 1428.4 7 屏山县 四川省宜宾市 28.83°N、104.34°E 18.2 1209.8 8 江安县 四川省宜宾市 28.72°N、105.06°E 18.1 1132.0 9 筠连县 四川省宜宾市 28.16°N、104.50°E 17.6 / 10 泸县 四川省泸州市 29.15°N、105.37°E 18.1 1179.4 11 纳溪区 四川省泸州市 28.77°N、105.36°E 17.4 1150.8 12 叙永县 四川省泸州市 28.15°N、105.44°E 17.9 1172.6 13 凤鸣镇 四川省泸州市合江县 28.71°N、105.88°E 18.2
(合江县年均温)1134.9
(合江县年降水量)14 榕山镇 四川省泸州市合江县 28.86°N、105.91°E 18.2
(合江县年均温)1134.9
(合江县年降水量)15 大桥镇 四川省泸州市合江县 28.83°N、105.71°E 18.2
(合江县年均温)1134.9
(合江县年降水量)16 福宝镇 四川省泸州市合江县 28.78°N、106.08°E 18.2
(合江县年均温)1134.9
(合江县年降水量)17 盐边县 四川省攀枝花市 26.68°N、101.85°E 19.2 1065.6 18 米易县 四川省攀枝花市 26.90°N、102.10°E 19.7 1112.6 19 德昌县 四川省凉山彝族自治州 27.40°N、102.17°E 17.7 1049.0 20 宁南县 四川省凉山彝族自治州 27.06°N、102.75°E 19.3 960.0 21 武陵镇 重庆市万州区 30.50°N、108.25°E 18.2 1155.8 注:表中位置来源于百度地图拾取坐标系统https://api.map.baidu.com/lbsapi/getpoint/,并已经将百度地图的BD09坐标系转换为WGS84坐标系;表中的年均温、年降水量数据来源于各地政府网及四川省公共气象服务网https://www.scggqx.com。 表 4 现代荔枝北界所处地2016—2020年年均温和年降水量
Table 4 Averages of annual temperature and annual precipitation of the sites where the modern northern boundary of lychee are located during 2016—2020
年份 年均温/℃ 年降水量/mm 泸州市 乐山市 宜宾市 眉山市 泸州市 乐山市 宜宾市 眉山市 2020 18.0 18.3 17.7 18.2 1396.8 1555.3 1746.0 1292.6 2019 18.1 18.2 17.5 17.7 1036.3 1483.1 1130.5 1167.9 2018 18.2 18.5 17.7 18.1 1388.4 1511.9 1190.0 924.6 2017 18.2 18.6 19.1 18.2 975.2 849.5 823.5 1083.1 2016 18.3 18.6 19.0 18.4 1443.7 1187.8 1482.5 990.2 平均值 18.2 1232.9 标准差 (σ) 0.4 260 表 5 三种“木棉”区别
Table 5 Differences among the three types of "Mumian" (cotton-like plants)
科属名 种名 古称或别名 生长条件 分布地区 锦葵科 棉属
(Gossypium)草棉 (herbaceum) 棉花、吉贝、木棉 年均温18~24 ℃为宜* 现主要分布于广东、云南、四川、
甘肃和新疆等省区**树棉 (arboreum) 橦树、梧桐木、娑罗木、木棉、吉贝树 年均温约17.8 ℃[46] 云南禄丰等县低海拔温暖地区尚有栽培[46] 木棉科 木棉属 (Gossampinus) 木棉 (malabarica) 攀枝花、英雄树 生长最适温度为23~
31 ℃,喜高温湿润气候,不耐寒[57]主要分布于广西、广东、贵州、
云南和四川南部等省区[57-58]注:参考郭声波[46]。*来源于https://www.cma.gov.cn/kppd/kppdqxyr/kppdnyqx/201212/t20121215_197016.html;**来源于https://baike.baidu.com/item/草棉/1914485。 表 6 学者对“西岭”的认识
Table 6 Different explanations to the “Western Mountains” by modern researchers
学者/官修 学者/著作
年代泛指/
确指山脉/山峰 出处 胡三省 宋末元初 泛指 西山在成都西 《资治通鉴音注》疏陈子昂上书云:“西山在成都西,松、茂二州都督府所统诸羌州,皆西山羌也。”[62] 仇兆鳌 明末清初 确指 岷山 《杜诗详注》引《元和郡县志》注《西山三首》云:“岷山,即汶山,南去青城山百里,天色晴朗,望见成都。山顶停雪,常深百丈,夏月融泮,江川为之洪溢,即陇之南首也。”仇注又引李宗谔《图经》云:“维州,南界,江城,岷山连岭而西,不知其极,北望高山,积雪如玉,东望成都若井底,一面孤峰,三面临江,是西蜀控吐蕃之要冲。”[60] 浦起龙 清 泛指 松、维等州诸山 《读杜心解》注云:“西山,即松、维等州诸山”[63] 杨伦 清 确指 岷山 《杜诗镜诠》引注与仇注同[18] 中国清朝官修地理总志 清 确指 九顶山 《嘉庆重修一统志》:在茂州列鹅村,去州四十里,实威茂彭灌之中。其高六十里,山有九峰,四时积雪,经暑不消,晨光射之,烂若红玉,去成都五百里。“西望之若在户牖,居人呼为九顶山,杜子美所咏西山即此也。”[61] 李绍明 1980年 泛指 唐代西山为成都平原以西岷江上游诸山之泛称 《唐代西山诸羌考略》引胡三省注[64] 高文德 1995年 泛指 西山为唐代对成都平原以西,岷江上游诸山的泛称 《中国少数民族史大辞典》[65] 张天健 1996年 确指 大雪塘 (又名庙基岭) 《杜甫“窗含西岭千秋”“西岭”考实》[66] 山东大学古典
文学教研室1998年 确指 松潘县南之雪栏山 《杜甫诗选》“西岭”即“雪岭”,在松潘县南之雪栏山,积雪终年不消,故名。[67] 郭声波 2002年 泛指 唐宋剑南道之“西山”,通常认为是对成都平原以西岷江上游的岷山、龙门山、甘松岭、羊拱山、鹧鸪山、大白山、邛崃山、夹金山等山脉的统称 《川西北羌族探源---唐宋岷江西山羁縻州部族研究》[20] 王军 2013年 确指 西岭,即西山,为岷山主峰 (雪宝顶) 《诗心:从<诗经>到<红楼梦>》西岭,即西山,因在成都西得名,为岷山主峰(雪宝顶)。因山顶终年积雪,故又称雪岭、雪山。“千秋雪”,岭上终年不化的积雪,言时间之久。[68] 赵华 2013年 泛指 西岭是泛指,千秋雪是四姑娘山幺妹峰 《由两幅照片看杜甫“窗含西岭千秋雪”的指向——成都的自然人文之旅[缩略版]》① 田峰 2016年 泛指 西山包含了松、茂、维等州的诸多山脉,即今之川西北的岷山 《杜甫从秦州到巴蜀荆湘到地理感知和文化体验》[69] 马玮 2017年 泛指 西岭指岷山,岷山在成都西 《杜甫诗歌赏析》[70] 吴石玉 2020年 确指 雪岭、西山,也称西岭、雪山,今名雪宝顶,为岷山主峰,在今四川松潘县东 嘉庆《四川通志》中的成都杜甫草堂历史资料综述[71] -
[1] 侯甬坚, 祝一志. 历史记录提取的近5~2.7ka黄河中下游平原重要气候事件及其环境意义[J]. 海洋地质与第四纪地质, 2000, 20(4):23-29 HOU Yongjian, ZHU Yizhi. Important climatic events showed by historical records from middle and lower reach plain of the Yellow River during 5~2.7 Ka and their environmental significance[J]. Marine Geology & Quaternary Geology, 2000, 20(4):23-29.]
[2] 费杰, 侯甬坚, 刘晓东, 等. 基于黄土高原南部地区历史文献记录的唐代气候冷暖波动特征研究[J]. 中国历史地理论丛, 2001, 16(4):74-81 doi: 10.3969/j.issn.1001-5205.2001.04.013 FEI Jie, HOU Yongjian, LIU Xiaodong, et al. Fluctuation characteristics of climatic change in temperature of Tang Dynasty based on historical document records in South Loess Plateau[J]. Collections of Essays on Chinese Historical Geography, 2001, 16(4):74-81.] doi: 10.3969/j.issn.1001-5205.2001.04.013
[3] 刘炳涛, 满志敏. 古代诗歌中的气候信息及其运用[J]. 中国历史地理论丛, 2010, 25(4):5-14 LIU Bingtao, MAN Zhimin. Climate information of ancient poetry and its application[J]. Journal of Chinese Historical Geography, 2010, 25(4):5-14.]
[4] 刘亚辰, 方修琦, 陶泽兴, 等. 诗歌中物候记录的基本特征及用于历史气候重建的处理方法[J]. 地理科学进展, 2017, 36(4):483-490 doi: 10.18306/dlkxjz.2017.04.009 LIU Yachen, FANG Xiuqi, TAO Zexing, et al. Basic features of phenological records in poetry and their usage for reconstructing past climate change[J]. Progress in Geography, 2017, 36(4):483-490.] doi: 10.18306/dlkxjz.2017.04.009
[5] 满志敏. 关于唐代气候冷暖问题的讨论[J]. 第四纪研究, 1998, 18(1):20-30 doi: 10.3321/j.issn:1001-7410.1998.01.003 MAN Zhimin. Climate in Tang Dynasty of China: Discussion for its evidence[J]. Quaternary Sciences, 1998, 18(1):20-30.] doi: 10.3321/j.issn:1001-7410.1998.01.003
[6] 马亚玲, 崔玉娟, 方修琦, 等. 杜诗记载的唐代荆湘地区寒冬及其古气候意义[J]. 古地理学报, 2015, 17(1):137-142 doi: 10.7605/gdlxb.2015.01.012 MA Yaling, CUI Yujuan, FANG Xiuqi, et al. Cold winters of Jing-Xiang region in Tang Dynasty recorded in Du Fu's poems and their palaeoclimatic significance[J]. Journal of Palaeogeography, 2015, 17(1):137-142.] doi: 10.7605/gdlxb.2015.01.012
[7] 蓝勇. 中国西南历史气候初步研究[J]. 中国历史地理论丛, 1993, 8(2):13-39 LAN Yong. A preliminary study of the historical climate in South-West of China[J]. Collections of Essays on Chinese Historical Geography, 1993, 8(2):13-39.]
[8] 费杰, 周杰, 安芷生. 历史文献记录的唐五代时期(618-959AD)气候冷暖变化[J]. 海洋地质与第四纪地质, 2004, 24(2):109-118 FEI Jie, ZHOU Jie, AN Zhisheng. Temperature change in China over 618-959 AD: Based on historical records[J]. Marine Geology & Quaternary Geology, 2004, 24(2):109-118.]
[9] 葛全胜, 刘浩龙, 郑景云, 等. 隋唐时期东中部地区温度变化的重建(601~920年)[J]. 科学通报, 2010, 55(31): 3048-3055 GE Quansheng, LIU Haolong, ZHENG Jingyun, et al. Reconstructing temperature change in Central East China during 601-920 AD[J]. Chinese Science Bulletin, 2010, 55(34): 3944-3949.]
[10] 郭忠明, 顾祝军, 吴红波, 等. 冰川雪线高度研究进展[J]. 遥感技术与应用, 2016, 31(4):645-652 GUO Zhongming, GU Zhujun, WU Hongbo, et al. Research progress of glacier snowline altitude[J]. Remote Sensing Technology and Application, 2016, 31(4):645-652.]
[11] 于希贤. 苍山雪与历史气候冷期变迁研究[J]. 中国历史地理论丛, 1996, 11(2):25-39 YU Xixian. Snow in Cangshan Mountain and the cold period variation of historical climate[J]. Collections of Essays on Chinese Historical Geography, 1996, 11(2):25-39.]
[12] 璩向宁, 汪一鸣. 近一千年来贺兰山积雪和气候变化[J]. 地理研究, 2006, 25(1):35-42 doi: 10.3321/j.issn:1000-0585.2006.01.005 QU Xiangning, WANG Yiming. Investigation on changes of snow cover and climate in Helanshan mountains for the past approximate millennium[J]. Geographical Research, 2006, 25(1):35-42.] doi: 10.3321/j.issn:1000-0585.2006.01.005
[13] 杨景春, 李有利. 地貌学原理(4版)[M]. 北京: 北京大学出版社, 2017: 78-79 YANG Jingchun, LI Youli. Principles of Geomorphology(4th ed.)[M]. Beijing: Peking University Press, 2017: 78-79.]
[14] 张籍. 张籍诗集[M]. 北京: 中华书局, 1959: 80, 67 ZHANG Ji. The Collected Poems of Zhang Ji[M]. Beijing: Zhonghua Book Company, 1959: 80, 67.]
[15] 吴淑娴. 中国果树志 荔枝卷[M]. 北京: 中国林业出版社, 1998: 20-39 WU Shuxian. Chinese Fruit Trees: Litchi Volume[M]. Beijing: China Forestry Publishing House, 1998: 20-39.]
[16] 白慧卿, 吴建国, 潘学标. 影响我国荔枝分布的关键气候要素分析[J]. 果树学报, 2016, 33(4):436-443 BAI Huiqing, WU Jianguo, PAN Xuebiao. Key climatic factors affecting the distribution of litchi in China[J]. Journal of Fruit Science, 2016, 33(4):436-443.]
[17] 吴怀东. 杜甫的美学——《绝句四首》之三“两个黄鹂鸣翠柳”解读[J]. 杜甫研究学刊, 2018, 38(3):37-44 WU Huaidong. On Du Fu's Aesthetics——Interpretation of "Two Golden Orioles Sing amid the Willows Green 'in' Four Quatrains"[J]. Journal of Dufu Studies, 2018, 38(3):37-44.]
[18] 杜甫. 杜诗镜铨[M]. 上海: 上海古籍出版社, 1998: 559-560, 317, 334, 575, 472 DU Fu. Du Fu's Poems[M]. Shanghai: Shanghai Classics Publishing House, 1998: 559-560, 317, 334, 575, 472.]
[19] 严钦尚, 曾昭璇. 地貌学[M]. 北京: 高等教育出版社, 1985: 124-125 YAN Qingshang, ZENG Zhaoxuan. Geomorphology[M]. Beijing: Higher Education Press, 1985: 124-125.]
[20] 郭声波. 川西北羌族探源——唐宋岷江西山羁縻州部族研究[J]. 中南民族大学学报: 人文社会科学版, 2002, 22(4):74-79 GUO Shengbo. The origin of the Ch'iang People in the Northwest Sichuan——A research into Xishan Jimi prefectures[J]. Journal of South-Central University for Nationalities: Humanities and Social Science, 2002, 22(4):74-79.]
[21] 刘亚辰, 陶泽兴, 呼雪梅. 唐宋诗歌中物候记录的规律总结及其对环境的指示意义[J]. 中国历史地理论丛, 2022, 37(4):12-21 LIU Yachen, TAO Zexing, HU Xuemei. Phenological records’ regularity in poems of the Tang and Song Dynasties and its indications of the environment[J]. Journal of Chinese Historical Geography, 2022, 37(4):12-21.]
[22] 卡列斯尼克 C B. 普通冰川学[M]. 兰州: 中国科学院地理研究所冰川冻土研究室, 1965: 32-33 Kalesnik C B. General Glaciology[M]. Lanzhou: Lanzhou Institute of Glaciology and Cryopedology, Chinese Academy of Sciences, 1965: 32-33.]
[23] 姚永慧, 张百平. 青藏高原气温空间分布规律及其生态意义[J]. 地理研究, 2015, 34(11):2084-2094 YAO Yonghui, ZHANG Baiping. The spatial pattern of monthly air temperature of the Tibetan Plateau and its implications for the geo-ecology pattern of the Plateau[J]. Geographical Research, 2015, 34(11):2084-2094.]
[24] 曾寒梅. 成渝两地城市形态特征形成与演变研究[D]. 重庆大学硕士学位论文, 2015 ZENG Hanmei. A study of the formation and evolution of urban morphology features in Chengdu and Chongqing[D]. Master Dissertation of Chongqing University, 2015.]
[25] 周斌, 邹强, 蒋虎, 等. 川西高原气候变化特征及泥石流动态危险性响应研究[J]. 自然灾害学报, 2022, 31(4):241-255 ZHOU Bin, ZOU Qiang, JIANG Hu, et al. Research on climate change characteristics and change of debris flow hazard in the Chuanxi plateau[J]. Journal of Natural Disasters, 2022, 31(4):241-255.]
[26] 星球研究所, 中国青藏高原研究会. 这里是中国[M]. 北京: 中信出版集团股份有限公司, 2019: 90 Institute for Planets, The China Society on Tibet Plateau. HI I'M China[M]. Beijing: CITIC Press Corporation, 2019: 90.]
[27] 雷传扬, 王波, 刘兆鑫, 等. 成都平原河流阶地的发育及其对古气候和新构造运动的指示[J]. 沉积与特提斯地质, 2024, 44(1):20-33 LEI Chuanyang, WANG Bo, LIU Zhaoxin, et al. Development of fluvial terraces in Chengdu Plain: Implications for the paleoclimate and neotectonic movement[J]. Sedimentary Geology and Tethyan Geology, 2024, 44(1):20-33.]
[28] 王羽珂, 陈浩, 冯兴雷, 等. 成都平原东郊台地中更新统合江组沉积特征及工程地质意义[J]. 沉积与特提斯地质, 2019, 39(3):33-39 doi: 10.3969/j.issn.1009-3850.2019.03.004 WANG Yuke, CHEN Hao, FENG Xinglei, et al. Sedimentary characteristics and engineering geological significance of the middle Pleistocene Hejiang Formation in the eastern suburb platform on the Chengdu Plain, Sichuan[J]. Sedimentary Geology and Tethyan Geology, 2019, 39(3):33-39.] doi: 10.3969/j.issn.1009-3850.2019.03.004
[29] 马志刚. 成都平原卵砾石层地震响应研究[D]. 成都理工大学硕士学位论文, 2009 MA Zhigang. The study on the seismic response of gravel formation on Chengdu Plain[D]. Master Dissertation of Chengdu University of Technology, 2009.]
[30] 蓝勇. 四川荔枝种植公布的历史考证[J]. 西南师范大学学报: 自然科学版, 1985, 10(4):86-99 LAN Yong. An investigation of the history and distribution of litchi planting in Sichuan[J]. Journal of Southwest China Normal University: Natural Science Edition, 1985, 10(4):86-99.]
[31] 四川省统计局, 国家统计局四川调查总队. 四川统计年鉴2017[M]. 北京: 中国统计出版社, 2017: 187-188 Sichuan Provincial Bureau of Statistics, NBS Survey Office in Sichuan. Sichuan Statistical Yearbook 2017[M]. Beijing: China Statistics Press, 2017: 187-188.]
[32] 四川省统计局, 国家统计局四川调查总队. 四川统计年鉴2018[M]. 北京: 中国统计出版社, 2018: 187-188 Sichuan Provincial Bureau of Statistics, NBS Survey Office in Sichuan. Sichuan Statistical Yearbook 2018[M]. Beijing: China Statistics Press, 2018: 187-188.]
[33] 四川省统计局, 国家统计局四川调查总队. 四川统计年鉴2019[M]. 北京: 中国统计出版社, 2019: 179-180 Sichuan Provincial Bureau of Statistics, NBS Survey Office in Sichuan. Sichuan Statistical Yearbook 2019[M]. Beijing: China Statistics Press, 2019: 179-180.]
[34] 四川省统计局, 国家统计局四川调查总队. 四川统计年鉴2020[M]. 北京: 中国统计出版社, 2020: 143-144 Sichuan Provincial Bureau of Statistics, NBS Survey Office in Sichuan. Sichuan Statistical Yearbook 2020[M]. Beijing: China Statistics Press, 2020: 143-144.]
[35] 四川省统计局, 国家统计局四川调查总队. 四川统计年鉴2021[M]. 北京: 中国统计出版社, 2021: 137-138 Sichuan Provincial Bureau of Statistics, NBS Survey Office in Sichuan. Sichuan Statistical Yearbook 2021[M]. Beijing: China Statistics Press, 2021: 137-138.]
[36] 徐新良. 中国多年度地市行政区划边界数据[DB/OL]. 资源环境科学数据注册与出版系统, 2023. 2024-03-21]. http: //www. resdc. cn/. [XU Xinliang. Multi-year administrative boundary data for municipalities in China[DB/OL]. Resource and Environmental Science Data Registration and Publication System, 2023. [2024-03-21]. http://www.resdc.cn/.
[37] 徐新良. 中国多年度省级行政区划边界数据[DB/OL]. 资源环境科学数据注册与出版系统, 2023. [2024-03-21]. http://www.resdc.cn/. XU Xinliang. Multi-year administrative boundary data for provinces in China[DB/OL]. Resource and Environmental Science Data Registration and Publication System, 2023. [2024-03-21]. http://www.resdc.cn/.
[38] 张立汉. 中国山河全书(上)[M]. 青岛: 青岛出版社, 2005: 1249 ZHANG Lihan. Mountains and Rivers of China( Volume One) [M]. Qingdao: Qingdao Publishing House, 2005: 1249.]
[39] 陈梦蝶. 青藏高原地区雪线时空变化动态研究[D]. 兰州大学硕士学位论文, 2014 CHEN Mengdie. Spatio-temporal change dynamic of snowline on Tibetan Plateau[D]. Master Dissertation of Lanzhou University, 2014.]
[40] 吴锡浩, 赵文龙. 固态降水率及其与气候雪线的关系[C]//第四纪冰川与第四纪地质文集 第五集. 北京: 地质出版社, 1988: 174-182 WU Xihao, ZHAO Wenlong. Solid precipitation rate and its relation to climatic snow line[C]//Quaternary Glaciation and Quaternary Geology 5th Collection. Beijing: Geology Press, 1988: 174-182.]
[41] 蒋复初, 吴锡浩, 王书兵, 等. 中国气候雪线空间分布特征[J]. 地质力学学报, 2002, 8(4):289-296 doi: 10.3969/j.issn.1006-6616.2002.04.001 JIANG Fuchu, WU Xihao, WANG Shubing, et al. Characteristics of space distribution of the climatic snowline in China[J]. Journal of Geomechanics, 2002, 8(4):289-296.] doi: 10.3969/j.issn.1006-6616.2002.04.001
[42] 周维扬, 丁浩. 杜甫草堂史话[M]. 成都: 天地出版社, 2009: 42-45 ZHOU Weiyang, DING Hao. The history of Du Fu thatched[M]. Chengdu: Tiandi Press, 2009: 42-45.]
[43] 竺可桢. 中国近五千年来气候变迁的初步研究[J]. 考古学报, 1972, 17(1):15-38 ZHU Kezhen. A preliminary study on the climatic fluctuations during the last 5000 years in China[J]. Acta Archaeologica Sinica, 1972, 17(1):15-38.]
[44] 蓝勇. 近2000年来长江上游荔枝分布北界的推移与气温波动[J]. 第四纪研究, 1998, 14(1):39-45 doi: 10.3321/j.issn:1001-7410.1998.01.005 LAN Yong. The movements of the northern boundary of Litchi distribution and fluctuations of temperature in the upper reaches of the Yangtze River in the past 2000 years[J]. Quaternary Sciences, 1998, 14(1):39-45.] doi: 10.3321/j.issn:1001-7410.1998.01.005
[45] 郭声波. 成都荔枝与十二世纪的寒冷气候[J]. 中国历史地理论丛, 1989, 4(3):38 GUO Shengbo. Chengdu Lychee and the cold climate of the twelfth century[J]. Collections of Essays on Chinese Historical Geography, 1989, 4(3):38.]
[46] 郭声波. 历史时期四川手工业原料作物的分布[J]. 中国历史地理论丛, 1990, 5(1):67-88 GUO Shengbo. Distribution of raw material crops for handicrafts in Sichuan during the historical period[J]. Collections of Essays on Chinese Historical Geography, 1990, 5(1):67-88.]
[47] 葛全胜. 中国历朝气候变化[M]. 北京: 科学出版社, 2011: 308 GE Quansheng. Climate Change in Chinese Dynasties[M]. Beijing: Science Press, 2011: 308.]
[48] 娄雨亭. 薛涛与唐代成都的荔枝及气候冷暖问题[J]. 中国史研究, 2001, 23(3):38 LOU Yuting. Xue Tao and Lychee in Chengdu in the Tang Dynasty and the Problem of Warm and Cold Climate[J]. Journal of Chinese Historical Studies, 2001, 23(3):38.]
[49] 聂顺新. 再论唐代长江上游地区的荔枝分布北界及其与气温波动的关系[J]. 中国历史地理论丛, 2011, 26(1):139-144,158 NIE Shunxin. Northern boundary of litchi's distribution in the upper Yangtze of Tang Dynasty and it's relationship with the temperature fluctuations[J]. Journal of Chinese Historical Geography, 2011, 26(1):139-144,158.]
[50] 满志敏. 唐代气候冷暖分期及各期气候冷暖特征的研究[J]. 历史地理, 1990, 7(2):1-15 MAN Zhimin. The division of cold and warm climatic periods of the Tang Dynasty Time and the clmracteristies of the respective term[J]. Historical Geography, 1990, 7(2):1-15.]
[51] 蓝勇. 采用物候学研究历史气候方法问题的讨论——答《再论唐代长江上游地区的荔枝分布北界及其与气温波动的关系》一文[J]. 中国历史地理论丛, 2011, 26(2):14-20 LAN Yong. A discussion on research's method of historical climate with the Phenology: answers on the northern boundary of Litchi's distribution in the upper Yangtze River in Tang Dynasty and the relationship with the temperature fluctuations[J]. Journal of Chinese Historical Geography, 2011, 26(2):14-20.]
[52] 唐圭璋. 全唐诗[M]. 北京: 中华书局, 2018 TANG Guizhang. Tang Poems[M]. Beijing: Zhonghua Book Company, 2018.]
[53] 张晓婷, 庄赟, 董嘉辉, 等. 荔枝种质资源抗寒性综合评价[J]. 果树学报, 2024, 41(3):403-425 ZHANG Xiaoting, ZHUANG Yun, DONG Jiahui, et al. Comprehensive evaluation of cold tolerance in Litchi Germplasm resources[J]. Journal of Fruit Science, 2024, 41(3):403-425.]
[54] 邓乐群. 杜甫诗歌所叙唐代陇蜀荆湘气候特征[J]. 南通大学学报: 社会科学版, 2009, 25(6):58-64 DENG Lequn. Climatic features of Long, Shu, Jin and Xiang in the Tang Dynasty revealed in DU Fu’s Poems[J]. Journal of Nantong University: Social Sciences Edition, 2009, 25(6):58-64.]
[55] 史学通. 我国历史上的木棉问题[J]. 中国史研究, 1981, 3(2):85-91 SHI Xutong. The problem of "Mumian" in China's history[J]. Journal of Chinese Historical Studies, 1981, 3(2):85-91.]
[56] 辛如如. 木棉栽培技术综述[J]. 防护林科技, 2014, 32(11):120-121 XIN Ruru. An overview of Bombax malabaricum cultivation technology[J]. Protection Forest Science and Technology, 2014, 32(11):120-121.]
[57] 高平, 谌振, 林忠, 等. 木棉名称考据与应用研究[J]. 南方农业, 2016, 10(22):24-27 GAO Ping, CHEN Zhen, LIN Zhong, et al. Research on the name evidence and application of Bombax malabaricum[J]. South China Agriculture, 2016, 10(22):24-27.]
[58] 韩茂莉. 中国历史农业地理[M]. 北京: 北京大学出版社, 2012: 626-627 HAN Maoli. Historical Agricultural Geography of China[M]. Beijing: Peking University Press, 2012: 626-627.]
[59] 吴锡浩. 青藏高原东南部现代雪线和林线及其关系的初步研究[J]. 冰川冻土, 1989, 11(2):113-124 WU Xihao. A preliminary study on existing snowline tember line and their relations in southeastern part of the Qinghai-Xizang Plateau[J]. Journal of Glaciology and Geocryology, 1989, 11(2):113-124.]
[60] 杜甫. 杜诗详注[M]. 上海: 上海古籍出版社, 1992: 412 DU Fu. The Detailed Annotations on Du Fu's Poems[M]. Shanghai: Shanghai Classics Publishing House, 1992: 412.]
[61] 清)仁宗敕撰. 四部丛刊续编 史部 嘉庆重修一统志 23 卷三百八十四[M]. 上海: 上海书店出版社, 1984 Qing dynasty) by Emperor Renzong. Si Bu Cong Kan Sequels, sequel to the history part, Yi Tong Zhi revision by Jiaqing 23 vol. 384[M]. Shanghai: Shanghai Bookstore Publishing House, 1984.]
[62] 胡三省. 资治通鉴音注 卷二百四[M]. 胡克家, 清嘉庆21 HU Sanxing. Zizhi Tongjian Yinzhu, Volume 244[M]. Hu Kejia, Qing Jiaqing 21.]
[63] 浦起龙. 读杜心解[M]. 北京: 中华书局, 1981: 462 PU Qilong. Explanation of Dufu's Poems[M]. Beijing: Zhonghua Book Company, 1981: 462.]
[64] 李绍明. 唐代西山诸羌考略[J]. 四川大学学报: 哲学社会科学版, 1980, 13(1):83-95 LI Shaoming. A study of the Qiang in the Western Mountains during the Tang Dynasty[J]. Journal of Sichuan University: Philosophy and Social Science Edition, 1980, 13(1):83-95.]
[65] 高文德. 中国少数民族史大辞典[M]. 长春: 吉林教育出版社, 1995: 700 GAO Wende. Dictionary of Chinese Minority History[M]. Changchun: Jilin Education Publishing House, 1995: 700.]
[66] 张天健. 杜甫“窗含西岭千秋雪”“西岭”考实[J]. 中国地名, 1996, 1(4):23 ZHANG Tianjian. The textual research of the "Western Mountain" in Du Fu's "The Window Contains Thousands of Autumn Snow in the Western Mountain"[J]. China Place Name, 1996, 1(4):23.]
[67] 山东大学中文系古典文学教研室选注, 袁世硕修订. 杜甫诗选[M]. 北京: 人民文学出版社, 1998: 258 Selected and annotated by the Department of Classical Literature, Shandong University, revised by YUAN Shishuo. Selected Poems of Du Fu[M]. Beijing: People's Literature Publishing House, 1998: 258.]
[68] 王军. 诗心: 从《诗经》到《红楼梦》[M]. 北京: 东方出版社, 2013: 76 WANG Jun. Poetic Heart: From the Book of Songs to The Story of the Stone[M]. Beijing: Oriental Publishing House, 2013: 76.]
[69] 田峰. 杜甫从秦州到巴蜀荆湘的地理感知与文化体验[J]. 中国韵文学刊, 2016, 30(1):11-17 doi: 10.3969/j.issn.1006-2491.2016.01.003 TIAN Feng. Du Fu's Geographical Perception and Cultural Experience from Qinzhou to Bashu Jingxiang[J]. Journal of Chinese Verse Studies, 2016, 30(1):11-17.] doi: 10.3969/j.issn.1006-2491.2016.01.003
[70] 马玮. 杜甫诗歌赏析[M]. 北京: 商务印书馆国际有限公司, 2017: 214 MA Wei. An Appreciation of Du Fu's Poetry[M]. Beijing: Business Book Printing Place International Co. , Ltd. , 2017: 214.]
[71] 吴石玉. 嘉庆《四川通志》中的成都杜甫草堂历史资料综述[J]. 杜甫研究学刊, 2020, 40(4):63-81 WU Shiyu. A textual research on historical documents of Du Fu Thatched Cottage in Chengdu in general history of Sichuan in Jiaqing Period of the Qing Dynasty[J]. Journal of Dufu Studies, 2020, 40(4):63-81.]
[72] 汤国安, 杨昕, 张海平. ArcGIS地理信息系统空间分析实验教程[M]. 3版. 北京: 科学出版社, 2021: 381-382, 385 TANG Guoan, YANG Xin, ZHANG Haiping. ArcGIS Geographic Information System Spatial Analysis Experimental Tutorial[M]. 3rd ed. Beijing: Science Press, 2021: 381-382, 385.]
[73] 肖时珍, 肖华, 吴宇辉. 基于GIS视域分析的项目建设对世界遗产景观美学价值的影响评价——以武陵源世界自然遗产地为例[J]. 桂林理工大学学报, 2020, 40(3):516-522 doi: 10.3969/j.issn.1674-9057.2020.03.008 XIAO Shizhen, XIAO Hua, WU Yuhui. Assessment of construction project on the aesthetic values of world heritage landscape based on GIS viewshed analysis: a case study of Wulingyuan World Natural Heritage Site[J]. Journal of Guilin University of Technology, 2020, 40(3):516-522.] doi: 10.3969/j.issn.1674-9057.2020.03.008
[74] 四川省阿坝藏族羌族自治州茂汶羌族自治县地方志编纂委员会. 茂汶羌族自治县志[M]. 成都: 四川辞书出版社, 1997: 80 Committee of Chorography Mao-wen Qiang Nationality Autonomous County, Aba Zang and Qiang Nationality Autonomous District, Sichuan Province. Chorography of Mao-wen Qiang Nationality Autonomous County[M]. Chengdu: Dictionary Press of Sichuan, 1997: 80.]
[75] 张宗福, 张晓英. 略论杜甫《西山三首》[J]. 杜甫研究学刊, 2012, 32(1):9-15 ZHANG Zongfu, ZHANG Xiaoying. A brief discussion of Du Fu's Three Songs from the Western Mountains[J]. Journal of Dufu Studies, 2012, 32(1):9-15.]
-
期刊类型引用(9)
1. 吴能友,李彦龙,蒋宇静,孙金声,谢文卫,胡高伟,王韧,于彦江,王金堂,陈强,申凯翔,孙志文,陈明涛. 海洋天然气水合物工程地质学的提出、学科内涵与展望. 地学前缘. 2025(02): 216-229 . 百度学术
2. 蒋宇静,颜鹏,栾恒杰,刘明康,梁伟,杜晓宇,马先壮,史一辰. 不同降压速率下天然气水合物垂直井分段降压开采特性试验研究. 岩土力学. 2024(09): 2682-2694 . 百度学术
3. 纪云开,王佳贤,邱晓倩,孟庆国,张永超,胡高伟,蔡峰,刘昌岭. 盐度影响下松散沉积物中甲烷水合物生成特征研究. 海洋地质与第四纪地质. 2024(06): 71-81 . 本站查看
4. 崔玉东,陆程,关子越,罗万静,滕柏路,孟凡璞,彭越. 南海海域天然气水合物降压开采储层蠕变对气井产能影响. 油气藏评价与开发. 2023(06): 809-818 . 百度学术
5. 李辉,张旭辉,陆程,谢鹏飞,鲁晓兵. 压缩加载条件下含水合物沉积物蠕变特性分析. 海洋地质与第四纪地质. 2023(06): 217-225 . 本站查看
6. Xu-wen Qin,Cheng Lu,Ping-kang Wang,Qian-yong Liang. Hydrate phase transition and seepage mechanism during natural gas hydrates production tests in the South China Sea: A review and prospect. China Geology. 2022(02): 201-217 . 必应学术
7. 董刚,蔡峰,孙治雷,闫桂京,梁杰,李清,孙运宝,李昂,骆迪,翟滨,郭建卫,窦振亚. 海洋浅表层天然气水合物地质取样技术及样品现场处置方法. 海洋地质前沿. 2022(07): 1-9 . 百度学术
8. 秦绪文,陆程,王平康,梁前勇. 中国南海天然气水合物开采储层水合物相变与渗流机理:综述与展望. 中国地质. 2022(03): 749-769 . 百度学术
9. 刘乐乐,李彦龙,张旭辉,罗大双,刘昌岭. 降压开采导致天然气水合物系统状态演化模拟实验. 工程地质学报. 2021(06): 1916-1925 . 百度学术
其他类型引用(5)