首页 | 官方网站   微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 343 毫秒
1.
研究泉水氢氧稳定同位素时空变化特征,探究泉水的补给来源,对水资源管理及旅游业可持续发展具有重要指导意义。以百脉泉泉群为研究区,通过分析降水及泉水中的氢氧稳定同位素特征,获取了本地大气降水线、泉水同位素时空变化特征及蒸发线,并探究了百脉泉泉水的补给来源。结果表明:大气降水中δ2H的波动范围为-151.55‰δ18O—-4.28‰,δ18O的波动范围为-19.93‰δ18O—-0.52‰,δ2H和δ18O在时间上的分布规律相似,均呈现双峰状(“M”型),本地大气降水线(LMWL)为:δ2H=7.58×δ18O+4.17。百脉泉泉水的δ2H的波动范围为-65.23‰δ18O—-37.01‰,δ18O的波动范围为-8.99‰δ18O—-4.04‰。泉水δ2H—δ18O均落在本地大气降水线上或...  相似文献   

2.
大气降水同位素的研究对于深入理解地区水循环特征具有重要意义.而基于降水样品实测数据在高纬度寒冷地区的研究尚不多见.本文基于2020年7月—2021年6月一个水文年周期内68个大气降水样品,测定了降水中的δD和δ18O,分析了哈尔滨降水稳定同位素特征及水汽来源.结果表明:(1)哈尔滨大气降水中δD和δ18O平均值分别为-89.33‰和-12.46‰,具有明显的冬低夏高的特点,年内变化范围均较大,分别为-275.50‰~-3.67‰和-35.68‰~-1.10‰,较全球平均值低67.33‰和8.46‰,较中国平均值低39.33‰和4.46‰.(2)哈尔滨大气降水线方程为:δD=7.66δ18O+6.29(R2=0.98),斜率和截距与全球的8和10和中国大气降水线的7.9和8.2接近,与东北地区大气降水线的7.20和-2.39相比,斜率差异较小而截距差异明显.(3)氢氧同位素与气温之间的拟合方程为:δD=0.34T-15.40(R2=0.41,p<0.01),δ18O=2.5...  相似文献   

3.
综合分析西北内陆区97个研究站点的降水稳定同位素数据,并结合相关气象资料,揭示了西北内陆区降水稳定同位素δ18O、δD和d-excess时空分布特征,明确了海拔、经纬度、温度和降雨量对降水δ18O的影响;利用水汽通量和HYSPLIT模型追踪了大气降水的水汽来源.同时,根据关键自然地理要素的空间差异将西北内陆区划分为4个子区域,对上述内容分区域进行系统分析和对比.结果表明:(1)西北内陆区降水δ18O和δD的变化范围分别为-21.20‰~1.70‰和-144.20‰~5.21‰,d-excess波动范围为-20.37‰~46.48‰,δ18O、δD和d-excess均存在显著的空间变化和季节变化特征,河西内陆区和塔里木地区δ18O和δD值相对偏正,柴达木-青海湖区和准噶尔-吐哈地区δ18O和δD偏负.(2)西北内陆区大气降水线方程的斜率和截距均小于中国和全球大气降水线,4个子区域亦低于全球水线,其中塔里木地区斜率最低.(3)西北内陆区海拔效应为-0.04‰·...  相似文献   

4.
刘佳驹  郭怀成 《中国环境科学》2021,41(11):5194-5200
为揭示流域内水汽循环过程变化规律、揭示河水稳定同位素与气象要素及海拔之间的关系,2018年先后采集雅鲁藏布江干流丰、平、枯,水期水样90个,对河水的氢氧同位素组成进行了分析.结果表明:利用δ18O与纬度和海拔定量关系模型研究了雅江流域降水氢氧同位素空间变化特征,发现降水中δ18O值变化呈现东西部(雅江上、下游)高中部(雅江中游)低的趋势;建立了雅江丰、平、枯水期河水氢氧同位素方程,研究发现δD与δ18O的线性方程斜率低于中国大气降水线,河水δD与δ18O有一定的季节性特征,其受大气降水的影响较大,且丰水期较平水期以及枯水期相比更为明显;雅江河水中δ18O空间特征表现为,从上游至下游整体表现出先降低再升高的变化趋势.南北空间上表现为河水δD-δ18O从南至北逐步增高的趋势.δD-δ18O值与海拔高度具有较弱的负相关关系;受地表蒸发分馏、降水云团迁移以及流域岩石风化的影响,雅江河水d值季节性、空间性变化明显.  相似文献   

