CN105092547B - 基于图谱形态的混合样品荧光分析方法 - Google Patents
基于图谱形态的混合样品荧光分析方法 Download PDFInfo
- Publication number
- CN105092547B CN105092547B CN201410400793.6A CN201410400793A CN105092547B CN 105092547 B CN105092547 B CN 105092547B CN 201410400793 A CN201410400793 A CN 201410400793A CN 105092547 B CN105092547 B CN 105092547B
- Authority
- CN
- China
- Prior art keywords
- collection
- illustrative plates
- oil
- mrow
- mixing
- 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
Links
Landscapes
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
Abstract
一种基于图谱形态的混合样品荧光分析方法,在对油田现场录井油基泥浆影响的荧光分析中,通过油基泥浆图谱与原油样品图谱建立原始图谱库,以及将不同比例的油基泥浆子图谱与原油样品子图谱相混合以得出混合图谱库,通过待分析的含油样品图谱与混合图谱库中的各个混合图谱进行相似度比对,并消除油基泥浆的影响,得到含油样品中的原油原始图谱。本发明提高了样品含油性分析的准确性,实施本方法过程不需要人工判断图谱中油基泥浆影响的程度,降低了人工判断的误差,便于在油田现场录井作业中推广和应用。
Description
技术领域
本发明涉及油田样品采样分析的技术领域,具体说是一种基于图谱形态的混合样品荧光分析方法。
背景技术
油田现场对荧光钻井液背景扣除法是使用油基泥浆浸泡干砂岩作为背景值,用含油岩屑图谱完全减去背景值图谱,来消除油基的影响。这里面假设含油岩屑中所含的油基与干砂岩泡出的油基含量相同。但实际情况是含油岩屑的孔隙中已经含有原油,所以油基的侵入应较干砂岩少,应用传统方法会产生扣过量的情况。因此,目前的方法对样品荧光分析中油基泥浆影响的排除并不适用。
发明内容
本发明要解决的技术问题是提供一种基于图谱形态的混合样品荧光分析方法。
本发明为解决公知技术中存在的技术问题所采取的技术方案是:
本发明的基于图谱形态的混合样品荧光分析方法,包括以下步骤:
A、分别采集油基泥浆图谱与原油样品图谱,作为原始图谱库;
B、保持波长不变,将油基泥浆图谱中的光强峰高分别按原荧光强度的1%至100%缩放,公差为1%,依次得到100个不同光强峰高的油基泥浆子图谱;保持波长不变,将原油样品图谱中的光强峰高分别按原荧光强度的1%至100%缩放,公差为1%,依次得到100个不同光强峰高的原油样品子图谱;
C、将上述100个油基泥浆子图谱与100个原油样品子图谱分别相叠加,得到10000个不同混合比例的混合图谱,作为混合图谱库;
D、对受油基泥浆污染的含油样品进行荧光图谱采集,得到含油样品图谱;
E、将含油样品图谱与混合图谱库中的所有图谱进行比对,得出混合图谱库中与此含油样品图谱相似度最高的混合图谱,并依此混合图谱中油基泥浆与原油的混合比例认定为含油样品中油基泥浆与原油的混合比例;
F、从含油样品图谱中减去相似度最高的混合图谱中所含的油基泥浆子图谱,即得到含油样品中的原油图谱;
G、根据步骤F中得到的含油样品中的原油图谱进行参数计算,得出含油样品中的含油状况。
本发明还可以采用以下技术措施:
原始图谱库中的图谱荧光强度都为不产生猝灭时的最高荧光强度。
从含油样品采集到的荧光强度不高于原始图谱库中图谱的荧光强度。
在含油样品图谱和某一混合图谱对比时,分别在两个图谱上的极大值点的前后两侧分别取多个特征点,根据特征点与极大值点之间的距离来判断两个图谱间的相似度。
含油样品图谱与混合图谱间的相似度判断步骤如下:
a、计算特征点与极大值点间波长的差值,方式为:
Δwl(y)=|wl(y)-wlpeak|,
其中Δwl(y)为特征点与极大值点波长的差值,
y为特征点与极大值点光强的比例,
wl(y)为特征点的波长,
wlpeak为极大值点的波长;
b、计算含油样品图谱在某一峰值处与混合图谱的相似度,方式为:,
其中,Sm为第m个峰值处含油样品图谱与混合图谱的相似度,
N为特征点个数,
Δpwlyn为含油样品图谱第n个前向特征点与极大值点波长的差值,
Δpwlsn为混合图谱第n个前向特征点与极大值点波长的差值,
Δbwlyn为含油样品图谱第n个后向特征点与极大值点波长的差值,
Δbwlsn为混合图谱第n个后向特征点与极大值点波长的差值;
c、计算含油样品图谱与混合图谱间的相似度,方式为:
其中S为含油样品图谱与混合图谱的相似度,
M为峰的个数,
Sm为第m个峰值处含油样品图谱与混合图谱的相似度。
本发明具有的优点和积极效果是:
本发明的基于图谱形态的混合样品荧光分析方法在对油田现场录井油基泥浆影响的荧光分析中,通过建立原始图谱库以及混合图谱库,将待分析的含油样品图谱与混合图谱库中的各个混合图谱进行相似度比对,并消除油基泥浆的影响,得到含油样品中的原油原始图谱。本发明提高了样品含油性分析的准确性,实施本方法过程不需要人工判断图谱中油基泥浆影响的程度,降低了人工判断的误差,便于在油田现场录井作业中推广和应用。
附图说明
图1是本发明的基于图谱形态的混合样品荧光分析方法的步骤流程图;
图2是实施例中原始图谱库中的原油样品荧光图谱;
图3是实施例中原始图谱库中的油基泥浆荧光图谱;
图4是混合图谱库中30%原油与70%油基泥浆混合比例时的混合图谱;
图5是混合图谱库中70%原油与30%油基泥浆混合比例时的混合图谱;
图6是待分析的含油样品的荧光图谱;
图7是混合图谱库中46%原油与54%油基泥浆混合比例时的混合图谱;
图8是含油样品图谱中扣除54%油基泥浆后的得到的原油图谱;
图9是使用传统的扣除法得到含油样品中的原油图谱。
具体实施方式
以下参照附图及实施例对本发明进行详细的说明。
如图1所示,本发明的基于图谱形态的混合样品荧光分析方法,包括以下步骤:
A、分别采集油基泥浆图谱与原油样品图谱,作为原始图谱库;
B、保持波长不变,将油基泥浆图谱中的光强峰高分别按原荧光强度的1%至100%缩放,公差为1%,依次得到100个不同光强峰高的油基泥浆子图谱,每个油基泥浆子图谱对应于1%、2%、3%直至100%比例的油基泥浆含量;保持波长不变,将原油样品图谱中的光强峰高分别按原荧光强度的1%至100%缩放,公差为1%,依次得到100个不同光强峰高的原油样品子图谱,每个原油样本子图谱对应于1%、2%、3%直至100%比例的原油含量;
C、将上述100个油基泥浆子图谱与100个原油样品子图谱分别相叠加,即将每个油基泥浆子图谱都与各原油样品子图谱分别叠加过一次,而每个原有样品子图谱也都与各油基泥浆子图谱分别叠加过一次,从而得到10000个不同混合比例的混合图谱,作为混合图谱库;
D、对受油基泥浆污染的含油样品进行荧光图谱采集,得到含油样品图谱;
E、将含油样品图谱与混合图谱库中的所有图谱进行比对,得出混合图谱库中与此含油样品图谱相似度最高的混合图谱,并依此混合图谱中油基泥浆与原油的混合比例认定为含油样品中油基泥浆与原油的混合比例;
F、从含油样品图谱中减去相似度最高的混合图谱中所含的油基泥浆子图谱,即得到含油样品中的原油图谱;
G、根据步骤F中得到的含油样品中的原油图谱进行参数计算,得出含油样品中的含油状况。
本方法适用于油基泥浆的荧光图谱形态与原油的荧光图谱形态不完全一致的情况,即对于油基泥浆与原油样品荧光图谱形态一致的情况本方法一般不适用。
由于在低浓度下样品的荧光强度与样品浓度成正比,当含量较高时会产生猝灭现象,所以原始图谱库中的图谱荧光强度都为不产生猝灭时的最高荧光强度。
从含油样品采集到的荧光强度不高于原始图谱库中图谱的荧光强度。
在含油样品图谱和某一混合图谱对比时,分别在两个图谱上的极大值点的前后两侧分别取多个特征点,根据特征点与极大值点之间的距离来判断两个图谱间的相似度。
含油样品图谱与混合图谱间的相似度判断步骤如下:
a、计算特征点与极大值点间波长的差值,方式为:
Δwl(y)=|wl(y)-wlpeak|,
其中Δwl(y)为特征点与极大值点波长的差值,
y为特征点与极大值点光强的比例,y可取0.5,0.6,0.7,0.8,0.9等,
wl(y)为特征点的波长,
wlpeak为极大值点的波长;
b、计算含油样品图谱在某一峰值处与混合图谱的相似度,方式为:
其中,Sm为第m个峰值处含油样品图谱与混合图谱的相似度,
N为特征点个数,
Δpwlyn为含油样品图谱第n个前向特征点与极大值点波长的差值,
Δpwlsn为混合图谱第n个前向特征点与极大值点波长的差值,
Δbwlyn为含油样品图谱第n个后向特征点与极大值点波长的差值,
Δbwlsn为混合图谱第n个后向特征点与极大值点波长的差值;
c、计算含油样品图谱与混合图谱间的相似度,方式为:
其中S为含油样品图谱与混合图谱的相似度,
M为峰的个数,
Sm为第m个峰值处含油样品图谱与混合图谱的相似度。
实施例:
1、分别采集油基泥浆与原油样品的单组份图谱,作为原始图谱库,如图2、图3分别所示;
2、将各油基泥浆子图谱和原有样品子图谱相混合,得到混合图谱库,限于篇幅,附图中仅给出两个混合图谱的示例,其中图4是混合图谱库中30%原油与70%油基泥浆混合比例时的混合图谱,图5是混合图谱库中70%原油与30%油基泥浆混合比例时的混合图谱;
3、对待分析的未知含油样品进行荧光图谱采集,得到含油样品图谱,如图6所示;
4、将未知的含油样品图谱与混合图谱库中的混合图谱进行相似度比对,找出相似度最高的图谱,如图7是混合图谱库中46%原油与54%油基泥浆混合图谱,两者相似度为82.31%,在混合图谱库中相似度最高;
5、用含油样品图谱减去相似度最高的混合图谱中所含的54%油基泥浆图谱,从而获得原油的原始图谱,结果如图8所示。
针对同一含油样品,图9是使用传统的扣除法得到含油样品中的原油图谱,对比图9和图8中的两个图谱可以明显看出由本发明减少了传统方法中背景值扣过量的情况,通过图谱形态分析得到的图谱更加精确。
在实际的荧光分析过程中,根据情况可以相应增大或者缩小图谱的公差,而图谱库也将相应减少或增大,以得到相对较为宽松或更精确的对比结果。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例公开如上,然而,并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当然会利用揭示的技术内容作出些许更动或修饰,成为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均属于本发明技术方案的范围内。
Claims (5)
1.一种基于图谱形态的混合样品荧光分析方法,其特征在于,包括以下步骤:
A、分别采集油基泥浆图谱与原油样品图谱,作为原始图谱库;
B、保持波长不变,将油基泥浆图谱中的光强峰高分别按原荧光强度的1%至100%缩放,公差为1%,依次得到100个不同光强峰高的油基泥浆子图谱;保持波长不变,将原油样品图谱中的光强峰高分别按原荧光强度的1%至100%缩放,公差为1%,依次得到100个不同光强峰高的原油样品子图谱;
C、将上述100个油基泥浆子图谱与100个原油样品子图谱分别相叠加,得到10000个不同混合比例的混合图谱,作为混合图谱库;
D、对受油基泥浆污染的含油样品进行荧光图谱采集,得到含油样品图谱;
E、将含油样品图谱与混合图谱库中的所有图谱进行比对,得出混合图谱库中与此含油样品图谱相似度最高的混合图谱,并依此混合图谱中油基泥浆与原油的混合比例认定为含油样品中油基泥浆与原油的混合比例;
F、从含油样品图谱中减去相似度最高的混合图谱中所含的油基泥浆子图谱,即得到含油样品中的原油图谱;
G、根据步骤F中得到的含油样品中的原油图谱进行参数计算,得出含油样品中的含油状况。
2.根据权利要求1所述的基于图谱形态的混合样品荧光分析方法,其特征在于:原始图谱库中的图谱荧光强度都为不产生猝灭时的最高荧光强度。
3.根据权利要求1或2所述的基于图谱形态的混合样品荧光分析方法,其特征在于:从含油样品采集到的荧光强度不高于原始图谱库中图谱的荧光强度。
4.根据权利要求3所述的基于图谱形态的混合样品荧光分析方法,其特征在于:在含油样品图谱和某一混合图谱对比时,分别在两个图谱上的极大值点的前后两侧分别取多个特征点,根据特征点与极大值点之间的距离来判断两个图谱间的相似度。
5.根据权利要求4所述的基于图谱形态的混合样品荧光分析方法,其特征在于:含油样品图谱与混合图谱间的相似度判断步骤如下:
a、计算特征点与极大值点间波长的差值,方式为:
Δwl(y)=|wl(y)-wlpeak|,
其中Δwl(y)为特征点与极大值点波长的差值,
y为特征点与极大值点光强的比例,
wl(y)为特征点的波长,
wlpeak为极大值点的波长;
b、计算待测的含油样品图谱在某一峰值处与混合图谱的相似度,方式为:
<mrow>
<mi>S</mi>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<mo>|</mo>
<msub>
<mi>&Delta;pwly</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;pwls</mi>
<mi>n</mi>
</msub>
<mo>|</mo>
<mo>+</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<mo>|</mo>
<msub>
<mi>&Delta;bwly</mi>
<mi>n</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;bwls</mi>
<mi>n</mi>
</msub>
<mo>|</mo>
</mrow>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msub>
<mi>&Delta;pwls</mi>
<mi>n</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</msubsup>
<msub>
<mi>&Delta;bwls</mi>
<mi>n</mi>
</msub>
</mrow>
</mfrac>
<mo>,</mo>
</mrow>
其中,Sm为第m个峰值处含油样品图谱与混合图谱的相似度,
N为特征点个数,
Δpwlyn为含油样品图谱第n个前向特征点与极大值点波长的差值,
Δpwlsn为混合图谱第n个前向特征点与极大值点波长的差值,
Δbwlyn为含油样品图谱第n个后向特征点与极大值点波长的差值,
Δbwlsn为混合图谱第n个后向特征点与极大值点波长的差值;
c、计算待测的含油样品图谱与混合图谱间的相似度,方式为:
<mrow>
<mi>S</mi>
<mo>=</mo>
<mfrac>
<mrow>
<msubsup>
<mi>&Sigma;</mi>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</msubsup>
<msub>
<mi>S</mi>
<mi>m</mi>
</msub>
</mrow>
<mi>M</mi>
</mfrac>
<mo>,</mo>
</mrow>
其中S为含油样品图谱与混合图谱的相似度,
M为峰的个数,
Sm为第m个峰值处含油样品图谱与混合图谱的相似度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410400793.6A CN105092547B (zh) | 2014-08-15 | 2014-08-15 | 基于图谱形态的混合样品荧光分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410400793.6A CN105092547B (zh) | 2014-08-15 | 2014-08-15 | 基于图谱形态的混合样品荧光分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105092547A CN105092547A (zh) | 2015-11-25 |
CN105092547B true CN105092547B (zh) | 2018-02-27 |
Family
ID=54573463
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410400793.6A Active CN105092547B (zh) | 2014-08-15 | 2014-08-15 | 基于图谱形态的混合样品荧光分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105092547B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110488710B (zh) * | 2019-08-26 | 2021-02-12 | 武汉和时利自动化系统工程有限公司 | 基于污水指纹图谱变化的河道排污点排查系统及方法 |
CN113607712B (zh) * | 2021-08-23 | 2023-12-19 | 天津陆海石油设备系统工程有限责任公司 | 一种排除油基泥浆污染的原油荧光光谱分析方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1967213A (zh) * | 2006-11-08 | 2007-05-23 | 北京中医药大学中药学院 | 一种检测不同物质混合均匀度的方法 |
JP2011033346A (ja) * | 2009-07-29 | 2011-02-17 | Shimadzu Corp | クロマトグラフ質量分析データ解析方法及び装置 |
CN101988895A (zh) * | 2009-07-30 | 2011-03-23 | 中国石油化工股份有限公司石油化工科学研究院 | 由近红外光谱预测混合原油中单种原油含量的方法 |
CN102590175A (zh) * | 2012-02-21 | 2012-07-18 | 浙江大学 | 基于拉曼光谱叠加的甲醇汽油甲醇含量快速测定方法 |
CN102928375A (zh) * | 2012-07-25 | 2013-02-13 | 苏州派尔精密仪器有限公司 | 采用紫外光谱分析方法进行地沟油检测的方法 |
CN103645171A (zh) * | 2013-12-03 | 2014-03-19 | 深圳市疾病预防控制中心 | 一种基于拉曼光谱法检测地沟油的方法 |
CN103983603A (zh) * | 2014-05-29 | 2014-08-13 | 南京中医药大学 | 矿物药禹余粮的傅里叶变化红外光谱指纹图谱的建立方法、标准指纹图谱及其应用 |
-
2014
- 2014-08-15 CN CN201410400793.6A patent/CN105092547B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1967213A (zh) * | 2006-11-08 | 2007-05-23 | 北京中医药大学中药学院 | 一种检测不同物质混合均匀度的方法 |
JP2011033346A (ja) * | 2009-07-29 | 2011-02-17 | Shimadzu Corp | クロマトグラフ質量分析データ解析方法及び装置 |
CN101988895A (zh) * | 2009-07-30 | 2011-03-23 | 中国石油化工股份有限公司石油化工科学研究院 | 由近红外光谱预测混合原油中单种原油含量的方法 |
CN102590175A (zh) * | 2012-02-21 | 2012-07-18 | 浙江大学 | 基于拉曼光谱叠加的甲醇汽油甲醇含量快速测定方法 |
CN102928375A (zh) * | 2012-07-25 | 2013-02-13 | 苏州派尔精密仪器有限公司 | 采用紫外光谱分析方法进行地沟油检测的方法 |
CN103645171A (zh) * | 2013-12-03 | 2014-03-19 | 深圳市疾病预防控制中心 | 一种基于拉曼光谱法检测地沟油的方法 |
CN103983603A (zh) * | 2014-05-29 | 2014-08-13 | 南京中医药大学 | 矿物药禹余粮的傅里叶变化红外光谱指纹图谱的建立方法、标准指纹图谱及其应用 |
Non-Patent Citations (1)
Title |
---|
遗传算法用于多组分拉曼光谱定量分析的研究;陈柳等;《光学与光电技术》;20080630;第6卷(第3期);第57-60页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105092547A (zh) | 2015-11-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109838230B (zh) | 油藏水淹层的定量评价方法 | |
CN104806231B (zh) | 稠油蒸汽吞吐水淹层的定量评价方法 | |
US10316655B2 (en) | Method and apparatus for consistent and robust fitting in oil based mud filtrate contamination monitoring from multiple downhole sensors | |
No et al. | Experimental design and optimisation (4): Plackett–Burman designs | |
RU2685600C1 (ru) | Способ определения внутрискважинных притоков флюида при многоступенчатом гидроразрыве пласта | |
CN106442600A (zh) | 页岩束缚水含量确定方法 | |
CN105201493B (zh) | 一种识别水驱油藏储层非均质性的双示踪剂方法 | |
CN109356564B (zh) | 一种基于三元动态参数的水驱油相渗时变规律表征方法 | |
CN106703794B (zh) | 一种碳酸盐岩岩性识别方法 | |
CN105092547B (zh) | 基于图谱形态的混合样品荧光分析方法 | |
CN106897531B (zh) | 一种低渗透石灰岩储层渗透率的定量评价方法 | |
CN104912547A (zh) | 应用电阻率成像测井资料连续定量评价储层非均质特征的方法 | |
EP3550300A1 (en) | Method to determine trace amounts of crude oil by spectroscopic absorption | |
CN106204302A (zh) | 一种计算原始含水饱和度的方法与应用 | |
CN106771071B (zh) | 一种基于油水相渗的密闭取心饱和度校正方法 | |
Velasco-Lozano et al. | Modeling of early-and late-time countercurrent spontaneous imbibition in porous media: A semi-analytical approach | |
CN110909466A (zh) | 提高差孔隙分选储层平均毛管压力曲线计算精度的方法 | |
CN108627490A (zh) | 一种基于三维定量荧光光谱积分计算含油浓度的方法 | |
CN106446587A (zh) | 一种气井压裂返排率校正方法 | |
Malik et al. | Tutorial: maximizing value from mudlogs-integrated approach to determine net pay | |
CN112324421B (zh) | 一种低电阻率稠油油层水淹前后饱和度计算方法 | |
CN106383215B (zh) | 一种高热演化盆地混源天然气的混合比的确定方法 | |
Saraf et al. | Analysis of the effect of residual oil on particle trapping during produced-water reinjection using X-ray tomography | |
CN108254505A (zh) | 一种利用氧同位素计算返排液中地层水含量的方法 | |
Holmes et al. | Mixed reservoir wetting in unconventional reservoirs and interpretation of porosity/resistivity cross plots, derived from triple-combo log data |
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 | ||
CP02 | Change in the address of a patent holder | ||
CP02 | Change in the address of a patent holder |
Address after: 300450 1st floor, no.6, Haitai Huake 10th Road, Huayuan Industrial Zone, Binhai New Area, Tianjin Patentee after: TIANJIN LAND & OCEAN PETROLEUM EQUIPMENT SYSTEMS ENGINEERING Co.,Ltd. Address before: No.16, kejiyan West Road, Nankai District, Tianjin 300192 Patentee before: TIANJIN LAND & OCEAN PETROLEUM EQUIPMENT SYSTEMS ENGINEERING Co.,Ltd. |