CN102289002A - 一种获取地球背景场重力加速度的方法 - Google Patents
一种获取地球背景场重力加速度的方法 Download PDFInfo
- Publication number
- CN102289002A CN102289002A CN2011102032959A CN201110203295A CN102289002A CN 102289002 A CN102289002 A CN 102289002A CN 2011102032959 A CN2011102032959 A CN 2011102032959A CN 201110203295 A CN201110203295 A CN 201110203295A CN 102289002 A CN102289002 A CN 102289002A
- Authority
- CN
- China
- Prior art keywords
- point
- time
- falling
- straight line
- speed
- 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.)
- Granted
Links
Images
Landscapes
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明涉及一种获取地球背景场高精度重力加速度的方法,包括:步骤1:同步记录绝对重力仪中落体自由下落过程中包含落体时间位移信息的干涉带信号和参考棱镜振动的速度信息;步骤2:提取落体下落轨迹上的时间位移坐标,即落体自由下落轨迹的重建点;步骤3:利用一簇斜率相同的平行直线,与参考棱镜速度时间坐标曲线相交,获得这些相交的交点所对应的时间点;步骤4:在落体时间位移坐标曲线上,选出与步骤2中获得的时间点所对应的重建点;步骤5:根据步骤4中所挑选出的重建点计算重力加速度。本发明的方法可以有效的去除地面振动的影响,整体测量结构避免了庞大的隔振系统,为绝对重力仪的小型化和流动测量提供了方法支持。
Description
技术领域
本发明涉及一种利用绝对重力仪获取重力加速度的方法,采用画线选点去除测量结果的振动误差。
背景技术
重力场是反映地球内部物质结构及其变迁的地球物理基本场,高精度重力观测资料是地震监测预报、地球科学研究、资源勘探等领域的基础,同时也是国家精密计量、飞行器导航和导弹制导等不可或缺的战略数据。激光干涉绝对重力仪是目前绝对重力仪中的一种主流设计方式,是获取微伽级(重力单位在厘米克秒制中的计量单位是厘米/秒2,为了纪念意大利人伽利略,又把厘米/秒2称为伽)绝对重力数据的一种重要的观测设备。国际计量局重力实验室与法国S.A.公司合作,于1983年研制了GA60型上抛法绝对重力仪,精度为10微伽,但是仪器大而笨重,没有得到推广。同时,中国、美国、俄罗斯、意大利、德国、乌克兰等国也相继开展了相关研究,目前唯一商用的绝对重力仪是美国Micro-G公司生产的FG5、FGL以及A10三种类型的绝对重力仪。1999年,美国斯坦福大学的诺贝尔奖获得者朱棣文教授领导的研究小组,以原子干涉技术,实现了绝对重力测量。将绝对重力观测由宏观引入微观层面,开创了一个新的领域。
激光干涉绝对重力仪主要包括落体控制机构、真空系统、干涉测量系统、地面振动隔离系统和数据采集与处理五大模块,图1为一种激光干涉绝对重力仪的实验样机的结构图。而绝对重力仪的主体包含了其中的落体控制机构和真空系统,落体控制机构在真空系统中控制落体做自由下落运动,其动作过程由计算机控制完成。激光干涉仪是干涉测量系统的具体实现,完成了测量光路的建立。地面振动隔离系统安装在干涉仪内部(例如中国的NIM系列激光干涉绝对重力仪)或者支撑整个干涉仪(例如FG5系列激光干涉绝对重力仪)。数据采集与处理包含了接收干涉信号的光电接收器、将干涉信号数字化的高速数字化仪。同时铷频标为高速数字化仪提供时间基准,而整个测量的长度基准由高稳频激光的波长(如碘吸收激光器)提供。数据处理算法主要完成从干涉信号中提取落体的时间-位移坐标以及重力加速度值的计算。
目前,激光干涉绝对重力仪是通过拟合落体下落时相对于参考棱镜的时间位移坐标,来提取测点的重力加速度值。但是由于来自落体控制机构产生的系统自振和来自自然环境中的大地脉动的影响,使得参考棱镜在落体下落的过程中呈现非线性运动。这样就使得提取的落体位移信息不仅包含了重力场对落体的作用,还包含了地面振动对参考棱镜的影响,这种影响由于其自身的非线性,带来的误差会传递到最终的测量结果。目前主要的处理方法是利用长周期弹簧将振动隔离。但是由于弹簧本身的制造工艺限制,使得隔振系统的体积变得庞大,不利于绝对重力仪的小型化和流动测量。
发明内容
本发明所要解决的技术问题是提供一种新的获取重力加速度的方法,以在不增大绝对重力仪样机体积的条件下,实现测量结果的振动误差剔除。
本发明解决上述技术问题的技术方案如下:
一种获取地球背景场重力加速度的方法,包括以下步骤:
步骤1:同步记录绝对重力仪中落体下落过程中包含落体时间位移信息干涉带信号和参考棱镜振动的速度时间信息;
步骤2:利用干涉带信号提取落体下落轨迹上的时间位移坐标,即落体自由下落轨迹的重建点;
步骤3:利用一簇斜率相同的平行直线,与所述参考棱镜振动的速度时间信息所对应的参考棱镜速度时间坐标曲线相交,获得这些相交的交点所对应的时间点;
步骤4:在所述落体的时间位移信息所对应的落体时间位移坐标曲线上,选出与步骤2中获得的时间点所对应的重建点;
步骤5:根据步骤4中所挑选出的重建点计算重力加速度。
本发明的有益效果是:利用直线的二次微分为零的基本原理,挑选落体时间位移坐标曲线上的重建点,在计算重力加速度时可以有效的去除地面振动的影响,整体测量结构避免了庞大的隔振系统,为绝对重力仪的小型化和流动测量提供了方法支持。
在上述技术方案的基础上,本发明还可以做如下改进。
进一步,步骤3中,所述斜率为所述参考棱镜速度时间坐标曲线上第一个采样点与最后一个采样点连成的直线的斜率。
采用上述进一步方案的有益效果是可以确定该速度时间曲线的线性变化趋势,确定下一步中所需的直线簇斜率。
进一步,步骤3中,所述平行直线满足以下条件:所述参考棱镜速度时间坐标曲线上的至少两个采样点均在同一条直线上;或者至少有一个采样点在该直线上,但不在该直线上的采样点速度坐标与和该直线对应与该采样点时刻的纵坐标之间的误差在所述参考棱镜速度时间坐标曲线峰峰值的1%~10%以内,其中所述交点为该同一条直线和所述参考棱镜速度时间坐标曲线的交点。
采用上述进一步方案的有益效果是,利用所确定的直线簇斜率,在误差范围内挑选落体自由下落轨迹的重建点,以提高下一步的计算精度。
进一步,步骤5中,对所述重建点采用最小二乘拟合的方法,计算重力加速度。
附图说明
图1为本发明方法中所使用的绝对重力仪结构示意图;
图2为本发明的画线选点获取重力加速度的方法的步骤流程图;
图3为本发明中的测量光路示意图;
图4为采用本发明的方法在H05测点的测量结果;
图5为采用本发明的方法在H08测点的测量结果;
图6为采用本发明的方法在H03测点的测量结果。
附图中,各标号所代表的部件列表如下:
LS、激光器,BS、分光镜,M1、落体,M2、参考棱镜,0、光束,L1、测量光束,L2、测量光束,L3、参考光束,1、重合后的光束
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
如图2所示,本发明的获取地球背景场重力加速度的方法,主要包括以下步骤:
步骤1:同步记录绝对重力仪中落体自由下落过程中包含落体的时间位移信息的干涉带信号和参考棱镜振动的速度时间信息;
步骤2:利用干涉带信号提取落体下落轨迹上的时间位移坐标,即落体自由下落轨迹的重建点;
步骤3:利用一簇斜率相同的平行直线,与所述参考棱镜振动的速度时间信息所对应的参考棱镜速度时间坐标曲线相交,获得这些相交的交点所对应的时间点;
步骤4:在所述落体的时间位移信息所对应的落体时间位移坐标曲线上,选出与步骤2中获得的时间点所对应的重建点;
步骤5:根据步骤3中所挑选出的重建点计算重力加速度。
以下对上述步骤进行细致阐述。
本发明的获取地球背景场重力加速度的方法目的是要剔除振动给重力加速度的计算所带来的影响,以提高计算重力加速度的精确。因此从另外一个角度来说,本发明的方法也可以叫做一种震动处理方法。其整个过程如下。
将绝对重力仪中的参考棱镜安装在频带为120秒~50Hz的高灵敏度甚宽频带拾振器上。启动绝对重力仪后,如图1、图3所示,高稳频激光器LS发射的光束0经过处理,入射分光镜BS,产生测量光束L1和参考光束L3。其中,光束0从分光镜BS上表面反射成为测量光束L1,经落体M1反射成为测量光束L2并入射分光镜BS上表面;光束0从分光镜BS上表面入射,从分光镜BS下表面透射成参考光束L3,参考光束L3经过参考棱镜M2反射回分光镜BS的下表面。测量光束L2经过分光镜BS透射以及参考光束L3经过分光镜BS反射后,测量光束L2和参考光束L3重合并产生干涉,成为重合后的光束1,该光束1由光电接收器R接收后再转换为电信号,将之称为干涉带信号。同时,参考棱镜M2的振动由拾振器采集,将之称为振动信号。干涉带信号和振动信号,同时发送给高速数字化仪进行采集。
由于接收到的干涉带信号可以认为是一些列频率逐渐增加的调频正弦波,因此通过过零点算法(经过去均值后提取每个周期的正弦波信号过零点精确时刻的算法称为过零点算法)从干涉带信号中提取落体下落轨迹上的时间位移坐标。当得到干涉带的过零点后,可以得到每个过零点的时刻,然后均分这些过零点,将每N个过零点作为一段。落体相对于参考棱镜每下落半个激光波长λ,产生一条干涉带,即一个周期的调频正弦波,因此,每N个过零点就等于落体相对于参考境下落了 的距离。这样就可以计算落体下落过程中每下落S对应的时间点,得到落体下落过程中的时间位移坐标,即落体轨迹的重建点。
从振动信号中提取速度时间信息,生成速度时间曲线。由于参考棱镜和拾振器是直接安装在绝对重力仪中的,测量时绝对重力仪又是放置于地面的,所以所获得的振动信号也即是地面的震动信号。利用一簇斜率相同的平行直线与该速度时间曲线相交,挑选与这些相交的交点时间所对应的落体时间位移坐标作为重建点,最后对这些重建点进行最小二乘拟合,计算所测量地点的重力加速度值。其理论过程如下:
首先设测量出落体下落过程中的参考棱镜的振动曲线为:
(1)
对(2)式进行两次微分得到:
由于干涉带信号和振动信号的采集是同步的,采集参数完全相同,因此挑选出与该簇平行直线与该振动曲线的交点所处时间所对应的落体轨迹上的时间位移坐标进行最小二乘拟合。在这些时间点上的干涉带波形受到的振动影响可以认为是线性的,在计算重力加速度时,由(3)式可知,振动的影响会随着二次微分的过程自动置零,即可以消除参考棱镜的振动对计算结果的影响。
这里直线簇的斜率一般取振动曲线第一个采样点与最后一个采样点连成的直线的斜率。该簇平行直线需要满足以下条件:振动曲线上的至少两个采样点均在同一条直线上;或者至少有一个采样点在该直线上,但不在该直线上的采样点的速度坐标与该直线对应与该采样点时刻的纵坐标之间的误差在振动曲线的峰峰值的1%~10%以内,其中所述的交点为该同一条直线和该振动曲线的交点。具体来说,直线簇与振动曲线相交时,由于振动曲线是离散性的(因为振动曲线是经过高速数字化仪采集后的数字信号,是离散信号),这里设定一个阈值TH,认为该簇平行直线与振动曲线相交后,纵坐标与相邻采样点之间的误差小于阈值TH时,就认为是一条有效画线。因为振动曲线的离散性,不可能所选的直线簇中的每条直线都与振动信号精确相交,而是在一定误差范围内的相交,因此这里需要设定一个阈值,也就是说满足这个阈值条件的直线都认为是有效的直线。一般来说,阈值TH是根据实验时的测量位置的背景振动情况来确定的。在背景振动大的测量位置,阈值TH一般取振动曲线非线性部分峰峰值的10%,在振动小的测量位置,可以取振动曲线非线性部分峰峰值的1%~5%。
上述方法利用直线的二次微分为零的基本原理,挑选落体轨迹上的重建点,在计算重力加速度过程中能够有效的去除地面振动的影响。整体测量结构避免了庞大的隔振系统,为绝对重力仪的小型化和流动测量提供了方法支持。
本发明中使用的最小二乘方法是一种参数估计算法(具体可参见:方崇智、萧德云的《过程辨识》,2004年清华大学出版社出版,参考页码133-156页),具体过程如下:
考虑到落体自由下落时光速有限性和地球表面的重力梯度值,落体下落的位置可以用下式表示:
其中,??为地球表面重力梯度,为考虑光速有限性后的时间,可以用表达。观测点重力加速度值,为落体初始位置相对与参考棱镜的距离,为第i时刻落体相对于参考棱镜的距离,则对与每次测量中获取的n个时刻对应的落体位移,可以利用(4)式得到如下表达:
(5)式可以写成向量形式:
通过最小二乘估计,可以得到参数向量的最小二乘解:
本发明的方法是一种画线选点拟合的方法,在以下的实施例中,画线选点需要以下五个参数:落体自由下落过程中相对于起始点开始取数的距离L1;落体自由下落过程中相对于起始点结束取数的距离L2;直线簇中每条直线的最小长度L3(以计算的时间位移坐标表示);直线簇中每条直线相对于其他直线的不平行度L4(以计算的时间位移坐标表示);遍历整条振动曲线,寻找每条直线的两个端点时的扫描间隔L5(以计算的时间位移坐标表示)。
以下列出三个计算实例显示画线选点拟合的结果。实验地点的地基分别为较厚土层(H05测点),碎石层(H08测点)和基岩(H03),画线选点拟合法的五个参数取相同的值,分别取L1=1cm,L2=7cm,L3=10,L4=1,L5=5。
算例1
H05测点,测量结果如图4,平行直线簇选取70条直线,选取70对时间位移坐标,采用直接拟合算法计算结果与参考值偏差为4592微伽,而采用本发明的画线选点法结果与参考值偏差为-1451微伽。
算例2
H08测点,测量结果如图5,平行直线簇选取120条直线,选取120对时间位移坐标,采用直接拟合算法计算结果与参考值偏差为1336微伽,而采用本发明的画线选点法结果与参考值偏差为-167微伽。
算例3
H03测点,测量结果如图6,平行直线簇选取93条直线,选取93对时间位移坐标,采用直接拟合算法计算结果与参考值偏差为830微伽,而采用本发明的画线选点法结果与参考值偏差为313微伽。
上述3个算例中,采用本发明的画线选点法因为因为有效的去除地面振动的影响,所以其结果均好于只限拟合算法的结果。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种获取地球背景场重力加速度的方法,其特征在于,包括以下步骤:
步骤1:同步记录绝对重力仪中落体自由下落过程中包含落体时间位移信息的干涉带信号和参考棱镜振动的速度时间信息;
步骤2:利用干涉带信号提取落体下落轨迹上的时间位移坐标,即落体自由下落轨迹的重建点;
步骤3:利用一簇斜率相同的平行直线,与所述参考棱镜振动的速度时间信息所对应的参考棱镜速度时间坐标曲线相交,获得这些相交的交点所对应的时间点;
步骤4:在所述落体的时间位移信息所对应的落体时间位移坐标曲线上,选出与步骤2中获得的时间点所对应的重建点;
步骤5:根据步骤4中所挑选出的重建点计算重力加速度。
2.根据权利要求1所述的获取地球背景场重力加速度的方法,其特征在于:步骤3中,所述斜率为所述参考棱镜速度时间坐标曲线上第一个采样点与最后一个采样点连成的直线的斜率。
3.根据权利要求1或2所述的获取地球背景场重力加速度的方法,其特征在于:步骤3中,所述平行直线满足以下条件:所述参考棱镜速度时间坐标曲线上的至少两个采样点均在直线簇中某一条直线上;或者至少有一个采样点在该直线上,不在该直线上的采样点速度坐标与该直线对应与该采样点时刻的纵坐标之间的误差在所述参考棱镜速度时间坐标曲线的峰峰值的1%~10%以内,其中所述交点为该同一条直线和所述参考棱镜速度时间坐标曲线的交点。
4.根据权利要求1所述的获取地球背景场重力加速度的方法,其特征在于:步骤5中,对所述重建点采用最小二乘拟合的方法,计算重力加速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011102032959A CN102289002B (zh) | 2011-07-20 | 2011-07-20 | 一种获取地球背景场重力加速度的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011102032959A CN102289002B (zh) | 2011-07-20 | 2011-07-20 | 一种获取地球背景场重力加速度的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102289002A true CN102289002A (zh) | 2011-12-21 |
CN102289002B CN102289002B (zh) | 2013-05-22 |
Family
ID=45335562
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011102032959A Expired - Fee Related CN102289002B (zh) | 2011-07-20 | 2011-07-20 | 一种获取地球背景场重力加速度的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102289002B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104749650A (zh) * | 2015-03-25 | 2015-07-01 | 北京航天控制仪器研究所 | 一种基于hom干涉原理的新型重力仪 |
CN106383367A (zh) * | 2015-07-31 | 2017-02-08 | 中国计量科学研究院 | 绝对重力测量方法和装置 |
CN109709620A (zh) * | 2019-02-26 | 2019-05-03 | 中国科学院测量与地球物理研究所 | 一种绝对重力仪 |
CN110132127A (zh) * | 2019-05-23 | 2019-08-16 | 王志文 | 一种天体引力场的测量装置及方法 |
CN111679334A (zh) * | 2020-06-19 | 2020-09-18 | 中国地震局地震预测研究所 | 一种绝对重力加速度的补偿方法和设备 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030136190A1 (en) * | 2002-01-21 | 2003-07-24 | Akito Araya | Micro-distance toss-up type absolute gravimeter |
US20080046185A1 (en) * | 2006-04-10 | 2008-02-21 | Tim Niebauer | Method and apparatus for processing an under-sampled chirped sinusoidal waveform using a complex-heterodyne |
-
2011
- 2011-07-20 CN CN2011102032959A patent/CN102289002B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030136190A1 (en) * | 2002-01-21 | 2003-07-24 | Akito Araya | Micro-distance toss-up type absolute gravimeter |
US20080046185A1 (en) * | 2006-04-10 | 2008-02-21 | Tim Niebauer | Method and apparatus for processing an under-sampled chirped sinusoidal waveform using a complex-heterodyne |
Non-Patent Citations (4)
Title |
---|
《APPLIED OPTICS》 19960701 Eric L.Canuteson等 Fiber-optic extrinsic Fabry-Perot vibration-isolated interferometer for use in absolute gravity meters 第35卷, 第19期 * |
ERIC L.CANUTESON等: "Fiber-optic extrinsic Fabry–Perot vibration-isolated interferometer for use in absolute gravity meters", 《APPLIED OPTICS》 * |
郭允晟等: "符合原理在激光绝对重力仪中的应用", 《计量学报》 * |
郭有光等: "NIM-Ⅱ型可移激光绝对重力仪", 《地球物理学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104749650A (zh) * | 2015-03-25 | 2015-07-01 | 北京航天控制仪器研究所 | 一种基于hom干涉原理的新型重力仪 |
CN104749650B (zh) * | 2015-03-25 | 2017-08-29 | 北京航天控制仪器研究所 | 一种基于hom干涉原理的新型重力仪 |
CN106383367A (zh) * | 2015-07-31 | 2017-02-08 | 中国计量科学研究院 | 绝对重力测量方法和装置 |
CN106383367B (zh) * | 2015-07-31 | 2018-07-10 | 中国计量科学研究院 | 绝对重力测量方法和装置 |
CN109709620A (zh) * | 2019-02-26 | 2019-05-03 | 中国科学院测量与地球物理研究所 | 一种绝对重力仪 |
CN109709620B (zh) * | 2019-02-26 | 2024-05-17 | 中国科学院测量与地球物理研究所 | 一种绝对重力仪 |
CN110132127A (zh) * | 2019-05-23 | 2019-08-16 | 王志文 | 一种天体引力场的测量装置及方法 |
CN111679334A (zh) * | 2020-06-19 | 2020-09-18 | 中国地震局地震预测研究所 | 一种绝对重力加速度的补偿方法和设备 |
Also Published As
Publication number | Publication date |
---|---|
CN102289002B (zh) | 2013-05-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103472494B (zh) | 基于原子干涉效应的重力势三阶微商测量传感器及其方法 | |
CN102289002B (zh) | 一种获取地球背景场重力加速度的方法 | |
CN106814391B (zh) | 基于菲涅尔体层析反演的地面微地震事件定位方法 | |
CN106291725B (zh) | 一种快速反演地下地质体空间位置的方法 | |
CN110133715A (zh) | 一种基于初至时差和波形叠加的微地震震源定位方法 | |
CN203480055U (zh) | 一种基于原子干涉效应的重力势三阶微商测量传感器 | |
CN105242328B (zh) | 古热岩石圈厚度的确定方法及装置 | |
CN107817523B (zh) | 绕射波偏移速度的分析方法及装置 | |
CN106842347B (zh) | 一种阵列式原子干涉重力梯度张量全分量的测量系统 | |
BR112013023985B1 (pt) | Método de estimar um parâmetro de uma anomalia em uma formação de terra e sistema para análise gravitacional | |
BR112012022931B1 (pt) | Método para identificar discrepâncias de temporizações de relógio em pelo menos um relógio de interesse que esteja associado com um receptor sísmico | |
CN110261903B (zh) | 一种基于逆时能量聚焦的地下震源被动定位方法 | |
CN101636669A (zh) | 地形校正系统 | |
CN106842080B (zh) | 一种磁场测量装置姿态摆动干扰去除方法 | |
CN102901985A (zh) | 一种适用于起伏地表的深度域层速度修正方法 | |
CN103792589A (zh) | 重力加速度的测量装置及测量方法 | |
CN105573963B (zh) | 一种电离层水平不均匀结构重构方法 | |
CN104698466B (zh) | 远程动态目标测距装置及方法 | |
CN109100792A (zh) | 基于台站与三维地震联合采集资料的速度反演方法 | |
Qiao et al. | Hybrid seismic-electrical data acquisition station based on cloud technology and green IoT | |
US20170329043A1 (en) | Measuring gravity curvature for mapping subterranean formations | |
CN103364823A (zh) | 震动源实时定位与分析系统 | |
CN105353406B (zh) | 一种生成角道集的方法和装置 | |
Piana Agostinetti et al. | Distributed acoustic sensing as a tool for subsurface mapping and seismic event monitoring: A proof of concept | |
CN109782355A (zh) | Obs检波点漂移的检测方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130522 Termination date: 20200720 |