(56)--门座式起重机臂架及自重平衡系统的优化设计研究.pdf

《(56)--门座式起重机臂架及自重平衡系统的优化设计研究.pdf》由会员分享，可在线阅读，更多相关《(56)--门座式起重机臂架及自重平衡系统的优化设计研究.pdf（8页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、 收稿日期：2016-09-02 门座式起重机臂架及自重平衡系统的优化设计研究 摘 要：随着经济的发展和国际物流业的繁荣，港口物流对门座式起重机装卸货物的效率、使用的安全可靠性等工作性能提出了更高的要求。因此，对门座式起重机各方面的研究一直是国内外的热点。四连杆铰接组合臂架变幅系统是目前应用较广泛的门座式起重机臂架形式。优秀的变幅系统设计以尽可能降低变幅机构的驱动功率和提高机构的操作性能为目的，臂架系统自重平衡的设计优劣是影响门及整体工作性能、安全性能的重要因素。关键词：臂架系统自重平衡 Matlab 优化设计 UG Ansys 中图分类号：TH21 文献标志码：A 文章编号：1006-894

2、5(2016)10-0062-08 Optimized Study of Portal Cranes Boom and Dead Load Blanced SystemCAO Na(Tianjin Port No.5 Stevedoring Co.，Ltd.，Tianjin 300456，China)Abstract：Along with the development of economy and the prosperity of international logistics industry，efficiency，safety and reliability of port logis

3、tics portal crane regarding loading and unloading of goods are facing higher require-ments Therefore，research on various aspects of portal crane has become a hot topic both at home and abroad Luffing sys-tem of four connecting rod hinge arm is the most widely used form of portal crane jib The portal

4、 crane is used to change theamplitude of the agency The design of the luffing system to reduce the driving power of luffing mechanism and improve the operating performance of institutions is an important factor affecting the performance of the door and the whole work and safety performance.Key words

5、：self balance arm frame system；Metlab；optimization design；UG；Ansys 0 引 言 天津港第五港埠有限公司(以下简称“公司”)有 4 台门机在满载作业时均存在着不同程度的拉杆、冲杆现象，不能达到最大工作幅度。拉杆、冲杆，指的是起重机在进行变幅运动时臂架不能在变幅范围内平稳进行变幅运动，在变幅过程中的某点产生很大的冲击载荷，使臂架系统不能继续伸展。这会带来以下问题：起重机在这种情况下按照原有变幅行程作业极有可能发生变幅齿条的损坏，乃至更严重的门机倾覆事故。门机作业不能达到额定变幅，造成门机抓取范围缩小，需要大量机械辅助作业，加大作业

6、成本。本文以 4043,m 门座式起重机为研究对象，对其进行臂架系统设计方案和理性分析，并以其结构尺寸为依据建立门座式起重机臂架系统的优化数学模型。利用 Matlab 优化工具箱对其进行优化求解，并根据优化结果设计可行的改造方案，然后利用 UG 软件对改造后的臂架系统进行三维建模，利用 ANSYS 软件对其关键部位进行静力学强度校核，以验证优化的效果。1 门机臂架系统四连杆机构合理性分析 1.1 设计目标 4043,m 门机采用的臂架形式是四连杆组合形式，该形式的臂架系统设计目标为：四连杆各部位长度合理，使吊重(货物)在变幅过程中运行轨迹趋近于平行于水平面的直线；门机变幅时由货物产生的变幅阻力

7、尽量小；臂架系统各杆件、平衡系统以 DOI:10.14099/ki.tjkj.2016.10.018 2015 年 8 月 曹娜：门座式起重机臂架及自重平衡系统的优化设计研究 63 及货物对臂架下铰点产生力矩的矢量和尽量小，且在变幅过程中的某个位置达到力矩的平衡；平衡重尽量小。验证臂架四连杆机构设计的合理性首先要建立该系统的简化模型，列出各部位参数相互制约的关系式，验证四连杆机构的运动轨迹是否符合预期。1.2 简化模型 门机四连杆组合臂架系统可以简化为如图 1 所示的以点 O、O、A、B 为顶点的四连杆模型，而点A、B、C 组成的三角形刚性梁结构即为象鼻梁的简化形式，臂架四连杆机构的杆件尺寸的

