金属学报  2015 , 51 (1): 114-120 https://doi.org/10.11900/0412.1961.2014.00400

高熔点金属区域熔炼中晶体生长角和凝固速率对熔区稳定性的影响

李双明1, 耿振博1, 胡锐1, 刘毅2, 罗锡明2

1 西北工业大学凝固技术国家重点实验室, 西安 710072
2 昆明贵金属研究所, 昆明 650106

EFFECT OF GROWTH ANGLE AND SOLIDIFICATION RATE ON THE FLOATING ZONE STABILITY FOR PROCESSING OF HIGH-TEMPERATURE PURE METALS

LI Shuangming1, GENG Zhenbo1, HU Rui1, LIU Yi2, LUO Ximing2

1 State Key Laboratory of Solidification Processing, Northwestern Polytechnical University, Xi'an 710072
2 Kunming Institute of Precious Metals, Kunming 650106

展开

摘要

对高熔点金属Nb, W, Ta, Mo及Ir电子束区熔熔区高度进行了稳定性分析, 发现在区熔相同尺寸试样时, 能够稳定熔区高度大小排序依次为Nb>Mo>W>Ta>Ir. 计算获得了这5种金属的晶体生长角在8°~13°之间, 发现生长角不为零对大尺寸试样熔区高度起主导作用, 同时金属的实际晶体生长角与界面生长机制有关. 当为粗糙界面生长机制时, 生长角随区熔凝固速率增加变化不大; 而为位错生长机制时, 生长角随区熔凝固速率增加而减少; 如为小面生长机制时, 生长角在低速下会大幅度减小, 并随凝固速率增加而增大. 采用较大的凝固速率(约1 mm/min)有利于控制Ir和Mo晶体生长角变化和熔区高度, 这一点与Mo区熔单晶生长实验结果基本吻合.

关键词: 高熔点金属 ; 电子束区熔 ; 生长角 ; 界面生长机制

Abstract

The height of floating zone and molten zone instability for five pure metals including Nb, W, Ta, Mo, and Ir with high melting points is investigated using electron beam floating zone method (EBFZM). The results show that the height level of floating zone for these five metals are in order with the sequence of Nb>Mo>W>Ta>Ir. The crystal growth angles for these metals are in the range of 8°~13° and the sample in large size can be developed by EBFZM as the growth angle is found not to be zero. Meanwhile, the actual growth angles are related with the interface growth mechanism. For continuous growth mechanism, the growth angles vary slightly with the solidification rate for rough interface, and for dislocation growth mechanism, the growth angles decrease with increasing the solidification rate. If faceting growth mechanism prevails, the growth angles drop remarkably at a low solidification rate and further increase with increasing the solidification rate. Additionally, by employing EBFZM growth of Ir and Mo pure metals, a solidification rate approaching 1 mm/min is available for controlling the growth angle and the height of floating zone. These calculations fit well with the experimental results of Mo single crystal prepared by EBFZM.

Keywords: high-melting point metal ; electron beam floating zone method ; growth angle ; interface growth mechanism

0

PDF (935KB) 元数据 多维度评价 相关文章 收藏文章

金属W, Mo, Nb, Ta, Ir及其合金单晶, 因具有优良的综合性能而广泛应用于电子电气、航空、航天、医学等许多高科技领域[ 1- 3]. 其中, W和Mo单晶因优异的高温抗疲劳性能、低的韧脆转变温度、高真空功率函数, 小的热中子捕获截面、低电阻率与核材料相容性好等综合性能而成为理想的空间反应堆离子能量转换器发射极材料, 已成功应用于TOPAZ型空间反应堆中. Nb单晶正作为热离子燃料元件的接受材料进行开发和应用. 对于Ir单晶, 它是金刚石材料最好的外延生长基材[ 4], 而Ir单晶坩埚拉出的单晶钇铝石榴石(熔点接近2000 ℃)质量高, 使用寿命可达2500~3000 h[ 5].

