CN105115998B - 一种粉末衍射择优取向参数确定方法 - Google Patents

一种粉末衍射择优取向参数确定方法 Download PDF

Info

Publication number
CN105115998B
CN105115998B CN201510516202.6A CN201510516202A CN105115998B CN 105115998 B CN105115998 B CN 105115998B CN 201510516202 A CN201510516202 A CN 201510516202A CN 105115998 B CN105115998 B CN 105115998B
Authority
CN
China
Prior art keywords
preferred orientation
crystal
distributed constant
function
diffraction
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201510516202.6A
Other languages
English (en)
Other versions
CN105115998A (zh
Inventor
王慎文
周建文
陈磊
蒲洪果
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Exploration and Development Research Institute of Sinopec Henan Oilfield Branch Co filed Critical China Petroleum and Chemical Corp
Priority to CN201510516202.6A priority Critical patent/CN105115998B/zh
Publication of CN105115998A publication Critical patent/CN105115998A/zh
Application granted granted Critical
Publication of CN105115998B publication Critical patent/CN105115998B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Analysing Materials By The Use Of Radiation (AREA)
  • Developing Agents For Electrophotography (AREA)

Abstract

本发明涉及一种粉末衍射择优取向参数确定方法,本发明以粉末样品定向压实时产生择优取向为模型,确定择优取向函数,利用所确定的函数计算每个晶面的择优取向函数值;然后根据每个晶面的择优取向函数值得到晶体结构的计算谱图;比较计算谱图与实际谱图之间的误差,调整所设定的择优取向分布参数,直至误差小于设定值或不再减小,以调整后的择优取向分布参数作为该晶体的择优取向分布参数;最后利用该分布参数计算得到的择优取向值即为本发明所要确定的择优取向。本发明充分考虑了衍射强度与晶面取向分布概率成正比和受择优取向影响的总体衍射强度恒定个因素,所确定的择优取向函数更加符合粉末衍射的实际情况,所得到的择优取向也更加精准。

Description

一种粉末衍射择优取向参数确定方法
技术领域
本发明涉及一种粉末衍射择优取向参数确定方法,属于技术领域。
背景技术
在一般多晶体中,每个晶粒有不同于邻晶的结晶学取向,从整体看,所有晶粒的取向是任意分布的;某些情况下,晶体的晶粒在不同程度上围绕某些特殊的取向排列,就称为择优取向或简称织构。择优取向会影响粉末衍射谱图的形态,在有关粉末衍射谱图的晶体学反演计算中,正确地描述择优取向的影响是非常的必要,择优取向对晶体学反演的影响如下:
Ip(hkl)=P(hkl)*I(hkl)
其中I(hkl)为原始计算衍射强度,Ip(hkl)为考虑择优取向影响后的衍射强度,P(hkl)为择优取向函数,择优取向函数描述了择优取向对晶面hkl衍射强度的影响。
P(hkl)=f(a,G)
其中为晶面hkl与衍射面的夹角,为锐角;G为择优取向分布状态参数,其意义主要包括择优取向的集中程度或强度。目前择优函数具体形式有以下6种,分别为D1~D6.
D1=exp(-Ga^2)
D2=exp(G(pi/2-a)^2)
D3=exp(-Gsin^2(a))
D4=exp(-G(1-cos^2(a))
D5=c+(1-c)exp(-Ga^2)
D6=(G^2cos^2(a)+sin^2(a)/G)^-3/2
Bragg-Brentano几何衍射面就是平行于样品台的平面,根据研究可知,只有平行衍射面的晶面才可以产生有效衍射,因此,在研究粉末样品择优取向时可只考虑法线包含于光路所在平面的晶面。在无择优取向时,晶面均匀分布,
P(hkl)=f(a,G)=1,-pi/2≤a≤pi/2
对上式进行积分可得
只考虑择优取向的影响,对于一个无择优取向的样品,假设一种特定情况,原来取向为a1的部分晶面,在有择优取向时,取向改变为a2,那么,在a1处P(hkl)应该有所降低,而在a2处P(hkl)应该有等量的增加,即择优取向函数的积分值应该恒定,如图1所示。而目前给出的择优取向函数均不满足上述要求,如图2所示,该图中给出了W.A.DOLLASE总结的5个单参数函数a>0右半边图形。目前的给出的择优取向函数均没有考虑受择优取向影响的总体衍射强度恒定的因素,导致所确定的择优取向不够准确。
发明内容
本发明的目的是提供一种粉末衍射择优取向参数确定方法,以解决目前由于没有考虑受择优取向影响的总体衍射强度恒定的因素导致所确定的择优取向不准确的问题。
本发明为解决上述技术问题提供了一种粉末衍射择优取向确定方法,该确定方法包括以下步骤:
1)测量某种晶体的单相衍射谱图,根据该衍射谱图确定晶体的结构;
2)根据晶体的结构分别计算晶体每个晶面与衍射面的夹角;
3)以粉末样品定向压实时产生择优取向为模型,确定择优取向函数,利用该函数和每个晶面与衍射面的夹角以及设定的择优取向分布参数计算每个晶面的择优取向函数值P,
P=1/(Gcos2(a)+sin2(a)/G)
其中G为择优取向分布状态,0<G≤1,a为压实后晶面与衍射面的夹角;
4)根据每个晶面的择优取向函数值P得到晶体结构的计算谱图;
5)计算谱图与实际谱图之间的误差,判断误差是否小于设定值,若小于设定值,以此时所选定的择优取向分布参数作为该晶体的择优取向分布参数;
6)调整所设定的择优取向分布参数,重复步骤3)-5),直至误差小于设定值,以调整后的择优取向分布参数作为该晶体的择优取向分布参数;
7)以晶体的择优取向分布参数计算晶体衍射择优取向,该结果即为所要确定的择优取向值。
所述步骤3)中择优取向函数P的确定过程如下:
a)测量未压实时的晶面与衍射面的夹角a和定向压实时晶面与衍射面的夹角b;
b)以压实的比例参数H作为择优取向分布参数G,0<H≤1,H为0时表示最强压实,H为1时表示无压实;
c)根据夹角a和夹角b与比例参数H之间的位置关系计算夹角b概率密度函数,所计算出的概率密度函数即为所求的择优取向函数。
所述步骤4)中误差采用均方误差法计算。
本发明的有益效果是:本发明以粉末样品定向压实时产生择优取向为模型,重新确定择优取向函数,利用所确定的函数和每个晶面与衍射面的夹角以及设定的择优取向分布参数计算每个晶面的择优取向函数值P;然后根据每个晶面的择优取向函数值P得到晶体结构的计算谱图;比较计算谱图与实际谱图之间的误差,调整所设定的择优取向分布参数,直至误差小于设定值,以调整后的择优取向分布参数作为该晶体的择优取向分布参数;最后根据所确定的晶体择优取向分布参数计算该晶体的择优取向值即为本发明所要确定的择优取向。本发明充分考虑了受择优取向影响的总体衍射强度恒定的因素,所确定的择优取向函数更加符合粉末衍射的实际情况,所得到的择优取向更加精准。
附图说明
图1是局部择优取向改变示意图;
图2是现有技术中5个择优函数的分布示意图;
图3是定向施压产生择优取向示意图;
图4是本发明所确定的择优函数的分布示意图。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步的说明。
本发明以粉末样品定向压实产生择优取向为模型,得到一个新的择优取向函数,与以往的择优取向函数相比,具备函数与晶面取向分布概率成正比和函数积分值恒定这两个特征,使根据该函数计算得到择优取向值更加精准。下面以某种晶体为例来说明本发明的具体实施过程。
1.测量某种晶体的单相衍射谱图,根据该衍射谱图确定晶体的结构。
2.根据晶体的结构分别计算晶体每个晶面与衍射面的夹角。
3.以粉末样品定向压实时产生择优取向为模型,重新确定择优取向函数,利用该函数和每个晶面与衍射面的夹角以及设定的择优取向分布参数计算每个晶面的择优取向函数值。
众所周知,在有片状、针状晶粒,如云母存在时,粉末样品定向压实会产生择优取向,如图3所示,只考虑其法线包含于光路所在平面的晶面,径向的直线表示片状晶面,垂直方向压缩了以后,会产生水平方向上的择优取向,是一种常见的择优取向分布状态。本发明以上述择优取向的分布状态为模型,以片状择优取向为例来说明择优取向函数的确定过程。
首先测量未压实时的晶面hkl与衍射面的夹角a和定向压实时晶面hkl与衍射面的夹角b,a和b均为锐角;
然后以压实的比例参数H作为择优取向分布参数G,G=H,0<G≤1,G为0时表示最强压实,G为1时表示无压实。在本实施例中G表示颓垣短轴长度与长轴长度的比值。
根据夹角a和夹角b与比例参数G之间的位置关系计算夹角b概率密度函数,所计算出的概率密度函数即为所求的择优取向函数。
b=atan(G*tan(a))
已知a在[0,pi/2]区间均匀分布,求解b的概率密度函数,以求解得到的b的概率密度函数作为择优取向函数P
P=1/(Gcos2(a)+sin2(a)/G)
其中G为择优取向分布状态,0<G≤1,a为晶面与衍射面的夹角。
其图形如图4所示,由一个平台和叠加的一个对称峰形组成,G值增大,则平台升高,峰形降低,以蒙特-卡罗方法模拟100000个数据点统计验证函数,结论符合。
4.根据每个晶面的择优取向函数值P得到晶体结构的计算谱图,该谱图计算过程为行业内众所周知的,属于常规方法,这里不再详述。
5.比较计算谱图与实际谱图之间的误差,这里的误差计算采用均方误差,判断误差是否小于设定值,若小于设定值,以此时所选定的择优取向分布参数作为该晶体的择优取向分布参数;
6.调整所设定的择优取向分布参数,重复步骤3-5,直至误差小于设定值或不再减小,以调整后的择优取向分布参数作为该晶体的择优取向分布参数。以上实施过程针对的是片状择优取向,对于线状择优取向也可以采用上述方法,只需将线状择优取向理解为是关于线轴法平面的片状择优取向,只需将片状择优取向函数中的夹角a替换为其余角pi/2-a即可。

