一项石油运移方向的示踪方法

文档序号:10623816阅读:285来源:国知局
一项石油运移方向的示踪方法
【专利摘要】一项石油运移方向的示踪方法,包括如下步骤:1)采集来源于同一烃源岩的同源油的至少两份石油样品;2)样品分析,进行傅立叶变换离子回旋共振质谱分析;3)数据处理,提取出N1S1类化合物的碳数、丰度及DBE结构信息,得到DBE=14、DBE=15、DBE=17三类N1S1化合物的丰度强度;4)确定石油运移方向:根据DBE=14、DBE=15、DBE=17三类N1S1化合物的丰度强度,分别计算三类化合物的丰度强度百分比,然后将三类化合物的丰度百分比进行做图分析,DBE=15类N1S1化合物丰度百分比增加方向即是石油运移方向。
【专利说明】
-项石油运移方向的示踪方法
技术领域
[0001] 本发明设及一种示踪方法,具体设及一种石油领域的示踪方法。
【背景技术】
[0002] 利用各种地球化学参数研究油气运移方向、路径甚至计算油气运移距离国内外学 者已进行了许多卓有成效工作,其中含N、S、0杂原子非控化合物的研究得到迅速发展。国 内外学者先后开发了各类含N、含S和含0化合物来进行油气运移的研究,并成功提取了多 种含N、含S和含0化合物指标,并在油气勘探实践中得到了广泛应用,取得了明显效果。
[0003] 原油中的有机含氮化合物主要是W中性(非碱性)的化咯类芳香化合物形式出 现,而碱性的化晚类芳香化合物则相对较少。一般说来,碱性氮的量占总氮含量的30%左 右,也就是说,原油中非碱性氮的含量占总氮含量的立分之二W上。Richter等用正戊烧沉 淀多个原油残余物中的渐青质,发现渐青质馈分中氮化合物含量最高,即多环芳香结构含 氮化合物主要存在于渐青质中。从高沸点到低沸点,氮的含量急剧下降。研究表明,非控中 的化咯类含氮化合物具有显著的油气运移分馈效应,可作为油气运移的示踪指标。其中巧 挫类化合物在油气运移研究中无疑是重点关注对象。巧挫类化合物应用于研究石油运移的 机理是因为,由于含有氮原子杂环化的巧挫类分子具有较强的极性,故其可W通过氮原子 键合的氨原子与地层中的有机质或粘±矿物上的负电性氧原子构成氨键(对于化晚型含 氮化合物来说,认为是通过离子键和氨键同时与储层表面及岩石中的水发生作用的),使得 部分巧挫类分子滞留在油气运移路径上或在储集层中,从而在油气运移过程中出现巧挫类 的地层色层分馈效应,见图1A。油气运移过程中巧挫类化合物由于分子结构的不同而产生 运移分馈效应:①氮屏蔽型异构体(1-或8-烷基巧挫,图IB ;图IC型苯并巧挫)相对于氮 裸露型异构体富集(非1-或非8-烷基巧挫,图IB ;图ID型苯并巧挫);②烷基巧挫相对 于烷基苯并巧挫富集;③高碳数同系物相对于低碳数同系物富集(Li等,1995 barter等, 1996)O
[0004] 原油中的有机含氧化合物包括脂肪酸、环烧酸及酪类等,目前应用于油气运移研 究的主要是酪类化合物。现有的分析研究仅设及到一些简单的烷基酪,对于多环酪的分析 应用很少。烷基酪类化合物可显示与化咯类含氮化合物类似的运移分馈效应,但控制酪类 化合物分布的因素不完全等同于含氮化合物。活性径基官能团暴露的烷基异构体(无邻位 取代)更易从油相分离进入水相、矿物表层(Taylor, 1994)。在石油运移过程中,烷基酪可 能大量进入水相,并逐渐与岩石表面发生关系。运个过程取决于酪类化合物在油、水中的分 离性质和在岩石、流体中的分布性质,而运些性质又取决于分子的结构。烷基酪氨氧官能团 暴露异构体有较低的分离系数,比官能团屏蔽异构体更易分离进入水相。因氨键是控制酪 类在矿物表面吸附自由能大小的主要因素,官能团暴露的烷基酪很可能在岩石表面有相对 较高的亲合性,随着原油的运移,官能团立体屏蔽异构体可能相对富集(图1E)。 阳005] 原油中的有机含硫化合物主要有硫醇、硫酸、环硫酸、二硫化物、嚷吩类、酱体硫戊 烧与蕾体硫戊烧等。目前应用的较多的是嚷吩类化合物,如二苯并嚷吩分子参数,除了用来 作为成熟度参数和用于判识有机质沉积环境外,还被用于研究油气运移。二苯并嚷吩类借 助于嚷吩环硫原子上的孤对电子(图1F),与周缘介质中正电性的氨原子形成氨键,导致二 苯并嚷吩类分子被介质吸附,同样引起运移分馈效应。与烷基巧挫类的分馈效应相似,1-或 8-烷基二苯并嚷吩异构体属于屏蔽型分子,氨键相对不易形成,随着运移距离增长,屏蔽 型或半屏蔽型(1-MDBT,1,4-DMDBT)分子的数量也相对增加,也导致4-/1-MDBT,2,4-/1, 4-DMDBT 和 4,6-/1,4-DMDBT 值下降。
[0006] 现有的石油中含N含0含S油气运移示踪化合物主要是低分子量的弱极性的简单 化合物,对运些化合物的分离分析主要采用传统的色谱质谱技术。研究表明,运些低分子量 的简单化合物容易受到生物降解作用的影响(黄海平,Later),分析研究大分子极性化合 物,找出具有较强抗生物降解作用同时又能反映油气运移的示踪化合物是有机地球化学者 追求的目标之一。石油中的极性大分子通常具有稠环结构,同时具有极性杂原子,传统色谱 分离质谱检测技术已不适用。同时,目前已取得的油气运移研究成果,主要是基于低分子量 低极性的巧挫类、烷基酪类和苯并嚷吩类化合物的色谱质谱分析和运用。由于运些小分子 化合物容易受到生物降解和水洗作用的影响,同时样品分析需要进行复杂的前处理,处理 过程容易造成样品二次污染和目标化合物的损失。因此,本领域亟需一种新的简化样品前 处理的更加简便、可靠的更适合于示踪石油运移方向的方法。