金属W, Mo, Nb, Ta, Ir都是高熔点材料, 其中Ir熔点最低, 为2443 ℃, W的熔点高达3370 ℃, 因此在熔体法生长单晶时, 并无合适的坩埚以供熔炼, 只能采用无坩埚的电子束悬浮区域熔炼法(electron beam floating zone method, EBFZM), 该方法也是制备难熔金属单晶最有效的工艺方法[ 6]. 目前采用该技术, 国内外已成功制备出了直径超过30 mm的W和Mo难熔单晶[ 7, 8]. 考虑到EBFZM制备难熔金属单晶的关键是保持金属熔区的稳定性, 而熔区稳定性靠表面张力来维持, 并与合金密度密切相关. 这5种金属中, Ir的固态金属密度最大, 为22.42 g/cm3, 比W的密度(19.35 g/cm3)大17%, 较Mo (10.2 g/cm3)大一倍, 因此金属Ir单晶制备难度要远大于W和Mo单晶. 同时, 熔区稳定性还与晶体生长角密切相关[ 9], 尽管国内外已报道成功制备出了W和Mo单晶, 但有关这5种金属晶体生长角与熔区稳定性的关系仍然没见报道.

本工作首先从理论上来讨论EBFZM中熔区的稳定性, 然后计算出了上述5种高熔点金属的晶体生长角, 给出了晶体生长角及凝固速率对熔区稳定性的作用, 获得了难熔金属单晶生长中熔区稳定合理的高度以及凝固速率范围, 上述研究不仅对Nb, Ta及Ir单晶的生长有很好的指导作用, 而且对其它难熔金属晶体生长也有较好的参考价值.

1 理论模型

图1是难熔金属电子束悬浮区域熔炼法的示意图, 其中图1a是熔区向上凝固的形状, 而图1b是熔区向下凝固的形状, 图中a为晶体生长角. 不同材料的a不同, 如对于半导体Si来说, a大约为10°~11°[ 10]; 而InSb的a为25°~30°[ 11], a的大小会严重影响熔区的形状及其稳定性. 对于图1中2种形状的熔区来说, 当熔区高度超过一定值时, 由于液态金属表面张力不能平衡熔区高度产生的静压力, 结果会导致熔区塌漏而使晶体生长实验失败, 为此研究晶体生长角对熔区稳定性影响对区熔晶体生长十分重要.

显示原图| 下载原图ZIP| 生成PPT

图1   电子束悬浮区域熔炼(EBFZM)晶体生长中向上凝固和向下凝固的熔区形状示意图

Fig.1   Schematic shapes of a molten floating zone growth of a cylindrical crystal by electron beam floating zone method (EBFZM) (M—melting interface, S—solidification interface, L—liquid metal, a—growth angle, rc—growing crystal radius, h—height of molten zone, R1 and R2—radii of principal curvature, r—axial coordinate, z—vertical coordinate)(a) upward growth(b) pulling down growth

在图1电子束区熔晶体稳态生长中, 如果不考虑晶体旋转造成的离心力(晶体生长时转速一般为1~2 r/min, 其产生的离心力与表面张力相比, 至少小1个数量级以上[ 9]), 则熔区高度产生的静压力需被金属熔体表面张力所支撑, 为此熔区上存在如下的力平衡关系:

γlv(1R1+1R2)=ρgz

式中, glv为区熔金属的表面张力(其值与金属液/气表面能数值相等); R1和R2是熔区的2个主曲率半径, 见图1b中截面所示; r为液态金属密度; g为重力加速度; z为熔区高度对应的纵坐标.

2 结果分析与讨论

简化条件下, 如果熔区形状为同一直径的圆柱(R1=R2=R), 且熔区侧面保持垂直状态, 即为晶体生长角a为零的情况, 则有:

a=0

此时式(1)可变为:

1R=ρghγlv

式中, R为熔区的半径, h为熔区的高度(见图1所示), 因此液态金属表面张力可支撑的熔区高度为:

h=γlvρgR