5.
干旱区强烈的蒸发作用导致地表水体和土壤水、浅层地下水中盐度不断升高,而盐分亦会导致蒸发强度降低,水体蒸发导致的氢氧稳定同位素分馏特征是确定蒸发量的有效手段,但目前水体蒸发过程中氢氧同位素分馏特征的研究主要集中于淡水和低矿化度盐水,为探求不同盐度对水体蒸发过程中的氢氧稳定同位素分馏的影响,开展浓度分别为0、10、20、50、100和200 g/L盐溶液的室内蒸发实验,比较纯水与不同浓度盐水蒸发过程中氢氧稳定同位素的差异,探究其分馏过程,定量评价盐度对水分蒸发分馏过程的影响。实验结果表明:蒸发量随着水体盐度的升高而呈指数形式降低;蒸发过程中,重同位素在水体中不断富集,0~200 g/L盐水中2H、18O和17O的瑞利分馏系数分别在1.089~1.155、1.019~1.033和1.009~1.018之间,并呈现随着盐水浓度的增大瑞利分馏系数逐渐增加的趋势;纯水和盐水中的δ2H-δ18O、δ2H-δ17O和δ17O-δ...  相似文献   

6.
不同水汽来源对湖南长沙地区降水中 δD、δ18O的影响   总被引:4,自引:1,他引:3  
根据2010年在长沙地区进行降水收集和气象要素观测的资料,分析了该地区降水中δ18O与温度、降水量之间的关系,揭示了降水中δD、δ18O的变化特征,讨论了水汽输送对降水中δ18O变化的影响。结果表明,在天气尺度下,长沙地区大气降水中δ18O与降水量、温度之间存在显著的负相关关系,即该地区降水中δ18O的变化具有显著的降水量效应及反温度效应。对长沙地区的降雪样和降雨样进行线性回归,得出大降水事件和降雪的大气降水线具有较大斜率和截距。随着降水量的减小,大气降水线方程的斜率和截距也逐渐减小,这主要由于小降水事件的雨滴在降落过程中受到二次蒸发强烈,同位素分馏强烈。利用HYSPLIT模式追踪该地区气流的轨迹发现,在季风降水期间(5—9月),δ18O值偏低的水汽主要来自孟加拉湾、南海洋面与西太平洋海区;在非季风降水期间(10—4月),δ18O值偏高的水汽主要来自西风带携带的水汽和局地水汽环流。  相似文献   

7.
水汽源区变化与ENSO事件显著影响季风区水循环过程。基于珠江中下游地区4个GNIP站点(中国香港、广州、桂林、柳州)的降水同位素及OLR(向外长波辐射)数据,研究了ENSO背景下δ18O的时空分布特征及ENSO事件对降水中氢氧同位素特征的影响机制。结果表明:ENSO事件是影响稳定同位素年际差异的主要因素,通过影响雨季降水的年内分配而实现;正常年大气降水线方程的斜率与截距均大于厄尔尼诺年而小于拉尼娜年;拉尼娜年加强 δ18O的反温度效应,厄尔尼诺年减弱反温度效应;厄尔尼诺与拉尼娜年热带辐合带的变化规律呈现出相反的趋势,其OLR场的变化与研究区稳定同位素特征有着较强的对应关系;ENSO事件年不同的水汽源区相对湿度特征是造成d值年际差异的主要因素。  相似文献   

8.
选取鄱阳湖典型洪泛湿地为研究对象,分析了2018年4~10月降水、湖水、河水和湿地地下水的氢氧同位素变化特征,利用δ18O~δD关系确定了不同水文时期湿地各类水体的转化关系,并结合同位素端元混合模型估算了不同水源对湿地地下水的贡献分量.结果表明,研究区降雨δ18O和δD值在6~7月份偏小,其余月份较高,存在明显季节变化和雨量效应.河水、湖水同位素与降水同位素的季节变化规律基本一致,但受蒸发分馏影响,重同位素更为富集,且变化幅度远小于降水同位素.湿地地下水同位素的季节变化较小,δ18O、δD均值(-5.26‰,-31.1‰)高于大气降水(-6.32‰,-40.1‰)、低于湖水(-3.60‰,-26.4‰),与河水同位素(-5.09‰,-34.4‰)较为接近,表明湿地地下水受降水、湖水和河水的共同影响.涨水期(4~5月)河水的补给源为降雨和流域内地下径流,湖水主要受河水和降水共同补给,湿地地下水主要受前期降水和河水补给的滞后影响,河水的贡献比重更大.丰水期(6~8月)地下水主要接受湖水和河水共同补给,湖水的补给贡献比例超过50%,退水期(9~10月)湿地地下水向河道和湖泊等地表水体排泄.  相似文献   

9.
北京西山大气降水中D和18O组成变化及水汽来源   总被引:2,自引:1,他引:1  
论文利用北京西山鹫峰地区2011年7月—2012年7月观测站点内大气降水δD和δ18O数据,研究了北京西山山前丘陵区降水稳定同位素的变化特征。结果表明:研究区的大气降雨δD和δ18O的平均值要低于中国和全球降水同位素的平均值,并与北京地区的降水同位素组成平均值之间显示出较大的差异,且同一个月份内降水的同位素特征也具有较大的变化范围。由于季风气候以及局地水汽循环的共同影响,使得夏季降水中的δD和δ18O值较小,并且在4—7月随着月降水量的增加和温度的升高,在温度效应和降雨量效应的交互作用中,降雨量效应此时起主导作用,月均δD和δ18O值呈现减小的趋势,而在7—9月,温度效应呈现出主导作用,δD和 δ18O值呈现增加的趋势,而冬季来自偏北风的大陆性气团形成的降水δD和δ18O值较大。同时局地水汽循环过程中降水的二次蒸发,也导致在典型长历时场降水过程中δD和δ18O值呈现出分段波动式下降的特征。研究区大气降水线方程为 δD = 7.17 δ18O + 1.46,斜率和截距小于全球和中国大气降水线方程,而与北京地区大气降水线方程的对比发现,斜率差异较小而截距差异明显,说明海洋水汽不是研究区唯一的水汽来源,同时受到大陆性气团以及局地较强的水汽循环和北京地区近年的干旱化过程的影响。而氘盈余的年内变化的无规律则间接反映了研究区水汽来源的复杂性。  相似文献   

10.
厦门地区大气降水氢氧同位素组成特征及水汽来源探讨   总被引:7,自引:0,他引:7  
采集厦门地区6个站位春、夏和冬季的大气降水样品,并用稳定同位素质谱仪分析降水样品中的氢氧同位素值(δD和δ18O).结果表明:厦门地区大气降水中δD和δ18O值春季最高(-7.86‰±8.07‰和-2.18‰±0.80‰),夏季最低(-61.17‰±4.85‰和-8.42‰±0.62‰).本文同时利用HYSPLIT模型对不同季节厦门地区水汽来源及输送路径进行追踪,发现厦门地区夏季降水主要受到来自南海及西太平洋气团的影响,期间降水量大,δD和δ18O值较低.厦门地区大气降水线方程为δD=8.35δ18O+12.52(R2=0.906),与全球降水线方程(δD=8.17δ18O+10.56)相比,截距及斜率略有偏高.厦门地区氘剩余值(d值)波动范围较大(-5.13‰~32.25‰),说明厦门地区降水的水汽来源较为多样,降雨条件较为复杂.厦门地区降水中d值表现为冬季高,春季次之,夏季低的季节性变化特征.年尺度下,厦门地区氢氧同位素与降水量在呈显著的负相关关系(r分别为-0.477和-0.369,p0.01).  相似文献   

11.
通过采集闪电河流域2020年2月至2022年2月的降水与2021年的丰(8月)和枯(10月)水期的地表水和地下水,运用稳定同位素技术,对流域“三水”的氢氧稳定同位素进行时空变化分析,探讨水体同位素与环境因子的关系,结合HYSPLIT模型追踪大气水汽来源,利用端元混合模型揭示水体转换关系.结果发现,当地降水线的斜率和截距均小于当地大气降水线,水汽主要来自西风水汽、极地气团和东亚季风环流,降水同位素有显著的温度效应;时间上,地表水与地下水同位素在季节变化上均表现出枯水期较丰水期更富集,地表水与地下水d-excess值均低于全球平均值,显现出当地强烈的蒸发作用;空间上,地表水δ18O值丰枯季节在空间上具有相同的变化特征,呈现上游至下游逐渐富集,地下水δ18O高值区空间分布不均,地下水δ18O值随埋深的增加更加贫化;地下水水线斜率最高在丰水期为7.87,与当地大气降水线和地表水水线斜率十分接近,表明丰水期“三水”存在复杂的水力联系.研究区在丰水期地表水主要接受降水的补给,其次是接受地下水的径流补给.研究结果有助于明确闪电河流域同...  相似文献   

12.
我国南部夏季季风降水水汽来源的稳定同位素证据   总被引:10,自引:1,他引:9  
我国南部地区夏季降水多受季风影响。不同的季风将来自不同通道的水汽带入我国境内,控制降水的时空分布。论文利用CHNIP(中国大气降水同位素网络)中位于南部地区的观测站点,在2005年7月间,收集了月大气降水样及同步观测的气象数据。分析表明,降水中稳定氢氧同位素的空间分布可以很好地示踪和反演该地区夏季季风降水的3个主要水汽来源以及传输路径--体现南亚季风的西南水汽通道、体现南海季风的南海水汽通道及体现副热带季风的东南水汽通道。得到的大气降水线方程:δD=5.15δ18O-15.5反映了我国南部地区的降水过程历经了一定的蒸发。对δ18O与各环境因子的关系进行探讨时发现,δ18O与降雨量和高程存在对数关系,而与温度和相对湿度间存在显著的二次函数关系。综合考虑各环境因子对δ18O的影响,给出多元线性回归方程:δ18O(‰)=0.007H(m)+1.47T(℃)-0.02P(mm)+0.24RH(%)-66.3。  相似文献   

13.
为探明内蒙古地区夏季大气降水中δD和δ18O组成特征及其对气象因子变化的响应关系,于2017~2019年夏季采集了内蒙古阿拉善左旗、呼和浩特市市区、正蓝旗、克什克腾旗(达里湖)、通辽科尔沁左翼中旗和呼伦贝尔新巴尔虎右旗(呼伦湖)等6个区域共计82次大气降水样品,结合来自全球降水同位素观测网(GNIP)的包头、张掖等6个区域大气降水样品中δD和δ18O数据,分析了内蒙古地区大气降水中δD和δ18O变化的区域差异及其主要影响因素,结果表明:内蒙古局地大气降水中δD和δ18O值存在自西向东不断偏负的趋势,其中呼伦贝尔新巴尔虎右旗大气降水中δD和δ18O值最偏负;相对地,西部阿拉善左旗大气降水中δD和δ18O值最偏正;内蒙古地区局地大气降水线斜率和截距同样表现出自西向东逐渐偏正的变化趋势,显示二次蒸发作用的影响逐渐下降,且局地蒸发水汽团对西部地区大气降水影响明显,如位于西风环流影响区的阿拉善左旗等区域局地蒸发气团占到8月部分单次大气降水来源水汽团的100%,而7月份,东亚夏季风环流对内蒙古局地大气降水的影响最为明显;整体上,虽然区域大气湿度变化引起的二次蒸发是影响内蒙古局地大气降水过程的一个关键因素,不过大气降水量对夏季大气降水稳定同位素组成的影响最明显.即西风环流影响区大气湿度变化对氘盈余指数d值的影响程度明显强于东亚夏季风区大气湿度变化的影响,而东亚夏季风环流影响区大气降水量对d值的影响程度则相对明显.  相似文献   

14.
长沙大气降水中稳定同位素变化及过量氘指示水汽来源   总被引:3,自引:1,他引:2  
基于2010 年1 月1 日至2012 年5 月31 日长沙日降水同位素资料,对长沙降水中稳定同位素、过量氘(记为d) 的变化特征以及它们与降水、温度和湿度的关系进行了分析。研究发现:①在季风系统下, 降水中稳定同位素、d 均具有明显的季节变化,表现出夏半年的低值与冬半年的高值交替变化的特点;②降水中稳定同位素在不同时段均存在降水量效应和湿度效应,另外,夏半年中表现为反温度效应,冬半年则表现出温度效应;③综合考虑降水中d 和δ18O 与大气湿度的关系, 可认为影响长沙降水中稳定同位素变化的主要原因与降水的气团性质有关。夏半年中,夏季风盛行,降水的水汽主要来源于西南季风、东南季风所携带的海洋水汽,空气湿度大,d 小,重同位素在水汽输送过程中因不断冷凝而大大贫化,从而降水中δ18O(δD) 较小;冬半年中,长沙受大陆性气团控制,降水的水汽则主要来源于西风带及当地蒸发,空气湿度小,降水中δ18O(δD) 以及d 均较高。  相似文献   