Claims (3)

1.一种粉末衍射择优取向确定方法,其特征在于,该确定方法包括以下步骤:
1)测量某种晶体的单相衍射谱图,根据该衍射谱图确定晶体的结构;
2)根据晶体的结构分别计算晶体每个晶面与衍射面的夹角;
3)以粉末样品定向压实时产生择优取向为模型,确定择优取向函数,利用该函数和每个晶面与衍射面的夹角以及设定的择优取向分布参数计算每个晶面的择优取向函数值P,
P=1/(Gcos2(a)+sin2(a)/G)
其中G为择优取向分布状态,0<G≤1,a为压实后晶面与衍射面的夹角;
4)根据每个晶面的择优取向函数值P得到晶体结构的计算谱图;
5)计算谱图与实际谱图之间的误差,判断误差是否小于设定值,若小于设定值,以此时所选定的择优取向分布参数作为该晶体的择优取向分布参数;
6)调整所设定的择优取向分布参数,重复步骤3)-5),直至误差小于设定值,以调整后的择优取向分布参数作为该晶体的择优取向分布参数;
7)以晶体的择优取向分布参数计算晶体衍射择优取向,该结果即为所要确定的择优取向值。
2.根据权利要求1所述粉末衍射择优取向确定方法,其特征在于,所述步骤3)中择优取向函数P的确定过程如下:
a)测量未压实时的晶面与衍射面的夹角a和定向压实时晶面与衍射面的夹角b;
b)以压实的比例参数H作为择优取向分布参数G,0<H≤1,H为0时表示最强压实,H为1时表示无压实;
c)根据夹角a和夹角b与比例参数H之间的位置关系计算夹角b概率密度函数,所计算出的概率密度函数即为所求的择优取向函数。
3.根据权利要求2所述的粉末衍射择优取向参数确定方法,其特征在于,所述步骤4)中误差采用均方误差法计算得到。
CN201510516202.6A 2015-08-20 2015-08-20 一种粉末衍射择优取向参数确定方法 Active CN105115998B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510516202.6A CN105115998B (zh) 2015-08-20 2015-08-20 一种粉末衍射择优取向参数确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510516202.6A CN105115998B (zh) 2015-08-20 2015-08-20 一种粉末衍射择优取向参数确定方法