【发明内容】

[0007] 为了解决上述现有技术中存在的问题,本发明提供了一种全新的石油运移方向的 示踪方法。本发明的方法样品简化了前处理,避免样品污染和目标化合物的后期损失,同 时,选取有代表性的原油样品进行分析和图谱解析,结合地质背景分析,筛选出有效的反映 油气运移的极性大分子示踪化合物指标,并形成研究石油运移方向的示踪方法,为油气勘 探开发提供新的科研支撑。
[0008] 本发明的石油运移方向的示踪方法利用了傅立叶变换离子回旋共振质谱(FT-ICR M巧技术。
[0009] 傅立叶变换离子回旋共振质谱(FT-ICR M巧技术的发展为研究石油中的极性稠环 大分子提供了新的科研手段。
[0010] FT-ICR MS的核屯、是一个处于高强磁场的大约6厘米长的圆筒状(也可W是其它 形状)的分析池,它是由=对电极,即一对捕陷电极,一对激发电极和一对检测电极构成。 捕陷电极用于捕获进入分析池中的带电离子,激发电极和检测电极是在垂直磁场方向上设 置的两组互相垂直的电极,一组电极用于激发带电离子,使其能够W较大半径产生旋转运 动;而另一组电极用W接收由周期性运动于两极之间带电离子产生的交变电流,检测电极 间接收的高频电流周期与离子在池中的回旋运动的周期相同,由于不同质荷比离子的回旋 周期不同,根据运一原理可W通过检测电流信号的频率计算出离子的质荷比,而信号的强 度反映离子的丰度。FT-ICR MS -般配置S种离子源,包括电喷雾电离源巧SI)、大气压光 电离源(APPI)和大气压化学电离源(APCI)。FT-ICR MS的数据处理有自带的数据处理软 件]。简言之,即将所有信噪比大于6的质谱峰导出到Excel表中,将质谱仪器所测并经内部 校正后的 IUPACdnternational Union of Pure and Applied Qiemistry)质量数(IUPAC Mass)通过下式转换为Ken化ick质量数化en化ick Mass): W11] Kemlrick 质量数=IUPAC 质量数 X (14/14. 01565)
[0012] 转换后的Kemlrick质量数与其最接近的整数质量的差值定义为质量偏差 Ofemlrick Mass Defect, KMD)。Kemlrick质量数的实质是将邸2的相对分子质量14. 01565 定义为整数质量单位,即14. 00000,运样转换后的质量表中所有相差14的整数质量单位所 对应的化合物即具有相同的母体结构单元,但具有不同的亚甲基数,也就是取代基不同的 同类型化合物具有相同的KMD数值。通过KMD值大小可W快速鉴定同类型化合物;通过分 子量计算程序计算出各个化合物分子中C、H、S、N、0等原子的组合方式,得到各质谱峰对应 的分子式(C;HhSgN"0。,C、h、n分别为分子中碳、氨、氮原子个数),最终能得到样品中所有类 型化合物的分子组成信息及其对应的等效双键数值ouble band equivalence, D邸),也即 分子结构中环烧环数和双键个数之和:DBE = c-h/2-n/2+l。
[0013] 本发明发现,通过傅立叶变换离子回旋共振质谱(FT-ICR M巧技术开发,对原油样 品进行分析和图谱解析,能够提取有效的反映油气运移的极性大分子示踪化合物指标。
[0014] 本发明的示踪方法总体上包括样品采集、样品分析、数据处理和确定石油运移方 向的步骤。
[0015] 本发明的方法示踪方法具体包括如下步骤:
[0016] 1)采集来源于同一控源岩的同源油的至少两份石油样品;
[0017] 。样品分析: 阳01引 2a)分别将每一份石油样品按5-20mg石油样品/Iml甲苯、优选IOmg石油样品 /Iml甲苯的浓度溶解于0. 2-5ml、优选Iml甲苯中,充分混合得到第一混合液,
[0019] 2b)取5-50 y L优选20 y L第一混合液按照每5-50 y L第一混合液/Iml第二混 合液的浓度、优选按照每20 y L第一混合液/Iml第二混合液的浓度溶于第二混合液,充分 混合得到第S混合液,所述第二混合液为体积比为0. 8-1. 2:1、优选为1:1的甲苯:甲醇的 混合溶液, 阳020] 2c)向第S混合液中加入1-20 y L、优选10 y L氨水,所述氨水的浓度为 20-50wt%、优选为28wt%,充分混合得到第四混合液,
[0021] 2d)对第四混合液进行傅立叶变换离子回旋共振质谱分析; 阳0巧扣数据处理,提取出NlSl类化合物的碳数、丰度及DBE结构信息,得到D邸=14、 D邸=15、D邸=17 S类NlSl化合物的丰度强度;
[0023] 4)确定石油运移方向:
[0024] 根据0邸=14、0邸=15、0邸=17^类化51化合物的丰度强度,分别计算^类化 合物的丰度强度百分比,然后将S类化合物的丰度百分比进行做图分析,DBE = 15类NlSl 化合物丰度百分比增加方向即是石油运移方向。
[00巧]在本发明的一个优选的实施方式中,在步骤1)所述的采集的过程中,在研究区域 尽可能多布点采集。
[00%] 在本发明的一个优选的实施方式中,在步骤1)所述的采集的过程中,避免原油样 品受到次生污染。
[0027] 在本发明的一个优选的实施方式中,在步骤1)所述的采集来源于同一控源岩的 同源油的至少两份石油样品之后,在进行步骤2)之前,将所述石油样品加盖密封在O-ICTC 下冷藏保存。
[0028] 在本发明的一个优选的实施方式中,步骤2d)所述的傅立叶变换离子回旋共振质 谱分析使用傅立叶变换离子回旋共振质谱仪,其采用ESI电离源,负离子模式。
[0029] 在本发明的一个优选的实施方式中,步骤2d)所述的傅立叶变换离子回旋共 振质谱分析的条件参数为:进样速度为IOO-SOOii LA、优选为180 iiL/h,极化电压为 2500-6000V、优选为4000V,毛细管入口电压为3000-6000V、优选为4500V,毛细管出口电压 为-600 至-200V、优选为-320V。
[0030] 在本发明的一个优选的实施方式中,步骤2d)所述的傅立叶变换离子回旋共振质 谱分析的条件参数为:离子源六极杆累积时间为0. 001-0. 05s、优选为0.0 ls,离子源六极 杆直流电压为1-5V、优选为2. 4V,射频电压为100-800化-P、优选为300化-P ;四极杆Qlm/z 为100-500、优选为300,射频为200-600化-P、优选为400化-P。
[0031] 在本发明的一个优选的实施方式中,步骤2d)所述的傅立叶变换离子回旋共振质 谱分析的条件参数为:碰撞池氣气流量为0. 1-0.化? S-I、优选为0.化? S-I,碰撞能量为-3 至-IV、优选为-1. 5V,胆集时间为0. 01-0. 5s、优选为0. 2s。
[0032] 在本发明的一个优选的实施方式中,步骤2d)所述的傅立叶变换离子回旋共振质 谱分析的条件参数为:激发衰减为10-15地、优选为11. 75地。
[0033] 在本发明的一个优选的实施方式中,步骤2d)所述的傅立叶变换离子回旋共振质 谱分析的条件参数为:采集质量范围为50-1500Da、优选为200-900Da,采样点数为2-8M、优 选为4M。
[0034] 在本发明的一个优选的实施方式中,步骤2d)所述的傅立叶变换离子回旋共振质 谱分析的条件参数为:扫描谱图叠加10-500次、优选64次W提高信噪比。
[0035] 在本发明的一个优选的实施方式中,在本发明的步骤3)所述的数据处理是在分 析结束之后利用仪器自带专用软件进行数据处理,来提取出NlSl类化合物的碳数、丰度及 D邸结构信息,得到DBE = 14、D邸=15、D邸=17 S类NlSl化合物的丰度强度。
[0036] 本发明的发明人出乎意料地发现,原油中含有丰富的DBE = 14、D邸=15、D邸= 17 =类NlSl化合物。通过结构分析,=类NlSl化合物是巧挫与嚷吩类化合物的稠合体,结 构中都含有一个极性N原子子和带有一对孤电子的S原子,N原子容易通过氨键与地质接 触体结合发生附着,同时分子上S原子上有对孤电子,易与周缘介质中的H原子发生作用形 成氨键,因此NlSl极性大分子化合物在油气运移过程中发生分馈作用,是很好的油气运移 示踪化合物,同时具有很好地抗生物降解作用。D邸=14、DBE = 15、DBE = 17 S类NlSl 化合物结构相似,但也存在差异。DBE = 15的NlSl类化合物主要有一个六元环处于极性 N原子的屏蔽位置上,降低了化合物的氨键作用力,因而其与周缘介质的附着能力相对DBE =14和DBE = 17的NlSl化合物而言要弱,不易被吸附,随油气运移距离的增加,其化合物 丰度相对富集。
[0037] 分析表明,原油中富含的DBE = 14、D邸=15、DBE = 17 S类NlSl化合物根据结 构推断具有图2的结构类型。
[0038] 在本发明的步骤4)中,根据样品分析处理后DBE = 14、D邸=15、D邸=17 S类 NlSl化合物的丰度强度,分别计算=类化合物的丰度强度百分比,然后将=类化合物的丰 度百分比进行做图分析,DBE = 15类NlSl化合物丰度百分比增加方向既是原油运移方向。
[0039] 本发明的有益效果在于:本发明采用高分辨质谱分析技术,可W分析检测传统色 谱质谱技术无法检测到的极性大分子,实验分析无需进行前处理,没有后期污染和损失,使 用的示踪化合物为极性Nisi类大分子,不受后期生物降解等次生作用的影响,根据DBE = 14、15、17 S类NlSl化合物之间的相对丰度的变化关系,W DBE = 15的NlSl类化合物相对 丰度增加方向示踪原油运移方向,判识方法简单明了,更加适用于示踪石油的运移方向。
【附图说明】
[0040] 图1为巧挫类含N化合物结构与运移分馈效应机理的示意图。
[0041] 图2为Le化C组原油中NlSl化合物的DBE、丰度、碳数分布及关键化合物结构图。
[0042] 图3为Le化C组原油运移方向的确定的示意图。
【具体实施方式】
[0043] W下通过具体的实施例来对本发明进行进一步的说明。需要注意的是,W下实施 例均为阐释性的,而并不对本发明要求保护的范围构成任何限定。本领域技术人员应理解 的是,对于W下【具体实施方式】的任何改进或变化,在没有脱离本发明的精神和范围的情况 下,均在本发明要求保护的范围内。
[0044] 实施例1 W45] 1)石油样品采集:
[0046] 多布点采集西加拿大盆地上泥盆统的化mb巧-Meadowbrook生物礁链油气藏的石 油样品。采集过程中确保原油样品不受到次生污染,样品采集后加盖密封低溫冷藏保存。
[0047] 。样品分析: W48] 取大约IOmg原油样品溶于1血甲苯中,再取其中20 y L样品溶液再溶于1血甲 苯:甲醇(1:1/V: V)混合溶液,向所得溶液中加入10 y L 28wt%氨水,轻轻振荡使其混合均 匀,然后进行傅立叶变换离子回旋共振质谱分析。分析仪为美国布鲁克SOlariX型FT-ICR MS,磁场强度12. 0T。采用ESI电离源,负离子模式。FT-ICR MS主要仪器参数:进样速度 180 y LA,极化电压4000V,毛细管入口电压4500V,毛细管出口电压-320V,离子源六极杆 累积时间0.0 ls,离子源六极杆直流电压2. 4V,射频电压300化-P ;四极杆Qlm/z 300,射频 400化-P ;碰撞池氣气流量0.化.S-I,碰撞能量-1. 5V,胆集时间0. 2s,激发衰减11. 75地, 采集质量范围200-900Da,采样点数4M,扫描谱图叠加64次W提高信噪比。 W例扣数据处理:
[0050] 分析结束后利用仪器自带专用软件进行数据处理,提取出NlSl类化合物的碳数、 丰度及DBE结构信息。参考图2。从分析结果可W看出,原油中含有丰富的DBE = 14、DBE =15、DBE = 17 S类NlSl化合物。通过结构分析,S类NlSl化合物是巧挫与嚷吩类化合 物的稠合体,结构中都含有一个极性N原子和带有一对孤电子的S原子,N原子容易通过氨 键与地质接触体结合发生附着,同时分子上S原子上有对孤电子,易与周缘介质中的H原子 发生作用形成氨键,因此NlSl极性大分子化合物在油气运移过程中发生分馈作用,是很好 的油气运移示踪化合物,同时具有很好地抗生物降解作用。D邸=14、D邸=15、D邸=17 S类NlSl化合物结构相似,但也存在差异。D邸=15的NlSl类化合物主要有一个六元环 处于极性N原子的屏蔽位置上,降低了化合物的氨键作用力,因而其与周缘介质的附着能 力相对DBE = 14和DBE = 17的NlSl化合物而言要弱,不易被吸附,随油气运移距离的增 加,其化合物丰度相对富集。
[0051] 4)原油运移方向的确定:
[0052] 分析表明,原油中富含08£=14、0邸=15、0邸=17^类化51化合物,运^类化 合物根据结构推断具有图2的结构类型。
[0053] 根据样品分析处理后DBE = 14、D邸=15、D邸=17 S类NlSl化合物的丰度强 度,分别计算=类化合物的丰度强度百分比,然后将=类化合物的丰度百分比进行做图分 析,D邸=15类NlSl化合物丰度百分比增加方向即是原油运移方向。参考图3和表1。
[0054] 表1 Le化C组原油中不同DBE的NlSl类化合物的相对丰度(% ) 阳化引
[0056] 本实施例所用的代表性原油取自西加拿大盆地上泥盆统的化mb巧-Meadowbrook 生物礁链油气藏。由于化mb巧-Meadowbrook中Le化C生物礁群镶边的白云岩台地形成了良 好的油气运移的通道,盆地中南部原油沿着化mbey-Meadowbrook链状生物礁进行了长达 数百千米的运移。样品分别取自化mbey油气藏中经历不同油气运移距离的原油样品。原 油样品经ESI电离源负离子模式FT-ICR MS分析后,在数据处理基础上,获取D邸=14、15、 17的NlSl类化合物丰度,并计算它们的相对丰度百分比,DBE = 15类NlSl化合物丰度增 加方向即是原油运移方向,证实了所述示踪方法的可靠性。
【主权项】
1. 一项石油运移方向的示踪方法,包括如下步骤: 1) 采集来源于同一烃源岩的同源油至少两份石油样品; 2) 样品分析: 2a)分别将每一份石油样品按5-20mg石油样品/lml甲苯、优选10mg石油样品/lml甲 苯的浓度溶解于〇. 2-5ml、优选lml甲苯中,充分混合得到第一混合液, 2b)取5-50 μ L、优选20 μ L第一混合液按照每5-50 μ L第一混合液/lml第二混合液 的浓度、优选按照每20 μ L第一混合液/lml第二混合液的浓度溶于第二混合液,充分混合 得到第三混合液,所述第二混合液为体积比为〇. 8-1. 2:1、优选为1:1的甲苯:甲醇的混合 溶液, 2c)向第三混合液中加入1-20 μ L、优选10 μ L氨水,所述氨水的浓度为20-50wt%、优 选为28wt %,充分混合得到第四混合液, 2d)对第四混合液进行傅立叶变换离子回旋共振质谱分析; 3) 数据处理,提取出N1S1类化合物的碳数、丰度及DBE结构信息,得到DBE = 14、DBE =15、DBE = 17三类N1S1化合物的丰度强度; 4) 确定石油运移方向: 根据DBE = 14、DBE = 15、DBE = 17三类N1S1化合物的丰度强度,分别计算三类化合 物的丰度强度百分比,然后将三类化合物的丰度百分比进行做图分析,DBE = 15类N1S1化 合物丰度百分比增加方向即是石油运移方向。2. 根据权利要求1所述的方法,其特征在于,在步骤1)所述的采集的过程中,避免原油 样品受到次生污染。3. 根据权利要求1或2所述的方法,其特征在于,在步骤1)所述的采集来源于同一烃 源岩的同源油的至少两份石油样品之后,在进行步骤2)之前,将所述石油样品加盖密封在 0-10 °C下冷藏保存。4. 根据权利要求1-3中任意一项所述的方法,其特征在于,步骤2d)所述傅立叶变换离 子回旋共振质谱分析使用傅立叶变换离子回旋共振质谱仪,其采用ESI电离源,负离子模 式。5. 根据权利要求1-4中任意一项所述的方法,其特征在于,步骤2d)所述的傅立叶变换 离子回旋共振质谱分析的条件参数为:进样速度为100-500 y L/h、优选为180 μ L/h,极化 电压为2500-6000V、优选为4000V,毛细管入口电压为3000-6000V、优选为4500V,毛细管出 口电压为-600至-200V、优选为-320V。6. 根据权利要求1-5中任意一项所述的方法,其特征在于,步骤2d)所述的傅立叶 变换离子回旋共振质谱分析的条件参数为:离子源六极杆累积时间为〇. 001-0. 〇5s、优选 为0. 01s,离子源六极杆直流电压为1-5V、优选为2. 4V,射频电压为100-800Vp-p、优选为 300Vp-p ;四极杆 Qlm/z 为 100-500、优选为 300,射频为 200-600Vp-p、优选为 400Vp-p。7. 根据权利要求1-6中任意一项所述的方法,其特征在于,步骤2d)所述的傅立叶变换 离子回旋共振质谱分析的条件参数为:碰撞池氩气流量为〇. 1-0. 5L *s \优选为0. 3L *s \ 碰撞能量为-3至-IV、优选为-1. 5V,贮集时间为0. 01-0. 5s、优选为0. 2s。8. 根据权利要求1-7中任意一项所述的方法,其特征在于,步骤2d)所述的傅立叶变换 离子回旋共振质谱分析的条件参数为:激发衰减为10_15dB、优选为11. 75dB。9. 根据权利要求1-8中任意一项所述的方法,其特征在于,步骤2d)所述的傅立叶变换 离子回旋共振质谱分析的条件参数为:采集质量范围为50-1500Da、优选为200-900Da,采 样点数为2-8M、优选为4M。10. 根据权利要求1-9中任意一项所述的方法,其特征在于,步骤2d)所述的傅立叶变 换离子回旋共振质谱分析的条件参数为:扫描谱图叠加10-500次、优选64次以提高信噪 比。
【文档编号】G01N27/62GK105987946SQ201510095535
【公开日】2016年10月5日
【申请日】2015年3月3日
【发明人】蒋启贵, 刘鹏, 黎茂稳, 申宝剑, 陶国亮, 马媛媛, 宋晓莹
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司石油勘探开发研究院
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1