切片采样算法在光释光年龄模型参数估计中的应用-地理科学进展.pdfVIP

切片采样算法在光释光年龄模型参数估计中的应用-地理科学进展.pdf

  1. 1、本文档共11页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
切片采样算法在光释光年龄模型参数估计中的应用-地理科学进展

第35 卷 第1期 地理科学进展 Vol.35, No.1 2016 年1月 ProgressinGeography Jan. 2016 切片采样算法在光释光年龄模型参数 估计中的应用 1 1 2 彭 俊 ,董治宝 ,韩凤清 (1. 中国科学院寒区旱区环境与工程研究所,兰州730000 ;2. 中国科学院青海盐湖研究所,西宁810008) 摘 要:光释光(OSL)年代学模型是基于数理统计学的一类概率密度模型,它根据特定的假设条件对样品等效剂量 (De)分布进行数学解释,由此估计具有不同沉积历史或者能够代表样品实际埋藏年龄的De 组分。年龄模型参数 估计常通过极大似然估计(MLE)算法实现,本文尝试了切片采样算法在年龄模型参数优化中的应用。切片采样属 于一种马尔科夫链蒙特卡罗采样(MCMC)算法,能根据测量数据与模型的联合似然函数进行随机采样,由此获得 参数的采样分布。本文编写了实现年龄模型切片采样算法的应用程序,并使用模拟及实测De 数据验证了该算法 估计的可靠性。相对于MLE 算法,MCMC 算法具有对参数初值依赖性低、误差估计更准确的特点,切片采样算法 提供了实现释光年龄模型参数估计的一种新方法。 关键词:OSL 测年;等效剂量;年龄模型;参数优化; 随机采样 1 引言 技术而是包含几十到上百粒样品的多颗粒单片 (Rodnight et al, 2006; Schmidt et al, 2012; Keen-Ze- 单片再生剂量法(SAR)(Murray et al, 2000)对样 bert et al, 2013) 。对于De 分布分散的样品,一般不 品不同单片的测量可获得一系列De 值及其误差估 能通过均值计算其埋藏剂量(Galbraith et al, 2012) 。 计。自然样品的De 分布常表现出一定的分散性, 有限组分混合年龄模型(FMM)(Galbraith et al, 这与样品颗粒的侵蚀、搬运、沉积历史与埋藏环境 1990; Roberts et al, 2000) 和最小年龄模型(MAM) 有关。导致样品De 分散的不确定性因素包括内部 (Van der Touw et al, 1997; Galbraith et al, 1999)是两 和外部误差(Duller, 2008) ,前者主要来源于实验操 种应用广泛的统计学年龄模型,其模型参数常使用 作、数据分析处理及仪器精度,包括光子计数统计 MLE 算法来优化。但MLE 算法的缺点在于优化结 不确定性、曲线拟合误差、仪器可重复率、人工贝塔 果对参数的初值敏感,模型容易陷入局部最优解。 源辐射剂量率的均一程度、实验室与自然条件下加 对于多参数模型或高噪声观测数据,MLE 算法需比 热与辐射过程的差异、样品不同颗粒释光特性的差 较以不同初值初始化的优化结果才能确定模型的 异等;后者指样品在搬运与埋藏过程中不同颗粒吸 全局最优解。此外,MLE 算法不能直接估计参数的 收辐射剂量的变化与不确定性,如埋藏前晒退程度 标准误差,参数误差常通过费歇尔信息矩阵来估 不同造成不同颗粒残余剂量的差异、埋藏环境中不 计。对于样本量无限大的样本数据,假定模型的对 同颗粒贝塔剂量率的空间差异、不同埋藏历史与年 数似然函数为二次型形式,则模型参数的方差—协 龄的颗粒发生混合等(Thomsen et al, 2007) 。这些 方差矩阵等价于费歇尔

您可能关注的文档

文档评论(0)

ailuojue + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档