万花镜
    首页社会国际娱乐科技时尚军事汽车探索美食旅游历史健康育儿
    购物旅游历史

    北斗专辑|范磊:多系统GNSS卫星轨道快速积分方法

    2017年6月22日 来源: 测绘学报

    《测绘学报》

    构建与学术的桥梁 拉近与权威的距离

    范磊, 李敏, 宋伟伟, 施闯, 王成

    武汉大学卫星导航定位技术研究中心, 武汉 430079

    收稿日期:2016-11-25; 修回日期:2016-12-20

    基金项目:国家重点研发计划(2016YFB0501800);国家自然科学基金(41231174;41325015;4157040195);湖北省自然科学重点基金(2015CFA057)

    第一作者简介:范磊(1989-), 男, 博士生, 研究方向为GNSS精密定轨理论与方法。

    E-mail:flei_sgg@whu.edu.cn

    通信作者:李敏

    E-mail: lim@whu.edu.cn

    摘要:快速高效且高精度的轨道数值积分算法是多系统GNSS卫星联合快速精密定轨的重要基础。本文从自适应变换Admas积分步长和多卫星同步积分两方面研究了多系统GNSS卫星轨道快速积分方法。为了验证该方法的精度和效率,利用武汉大学(WHU)与欧洲定轨中心(CODE)发布的事后精密星历进行轨道动力学拟合。试验结果表明:GPS/GLONASS/BDS/Galileo 4个系统卫星平均三维RMS均优于20 mm;在不损失传统方法精度的前提下,单颗卫星平均积分与拟合耗时仅需0.09 s,较传统逐颗卫星固定步长积分算法提升了14倍,并且随着卫星数的增加,效率提升越明显。

    关键词:全球卫星导航系统 快速轨道积分算法 轨道拟合 多系统

    A Rapid Orbit Integration Algorithm for Multi-GNSS Satellites

    FAN Lei, LI Min, SONG Weiwei, SHI Chuang, WANG Cheng

    Abstract: A rapid and efficient orbit numerical integration algorithm with high accuracy is needed in multi-GNSS rapid precise orbit determination. In order to improve the compute efficiency, an adaptive step-changed Admas integration method and a synchronous integration algoritm for multi-GNSS satellites are developed in this paper. To validate the precision and efficiency of the proposed method, the multi-GNSS precise orbit products calculated by Wuhan University (WHU) and Center for Orbit Determination in Europe (CODE) are used for orbit fitting. Results show that, the average 3DRMS of GPS, GLONASS, BDS and Galileo satellites are all below 20 mm. Comparing with the traditional step-fixed orbit integraion method applied for each satellite separately, the computational efficiency of the proposed method is improved significantly:without damaging the accuracy, it takes only 0.09 s for a single satellite, which is 14 times faster than the traditional method. Besides, further improvement can be achieved when the number of satellites is increased.

    Key words: GNSS rapid orbit integration algorithm orbit fitting multi-system

    高精度轨道数值积分算法是导航卫星动力学精密定轨的基础。现阶段常用的高精度数值积分算法有单步法(如Runge-Kutta-Fehlberg算法,RKF)和多步法(如Admas预估校正算法)两种[1-2]。为了充分利用这两种积分方法的优势,实际应用中常常先通过RKF计算出Admas起始状态点,再进行Admas积分(RKF-Admas方法),有效提高单颗卫星长弧段轨道数值积分效率[3-5]随着我国北斗导航卫星区域系统(BDS)建成运行[6]以及其他GNSS(global navigation satellite system)系统(美国GPS、俄罗斯GLONASS、欧洲Galileo)的快速发展[7],多GNSS系统兼容应用是卫星导航领域的发展趋势,其中包括多系统GNSS卫星联合精密定轨[8-11]。在可预见的未来,空间中将有多达90余颗GNSS卫星投入使用[12],而传统RKF-Admas方法将会显著影响多系统GNSS卫星数值积分效率,具体表现在:一方面,由于各颗卫星积分过程较为独立,传统方法采用逐颗卫星依次积分,这种方式导致了同一时刻各颗卫星摄动公共计算部分的过多冗余计算,并会随着卫星数的增加显著降低多卫星积分效率,这在效率要求较高的实时滤波定轨中尤为明显[13-14];另一方面,在长弧段轨道积分中,为了保证所有卫星的积分精度,传统方法在实际应用中通常采用较为保守的固定积分步长(如90 s),从而限制了计算效率。这是由于数值积分方法易受动力学模型异常扰动的影响(如卫星地影机动[15-17]),大步长会严重损害积分精度,并且步长对GNSS卫星的影响有待进行研究。鉴于此,文献[18—20]研究了变步长Admas积分算法,通过计算不同阶公式实时估计每一步的误差,并自适应确定每步积分所需步长。这种方法需要增加不同阶的计算量,适用于一般轨道类型卫星的变步长数值积分。然而GNSS卫星轨道类型较单一,均采用近圆形轨道[21](偏心率e≤0.01),运行速度和摄动力影响都较稳定,无需频繁变更积分步长。

    综上所述,研究一种适用于多系统GNSS卫星的快速积分方法具有重要的意义。本文在RKF-Admas方法基础上,利用Admas预估校正算法特点,提出在大步长积分过程中通过预估值与校正值的差异探测动力学模型异常扰动,自适应降低Admas积分步长以保证积分精度,同时充分利用各颗卫星积分过程中的相同计算部分实现多卫星联合同步积分,从而达到快速积分目的。

    1 多系统GNSS卫星快速积分算法1.1 自适应变换Admas积分步长

    Admas预估校正算法使用被积函数值的线型多项式代替被积函数,在一步积分过程中只需计算两次动力学右函数,并且可以选取较大的步长,因此效率较单步法RKF高,常用于卫星长弧段轨道数值积分。该算法由两步构成,首先利用Adams-Bashforth显式公式通过历史存储的同阶数状态点对下一个状态点的值进行预报,然后利用Admas-Moulton隐式公式计算出下一个状态点的改进值[2, 22] (1)式中,表示包含卫星位置、速度及动力学参数的状态向量;x0表示在积分起始时刻的状态向量。

    n阶Adams-Bashforth显性公式由式(2) 给出

    (2)式中,xi+1p表示下一个状态点的预估值;h表示积分步长;ψnpi, i+1)表示由被积函数值组成的n阶线型多项式,其多项式系数αnj可通过查表获得。

    n阶Adams-Moulton隐性公式由式(3) 给出

    (3)式中,x表示下一个状态点的校正值;ψn(i, i+1)包含了由式(2) 给出的预估值计算得到的被积函数组成的n阶线性多项式,其系数βnj亦可通过查表获得。

    由式(2) 与式(3) 可以看出,Adams-Bashforth显式公式计算得到的预估值只受历史状态被积函数值的影响,而Adams-Moulton隐式公式则受下一状态被积函数的影响。当被积函数受动力学模型异常扰动影响时(如进出地影或轨道机动),将会通过校正值与预估值之间差异反映出来。

    n阶Adams-Moulton隐式公式截断误差可表示为

    (4)γ表示后向差分系数[2]。由此可见,当被积函数受动力学模型异常扰动影响时,被积函数的n+1阶导数fi(n+1)将变得不连续,若步长h选取较大时,则会放大被积函数的不连续导致的误差,破坏了数值稳定性,导致长弧段积分误差积累不可忍受。

    因此,可以通过构造预估值与校正值的差异探测动力学模型异常扰动,自适应降低Admas积分步长以保证积分精度。判断条件式如下

    (5)式中,表示前i个差异值的平均值;σi表示前i个差异值的标准偏差。当被积函数光滑时,Δxi+1会接近于零且变化稳定,而被积函数连续性被破坏时,Δxi+1将会呈数量级增大,因此k值需要根据卫星的受力变化或Δx变化特征分析进行合理选取。当Δxi+1满足式(5) 时,可认为动力学模型有异常发生。为保证积分精度,此时终止原先的积分过程,根据Admas阶数选取先前某个未发生异常的状态点,利用RKF积分得出较小步长Admas起始点后再继续积分过程,由此可以完成自适应Admas变换步长积分算法。值得说明的是,由于被积函数fi与积分步长无关,因此积分步长的调整不会影响积分轨道的连续性。以11阶Admas算法为例,其变换步长过程如图 1所示。点击查看原图图 1 自适应变换Admas积分步长方法示意Fig. 1 Diagram of the adaptive step-changed Admas integration method图选项
    1.2 多系统GNSS卫星快速积分算法

    多系统GNSS卫星精密定轨对联合数据处理效率要求较高,其中就包括需要在较短时间内完成几十颗卫星的数值积分工作。由于各颗卫星积分过程较为独立,现阶段通常采用逐卫星依次积分,而同一时刻惯性系与地固系转换矩阵及其偏导数、行星星历的计算、固体潮海潮摄动和重力场摄动系数等对所有卫星影响相同,因此逐颗卫星依次积分存在过多的冗余计算。由于导航卫星轨道类型较单一,积分策略基本一致,这为多卫星同步积分的实现提供了便利。鉴于此,结合自适应变换Admas步长和多卫星联合同步积分两方面,提出一种多系统GNSS卫星快速积分算法:以RKF-Admas算法为基础,充分利用大步长及多卫星联合同步积分的效率优势,通过自适应变换Admas步长控制动力学模型扰动对卫星积分精度的损害,以达到快速精确完成多卫星积分目的。多系统GNSS卫星快速积分算法流程如图 2所示。

    点击查看原图图 2 多系统GNSS卫星快速积分算法流程Fig. 2 The flow chart of the rapid orbit integration algorithm for multi-GNSS satellites图选项
    2 试验与分析2.1 积分步长对轨道动力学拟合精度影响分析

    为了提高积分效率,需要选取较大的积分步长。然而大步长会影响积分精度,因此需要分析不同积分步长(尤其是大步长)对导航卫星积分精度的影响。由于积分产生的误差会最终反映至轨道动力学拟合RMS中,因此轨道拟合可以用来评估不同步长对积分精度的影响。

    利用武汉大学(Wuhan University, WHU)发布的GNSS卫星精密轨道产品[8],通过设置60 s、90 s、300 s和450 s 4种积分步长对GPS、GLONASS、BDS和Galileo卫星进行RKF-Admas(分别选取RKF8(9) 阶和11阶Admas)固定步长的单天数值积分与拟合计算(2014年,年积日275),其三维RMS由图 3给出。其中动力学模型选用IGS推荐的策略(http://acc.igs.org/reprocess2.html)。图中G03与G25两颗卫星因在精密星历中缺失未参与统计。从图中可以看出,有7颗GPS卫星(G05、G10、G12、G18、G20、G22和G32)、5颗BDS的GEO卫星(C01至C05) 以及1颗BDS的MEO卫星(C14) 对不同积分步长较敏感,其轨道拟合精度在不同步长间差异较大,尤其是当步长增加到450 s时,轨道拟合精度最大衰减至分米级别;而当步长在90 s以内时,积分仍能反映出较好的稳定性,其差异在2 mm以内。除去这13颗卫星,其余卫星轨道拟合精度对这4种积分步长不敏感,这说明即使将积分步长增加到450 s,大部分卫星的积分精度也并无损害(差异小于0.5 mm)。点击查看原图图 3 不同积分步长GPS、GLONASS、BDS和Galileo卫星轨道拟合三维RMSFig. 3 The 3D-RMS of 1-day orbit fitting applied for GPS, GLONASS, BDS and Galileo satellites based on different intergration step图选项图 4给出了450 s积分步长时BDS的C08卫星与C14卫星在Admas积分过程中三维坐标位置预估值与校正值差异的单天时间序列(总共192个积分点)。由于Admas阶数为11,其起始的11个积分点由RKF计算获得,因此差异数值为零。从图中可以看出,C08卫星预估值与校正值差异的时间序列稳定(STD约为1×10-8 m),并且其值接近于零;而C14卫星的时间序列产生了两次明显的扰动,并且扰动幅度呈数量级增大,远远超过平稳数据点的变化范围。点击查看原图图 4 Admas积分(450 s步长)卫星三维坐标位置预估值与校正值差异Fig. 4 The difference of Admas (450 s step) prediction and correction for the satellite's position图选项为了进一步探究扰动的原因,利用双锥体地影模型[15]统计当天各颗蚀卫星通过地影的起止时刻,其结果列于表 1中,单位为秒。从表 1中可以看出,C08卫星未受地影影响,被积函数光滑连续,因此Admas积分预报值与估计值差异时间序列较稳定且数值接近于零;C14卫星两次通过地影区间,并且进入地影时刻与图 4中的两次波动的起始时刻吻合,这是由于卫星进入地影时破坏了被积函数的连续性,采用大积分步长时在半影期间未被积分算法所采样,因此导致产生了明显的波动[2],从而影响积分精度。另一方面,所有检测到通过地影的卫星与图 3中不同积分步长结果差异较大的卫星一致,并且通过试验得知这些卫星具有与C14卫星相似的变化规律(GEO卫星略有不同,其时间序列只有一次波动,这与进出地影次数相符,由于篇幅有限,结果不再给出),而未通过地影的卫星具有与C08卫星相似的稳定性,因此进出地影为导致GNSS卫星被积函数不连续的主要原因。本文基于以上GNSS卫星Admas积分预报值与估计值差异变化特征,假定差异数值上两个数量级的跳跃即有可能发生了动力学异常扰动,即k取值100,则地影机动通过式(5) 同时结合双锥体地影模型可以准确检测出来。

    表 1 蚀卫星通过地影的起止时刻(2014年,年积日275)Tab. 1 The duration (second of day) of satellite eclipse on DOY=275, 2014

    s卫星起始时刻结束时刻G05660342043 62046 560G127320798050 28051 240G2028 62031 86071 70074 940G3224 90027 00067 92070 140C0357 48061 200C0569 78073 680G1038 16041 04081 24084 180G1812 12015 06055 20058 200G2216 02019 02059 10062 100C0150 28054 240C0264 62068 580C0445 54049 320C1426 94027 96073 14074 520

    表选项

    2.2 自适应变换Admas积分步长算法效果分析

    由于大步长会导致蚀卫星积分精度的极大衰减,因此可以利用本文提出的自适应变换Admas积分步长算法来控制这种精度的衰减。为了验证自适应变换Admas积分步长算法效果,分别用固定步长Admas积分算法和自适应变换步长Admas积分算法对年积日275当天蚀卫星进行轨道拟合,其三维RMS由图 5给出。其中自适应算法中小步长和大步长分别设置为90 s和450 s。从图中可以看出,自适应变换步长算法可以有效避免大步长对蚀卫星积分精度的损害,其轨道拟合精度与固定为小步长的Admas积分算法精度相同,证明了本文提出的自适应变换Admas积分步长算法的有效性。

    点击查看原图图 5 不同积分策略的蚀卫星轨道拟合三维RMSFig. 5 The 3D-RMS of 1-day orbit fitting applied for satellites across eclipse based on different integration strategies图选项
    2.3 多系统GNSS卫星快速积分算法精度及效率分析为了充分验证多系统GNSS卫星快速积分方法的精度,分别利用IGS多模GNSS试验跟踪网(multi-GNSS experiment, MGEX)[7]的两个分析中心武汉大学(WHU)[8]和欧洲定轨中心(Center for Orbit Determination in Europe, CODE)[23]发布的GNSS精密轨道产品进行一个月的单天轨道拟合试验(2014年,年积日245至275),统计平均三维RMS。结合上节的分析,小步长选取90 s,大步长选取450 s,通过自适应变换Admas步长算法保证蚀卫星的积分精度。

    图 6给出了一个月内四系统GNSS卫星单天轨道拟合平均三维RMS。从图中可以看出,算法具有较高的精度,其中WHU拟合三维RMS的平均值:GPS为11 mm,GLONASS为14 mm,BDS为16 mm,Galileo为12 mm;CODE拟合三维RMS的平均值:GPS为17 mm,GLONASS为18 mm,BDS为18 mm,Galileo为14 mm。由于本文采用的动力学模型与WHU产品一致,因此拟合精度略优于CODE产品。两个分析中心产品的平均三维拟合精度均优于20 mm,说明本文提出的快速积分算法能够满足精密定轨的需求。

    点击查看原图图 6 4系统GNSS卫星单天轨道拟合平均三维RMSFig. 6 The average 3D-RMS of 1-day orbit fitting for GNSS satellites图选项

    表 2给出了本文快速积分算法与传统逐卫星固定步长算法在一个月内单天轨道拟合RMS差异的相关统计数据。为了保证包括蚀卫星在内的所有卫星的积分精度,传统积分方法固定积分步长选取为90 s。从表中可以看出,采用两种方法计算得到的差异均值均小于0.5 mm,最大值小于2 mm,说明本文提出的算法与传统方法具有相同的积分精度。由于快速积分算法中非蚀卫星实质上采用的是450 s大步长,因此这种积分差异主要是由两种方法积分步长不同导致的。此外,蚀卫星拟合RMS差异较非蚀卫星稍大,这是由于蚀卫星动力学模型的不连续性对积分步长变化的影响稍大,不过其值本身依然与非蚀卫星差异较小(0.1 mm左右),说明新算法能够有效处理大步长对蚀卫星的影响。

    表 2 本文快速积分算法与传统逐卫星固定步长积分算法单天轨道拟合RMS差异统计Tab. 2 The statistics of 1-day orbit fitting RMS difference between the traditional step-fixed orbit integraion method applied for each satellite separately and the rapid integration algorithm developed in this paper

    mm
    WHUCODE非蚀卫星平均值0.40.3最大值1.31.2蚀卫星平均值0.50.4最大值1.81.4

    表选项

    为了进一步验证多系统GNSS卫星快速积分方法对计算效率的提升作用,分别利用传统逐卫星固定步长积分、多卫星联合固定步长积分和多卫星联合自适应变换步长积分(即本文提出的快速积分方法)3种积分策略进行试验,统计平均单颗卫星的耗时以及总耗时随卫星数的变化(CPU主频2.60 GHz,内存容量2 GB)。图 7给出了以上3种积分策略下平均单颗卫星耗时。从图中可以看出,同时采用固定步长时,多卫星联合积分算法较逐卫星积分算法效率有明显的提高,其平均单颗卫星耗时均值分别为0.21 s和1.38 s,提高了约6倍,这是由于不同卫星共用计算部分在同一时刻只计算一次,大幅度提高了计算效率;同时采用多卫星联合积分时,自适应变换步长较固定步长算法效率也有所提高,其平均单颗卫星耗时均值为0.09 s,提高了约1倍,这是由于大部分未受地影影响的卫星有效利用了大步长的效率优势,减少了动力学模型计算次数。与传统逐卫星固定步长积分算法相比,自适应变换步长多卫星联合快速积分算法效率提升了14倍,说明本文提出的新算法具有较高的效率。

    点击查看原图图 7 不同积分策略平均单颗卫星耗时Fig. 7 The average computing time of a single satellite based on different integration strategies图选项

    图 8给出了3种积分策略下计算总耗时与卫星数的关系。从图中可以看出随着卫星数的增加,3种积分策略的计算总耗时也增加,并且呈线性相关。然而,在线性变化趋势上不同积分策略具有不同的斜率,其中采用自适应变换步长的多卫星联合积分算法的斜率明显小于采用传统逐卫星固定步长积分算法的斜率,这说明参与积分计算卫星数目越多,效率提升越明显。可以预计在多达100颗卫星同时参与计算情况下,新算法仍能在10 s内完成单天轨道积分与拟合。

    点击查看原图图 8 不同积分策略总耗时与卫星数的关系Fig. 8 The relationship between the computing time and the number of satellites based on different integration strategies图选项
    3 结论

    随着多系统GNSS在轨卫星数的不断增加,传统逐卫星固定步长积分算法成为制约多系统GNSS卫星数值积分效率的重要因素。本文结合自适应变换Admas积分步长和多卫星联合同步积分两个方面实现了多系统GNSS卫星快速积分方法,并研究了不同积分步长对GNSS卫星积分精度的影响。试验结果表明:① 对于未受地影影响的卫星,积分步长增大到450 s时也不损失轨道拟合精度,而对于受地影影响的卫星,积分步长在90 s以内时轨道拟合精度差异小于2 mm;② 自适应变换Admas积分步长可以在充分利用大步长效率优势的同时,避免了大步长对蚀卫星积分精度的损害;③ 利用WHU和CODE发布的精密轨道产品,GPS/GLONASS/BDS/Galileo 4个系统轨道拟合平均三维RMS均优于20 mm,说明多系统GNSS卫星快速积分算法满足精密定轨的需求;④ 相比传统逐卫星固定步长积分算法,新算法在不损失传统算法积分精度前提下能显著提高计算效率,可用于精密定轨中变分方程的计算,其平均单颗卫星计算速度提升14倍,并且随着卫星数的增加,提升效果越明显。

    致谢:感谢武汉大学和欧洲定轨中心提供的GPS、GLONASS、BDS和Galileo卫星精密轨道产品。

    【引文格式】范磊,李敏,宋伟伟,等。 多系统GNSS卫星轨道快速积分方法[J]. 测绘学报,2016,45(S2):93-100. DOI: 10.11947/j.AGCS.2016.F030

    权威 | 专业 | 学术 | 前沿

    微信投稿邮箱 | song_qi_fan@163.com

    欢迎加入《测绘学报》作者QQ群:297834524

    进群请备注:姓名+单位+稿件编号

    位置:首页 > 探索
    加载更多评论...
    本类推荐
    分分钟消灭翻译官?揭秘搜狗语音实时翻译技术
    分分钟消灭翻译官?揭秘搜狗语音实时翻译技术

    英语没学好,怎么与老外交流?有钱人可以请翻译官,那屌丝怎么办?在不久前结束的乌镇互联网大会上,搜狗 CEO 王小川曾演示了搜狗最新的语音实时翻译技术,当时王小川的中文讲话内容被实时识别为文字并翻译成英文显示在了屏幕上,成功将王小川关于人工智能专业领域的报告进行机器同传。

    4H-SiCMESFET特性对比及仿真
    4H-SiCMESFET特性对比及仿真

    侯 斌,邢 鼎,张战国,臧继超,马 磊(航天科技集团九院七七一研究所,陕西 西安710000)通过对双凹栅结构和阶梯栅结构4H-SiC MESFET的直流特性对比,得出阶梯栅结构的直流特性优于双凹栅结构。

    LOL模型重做后非常好看的十款皮肤,琴女简直脱胎换骨!
    LOL模型重做后非常好看的十款皮肤,琴女简直脱胎换骨!

    LOL有些皮肤出一直被吐槽着,而且为了联盟新血液,偶尔也要模型重做下。下面我们就来盘点下模型重做后非常好看的十款皮肤:1、特种部队 普朗克这皮肤重做后整个人物比例接近完美,最大的受益就是动作流畅上的提升,还有技能特效在视觉上的提升。

    人体是如何调控血糖的?数学模型告诉你详细过程
    人体是如何调控血糖的?数学模型告诉你详细过程

    随着我们生活水平的提高,人们已经渐渐地不像过去那样担心温饱问题了,每天山珍海味的生活也不再稀罕——午餐来一顿正宗的羊杂火锅,晚餐煎两块半熟的神户牛排,再和三五好友撸几串烧烤作为夜宵,拍个照发发朋友圈,大快朵颐又格调十足的一天就这么过去了。

    能根据文字生成图片的GAN,深度学习领域的又一新星
    能根据文字生成图片的GAN,深度学习领域的又一新星

    2014 年 6 月,Ian Goodfellow等学者发表了论文《Generative Adversarial Nets》,题目即“生成对抗网络”,标志着 GAN 的诞生。2015 年还名不见经传的 GAN ...

    论文推荐|董杨:主成分分析的匹配点对提纯方法
    论文推荐|董杨:主成分分析的匹配点对提纯方法

    《测绘学报》构建与学术的桥梁 拉近与权威的距离测绘地理信息与导航高端论坛 ——《测绘学报》创刊60周年学术研讨会通知(第一号)董杨, 范大昭, 纪松, 雷蓉 信息工程大学, 河南 郑州 450000收稿日期:2016-05-24;

    教学使用的人骨模型居然是真人骨!?
    教学使用的人骨模型居然是真人骨!?

    学校里使用的人骨模型一般是用PVC(聚氯乙烯)材料制造的。没什么可怕,对不对!!然而,要是真的人骨呢?英国Merseyside 一高中(Haydock High School)有一具使用了40年的“人骨模型”,因传说这是真人骨。

    问答Goodfellow:没有样例能不能训练机器学习算法?
    问答Goodfellow:没有样例能不能训练机器学习算法?

    李杉 编译整理量子位 出品 | 公众号 QbitAI从图像和语音识别到自然语言分析,神经网络已经在很多领域大展身手。过去几年,它们的精确度已经几乎可以与人类媲美。但仍有很多神经网络无法完成的任务——例如,这种技术还无法取代人类的创造力。

    捐款公示20万只给2000元,你的无耻能值19万8000元?
    捐款公示20万只给2000元,你的无耻能值19万8000元?

    知道“无耻”两个字咋写吗?文/吴元中据央广网报道,河北邯郸涉县后池村民义务上山修路,靠双手改变贫穷落后面貌被誉为太行“新愚公”。近日在当地新建成的希望小学爱心拍卖会上,有单位将写有20万元的现金支票模型送到校领导手中,但事后拿到手的却是2000元。

    连锁便利店与传统小卖部差距在哪里
    连锁便利店与传统小卖部差距在哪里

    通常品牌连锁性质的便利店会为了满足消费者贪新厌旧的消费心理,不断的引进和创造新的消费产品,来吸引消费者的注意,形成独特的产品消费模式。而对于许多便利店经营者来说,都只会停留在缺货补货,收钱找钱的阶段上...

    中国古代刑具模型图,跟真的一样,看着都吓人
    中国古代刑具模型图,跟真的一样,看着都吓人

    站笼:又称为立枷,笼上有一个口,卡往囚犯的颈部,让犯人日夜站立,直到死去拶刑:用拶子套入手指,再用力紧收,古代对女犯施用的一种酷刑,大家在影视剧中,县衙老爷审问女犯时应该见过削足:这个不说了,非常残忍...

    古代刑具模型图,跟真的一样,看着都吓人
    古代刑具模型图,跟真的一样,看着都吓人

    站笼:又称为立枷,笼上有一个口,卡往囚犯的颈部,让犯人日夜站立,直到死去拶刑:用拶子套入手指,再用力紧收,古代对女犯施用的一种酷刑,大家在影视剧中,县衙老爷审问女犯时应该见过削足:这个不说了,非常残忍...

    延伸热词
    首页社会国际娱乐科技时尚军事汽车探索美食旅游历史健康育儿
    万花镜 版权所有 京ICP备14059027号
    值班QQ:3012642954
    邮箱:wanhuajingnews@qq.com