8、设计合理性首先要验证象鼻梁端点在变幅行程内的运动轨迹，也就是C 点的轨迹是否接近一条水平直线。图 1 臂架系统四连杆简化模型 Fig.1 Simplified model of luffing systems four-bar linkage 1.3 象鼻梁端点轨迹计算 为了方便计算，首先将臂架四连杆相关的尺寸参数用字母表示，如图 2 所示：图 2 臂架四连杆 Fig.2 Luffing systems four-bar linkage 图 2 中：1lOA=；2lAB=；31lO B=；4lAC=；5lAD=；6lAE=；7lOE=；81lO A=；911lOO=；101lOO=；111lO

9、O=；A 点坐标为(Xa，Ya)；C 点坐标为(Xc，Yc)；Smax为臂架变幅最大时 OC 的水平距离；Smin为臂架最小幅度时 OC 的水平距离，H 为变幅最大时的起升高度。另外，令1ABC=；2AOA=；3ACC=；4EOA=；51OOE=；61O AB=；7EAO=；81AEO=；9OAA=；10CAA=。则在最大幅度时有：12minmax43mincoscoslSl=(1)12minmax43minsinsinlHl=(2)将两式分别平方并求和得：2221max1max2min2coslSHl S+212min42sinl Hl=(3)令：a122221max4lSHl+；a21ma

10、x2l S；a312l H；则有：a1-a3sin2mina2cos2min a222min1 sin(4)展开并整理得：(2232aa+)22minsin 22132min122sin()0a aaa+=(5)解此一元二次方程可得最大幅度时主臂架摆角：2minarcsin ()()222222131332122223a aa aaaaaaa+)(6)同理可求得最小幅度时主臂架摆角：2maxarcsin ()()222222434335452235a aa aaaaaaa+)(7)式中：a422221min4lSHl+；a51min2l S 以主臂架摆角2作为自变量，2min到2max是其变化

11、范围，则有：12410cossincxll=+(8)12410sincoscyll=(9)利用 MATLAB 软件的绘图功能可以绘制出点C(Xc，Yc)的曲线。编写程序需要的相关数值根据几何关系可计算如下：55124arccosarccosllll=+；(10)3102=；(10)222476447arccos2lmml m+=；(11)952410arctanll=；(12)64 天 津 科 技 第 43 卷 第 10 期 2211910lll=+；(12)2281111 11452cos()lllll=+；(13)22228362 8arccos2llll l+=；(14)7746arcs