Publications (2)

Publication Number Publication Date
CN105115998A CN105115998A (zh) 2015-12-02
CN105115998B true CN105115998B (zh) 2017-11-28

Family

ID=54664038

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510516202.6A Active CN105115998B (zh) 2015-08-20 2015-08-20 一种粉末衍射择优取向参数确定方法

Country Status (1)

Country Link
CN (1) CN105115998B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107590296A (zh) * 2016-07-08 2018-01-16 深圳大学 一种小角x射线散射的全谱拟合方法及系统
CN108896591B (zh) * 2018-05-11 2020-10-16 深圳晶泰科技有限公司 一种晶体粉末衍射的指标化方法
CN109916928B (zh) * 2019-03-15 2020-11-17 浙江大学 一种判断各向异性纳米晶体择优取向性组装结果的方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA822213A (en) * 1969-09-02 Esso Research And Engineering Company Elimination of the preferred orientation distortions of an x-ray diffraction pattern
US3596092A (en) * 1968-05-10 1971-07-27 Rolls Royce X-ray diffraction method for determining the value of a mechanical property or properties of a fiber
EP1365230A2 (en) * 2002-05-23 2003-11-26 Rigaku Corporation Method for estimating preferred orientation of polycrystalline material
WO2011053692A1 (en) * 2009-10-29 2011-05-05 Glaxo Group Limited Sample plate for an x-ray powder diffraction apparatus
CN102721712A (zh) * 2012-04-27 2012-10-10 中国航空工业集团公司北京航空材料研究院 一种用于测量金属管材织构的方法和装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA822213A (en) * 1969-09-02 Esso Research And Engineering Company Elimination of the preferred orientation distortions of an x-ray diffraction pattern
US3596092A (en) * 1968-05-10 1971-07-27 Rolls Royce X-ray diffraction method for determining the value of a mechanical property or properties of a fiber
EP1365230A2 (en) * 2002-05-23 2003-11-26 Rigaku Corporation Method for estimating preferred orientation of polycrystalline material
WO2011053692A1 (en) * 2009-10-29 2011-05-05 Glaxo Group Limited Sample plate for an x-ray powder diffraction apparatus
CN102721712A (zh) * 2012-04-27 2012-10-10 中国航空工业集团公司北京航空材料研究院 一种用于测量金属管材织构的方法和装置