15.
大气降水的氢氧同位素不仅在长时间尺度和短时间尺度上对气候环境有一定的响应,而且在极端天气事件中也响应强烈。论文基于2014年7月23—24日10号台风麦德姆登陆前后福州降水的氢氧稳定同位素分析,发现台风麦德姆降水的δ18O变化范围为-7.2‰~-15.4‰,整个降雨过程分为3个阶段:阶段1(23日3:00至23日12:00),雨水δ18O相对偏正(加权平均值为-7.8‰),降水开始增加,主要受台风外围环流的影响;阶段2(23日12:00至24日8:00),雨水δ18O显著偏负(加权平均值为-13.9‰),采样点降水量达134.8 mm,约占观测时段降水量的78%,降雨量效应明显,主要受台风云雨区影响;阶段3(24日8:00至24日17:00),雨水δ18O偏正(加权平均值为-9.2‰),降水减少,台风远离福州。通过降水δ18O的阶段性特征以及过量氘,结合水汽输送通量与HYSPLIT-4模拟水汽轨迹分析,确定此次降水的水汽源主要来自西北太平洋水汽通道与孟加拉湾和南海水汽通道。  相似文献   

16.
青海湖流域降水和河水中δ18O和δD变化特征   总被引:1,自引:0,他引:1  
稳定同位素方法已被广泛用来研究水循环过程中的水汽来源分析、不同水体间补给关系及水量平衡。论文基于青海湖流域2012 年夏季所收集的河水和逐次大气降水中δ18O和δD 及实测的气象数据,分析了它们的氢氧稳定同位素时空变化特征。结果表明:降水中同位素值组成变化存在两个明显的变化阶段,在8 月中旬之前降水中稳定同位素值较低,而之后明显偏高,这与受夏季风和局地再循环水汽的影响有关;该流域大气降水线(LMWL)的斜率(8.69)和截距(17.5)明显大于全球大气降水线(GMWL),表明青海湖流域夏季大气降水在一定程度上受夏季风所携带湿润海洋性气团的影响;降水中稳定同位素与降水量及温度之间存在负相关关系,表明青海湖流域夏季降水中稳定同位素存在显著的降水量效应,但却不存在温度效应;该流域内河水受到大气降水、地下水及冰川融水的混合和调节作用,这使得河水中稳定同位素的波动范围比降水小;通过瑞利分馏模型模拟了降水中稳定同位素的变化,降水中δ18O、δD和过量氘 (d_excess)的值受青海湖湖面再循环水汽补充的影响。因此,这也将为今后展开青海湖流域水循环过程中大气降水-地下水-地表水的研究提供科学依据,掌握该流域不同水体间的转化关系及水资源利用、管理具有重要的意义。  相似文献   

17.
庐山地区大气降水中稳定同位素变化特征   总被引:3,自引:0,他引:3  
氢氧稳定同位素技术被广泛用来研究水循环过程中的水汽来源、水量平衡及不同水体间的补给关系。以2016年4月至2017年4月在庐山地区三个不同研究点(庐山西北面莲花镇,山顶牯岭镇和东南面海会镇)采集的102个次降水样品同位素资料为基础,应用线性回归分析和对比分析等方法,对庐山地区大气降水中氢氧稳定同位素和氘盈余的时空分布特征及大气水汽来源进行了研究。结果发现:庐山地区夏半年降水中的稳定同位素值δ18O平均值(-6.1‰)小于冬半年(-4.8‰);氢氧同位素特征和氘盈余呈现明显的季节差异;平均氘盈余值(10.6‰)大于全球大部分地区的评估值(10.0‰);当地大气降水线(LMWL)δD=7.45δ18O+8.36与全球大气降水线(GMWL)δD=8δ18O+10相比,其斜率和截距均偏小。结合HYSPLIT后向轨迹模型分析同位素特征发现,庐山地区大气水汽夏半年主要来源于低纬度南海和印度洋,冬半年来自于干燥的华北和西北内陆;局地水汽影响和地理位置差异导致了降雨同位素特征的空间差异性。本研究可为今后展开庐山地区水循环过程的研究提供科学依据。  相似文献   