12、in(sin)ll=；(15)1184578arccossin()ll=+；(16)9210,2=167892(17、18)带入已知 4043,m 门机的臂架系统参数，在Matlab 中编写绘图函数 M 文件(部分)如下：a134.2239.743216.87216.4232；a22*39.743*34.2；t1acos(0.854/16.423)acos(0.854/5.935)；t4 acos(34.2)2(8.031)2 (26.918)2)/(2*(34.2)*(8.031)得到 C 点轨迹图如图 3。图 3 象鼻梁端点运动轨迹 Fig.3 The trajectory of Trun

13、k beam end point 由图 3 可以看出，Yc 在整个变幅过程中虽然变化值不是很大，但在接近最大变幅时曲线上扬有一个明显拐点，这很可能就是臂架变幅在接近最大时发生失稳的重要原因。1.4 臂架失稳原因的确定 为进一步确定臂架在变幅过程中失稳的原因，我们采用力矩检测工具对门机变幅过程中变幅系统减速器高速轴的力矩进行了测量，得到结果如表 1 所示：表 1 变幅系统减速器高速轴扭矩测量结果 Tab.1 Torque measurement results of high speed shaft of reducer 幅度/m 41 39 36 33 30 27 24 21 18 15 高速

14、轴力矩/(kgm-1)-29-33-33-33-26.4-26.4-19.8-13.2-5.2-6.6 在门机变幅系统的设计中，臂架进行变幅运动时应存在一个力矩平衡点，即力矩零点，而由表 1 可以看出在整个变幅过程中减速器高速轴的扭矩方向始终没有变化(数值始终为负)，因此在变幅过程中不存在扭矩为零的变幅平衡点。这一点导致了门机臂架进行变幅的过程中稳定性不佳，进而不能达到最大幅度的结果。2 门机臂架及平衡系统优化数学模型的建立 门机臂架及平衡系统的优化数学模型包括 3 个要素，即：设计变量、目标函数和约束条件。2.1 设计变量的确定 以上文建立的臂架四连杆简化模型为基础，添加上平衡系统的相关参数

15、，即构成了臂架及平衡系统的整体简化模型(见图 4)。图 4 门机臂架平衡系统优化设计变量 Fig.4Designed variables of crane boom balance system optimization 如图 4 所示，由于点 O 和 O1位置固定，则 l9，l10，l11的长度固定，l5的值一般也取为常量，这里不做改变，需要优化的参数包括：l1主臂架的长度；l2象鼻梁后臂长度；l3大拉杆长度；l4象鼻梁前臂长度；l6象鼻梁铰点与小拉杆下铰点距离；l7主臂架下铰点到小拉杆下铰点距离；l12小拉杆长度；l13平衡 2015 年 8 月 曹娜：门座式起重机臂架及自重平衡系统的优化

16、设计研究 65 梁前臂长度；l14平衡梁后臂长度；l15配重中心点到小拉杆上铰点距离；Gp配重重量。这其中 l14平衡梁后臂长度由于受到机房位置尺寸的影响不宜变动，因此也不做变化。由此就得到了臂架及平衡系统的 10 个设计变量，用向量表示为：Xl1，l2，l3，l4，l6，l7，l12，l13，l15，GpT x1，x2，x3，x4，x5，x6，x7，x8，x9，x10T 根据经验值，设定设计变量上限：lb30；4；24；10；20；5；6；1；6；45,000；设定变量上限为：Ub40，8；35；20；30；10；12；5；13；55,000。2.2 目标函数的确定 2.2.1 象鼻梁端点运

17、动轨迹 根据 1.3 的具体运算，得到象鼻梁端点坐标函数方程：1241012410cossin,sincosccxllyll=+=在cx在最大幅度与最小幅度之间变化时，cy应该接近于cyH=的一条直线，则目标函数可表示为：112410()sincoscf XyHllH=(目标函数 1)2.2.2 变幅阻力矩 臂架的变幅阻力矩可以利用象鼻梁端点 C 轨迹表达式并根据功能原理进行计算。设臂架从位置一摆动到位置二时臂架与水平面夹角2角变化值为2，象鼻梁端点 C 轨迹的高度变化为yc，负载货物重量为 Q。则 C 点为取得克服货物重力升高yc所需的能量，需在臂架上作用力矩M，且有：2cyMQ=由此得：c

18、2yMQ=(19)因此当臂架摆角幅度足够小时，每相邻两个幅度位置的象鼻梁端点高度差值为单位重量货物所引起的在该微量摆角行程上的臂架力矩。2.2.3 臂架平衡系统不平衡力矩 若将臂架平衡系统中各杆件设为均质等截面杆件，则平衡系统重量可表示为：1 12243 3()bjGllll=+412131415()llll+(20)式中：1为主臂架长度重量，2为象鼻梁长度重量，3为大拉杆长度重量，4为平衡系统长度重量。对于本文引用的 4043,m 门机，其长度重量计算如表 2。表 2 臂架平衡系统长度重量计算 Tab.2 Boom balance system length weight calculati

19、on 杆件名称 杆件重量/kg 杆件简化 尺寸/m 长度重量/(kgm-1)主臂架 39,896 34.200 1,166.5 象鼻梁 17,382 22.358 777.4 大拉杆 6,667 27.941 238.6 平衡系统 58,794 26.232 2,241.3 定义了长度重量后，整个臂架平衡系统则完全简化为有长度、有重量，无宽度、无高度的杆件系统模型。由于在变幅过程中主臂架下铰点 O 为臂架系统的回转中心，根据瞬心回转功率法，可以计算各杆件自重对 O 点的力矩。对重应与臂架重力、象鼻梁重力，以及拉杆重力的一半(如不计杠杆系统自重重力)相平衡。2.2.3.1 象鼻梁自重力矩xblM

20、 由于整个象鼻梁长度和为24ll+，而长度重量的定义把所有杆件设为均质杆件，则重心位置应位于242ll+处，即 BC 上位于 B 点422ll处。则象鼻梁对主臂架下铰点 O 的力矩为：412MxblxblmGdd=(21)式中：224()xblGll=+象鼻梁自重(kg)。2.2.3.2 大拉杆自重力矩lgdM 设大拉杆自重的一半lgdG作用在 F 点，则：3lglg32ddlMGdd=(22)式中：3 3lg2dlG=大拉杆自重的一半(kg)。2.2.3.3 主臂架自重力矩zbjM 122coszbjzbjGlM=(23)式中：1 1zbjGl=主臂架自重(kg)。2.2.3.4 平衡系统自

21、重力矩phM 设平衡重系统(包括配重、平衡梁、小拉杆)的合成中心在 P 点，则：111341213141512()phPnMnGlllln=+(24)式中：PG平衡重重量。2.2.3.5 系统不平衡力矩 臂架系统自重力矩为主臂架、象鼻梁、大拉杆及变幅齿条自重力矩之和，即：lgbjzbixbldMMMM=+(25)整个臂架及平衡系统的不平衡力矩为臂架系统 66 天 津 科 技 第 43 卷 第 10 期 力矩与平衡系统力矩之差，即为：wphbjphMMM=(26)f2(X)Mwph(目标函数 2)2.2.4 配重重量 在满足平衡条件的情况下，配重的重量应该取尽量小的值，因此：f3(X)Gp(目标

22、函数 3)其中，Gp为配重重量。2.3 约束条件的确定 在臂架平衡系统的设计中，约束限制主要包括结构布局方面的限制、设计准则要求限制、经验数据限定的变量范围等方面，约束条件的限定既可以避免计算过程中可能发生的超界溢出现象，又可以提高优化搜索速度，在优化设计的过程中具有相当重要的作用。本文根据以下几个方面制定优化约束条件：2.3.1 象鼻梁前后臂的长度比 象鼻梁的前后臂的比值根据起重机设计手册为：象鼻梁前臂长度为后臂的 0.30.5 倍，即：1413140.30.5lll 这里可表现为两个约束关系式：11413()0.30g xll=21314()0.50gxll=2.3.2 臂架与水平面夹角

23、臂架于水平面的夹角根据经验数据其变化范围在 3080 之间，即：2469 这里也需要两个约束表达式：32()06gx=424()09gx=2.3.3 象鼻梁与水平线夹角 象鼻梁同水平线之间的夹角变动范围在 890,之间，即：32452 用约束表达式可以表达为：532()045gx=63()02gx=2.3.4 臂架系统四连杆机构约束 在臂架系统中，大拉杆、象鼻梁后臂、主臂架和1OO组成了一个四连杆机构，四杆的长度就必须满足四连杆的基本条件，也就是在最大幅度时，相邻的两边之和大于对角线；最小幅度时，相邻两边之差应小于对角线。即：78max23()0gxlll=8328min()0gxlll=2.

24、3.5 平衡系统四连杆机构约束 同样，在平衡系统中，平衡梁后臂、小拉杆以及EO和1OO组成了一个四连杆机构，则有：911213max()0gxO Ell=1012131min()0gxllO E=2.3.6 铰点 E 的设计 根据经验，EO的长度与1OO相近则比较可靠。即：11117max()0gxll=127min11()0gxll=2.3.7 平衡量前后臂夹角 由于平衡梁后臂应该设计得尽量长，所以平衡梁前后臂夹角应为钝角，设平衡量前后臂夹角为11，即有：1311()02gx=1411()0gx=还有，应要求臂架处于最小幅度时，平衡重不与机房碰撞。设平衡梁后臂与竖直线夹角为12，即：1514

25、()00gx=2.3.8 臂架不平衡力矩 在最大、最小幅度时，臂架系统的不平衡力矩应有利于臂架恢复正常位置。若设载荷对主臂下铰点 O 引起的力矩有使臂架系统向大幅度方向运动的趋势为正，有向小幅度方向运动的趋势为负，则在最大幅度时，臂架系统不平衡力矩应为负，在最小幅度时，臂架系统不平衡力矩应为正，则有约束：16max()0wphgxm=17wphmin()0gxm=式中：maxwphm、wphminm分别为载荷在最大幅度，最小幅度时的杆件自重不平衡力矩。此外，根据起重机设计手册，臂架平衡系统不平衡力矩应小于臂架四连杆对 O 点力矩之和的十分之一，则有：110wphbjMM 写成约束函数则为：20

26、15 年 8 月 曹娜：门座式起重机臂架及自重平衡系统的优化设计研究 67 18wph()010bjMgxM=3 臂架及平衡系统的优化 3.1 目标函数 f 文件的编写 根据上文对优化目标的定义，一共有 3 个目标函数，将目标函数 f 文件命名为 BJPH_3,mb_MB m，编写程序如下(部分)：function fBJPH_3,mb_MB(x)Smin8.76；Smax39.743；H16.87；l50.854；l914.4；l105；l140.6；已知条件定义(最小、最大幅值；起升高度；象鼻梁铰点偏移距离；主臂架下铰点与大拉杆下铰点水平、竖直距离；象鼻梁后臂距离)。a1x(1).2Sma

27、x.2H.2x(4).2；a22.*x(1).*Smax；a32.*x(1).*H；f(1)abs(ycH)；%,目标函数 1 f(2)MbjMph；%,目标函数 2 f(3)x(10)；%,目标函数 3 3.2 约束条件 M 文件的编写 根据上文对约束条件的定义，将目标函数 M 文件命名为 BJPH_3,mb_YS m，编写程序如下(部分)：functiong，ceqBJPH_3,mb_YS(x)l140.6；%,已知条件象鼻梁后臂长度 Smin12；Smax48；H16.87；L50.854；l914.4；L105；a1x(1).2Smax.2H.2-x(4).2；g(16)Mbjmax-

28、Mphmax；g(17)Mphmin-Mbjmin；g(18)Mbj-Mph-0.1.*Mbj；ceq；3.3 优化函数的调用及优化结果 调用优化工具箱的 fgoalattain 函数需要首先进入 MATLAB 软件主页面，开启优化工具箱后在弹出的优化工具箱界面中选择优化函数以及设定基本参数，优化结果显示(部分)：x1.0,e04*0.003,40.000,0,i 0.000,60.000,0,i 0.002,8 0.000,9 0.000,3 0.000,8 5.300,00.000,0,i fval 1.0,e04*0 0 5.300,00.000,0,i 根据 MATLAB 计算结果，4

29、043,m 门座式起重机优化前后相关参数变化如表 3。表 3 优化参数比较表 Tab.3 Comparison of optimization parameters 设计变量 优化前 优化后主臂架的长度 34.2,m 33.93,m象鼻梁后臂长度 5.9,m 5.97,m 大拉杆长度 27.9,m 28.7,m 象鼻梁前臂长度 16.4,m 16.29,m象鼻梁铰点与小拉杆下铰点距离 26.9,m 27.37,m主臂架下铰点到小拉杆下铰点距离 8.03,m 7.46,m 小拉杆长度 9.395,m9.14,m 平衡梁前臂长度 2.54,m 3.32,m 配重中心点到小拉杆上铰点距离 8.297

30、,m8.07,m 配重重量 45,000,kg53,023,kg 4 优化后模型的三维建模及静力分析 4.1 三维模型建立 根据优化设计结果，在调整主臂架、象鼻梁前后臂、大拉杆及小拉杆的长度后，在 UG 中建模如图59。图 5 主臂架三维模型 Fig.5 3D model of the main arm 图 6 象鼻梁三维模型 Fig.6 3D model of the Trunk beam 图 7 大拉杆三维模型 Fig.7 3D model of the big draw rod 68 天 津 科 技 第 43 卷 第 10 期 图 8 小拉杆三维模型 Fig.8 3D model of

31、the small draw rod 图 9 平衡梁三维模型 Fig.9 3D model of the balance beam 根据各杆件在各自铰点上相连接的相互约束的条件，将主臂架、象鼻梁、大拉杆进行装配，并设计一部件代替门机机身主体，将其设为固定杆件，则主臂架可以绕主臂架下铰点和大拉杆下铰点作变幅运动。臂架系统装配模型如图 10。图 10 臂架系统装配模型 Fig.10 3D model of the arm rack system assembly 4.2 平衡系统的 ANSYS有限元分析 将 UG 中建立的平衡梁三维模型导入 ANSYS，在 ANSYS 中进行参数设置以及网格划分。

32、模型的工程单位、单元类型和材料性质设置 如下：工程单位，长度，mm；质量，kg；重量，N；单元类型，二节点的三维梁单元；图 11 平衡梁及小拉杆有限元模型 Fig.11 Finite element modle of balance beam and small draw rod 材料性质，杨氏弹性模量 E2.06,e5,MPa；泊松比 0.3；建立完成的平衡梁有限元模型如图 11 所示。4.2.1 主要载荷 由结构自重引起的载荷 SG。4.2.2 静力计算结果 下面分别给出结构增加活配重前后平衡梁、小拉杆应力计算结果表 4 所示：表 4 1类工况应力计算结果 Tab.4Stress calc

33、ulation results of Category 1 working conditions 杆件名称 优化之前 优化之后 平衡梁 96.809,MPa 116.941,MPa 小拉杆 104.096,MPa 120.017,MPa 由 表 4 可 以 看 出，小 拉 杆 结 构 应 力 增 加15.921,MPa，平衡梁结构应力增加 20.132,MPa，其所受应力云图见图 12、13。图 12 平衡梁应力云图 Fig.12 Stress cloud of balance beam 图 13 小拉杆应力云图 Fig.13 Stress cloud of small draw rod 4.

34、2.3 弹性极限评估 在结构中，所用到的材料为 Q235C，其机械性能如下：弹性极限 E235,MPa 1、2、3 类工况的许用应力 如下：1 类工况 1E/1.5156.7,MPa 2 类工况 2E/1.33176.7,MPa 综上所述，该平衡系统在优化之后结构仍满足静力强度的要求。2015 年 8 月 曹娜：门座式起重机臂架及自重平衡系统的优化设计研究 69 5 平衡系统改造 根据表 3 比价系统杆件长度变化量相对于初始长度变化幅度很小，而配重重量变化幅度较大，可见该门机逼架平衡系统的主要问题在于配重不足。考虑到臂架系统改造难度大、成本高，最终决定仅对平衡系统进行改造。其方案是：增加的后配

35、重为两个圆柱形箱体，焊接于原配重两侧，以增加配重重量。改造后的优点为：提高门机的抓取范围，减少下舱机械的辅助作业时间，降低作业成本；变幅系统整体 载荷状况改善，降低了便服系统故障频次级门及维修成本。参考文献 1 刘剑波.单臂架门座式起重机臂架变幅系统的优化设计D.上海：上海交通大学，2011.2 孙翔.门座式起重机臂架四连杆结构设计J.工程技术：引文版，2016(5)：248.3 黄陈娣.浅析门座式起重机臂架系统的优化设计J.科技促进发展，2012(5)：128-130.上接第 61 页 入速度的不断增加，泡沫在地层中的封堵效果也在不断增加，但当其注入速度增加到一定值后，泡沫封堵效果逐渐降低，

36、主要是由于随着注入速度的不断增加，地层中的气量过高引起气体窜流，降低了其封堵效果。通过实验可知，该体系的最佳注入速度为2,mL/min，阻力因子为 195。3 结 论 通过起泡剂初筛选、耐油性能评价、界面张力评价和浓度筛选以及稳泡剂浓度筛选，最终筛选出JZ25-1S 油田 A22 井氮气泡沫调驱用氮气泡沫体系为：0.5%,QP18800,mg/LHPAM。通过阻力因子测定实验，优选出 JZ25-1S 油田A22 井氮气泡沫调驱体系的最佳注入气液比为31，最佳注入速度为 2,mL/min，阻力因子为 195。参考文献 1 朱明，姚凯，叶惠民，等.高温氮气泡沫调剖技术在Girasol 油田的应用J

37、.特种油气藏，2015(2)：137-139.2 张杰.栲胶复合氮气泡沫体系调堵剂性能评价及应 用J.断块油气田，2012(5)：642-645.3 蒋晓波.超稠油氮气泡沫凝胶调剖体系的研究与应 用J.中外能源，2012(3)：61-64.4 石海磊，王彦春.海上油田氮气泡沫压锥堵水技术应用J.化工进展，2014(S1)：299-302.5 王晓东，王利敏，张红华，等.河南稠油水驱油藏氮气泡沫调驱体系的研究及应用J.石油地质与工程，2014(4)：141-143.6 马跃，唐晓旭，刘晖，等.海上底水稠油油藏氮气泡沫压锥技术研究与应用J.中国海上油气田，2012(1)：45-50.7 周根荣.蒸汽氮气泡沫调驱实验研究J.科学技术与工程，2012(6)：1393-1396.8 金伟.氮气泡沫驱体系筛选与评价J.内蒙古石油化工，2011(15)：22-23.9 刘清栋，祝红爽，权莉，等.超深稠油泡沫驱泡沫体系优选J.油田化学，2014(2)：247-251.10 曹健.春光油田氮气泡沫堵水技术可行性实验研 究J.石油地质与工程，2015(2)：140-142.