一种用于岩石损伤特征辨析的声发射信号处理方法
一种用于岩石损伤特征辨析的声发射信号处理方法
【技术领域】
[0001] 本发明设及一种损伤辨析的声发射信号处理方法，特别是一种用于岩石损伤特征 辨析的声发射信号处理方法。
【背景技术】
[0002] 目前，一般的岩石力学室内试验系统在实验过程中只能测试到力-位移曲线，运些 数据只能提供很有限的特征性质，也不能检测到在整个变形过程中的不均匀的裂隙的形成 与演化。声发射信号在整个压裂过程中大量产生，它携带了在加载过程中有关裂纹扩展和 能量释放的信息，在检测岩石加载破坏运样一个动态的过程中是一种很有效的方法。声发 射技术作为一种有效的无损评价的工具，声发射技术已被广泛用来定量研究岩石破裂过 程，分析运些信号可W捕获内部结构的演变和检测变形阶段。同时，声发射信号包括大量的 时频信息，如果采用适当的数学方法来处理运些声发射信号，会有助于对岩石破裂过程的 理解，信号是数字信号，可W通过信号处理方法研究。部分研究人员研究了声发射信号的时 频特性，他们通过对声发射信号进行小波变换W及傅里叶变换来获取声发射信号的时频特 性，还研究了微振信号的光谱特性，上述所有的研究均不能解决光谱特征和岩石变形行为 之间的关系。
[0003] 从技术上讲，声发射信号是非线性和非稳定的，不遵循传统的光谱分析的限制。例 如，用傅里叶变换的信号必须是线性的和严格的周期性或相对静止的，否则，形成的光谱没 有什么物理意义。同样，基于小波方法的数据序列的谱特征分析取决于小波函数的选择，所 W运些结果高度依赖于分析者的经验去选取合适的函数。
【发明内容】
[0004] 本发明的目的是要提供一种用于岩石损伤特征辨析的声发射信号处理方法，解决 基于小波变换谱特征分析中小波函数选取不准的问题，使得加载下岩石的损伤与衰减可W 在宏观和微观下进行描述。
[000引本发明的目的是运样实现的：声发射信号处理方法，包括如下步骤：
[0006] 步骤一:在岩石力学特性的测试系统上加上声发射检测系统，采集岩石加载过程 中声发射信号的变化；
[0007] 步骤二:将声发射系统采集到的岩石试件破裂过程中声发射信号导出；
[0008] 步骤Ξ:对导出的声发射信号进行提取，提取其中的有效能量作为声发射参数信 号，W及提取声发射信号中比较振幅比较大的信号作为分析对象；
[0009] 步骤四：将步骤Ξ提取的信号分别进行希尔伯特黄变换包括两个部分:变换一为 经验模态分解或集合经验模分解变换;变换二为希尔伯特谱分析；
[0010] 步骤五:对变换结果进行分析，进而通过声发射信号从宏观和微观两个层面上来 说明岩石破坏过程中的损伤与衰减。
[0011] 有益效果，由于采用了上述方案，利用本发明的声发射信号处理办法，可W将非线 性非稳定的信号变换成不同的固有模态函数，然后获得它们的希尔伯特谱，运样就使得岩 石在加载下的损伤与衰减能够从宏观与微观两个层面上进行描述，达到了本发明的目的。
[0012]其优点在于：
[001引1.解决了基于小波变换谱特征分析中小波函数选取不准的问题；
[0014] 2.基于声发射信号参数的希尔伯特黄变换提供了一个更好的区分岩石破坏过程 不同阶段的方法，运是相对于那些只利用声发射信号分析方法的优势所在。因此，它可W用 来预测岩石试件的破坏。
【附图说明】
[0015] 图1为本发明岩石损伤特征辨析的声发射处理系统的实验流程图。
[0016] 图2为岩石试件在溫度-应力禪合作用下的应力应变曲线W及曲线上各点对应的 试件形态图。
[0017] 图3为本发明岩石加载在线性变形阶段声发射信号有效能量与它的希尔伯特谱 图。
[0018] 图4为本发明岩石加载在过度阶段的声发射信号有效能量与它的希尔伯特谱图。
[0019] 图5为本发明岩石加载过程在过度阶段声发射信号的希尔伯特谱特征图。
[0020] 图6为本发明岩石加载过程的软化阶段某点的典型的声发射信号与它的希尔伯特 谱图。
[0021 ]图7为本发明在岩石非线性压密阶段典型的声发射信号与它的希尔伯特谱图。
【具体实施方式】
[0022] 下面结合附图和具体的实例对本发明提供的一种岩石损伤特征辨析的声发射信 号处理方法进行详细的说明。
[0023] 实施例1:本发明提供的一种岩石损伤特征辨析的声发射信号处理方法，
[0024] 图1为实验流程图，声发射信号处理具体的实施方法如下：
[0025] 步骤一:在岩石力学特性的测试系统上加上声发射检测系统，采集岩石加载过程 中声发射信号的变化；
[0026] 步骤二:将声发射系统采集到的岩石试件破裂过程中声发射信号导出；
[0027] 步骤对导出的声发射信号进行提取，提取其中的有效能量作为声发射参数信 号，W及提取声发射信号中比较振幅比较大的信号作为分析对象；
[0028] 步骤四：将步骤Ξ提取的信号分别进行希尔伯特黄变换包括两个部分:变换一为 经验模态分解或集合经验模分解变换;变换二为希尔伯特谱分析；
[0029] 步骤五:对变换结果进行分析，进而通过声发射信号从宏观和微观两个层面上来 说明岩石破坏过程中的损伤与衰减。
[0030] 经验模态分解和集合经验模分解变换的具体实施方法:利用经验模态分解和集合 经验模态分解处理工具，可W将信号分解成不同的固有模态函数：
[0031 ]经验模态分解数据预处理遵循下面的步骤：
[0032] 1.通过差值法计算输入信号x〇(t)在每一点上的局部最大值emax(t)和局部最小值 emin(t)，每一点上信封运些信号，运用Ξ次样条方法形成上包络线和下包络线；
[003引 2.计算局部最大值emax(t)和局部最小值emin(t)的平均值m(t)，然后计算xo(t)和m (t)之间的差异，得到h(t);
[0034] 3.检查h(t)是否是固有模态函数，如果不是，将h(t)设为输入信号继续进行上述 的步骤，直到第一个固有模态函数的出现；
[0035] 4.如果一旦发现第一个固有模态函数，则从原始信号减去它;运时xi(t)=xo(t)-Cl(t)，把Xl(t)作为输入信号重复上述过程获得固有模态函数，重复上述过程，得到第i个 固 有模态函数；
[0036] 5.继续上面的步骤，直到不能分离出固有模态函数，剩余的部分被称为一个单调 函数;然后，得到了 η个ci(t)，i = l-n和一个单调函数r(t);
[0037 ]原始信号经过变换W后就可W写作：
[0038] 集合经验模态分解的具体实施方法：
[0039] 岩石破裂发生时声发射信号的振幅会突然增加;然而，当面临破裂时会出现模态 混杂的问题，它会严重混淆时频分析和隐藏固有模态函数的物理意义;所W提出了集合经 验模态分解：
[0040] 1、集合经验模态分解，首先添加白噪音序列到原始信号;第j个原始信号和白噪声 序列合并用xW(t)表示；
[0041] 2、对于每个混合信号，
可W通过使用经验模态分解方法获得;注意， 意味着第j个组合序列的第i个固有模态函数；
[0042] 3、原始信号的第i个固有模态函数ci(t)是所有之前所获的的第i个固有模态函数 的总平均值。
[0043] 希尔伯特变换用于计算瞬时频率和振幅，其具体的实施方法为：
[0044] (1)
[0045]
[0046] 通过运种变换，任何时间序列代表一个分析信号：
[0050] (2)本方法选择有效能量的声发射率作为声发射参数信号。声发射有效能量Eae被 定义为：
[0051 ]
代表的是在t时刻的电压值；
[0052] (3)岩石破裂时会释放能量，声发射信号的产生伴随着弹性能的释放，声发射有效 能量和断裂能量可能存在一定得关系;声发射有效能量分析的目的是测量岩石的断裂瞬间 释放的能量;岩石破裂事件造成的声发射信号持续时间只有10毫秒，而且信号是瞬间的，很 快就会降低到阀值W下，我们假设每一个时间点只有一个岩石的破裂事件产生，所W每个 时刻的声发射信号都是独立的。因此，第i次岩石破裂产生的声发射有效能量为：
[0053]
[0054] (4)如果在一个单位时间内发生了 η次的破坏事件，那么声发射有效能量率馬f为：
[0055]
[0056] (5)经过希尔伯特黄变换，声发射有效能量率可W表示为：
[0057]
[005引在上述方程中：
[0059] xi(t)第i个输入的声发射信号；
[0060] x〇(t) 一个典型的声发射输入信号；
[0061] emax(t)，emin(t)，m(t)输入信号的局部最大值与局部最小值，W及它们的平均值；
[0062] h(t)输入信号与平均值之间的差值；
[0063] Ci(t)第i个固有模态函数；
[0064] r(t) -个单调函数；
[0065] xW(t)第j个典型声发射信号与白噪音序列的混合信号；
[0066] 的第j个混合信号的第i个固有模态函数；
[0067] y(t)x(t)的希尔伯特变换；
[0068] ω (t)，a(t)瞬时频率与瞬时振幅；
[0069] Eae声发射信号有效能量；
[0070] 巧》C第i次岩石破裂事件发生时的声发射有效能量；
[0071] 馬声发射有效能率；
[0072] Wt)与声发射信号强度有关的电压；
[0073] Δ ti每一个岩石破裂事件持续的时间；
[0074] 〇如)第j个固有模态函数的瞬时频率；
[007引 w(t)，a/ (t)第j个固有模态函数的瞬时振幅和它的对数。
[0076]岩石材料的室内实验是一个开放系统热力学过程，能量从岩石实验系统到岩石试 件连续不断的转化成Ξ种能量:弹性能量、塑料能和消失的能量。弹性变形产生弹性能量可 W在卸载时完全释放，相反，塑性变形产生的能量无法释放。此外，能量的消失与福射现象， 相变和化学能源、压裂能等有关。为了简化分析方法，我们假设消失的能量等于从弹性能中 生成的破裂能。除此之外，先前的研究表明，声发射有效能量率与断裂能量释放率在整个故 障过程中是线性相关，无论材料是金属或非金属，宏观或微观尺度，同时声发射有效能量率 可W反映在上述方程瞬时振幅中。因此，基于声发射信号参数的希尔伯特黄变换提供了一 个更好的区分岩石破坏过程不同阶段的方法，运是相对于那些只利用声发射信号分析方法 的优势所在。由此可得，该方法可W用来预测岩石试件的破坏。
【主权项】
1. 一种基于用于岩石损伤特征辨析的声发射信号处理方法，其特征是：声发射信号处 理方法，包括如下步骤： 步骤一:在岩石力学特性的测试系统上加上声发射检测系统，采集岩石加载过程中声 发射信号的变化； 步骤二:将声发射系统采集到的岩石试件破裂过程中声发射信号导出； 步骤三:对导出的声发射信号进行提取，提取其中的有效能量作为声发射参数信号，以 及提取声发射信号中振幅比较大的信号作为分析对象； 步骤四：将步骤三提取的信号分别进行希尔伯特黄变换包括两个部分:变换一为经验 模态分解或集合经验模分解变换;变换二为希尔伯特谱分析； 步骤五:对变换结果进行分析，进而通过声发射信号从宏观和微观两个层面上来说明 岩石破坏过程中的损伤与衰减。
【专利摘要】一种用于岩石损伤特征辨析的声发射信号处理方法，属于损伤辨析的声发射信号处理方法。该方法通过对岩石力学特性测试系统中采集到的声发射信号进行变换，对声发射信号进行提取，选取声发射信号中的有效能量参数，通过电压对时间的积分得到有关于能量的声发射参数信号，然后分别对比较强的声发射信号和声发射参数信号进行希尔伯特黄变换，通过变换提取信号中的瞬时频率与瞬时振幅，瞬时振幅代表从岩石试件在实验过程中释放的能量，瞬时频率与岩石试件在加载过程中的损伤的衰减相关。将原始的声发射信号通过希尔伯特黄变换分解成不同的固有模式函数，然后获得希尔伯特光谱，这种变换方法使得加载下岩石的损伤描述更加清楚，更加直观。
【IPC分类】G01N29/44
【公开号】CN105486756
【申请号】CN201510863621
【发明人】彭维红, 刘峰宇, 张骥, 张海翔, 杜梦琳, 曹国华, 程红梅, 张志镇
【申请人】中国矿业大学
【公开日】2016年4月13日
【申请日】2015年12月1日
【技术领域】
[0001] 本发明设及一种损伤辨析的声发射信号处理方法，特别是一种用于岩石损伤特征 辨析的声发射信号处理方法。
【背景技术】
[0002] 目前，一般的岩石力学室内试验系统在实验过程中只能测试到力-位移曲线，运些 数据只能提供很有限的特征性质，也不能检测到在整个变形过程中的不均匀的裂隙的形成 与演化。声发射信号在整个压裂过程中大量产生，它携带了在加载过程中有关裂纹扩展和 能量释放的信息，在检测岩石加载破坏运样一个动态的过程中是一种很有效的方法。声发 射技术作为一种有效的无损评价的工具，声发射技术已被广泛用来定量研究岩石破裂过 程，分析运些信号可W捕获内部结构的演变和检测变形阶段。同时，声发射信号包括大量的 时频信息，如果采用适当的数学方法来处理运些声发射信号，会有助于对岩石破裂过程的 理解，信号是数字信号，可W通过信号处理方法研究。部分研究人员研究了声发射信号的时 频特性，他们通过对声发射信号进行小波变换W及傅里叶变换来获取声发射信号的时频特 性，还研究了微振信号的光谱特性，上述所有的研究均不能解决光谱特征和岩石变形行为 之间的关系。
[0003] 从技术上讲，声发射信号是非线性和非稳定的，不遵循传统的光谱分析的限制。例 如，用傅里叶变换的信号必须是线性的和严格的周期性或相对静止的，否则，形成的光谱没 有什么物理意义。同样，基于小波方法的数据序列的谱特征分析取决于小波函数的选择，所 W运些结果高度依赖于分析者的经验去选取合适的函数。
【发明内容】
[0004] 本发明的目的是要提供一种用于岩石损伤特征辨析的声发射信号处理方法，解决 基于小波变换谱特征分析中小波函数选取不准的问题，使得加载下岩石的损伤与衰减可W 在宏观和微观下进行描述。
[000引本发明的目的是运样实现的：声发射信号处理方法，包括如下步骤：
[0006] 步骤一:在岩石力学特性的测试系统上加上声发射检测系统，采集岩石加载过程 中声发射信号的变化；
[0007] 步骤二:将声发射系统采集到的岩石试件破裂过程中声发射信号导出；
[0008] 步骤Ξ:对导出的声发射信号进行提取，提取其中的有效能量作为声发射参数信 号，W及提取声发射信号中比较振幅比较大的信号作为分析对象；
[0009] 步骤四：将步骤Ξ提取的信号分别进行希尔伯特黄变换包括两个部分:变换一为 经验模态分解或集合经验模分解变换;变换二为希尔伯特谱分析；
[0010] 步骤五:对变换结果进行分析，进而通过声发射信号从宏观和微观两个层面上来 说明岩石破坏过程中的损伤与衰减。
[0011] 有益效果，由于采用了上述方案，利用本发明的声发射信号处理办法，可W将非线 性非稳定的信号变换成不同的固有模态函数，然后获得它们的希尔伯特谱，运样就使得岩 石在加载下的损伤与衰减能够从宏观与微观两个层面上进行描述，达到了本发明的目的。
[0012]其优点在于：
[001引1.解决了基于小波变换谱特征分析中小波函数选取不准的问题；
[0014] 2.基于声发射信号参数的希尔伯特黄变换提供了一个更好的区分岩石破坏过程 不同阶段的方法，运是相对于那些只利用声发射信号分析方法的优势所在。因此，它可W用 来预测岩石试件的破坏。
【附图说明】
[0015] 图1为本发明岩石损伤特征辨析的声发射处理系统的实验流程图。
[0016] 图2为岩石试件在溫度-应力禪合作用下的应力应变曲线W及曲线上各点对应的 试件形态图。
[0017] 图3为本发明岩石加载在线性变形阶段声发射信号有效能量与它的希尔伯特谱 图。
[0018] 图4为本发明岩石加载在过度阶段的声发射信号有效能量与它的希尔伯特谱图。
[0019] 图5为本发明岩石加载过程在过度阶段声发射信号的希尔伯特谱特征图。
[0020] 图6为本发明岩石加载过程的软化阶段某点的典型的声发射信号与它的希尔伯特 谱图。
[0021 ]图7为本发明在岩石非线性压密阶段典型的声发射信号与它的希尔伯特谱图。
【具体实施方式】
[0022] 下面结合附图和具体的实例对本发明提供的一种岩石损伤特征辨析的声发射信 号处理方法进行详细的说明。
[0023] 实施例1:本发明提供的一种岩石损伤特征辨析的声发射信号处理方法，
[0024] 图1为实验流程图，声发射信号处理具体的实施方法如下：
[0025] 步骤一:在岩石力学特性的测试系统上加上声发射检测系统，采集岩石加载过程 中声发射信号的变化；
[0026] 步骤二:将声发射系统采集到的岩石试件破裂过程中声发射信号导出；
[0027] 步骤对导出的声发射信号进行提取，提取其中的有效能量作为声发射参数信 号，W及提取声发射信号中比较振幅比较大的信号作为分析对象；
[0028] 步骤四：将步骤Ξ提取的信号分别进行希尔伯特黄变换包括两个部分:变换一为 经验模态分解或集合经验模分解变换;变换二为希尔伯特谱分析；
[0029] 步骤五:对变换结果进行分析，进而通过声发射信号从宏观和微观两个层面上来 说明岩石破坏过程中的损伤与衰减。
[0030] 经验模态分解和集合经验模分解变换的具体实施方法:利用经验模态分解和集合 经验模态分解处理工具，可W将信号分解成不同的固有模态函数：
[0031 ]经验模态分解数据预处理遵循下面的步骤：
[0032] 1.通过差值法计算输入信号x〇(t)在每一点上的局部最大值emax(t)和局部最小值 emin(t)，每一点上信封运些信号，运用Ξ次样条方法形成上包络线和下包络线；
[003引 2.计算局部最大值emax(t)和局部最小值emin(t)的平均值m(t)，然后计算xo(t)和m (t)之间的差异，得到h(t);
[0034] 3.检查h(t)是否是固有模态函数，如果不是，将h(t)设为输入信号继续进行上述 的步骤，直到第一个固有模态函数的出现；
[0035] 4.如果一旦发现第一个固有模态函数，则从原始信号减去它;运时xi(t)=xo(t)-Cl(t)，把Xl(t)作为输入信号重复上述过程获得固有模态函数，重复上述过程，得到第i个 固 有模态函数；
[0036] 5.继续上面的步骤，直到不能分离出固有模态函数，剩余的部分被称为一个单调 函数;然后，得到了 η个ci(t)，i = l-n和一个单调函数r(t);
[0037 ]原始信号经过变换W后就可W写作：
[0038] 集合经验模态分解的具体实施方法：
[0039] 岩石破裂发生时声发射信号的振幅会突然增加;然而，当面临破裂时会出现模态 混杂的问题，它会严重混淆时频分析和隐藏固有模态函数的物理意义;所W提出了集合经 验模态分解：
[0040] 1、集合经验模态分解，首先添加白噪音序列到原始信号;第j个原始信号和白噪声 序列合并用xW(t)表示；
[0041] 2、对于每个混合信号，
可W通过使用经验模态分解方法获得;注意， 意味着第j个组合序列的第i个固有模态函数；
[0042] 3、原始信号的第i个固有模态函数ci(t)是所有之前所获的的第i个固有模态函数 的总平均值。
[0043] 希尔伯特变换用于计算瞬时频率和振幅，其具体的实施方法为：
[0044] (1)
[0045]
[0046] 通过运种变换，任何时间序列代表一个分析信号：
[0050] (2)本方法选择有效能量的声发射率作为声发射参数信号。声发射有效能量Eae被 定义为：
[0051 ]
代表的是在t时刻的电压值；
[0052] (3)岩石破裂时会释放能量，声发射信号的产生伴随着弹性能的释放，声发射有效 能量和断裂能量可能存在一定得关系;声发射有效能量分析的目的是测量岩石的断裂瞬间 释放的能量;岩石破裂事件造成的声发射信号持续时间只有10毫秒，而且信号是瞬间的，很 快就会降低到阀值W下，我们假设每一个时间点只有一个岩石的破裂事件产生，所W每个 时刻的声发射信号都是独立的。因此，第i次岩石破裂产生的声发射有效能量为：
[0053]
[0054] (4)如果在一个单位时间内发生了 η次的破坏事件，那么声发射有效能量率馬f为：
[0055]
[0056] (5)经过希尔伯特黄变换，声发射有效能量率可W表示为：
[0057]
[005引在上述方程中：
[0059] xi(t)第i个输入的声发射信号；
[0060] x〇(t) 一个典型的声发射输入信号；
[0061] emax(t)，emin(t)，m(t)输入信号的局部最大值与局部最小值，W及它们的平均值；
[0062] h(t)输入信号与平均值之间的差值；
[0063] Ci(t)第i个固有模态函数；
[0064] r(t) -个单调函数；
[0065] xW(t)第j个典型声发射信号与白噪音序列的混合信号；
[0066] 的第j个混合信号的第i个固有模态函数；
[0067] y(t)x(t)的希尔伯特变换；
[0068] ω (t)，a(t)瞬时频率与瞬时振幅；
[0069] Eae声发射信号有效能量；
[0070] 巧》C第i次岩石破裂事件发生时的声发射有效能量；
[0071] 馬声发射有效能率；
[0072] Wt)与声发射信号强度有关的电压；
[0073] Δ ti每一个岩石破裂事件持续的时间；
[0074] 〇如)第j个固有模态函数的瞬时频率；
[007引 w(t)，a/ (t)第j个固有模态函数的瞬时振幅和它的对数。
[0076]岩石材料的室内实验是一个开放系统热力学过程，能量从岩石实验系统到岩石试 件连续不断的转化成Ξ种能量:弹性能量、塑料能和消失的能量。弹性变形产生弹性能量可 W在卸载时完全释放，相反，塑性变形产生的能量无法释放。此外，能量的消失与福射现象， 相变和化学能源、压裂能等有关。为了简化分析方法，我们假设消失的能量等于从弹性能中 生成的破裂能。除此之外，先前的研究表明，声发射有效能量率与断裂能量释放率在整个故 障过程中是线性相关，无论材料是金属或非金属，宏观或微观尺度，同时声发射有效能量率 可W反映在上述方程瞬时振幅中。因此，基于声发射信号参数的希尔伯特黄变换提供了一 个更好的区分岩石破坏过程不同阶段的方法，运是相对于那些只利用声发射信号分析方法 的优势所在。由此可得，该方法可W用来预测岩石试件的破坏。
【主权项】
1. 一种基于用于岩石损伤特征辨析的声发射信号处理方法，其特征是：声发射信号处 理方法，包括如下步骤： 步骤一:在岩石力学特性的测试系统上加上声发射检测系统，采集岩石加载过程中声 发射信号的变化； 步骤二:将声发射系统采集到的岩石试件破裂过程中声发射信号导出； 步骤三:对导出的声发射信号进行提取，提取其中的有效能量作为声发射参数信号，以 及提取声发射信号中振幅比较大的信号作为分析对象； 步骤四：将步骤三提取的信号分别进行希尔伯特黄变换包括两个部分:变换一为经验 模态分解或集合经验模分解变换;变换二为希尔伯特谱分析； 步骤五:对变换结果进行分析，进而通过声发射信号从宏观和微观两个层面上来说明 岩石破坏过程中的损伤与衰减。
【专利摘要】一种用于岩石损伤特征辨析的声发射信号处理方法，属于损伤辨析的声发射信号处理方法。该方法通过对岩石力学特性测试系统中采集到的声发射信号进行变换，对声发射信号进行提取，选取声发射信号中的有效能量参数，通过电压对时间的积分得到有关于能量的声发射参数信号，然后分别对比较强的声发射信号和声发射参数信号进行希尔伯特黄变换，通过变换提取信号中的瞬时频率与瞬时振幅，瞬时振幅代表从岩石试件在实验过程中释放的能量，瞬时频率与岩石试件在加载过程中的损伤的衰减相关。将原始的声发射信号通过希尔伯特黄变换分解成不同的固有模式函数，然后获得希尔伯特光谱，这种变换方法使得加载下岩石的损伤描述更加清楚，更加直观。
【IPC分类】G01N29/44
【公开号】CN105486756
【申请号】CN201510863621
【发明人】彭维红, 刘峰宇, 张骥, 张海翔, 杜梦琳, 曹国华, 程红梅, 张志镇
【申请人】中国矿业大学
【公开日】2016年4月13日
【申请日】2015年12月1日
最新回复(0)