18.
以影响我国大陆干湿状况的孟加拉湾和南海两股水汽的疑似交界影响区域之一的云南省、广西壮族自治区为研究区,借助稳定性同位素质谱仪MAT253测定2014年16个点、239个有效样本的雨季大气降水的氢氧稳定同位素组成,完成了基于GIS平台的δD、δ18O空间格局分析,实现了中国夏季风西南水汽和东南水汽的交互区域界定。主要研究结果有三点:(1)2014年雨季及6月中下旬一次降水过程δD和δ18O空间变化格局基本一致,因为夏季风大气降水δD和δ18O均沿水汽输送路径不断衰减;(2)大气降水氢氧稳定同位素空间分布主要受降雨量效应和大陆效应的影响,哀牢山高大地形阻隔及云南高原正地形水汽截留作用可能是其数值发生突变的主要原因;(3)2014年夏季的西南水汽在越过哀牢山后与东南水汽在红河、个旧、蒙自附近交互影响使δD和δ18O发生显著变化,6月中下旬一次降水过程中的氢氧同位素数值在红河、个旧附近达到最低,两者互为验证红河、个旧应该是西南水汽和东南水汽影响区域分界。  相似文献   

19.
长沙地区近地面水汽中氢氧稳定同位素的变化特征   总被引:4,自引:0,他引:4  
基于在长沙地区2014年10月—2015年10月监测的近地面水汽中氢氧稳定同位素组成(δv)及相关气象要素,对水汽中氢氧稳定同位素的变化特征、影响因素及与降水中稳定同位素(δp)的相互关系进行了分析.结果表明:大气水汽中氢氧稳定同位素存在显著的季节变化和日变化,冬、春季值高于夏、秋季值,夜晚值高于白天值.δv~(18)O的季节变化与大尺度水汽输送的季节性变化有关,日变化则与地表蒸散发、大气湍流等局地气象条件有关.通过对δv~(18)O的平衡模拟发现,水汽中和降水中稳定同位素在暖季处于或接近平衡状态,在冷季处于非平衡状态.不同季节的大气水汽线和大气水线具有一定程度的相似性,两者的斜率均为暖季大于冷季;受下垫面新降水蒸发的影响,降水日大气水汽线的斜率和截距相对于无降水日均有增加,暖季分别增加0.11和3.52‰,冷季分别增加0.07和0.14‰.  相似文献   

20.
以漓江流域为例查明不同水体之间氢氧稳定同位素组成特征,并探讨氢氧稳定同位素对漓江流域的蒸发过程指示作用.结果表明:流域内不同水体之间,氢氧稳定同位素表现出不同的组成特征,地表水和地下水,在丰水期δD值和δ18O值要比枯水期更容易富集,地下水的δD值和δ18O值分布范围较地表水小;随着水温的升高,δ18O的变化趋势比d-excess明显.从漓江上游到下游高程逐渐下降,河水线的斜率和截距也在逐渐减小,其中漓江下游河水线的斜率和截距要低于当地大气降水线,表明下游受到蒸发作用较强烈;地下河水线、地表河水线在一定程度上偏离当地大气降水线,但偏离程度较小,表明三者之间有很好的水力联系.受温度和湿度的共同影响,漓江干流丰水期河水的蒸发量占最初水体总量的0.7%~9.1%,枯水期河水的蒸发量占最初水体总量的2.6%~9.7%,丰水期的蒸发比例低于枯水期,从上游到下游蒸发比例在逐渐上升.研究区蒸发量估算值为959.40mm,与多年实测值少43.11mm,相对误差4.70%.氢氧稳定同位素对研究区降水、地表水、地下水之间的转换规律具有重要的实际意义,在今后的漓江流域水文研究中有着更加广阔的空间.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司    京ICP备09084417号-23

京公网安备 11010802026262号