从式(3)可看出, 随熔区半径R减小, 区熔熔区稳定性高度增加. 当 R0时, 有 h, 这显然不合理, 因为随熔区高度增加, 熔区在表面张力的作用下, 会自身收缩, 破裂成球形, 以减少表面积, 使得整个体系的自由能减小, 即存在所谓的Rayleigh不稳定性(Rayleigh′s instability)[ 12], 根据Rayleigh不稳定性判据可知, 当熔区高度满足下述条件时, 熔区不会收缩成球形:

h2πR

式(4)中等号为熔区稳定最大高度, 当h大于式(4)给出的数值, 熔区本身发生分解.

表1   5种高熔点金属的物性参数[13,14]

Table 1   Physical parameters for five pure metals with high melting points[13,14]

Metalr / (g·m-3)Tm / ℃gsl / (mN·m-1)glv / (mN·m-1)a / (°)hmax / mm
Ir20.002443411224110.68.4
Mo9.342620464211012.711.9
Nb7.83246834723358.413.7
Ta15.00295041524679.710.1
W16.203370510267612.010.2

新窗口打开

显示原图| 下载原图ZIP| 生成PPT

图2   电子束悬浮区域熔炼下5种金属稳定生长的熔区高度与熔区半径之间的关系

Fig.2   Relationship of heights of floating zone of Ir, Mo, Nb, Ta, and W with the radius of the floating zone stable growth by EBFZM, and imposed with the Rayleigh′s instability results

表1给出了5种高熔点金属物性参数[ 13, 14]. 图2是采用表1中参数, 结合式(3)和(4)计算出的可稳定生长的熔区高度曲线, 其中可支撑熔区高度最大值(hmax)为:

hmax=2πγlvρg

由于式(4)中Rayleigh不稳定性造成的熔区高度与熔区半径成正比, 而式(3)中表面张力稳定的熔区高度与熔区半径成反比, 因此这2种曲线有1个交点, 见图2中2个虚线的交点所示. 在交点的左侧, 表面张力稳定的熔区高度远大于Rayleigh不稳定造成的熔区高度, 因此在该情况下, 区熔熔区高度主要取决于Rayleigh不稳定性造成的影响; 在交点的右侧, Rayleigh不稳定性造成的熔区高度逐渐增加, 而表面张力所支撑的熔区高度随熔区尺寸增加逐渐减小, 因此在交点的右侧, 稳定区熔熔区的高度取决于表面张力的大小, 而交点处对应的熔区半径 R=γlv2πρg, 所以对于a=0的情况, 稳定区熔熔区的高度(h)应该为:

h=2πR(R<γlv2πρg)γlvρgR(Rγlv2πρg)

另外从表1数据还可看出, 高熔点金属的表面张力(glv)相差不大, 而金属的密度相差很大, 因此密度小的金属, 能够在表面张力作用下支撑较大的熔区高度, 如Nb的液态金属密度为7.83 g/cm3, 而Ir为20 g/cm3, 在区熔相同尺寸试样时, 表面张力能支撑的Nb熔区高度就为Ir的2.5倍, 因此图2中表面张力能够稳定的区熔熔区高度排序依次为Nb>Mo>W>Ta>Ir.

同时根据式(5), 在a=0的前提下, 可分别获得高熔点金属Nb, Mo, Ta, W和Ir熔区高度最大值分别为13.7, 11.9, 10.1, 10.2和8.4 mm (表1), 该结果也为图2中曲线与直线交点的纵坐标数值, 而结果对应的横坐标熔区半径在2.5 mm以下. 由此看来, 稳定区熔大尺寸试样, 能够实现晶体生长的熔区高度会相应减小.

由于a=0时, 熔区侧面是垂直的, 因此熔区形状与生长的晶体形状是一致的, 但实验中, 获得的晶体尺寸并不一定与熔区形状相同, 这说明实际凝固中晶体生长角a可以不为零. 另外已有Mo区熔单晶生长实验结果表明, 区熔大尺寸试样时, 能够稳定生长的熔区高度可以超过5 mm[ 15], 比图2中结果要大, 说明晶体生长角能够在一定程度上增大熔区生长的稳定性. 对于 α0, 可按下列情况进行分析:

α0, 只考虑1个主曲率半径

R1主曲率半径的数学表达式代入式(1)中有:

d2zdr2[1+(dzdr)2]3/2=ρgzγlv

结合图1b的坐标, 则:

dzdr=-tan(π2-α)

d2zdr2=ddr(-tanβ)=-1cos2βdβdr

其中, β=π2-α, 将式(8)和(9)代入式(7), 有:

dβdr=-ρgzγlv1cosβ

将式(8)代入式(10), 即:

tanβdβdz=ρgzγlv1cosβ

结合图1b的坐标(由于熔区对称, 仅对一半熔区进行分析), 对式(11)进行积分:

β190osinβdβ=0h/2ρgzdzγlv

h=22γlvρgcosβ1

由于 β1=π2-α1, α1为区熔晶体生长中的实际生长角. 因此式(13)可变为:

h=22γlvρgsinα1

结合a1=0和 α10这2个部分, 则最终表面张力能支撑熔区高度(h)为:

h=22γlvρgsinα1α10+2πRα1=0(R<γlv2πρg)22γlvρgsinα1α10+γlvρgRα1=0(Rγlv2πρg)

另外, 有研究[ 16]表明实际晶体生长角与凝固速率存在如下关系:

α1=α0-kΔTV-1/3

式中, DT和V分别为金属凝固生长的过冷度与凝固速率; α0是凝固速率为零时晶体的生长角, 其与固/液界面能(gsl)和液/气表面能(glv)有关[ 17], 有:

γlvsinα0=γsl

根据表1中固/液界面能和液/气表面能数值, 计算获得的5种金属的晶体生长角 α0在8.4°~12.7°之间, 差别并不是很大. 这5种金属中, 目前仅有W有晶体生长角报道[ 8], 其值为20°±5°. 据表1计算出的W晶体生长角为12.0°, 比较接近测量的最小生长角数值15°, 而与最大值25°有较大的偏离.

在晶体生长中, 根据金属凝固理论可知, ΔT和V存在如下关系[ 18]:

V=μcΔTμsΔT2Aexp(-B/ΔT)

式中, μc为连续生长机制的动力学系数, μs为位错生长机制的动力学系数, 而系数A和B为小面生长动力学系数. 在高熔点金属区熔生长中, 这3种机制都有可能性, 如对Mo区熔单晶生长中, 粗糙界面生长可获得 (111)面取向的单晶, 而如为小面生长方式, 获得的则为多晶组织, 不同界面生长机制获得的晶体取向有较大的差别[ 19].

显示原图| 下载原图ZIP| 生成PPT

图3   生长角不为零电子束悬浮区域熔炼下Mo和Ir金属稳定生长的熔区高度与熔区半径之间的关系

Fig.3   Relationship of heights of floating zone of Ir and Mo with the radius of floating zone stable growth by EBFZM as the growth angle does not equal to zero

考虑到Mo区熔单晶生长有一定的实验报道, 而Ir单晶生长难度较大, 下面以这2种金属区熔生长为例进行分析. 图3是根据式(14)和(15), 在凝固速率为零的条件下, 2种金属熔区稳定性高度与熔区半径之间的关系. 从中可看出, 熔区半径较小的情况下, 式(15)计算获得的熔区高度远大于式(14), 说明小尺寸试样熔区区熔主要为侧面垂直状态(a=0); 而熔区半径较大时, 式(15)和(14)计算结果接近, 说明生长角不为零对熔区高度支撑作用占主导地位. 为此利用生长角不为零可以实现大尺寸试样区熔熔炼, 这一点也是needle-eye technique和Pedestal technique能够区熔生长大尺寸单晶的原因. 如对半导体Si来说, 采用needle-eye technique可以区熔生长直径200 mm的Si单晶, 区熔高度也有16 mm左右[ 17]. 具体对本工作中金属Mo和Ir来说, 在熔区半径达到20 mm时, 式(15)计算获得金属Mo的区熔高度为7.4 mm, Ir为4.6 mm, 而根据式(14)计算的Mo和Ir的结果分别为6.3和4.1 mm, 两者相差较小, 因此对大尺寸熔区, 式(15)和(14)都能为区熔生长提供熔区高度参数.

显示原图| 下载原图ZIP| 生成PPT

图4   Ir和Mo金属电子束悬浮区域熔炼下凝固速率对晶体生长角的影响

Fig.4   Crystal growth angles as a function of the solidification rate for Ir (a) and Mo (b) prepared by EBFZM

显示原图| 下载原图ZIP| 生成PPT

图5   不同晶体生长机制下直径30 mm的Ir和Mo金属电子束悬浮区域熔炼中凝固速率对熔区高度的影响

Fig.5   Heights of floating zone vs the solidification rate of Ir (a) and Mo (b) in diameter of 30 mm prepared by EBFZM at different growth mechanisms

进一步考虑金属Ir和Mo凝固速率对晶体生长角的影响, 如图4a和b所示. 其中式(18)参数取值参考文献[18], 对图4a中金属Ir来说, 取k=1, mc=1, ms=10-2, A=102, B=1; 而对金属Mo来说, 取k=1.5, mc=0.9, ms=6.94×10-3, A=103, B=2. 从图4给出的曲线来看, 尽管参数取值略有差别, 但金属Ir和Mo凝固速率对晶体生长角的影响基本上是一致的.

对于粗糙界面连续生长机制, 随凝固速率的增加, 金属Ir和Mo实际晶体生长角变化不大. 对位错生长机制, 随凝固速率的增加, 2种金属实际晶体生长角变小; 与粗糙界面和位错生长机制相反, 小面生长过程中, 2种金属实际生长角在低速下大幅度减小, 并随凝固速率增加而增大.

考虑到实际凝固过程中, 凝固速率一般处于0.1~1.0 mm/min之间[ 1], 如果凝固速率继续增加, 晶体生长机制会从相应的位错生长机制或小面生长机制转变为粗糙界面生长机制[ 20], 因此图4中位错生长机制和小面生长机制只能在低凝固速率下起作用. 从图4可看出, 该速率下粗糙界面和位错生长机制对实际晶体生长角的影响较小, 而小面生长中实际生长角会变化很大, 上述结论与实验观测结果是一致的, 因为在Mo区熔单晶生长中[ 7], 如果为小面生长机制, 则试样尺寸很不均匀(生长角发生变化), 而如果为取向一致的粗糙界面或位错生长单晶, 试样尺寸的一致性很好(生长角变化较小).

同样, 凝固速率对生长角的影响会造成熔区高度的变化, 图5是金属Ir和Mo在不同晶体生长机制下可稳定支撑的熔区高度. 从中可看出, 在凝固速率0.1~1.0 mm/min之间, 粗糙界面生长的金属Ir和Mo在表面张力的作用下可分别支撑4.8和7.8 mm, 而在小面生长机制下, 只能支撑3.6和4.6 mm. 随凝固速率的增加, 可支撑的熔区高度也相应增大, 因此要实现高熔点金属Ir和Mo区熔单晶生长, 较好的工艺是采用较大的凝固速率, 尽量避开很低的凝固速率. 在正常的凝固速率范围0.1~1.0 mm/min内, 金属Ir不同生长机制稍微有所变化, 如对位错生长机制来说, 可支撑的熔区高度从4.8 mm下降到4.5 mm, 粗糙界面生长机制几乎没有太大的变化, 而小面生长机制, 可支撑熔区高度从3.6 mm增加到了4.2 mm. 上述规律也适用于金属Mo的区熔生长, 因此采用较大的凝固速率(约1.0 mm/min)不仅使生长角变化小, 试样均匀性好, 同样对可支撑的熔区高度来说也是有利的.

如果生长角不为零, 同时存在2个曲率半径的情况下, 式(1)变为一个非线形方程, 不能获得解析解, 只能通过数值计算方法来获得有关熔区高度与生长角的关系, 这需要进一步加以详细讨论. 本计算中将生长角应用到高熔点金属的电子束区熔晶体生长工艺中, 并与熔区稳定性高度和凝固工艺参数(抽拉速率)联系在一起, 且与界面生长机制相关, 这方面的研究揭示出了高熔点金属单晶生长中杂晶出现的原因可能是多方面的, 体现出工艺控制的重要性. 考虑到高温熔体晶体生长实验中还没有方法检测晶体不同生长机制, 为此理论研究中有关熔区高度微小变化造成生长机制改变, 目前还不能通过实验加以验证. 在实验难以检测的前提下, 开展理论研究, 有助于对难熔及高熔点金属电子束区熔晶体生长有更好的认识, 这为后续准确控制晶体生长工艺提供了指导.

3 结论

(1) 计算获得了高熔点金属Nb, W, Ta, Mo及Ir的晶体生长角在8°~13°之间. 在区熔相同尺寸试样时, 能够稳定熔区高度大小排序依次为Nb>Mo>W>Ta>Ir. 大尺寸试样区熔中, 晶体生长角不为零对熔区高度稳定性作用占主导地位.

(2) 高熔点金属Nb, W, Ta, Mo及Ir的实际晶体生长角与界面生长机制相关. 当为粗糙界面生长机制时, 生长角随区熔凝固速率增加变化不大; 而为位错生长机制时, 生长角随区熔凝固速率增加而减少; 如为小面生长机制时, 生长角在低速下会大幅度减小, 并随凝固速率增加而增大. 采用较大的凝固速率(约1.0 mm/min)有利于控制晶体生长角变化和稳定熔区高度.


参考文献

[1] Zee R H, Xiao Z, Chin B A, Liu J. J Mater Proc Technol, 2001; 113: 75

[本文引用: 2]     

[2] Baars R E. Evolution of Refractory Metals and Alloys. Pennsyvania: The Minerals, Metals & Materials Society, 1994: 185

[3] Ohriner E K. Platinum Met Rev, 2008; 52: 186

[本文引用: 1]     

[4] Verstraete M J, Christpehe C J. Appl Phys Lett, 2005; 86: 191917

[本文引用: 1]     

[5] Glebovsky V G, Semonov V N. Int J Refractory Met Hard Mater, 1993-1994; 12: 295

[本文引用: 1]     

[6] Glebovsky V G, Semonov V N. High Temp Mater Proc, 1995; 14: 121

[本文引用: 1]     

[7] Hu Z W, Li Z K, Zhang Q, Zhang T J, Zhang J L, Yin T. Rare Met Mater Eng, 2007; 36: 367

[本文引用: 2]     

(胡忠武, 李中奎, 张 清, 张廷杰, 张军良, 殷 涛. 稀有金属材料与工程, 2007; 36: 367)

[本文引用: 2]     

[8] Glebovsky V G, Semonov V N, Lomeko V V. J Crystal Growth, 1989; 98: 487

[本文引用: 2]     

[9] Duffar T. Crystal Growth Processes Based on Capillarity. Chichester: John Wiley & Sons Ltd, 2010: 204

[本文引用: 2]     

[10] Virozub A, Rasin I G, Brandon S. J Cryst Growth, 2008; 310: 5416

[本文引用: 1]     

[11] Satunkin G A. J Cryst Growth, 2005; 255: 170

[本文引用: 1]     

[12] Johns L E, Narayanan R. Interfacial Instability. New York: Springer-Verlag, Inc., 2002: 8

[本文引用: 1]     

[13] Keene B J. Int Mater Rev, 1993; 38: 157

[本文引用: 1]     

[14] Jiang Q, Lu H M. Surf Sci Rep, 2008; 63: 427

[本文引用: 1]     

[15] Wang H. Master Thesis, Northwestern Polytechncial University, Xi′an, 2007: 23

[本文引用: 1]     

(王 红. 西北工业大学硕士学位论文, 西安, 2007: 23)

[本文引用: 1]     

[16] Weinstein O, Brandon S. J Cryst Growth, 2004; 268: 299

[本文引用: 1]     

[17] Ludge A, Riemann H, Wunscher M, Behr G, Loser W, Muiznieks A, Croll A. In: Duffar T ed., Crystal Growth Processes Based on Capillarity. Chichester: John Wiley & Sons Ltd, 2010: 203

[本文引用: 2]     

[18] Min N B. Physical Fundamentals of Crystal Growth. Shanghai: Shanghai Science & Technology Press, 1982: 398

[本文引用: 1]     

(闵乃本. 晶体生长的物理基础. 上海: 上海科学技术出版社, 1982: 398)

[本文引用: 1]     

[19] Zhang Q, Li Z K, Zhen Y Z, Zhang J L, Hu Z W, Zhang T J, Yin T, Ding X F, Ding X. Rare Met Mater Eng, 2005; 34: 1498

[本文引用: 1]     

(张 清, 李中奎, 郑玉柱, 张军良, 胡忠武, 张廷杰, 殷 涛, 丁学峰, 丁 旭. 稀有金属材料与工程, 2005; 34: 1498)

[本文引用: 1]     

[20] Tiller W A. The Science of Crystallization:Microscopic Interfacial Phenomena. Cambridge: Cambridge University Press, 1991: 87

[本文引用: 1]     

/

深圳SEO优化公司唐山网站建设报价罗湖高端网站设计报价阿里网站seo优化推荐临夏网站建设公司观澜网站优化按天扣费价格仙桃阿里店铺托管推荐拉萨百度标王迪庆网站改版多少钱邢台外贸网站建设公司临沂网站优化软件公司南平seo网站推广价格黑河关键词排名包年推广推荐福田网站优化推广公司玉溪网站推广方案哪家好岳阳网站推广系统公司张家界网站搜索优化价格黄石品牌网站设计多少钱光明seo网站推广公司淮北百度网站优化光明网站设计价格濮阳百度关键词包年推广哪家好光明百度竞价廊坊百度竞价报价丽江关键词排名包年推广报价鹰潭网站建设设计哪家好定西至尊标王哪家好盐田seo网站推广公司商洛建设网站多少钱来宾百度关键词包年推广价格安康关键词按天计费哪家好歼20紧急升空逼退外机英媒称团队夜以继日筹划王妃复出草木蔓发 春山在望成都发生巨响 当地回应60岁老人炒菠菜未焯水致肾病恶化男子涉嫌走私被判11年却一天牢没坐劳斯莱斯右转逼停直行车网传落水者说“没让你救”系谣言广东通报13岁男孩性侵女童不予立案贵州小伙回应在美国卖三蹦子火了淀粉肠小王子日销售额涨超10倍有个姐真把千机伞做出来了近3万元金手镯仅含足金十克呼北高速交通事故已致14人死亡杨洋拄拐现身医院国产伟哥去年销售近13亿男子给前妻转账 现任妻子起诉要回新基金只募集到26元还是员工自购男孩疑遭霸凌 家长讨说法被踢出群充个话费竟沦为间接洗钱工具新的一天从800个哈欠开始单亲妈妈陷入热恋 14岁儿子报警#春分立蛋大挑战#中国投资客涌入日本东京买房两大学生合买彩票中奖一人不认账新加坡主帅:唯一目标击败中国队月嫂回应掌掴婴儿是在赶虫子19岁小伙救下5人后溺亡 多方发声清明节放假3天调休1天张家界的山上“长”满了韩国人?开封王婆为何火了主播靠辱骂母亲走红被批捕封号代拍被何赛飞拿着魔杖追着打阿根廷将发行1万与2万面值的纸币库克现身上海为江西彩礼“减负”的“试婚人”因自嘲式简历走红的教授更新简介殡仪馆花卉高于市场价3倍还重复用网友称在豆瓣酱里吃出老鼠头315晚会后胖东来又人满为患了网友建议重庆地铁不准乘客携带菜筐特朗普谈“凯特王妃P图照”罗斯否认插足凯特王妃婚姻青海通报栏杆断裂小学生跌落住进ICU恒大被罚41.75亿到底怎么缴湖南一县政协主席疑涉刑案被控制茶百道就改标签日期致歉王树国3次鞠躬告别西交大师生张立群任西安交通大学校长杨倩无缘巴黎奥运

深圳SEO优化公司 XML地图 TXT地图 虚拟主机 SEO 网站制作 网站优化