Also Published As

Publication number Publication date
CN105115998A (zh) 2015-12-02

Similar Documents

Publication Publication Date Title
CN105115998B (zh) 一种粉末衍射择优取向参数确定方法
CN107085148B (zh) 一种并联电容器组的在线容限预警方法及装置
Tang et al. Vibration analysis based on empirical mode decomposition and partial least square
CN110685651B (zh) 一种多层合采气井产量劈分方法及系统
CN103617816A (zh) 反应堆堆芯功率分布的测量方法
TW200849345A (en) Dual-phase virtual metrology method
CN106786567B (zh) 一种基于pmu类噪声数据的在线负荷建模方法
CN103604821A (zh) 一种钢中奥氏体含量的测量方法
US11402341B2 (en) Quantitative phase analysis device for analyzing non-crystalline phases, quantitative phase analysis method for analyzing Non-Crystalline phases, and non-transitory computer-readable storage medium storing quantitative phase analysis program for analyzing Non-Crystalline Phases
CN105141253A (zh) 基于光伏出力曲线斜率的光伏窃电识别方法
CN106684866A (zh) 一种静态开环综合增益的计算方法及装置
CN109783934A (zh) 一种基于h-adcp的断面平均流速拟合率定方法
CN105696064B (zh) 一种图形电镀参数的获取方法
CN109444537B (zh) 一种计及带外干扰的自适应同步相量测量方法
CN104655496B (zh) 自重对离轴反射镜面形影响的测试方法
Hazebroek et al. The optimum adjustment of regulators
WO2012093400A9 (en) Method and system for use in measuring in complex patterned structures
CN104614144B (zh) 泄洪消能诱发场地振动的预测方法
Nadal et al. Statistical methods for analyzing color difference distributions
CN112766630A (zh) 一种评价低渗气藏气井无阻流量的方法
CN108984484B (zh) 基于回归分析评估双面双玻光伏组件辐照增益的方法
Saha et al. Sample shifting technique (SST) for estimation of harmonic power in polluted environment
CN109143126A (zh) 一种硬质合金矫顽磁力或磁饱和标准样品的制备方法
CN109145258A (zh) 基于非线性拟合的威布尔分布参数置信区间估计方法
CN110782100B (zh) 一种低渗透气藏产能快速预测方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant