完整版冻土带水合物有效开采方法研究毕业设计.docx
- 文档编号:12062950
- 上传时间:2023-06-04
- 格式:DOCX
- 页数:29
- 大小:3.92MB
完整版冻土带水合物有效开采方法研究毕业设计.docx
《完整版冻土带水合物有效开采方法研究毕业设计.docx》由会员分享,可在线阅读,更多相关《完整版冻土带水合物有效开采方法研究毕业设计.docx(29页珍藏版)》请在冰点文库上搜索。
完整版冻土带水合物有效开采方法研究毕业设计
中国石油大学(华东)毕业设计(论文)
冻土带水合物有效开采方法研究
摘要
天然气水合物作为未来化石燃料的一种重要替代品,以其分布广泛、资源量巨大、能量密度高、清洁无污染等诸多优点为世界各国重视。
而冻土带水合物比海底水合物开发难度小的特点。
本文通过分析冻土带的地质条件,阐述了冻土带与水合物稳定存在的关系,得到了天然气水合物在冻土带存在条件。
分析推测了常见几种冻土带水合物微观结构及其形成过程。
以麦索亚哈气田开采实例,分析了注热降压开采的原理,同时对水力压裂开采水合物压裂的基本原理进行分析,在简化模型下分析了水力压裂增产效果。
关键字:
天然气水合物;冻土带;开采方法
ABSTRAC
Naturalgashydrateisanimportantproductasanalternativetofossilfuelsinthefuture.Withitsadvantageslikewidelydistributed,ahugeamountofresources,highenergydensity,andpollution-free,manycountriesintheworldarehighlyconcernedaboutit.Ratherthanunderseagashydrateformation,gashydrateinthepermafrostareahassmalldevelopingdifficulty.Throughanalysisofthegeologicalconditionsofthepermafrostareas,expoundsitsstabilitybetweenhydrateandpermafrostareas,obtainedtherelationshipofnaturalgashydrateinthepermafrostareascondition.Speculatingthemicrocosmicstructureofhydrateformationprocessofpermafrostareas.AnalyseMaiSuoahabinminingfield,theprincipleofthermalinjectionanddepressurization,atthesametime,Analysethebasicprincipleofhydraulicfracturing,underthesimplifiedmodelisanalyzedeffectonincreasingproductioninhydraulicfracturing.
Keywords:
Gashydrate;permafrostareas;methodofexploitation
第1章前言
1.1研究目的及意义
随着常规油气资源开发利用,人类将面临着巨大能源危机。
目前已经发现的石油储备量还可用40年,天然气还可用70年,煤炭还可用190年,作为世界上最大的发展中大国,我国能源短缺十分突出。
因此急需开发新能源以满足中国经济的高速发展,而天然气水合物作为新型有效能源,2000年开始,天然气水合物的研究与勘探进入高峰期,世界上至少有30多个国家和地区参与其中[1]。
天然气体水合物广泛分布于多年冻土区、大陆架边缘的深海沉积物和深湖泊沉积物中,估计全球天然气水合物中的碳储量为2×1016m3,相当于全球已探明常规化石燃料总碳量的2倍以上[2]。
与海洋相比,陆地开采天然气水合物的难度要小得多,同时陆地上开展天然气水合物开采试验无疑为海洋大规模开发天然气水合物提供良好的基础[3]。
2009年9月,我国在青海发现“可燃冰”,使我国成为世界上第一次在中低纬度冻土区发现天然气水合物的国家,也是继加拿大1992年在北美麦肯齐三角洲、美国2007年在阿拉斯加北坡通过国家计划钻探发现天然气水合物之后,在陆域通过钻探获得天然气水合物样品的第三个国家。
初略的估算,远景资源量至少有350亿吨油当量[4]。
可见冻土带水合物作为一种新的替代能源越来越受到人们的重视,随着我国青藏高原水合物勘探的进步,冻土带水合物的动用必然会对我国油气资源做出巨大贡献。
但要对其进行有效地开采,需要对冻土带地质特征,水合物形成条件,成藏模式进行认识。
同时借助于冻土带水合物开采成功典型例子,麦索雅哈气田的开发,对其开采方式进行分析,以求一种新的开采方法。
1.2研究现状
前苏联二十世纪60年在卖索亚哈油田发现了陆地永久冻土带发现了天然气水合物,1998年,加拿大在北极附近的西北部马更些三角洲(MackenzieDelta)的Mallik地区实施了世界上第一个专门进行陆地永久冻土区天然气水合物调查研究的钻探工程,曾先后完成了5口天然气水合物调查研究探井(Mallik-L38,2L-38,3L-38,3L-3,5L-38),并进行了试开采。
这项世界上第一个对陆地永久冻土区进行天然气水合物全面调查评价的科学研究工程引起了许多国家的兴趣和关注,先后由加拿大、日本、美国、德国、印度、国际大陆钻探计划(ICDP)等国际和机构参与。
除了进行了钻探资源调查评价外,还进行了试开采研究[5]。
2008年3月10~16日,日本和加拿大科学家再次对Mallik2L238井进行了水合物开发测试。
此次开发试验由于安装了防沙装置,开发试验进展顺利,在井底和井口分别采集到了压力、温度、气体和液体流量等数据。
连续6d(139h)的开发,天然气产量达到(2000m3~4000m3)/d,累计产量约为13000m3。
初步数据显示,该区采用水合物降压开采是一种正确的开发方式。
另外,日本国家石油、天然气及金属公司FujiiK女士介绍日本在2009—2011年进行水合物海上开发试验,2012—2016年进行水合物商业开发准备的有观计划[6]。
中国地质调查局于2002年开始先后设立了4个调查研究项目,即“青藏高原多年冻土区天然气水合物地球化学勘查预研究”、“青藏铁路沿线天然气水合物遥感识别标志研究”、“我国陆域永久冻土带天然气水合物资源远景调查”、“我国陆域永久冻土带天然气水合物钻探技术研究”。
国家自然科学基金委员会也于2005年设立了“青藏高原多年冻土区天然气水合物形成条件探讨”的面上科研项目。
国土资源部的科研项目以资源类为主。
目前最重要的当属《我国海域天然气水合物资源调查与评价》专项中的3个配套科研项目。
每年都投入不少资金开展技术方法、环境效应和资源战略评价等方面的研究。
由广州海洋地质调查局、青岛海洋地质研究所、国家海洋局、中国地质科学院、中国科学院、中国地质大学、南京大学等单位承担。
国家自然科学基金委员会自1998年开始设立“气体水合物生成过程机理及其分形理论的研究”(中国科学院广州能源研究所舒碧芬承担)以来。
至2006年共设立了与天然气水合物有关的项目43个。
其中1998年1项。
1999年1项。
2000年4项,2001年6项。
2002年4项,2003年8项。
2004年1项。
2005年10项。
2006年8项。
除2004年外基本上呈逐年增加的趋势。
内容涉及到水合物的物性、资源、抑制、储运、开采、环境等方面的基础研究和应用基础研究。
中国科学院也大力重视天然气水合物的研究工作。
于2004年组织了跨研究所、跨学科的优势研究力量组建了“中国科学院广州天然气水合物研究中心”。
并先后设立了“天然气水合物勘探开采模拟研究”、“天然气水合物开采中若干关键问题的研究”、“大陆坡天然气水合物形成的地质条件与成藏机理研究”、“海底天然气渗漏系统水合物成藏机制及识别方法”等知识创新重要方向项目4个、百人项目3个以及其他项目多个。
研究内容涉及到天然气水合物合成、物性测试、开采模拟、成藏机理、资源评价等领域。
除上述四大部委外,三大油公司、教育部、中国工程院以及前述的中国大洋矿产资源研究开发协会等也先后设立过与天然气水合物有关的科研项目[7]。
最近,有关专家在西宁对中国科学院广州天然气水合物研究中心完成的“青海省天然气水合物勘探开发方案研究报告”进行了审查,指出青海省多年冻土区具备较好的“可燃冰”找矿前景,并建议在青海建立一个陆地冻土区天然气水合物勘探开发技术方法试验基地,以引领带动我国对天然气水合物的勘探开发[8]。
1.3主要研究内容
本课题旨在搞清冻土带地质特点,冻土带厚度,天然气水合物在冻土带形成条件,以及常见的冻土带天然气水合物矿藏模式。
对比常见冻土带天然气水合物开采方法,分析麦索亚哈水合物开采方法,以及找出一种新的适合冻土带水合物开采方法的。
为以后冻土带水合物开采提供依据。
研究内容包括以下几个方面:
通过分析给出冻土带地质特点;研究冻土带厚度;天然气水合物在冻土带稳定存在的条件;分析麦索亚哈气田水合物开采的具体方法;提出一种新的水合物开采方法,以及分析其可行性。
第2章冻土带天然气水合物形成以及存在形式
冻土带在地质学是指0℃以下,并含有冰的各种岩石和土壤。
一般可分为短时冻土(数小时、数日以至半月)、季节冻土(半月至数月)以及多年冻土(数年至数万年以上)。
天然气水合物,在一定地质条件,存在于冻土带。
2.1冻土带水合物的形成
2.1.1冻土带分布情况
世界冻土带分布冻土分布于高纬地带和高山垂直带上部,其中冰沼土广泛分布于北极圈以北的北冰洋沿岸地区,包括欧亚大陆和北美大陆的极北部分和北冰洋的许多岛屿,在这些地区的冰沼土东西延展呈带状分布,在南美洲无冰盖处亦有一些分布。
据估计,冰沼土的总面积约590万平方公里,占陆地总面积的5.5%。
其中北半球20%为冻土面积[9],如图2-1所示北半球冻土带分布。
图2-1北半球多年冻土带水合物
中国多年冻土区面积到215×10km2,占国土总面积的22.4%,是仅次于俄罗斯、加拿大的世界第三冻土大国。
冻土区主要分布在青藏高原、大兴安岭以及其他高山地区。
如图2-2所示为我国冻土地带分布和找矿前景图。
图2-2中国冻土区和天然水合物找矿远景图[10]
图2-3我国青藏高原冻土带分布[11]
2.1.2冻土带厚度
多年冻土厚度是多年冻土区天然气水合物形成的一个非常重要的热力学条件,多年冻土地板附近的地温梯度决定了多年冻土厚度大小。
土体在融化层和冻结层相界面满足能量守恒定律,在多年冻土地板的冻融相界面有:
(2-1)
式中分别是融化层和冻结层之间的相界面的热流,是单位土体的相变热,是相界面移动速率。
当时,说明多年冻土底板附近冻土温度升高,冻土厚度减小。
当时,多年冻土底板附近冻土温度降低,冻土厚度会增大。
如果=,融化层和冻结层相界面不发生热交换,此时相界面移动速率=0。
如果地表温度和地中热流的变化足够小,含水量不大的土(岩)体能够很快地融化和冻结,使得土(岩)体满足稳定状态条件。
此时多年冻土厚度有:
(2-2)
式中,H为多年冻土厚度(m);为年均地温(℃);为冻土(岩)的导热系数;为多年冻土底板附近融土的地温梯度;为多年冻土地板附近融土的导热系数;为地温年变化深度(m),通常取15m;导热系数取高原上土质的平均状态,如=1.31。
图给出了多年冻土厚度计算结果,表明多年冻土底板下融土对多年厚度有较大影响,但多年冻土厚度与冻土低温梯度有密切关系。
图2-4多年冻土厚度计算
综上所述对天然气水合物冻土带厚度决定性影响的客观条件是冻土层的年平均温度和融土层的地温梯度[12]。
长期观测所得数据显示,青藏高原绝大部分位于盆地地区多年冻土,其冻土层年平均温度为-0.5℃—3℃之间,多年冻土厚度介于25-200m之间。
多年冻土区天然气水合物形成,最重要的因素是地热梯度、气体组分、孔隙流体盐度、孔隙压力等,其中地热梯度和气体组分最为重要。
这2个因素决定了多年冻土区天然气水合物形成的温压条件和储藏量。
地热梯度越小,那么天然气水合物可以在更大深度的地质环境条件下形成;气体组分决定了天然气水合物的相平衡状态,比如甲烷水合。
物形成时的温压条件要比丙烷水合物形成条件苛刻得多。
多年冻土区天然气水合物的形成需要充足的烃类气源、地下水、适合的温度和压力4个基本条件。
烃类气源和地下水也许与多年冻土无关;但温度和压力条件与多年冻土有密切的关系,受到了多年冻土热状态、多年冻土厚度以及多年冻土层下地热梯度等的控制和影响。
天然气水合物形成的相平衡条件因为多年冻土区形成的天然气水合物基本为甲烷水合物,所以本文采用由SloanED和KamathVA。
提出的甲烷气水合物形成的相平衡条件统计方程来描述:
(2-3)
(2-4)
式中:
P为压力,kPa;T为温度。
多年冻土底板处温度为0℃,因此,公式(2-3)为多年冻土底板以上天然气水合物形成的相平衡条件;公式(2-4)为多年冻土底板以下沉积层天然气水合物形成的相平衡条件[11]。
冻土层内和冻土层之下的温度与深度可表示为
(2-5)
(2-6)
式中为冻土层内深度为的温度,为地表温度,为冻土层地温梯度(℃/m),为在冻土层之下深度为的温度,为冻土层之下沉积层的地温梯度。
为冻土层底界向下的深度(m),冻土层底界温度为0℃。
冻土层及其下沉积层的压力与深度的关系分别与静岩压力和静水压力计算
(2-7)
(2-8)
式中为地表大气压力(0.1MPa),g为重力加速度常数(取9.81m/s2),为冻土层密度,实验测定冻土密度为1500~2000kg/m3,计算时取1750kg/m3,为冻土层之下孔隙流体密度(1000kg/m3)。
水合物稳定带厚度为:
(2-9)
式中为水合物稳定带厚度,()表示水合物底界埋深,为水合物顶界埋深。
联解方程式获得水合物稳定带底界、顶界埋藏深度和厚度。
图2-5冻土带天然气水合物稳定的深度一温度带[13]
如图2-5给出了多年冻土区天然气水合物稳定带的示意图。
由图2-5可知,多年冻土地温梯度与多年冻土层下地热梯度与气水合物的相平衡边界共同构成了天然气水合物的稳定带,灰色区域内为天然气水合物的稳定带,理论上天然气水合物均可以在这一区域内形成。
多年冻土地温梯度越小,多年冻土厚度越大,温度和压力条件就越有利于形成天然气水合物,其稳定带厚度也越大。
由于多年冻土底板处温度为,所以多年冻土底板深度越大,多年冻土层下天然水合物稳定带的下界深度越大。
天然气水合物层可能的埋深估算回归所得生物成因的天然气水合物温度、压力平衡公式计算显示,最浅的水合物层顶界埋深可能为174m左右,最深的底界埋深有可能到上千米[14]。
2.2冻土带天然气水合物的存在形式
水合物所赋存的沉积物多是新生代沉积。
在沉积层中,水合物以分散状胶结尚未固结的沉积物颗粒,或者是以结核状、团块状、薄层状的集合体赋存于沉积物中,还可能以细脉状、网脉状充填于裂隙中。
冻土带水合物与常规油气藏关系密切。
包括汇聚和离散两种类型的沉积盆地,是比较有利的水合物形成区。
对水合物气藏形成,前人虽曾提出冷冻、渗流等模式,但对构造动力学、古气候演化、热史变异、气体运移、深部气源等问题研究不足,因而远落后于对常规油气藏的认识[15]。
2.2.1青藏高原地质运动分析
青藏高原多年冻土区主要分布在昆仑山北坡至唐古拉山南坡安多县北的广大地域,沿着青藏铁路线横跨约550km长。
其内分布着一系列中新生代残留盆地,如可可西里盆地、沱沱河盆地、羌塘盆地、唐古拉盆地等,它们特别是羌塘盆地,分布面积广、沉积厚度大、有机质含量丰富,且已发现有大量的油气显示或工业油流,这些条件可为研究区天然气水合物成藏提供良好物质基础。
上述盆地是在伴随着青藏高原多期脉动性急剧而又整体隆升过程中形成的构造残留盆地,在高原隆升过程中由于不同的构造演化,后期改造作用一方面形成相应的构造圈闭,使油气相对聚集成藏,另一方面构造运动太强烈,抬升和切割幅度太大,对油气的保存起着破坏作用,使大量的轻烃难以保存而逃逸。
这种构造破坏作用造成的烃类气体重新分配现象在上新世一更新世以来尤为明显。
此间,在上新世至更新世期间,深海有孔虫氧的同位素记录显示全球总体上为强烈的降温时期;在更新世末期至全新世期间(约在0.71Ma),青藏高原也整体进入冰冻期,形成了大面积的冻土,出现了大规模的山地冰川,除夏季以外,高原上稳定的积雪及大规模的冰川提高了地表反射率,增强了冬季高原上空的冷高压,促使高原进一步变冷,其中,晚更新世末期是末次冰期最盛期,青藏高原处于严寒、干旱的冰缘环境,在这种干冷的气候下,青藏高原开始形成大面积的多年冻土。
进入全新世,虽然大规模的冰期已经结束,但由于整个高原仍处于不断隆升状态,且高原海拔已上升到4000m以上,特殊的海拔高度决定着整个高原仍处在冰缘环境下。
上新世至更新世期间的降温环境,特别是更新世末期至全新世期间的冰川环境为同期构造作用破坏的油气藏烃类气体进一步被圈闭提供了一种耦合机制,即游离的烃类气体在低温(冻土层)环境和地层压力条件作用下重新被圈闭形成天然气水合物藏[16]。
2.2.2几种常见成矿模式分析
冻土带水合物形成之前,一般已经有天然气藏,随着温度的变化,达到水合物形成条件以下,就会形成水合物,并且由于陆地上的温度的循环交替变化,水合物也在不断形成和分解,所以水合物矿藏也伴随有天然气藏[17]。
麦索亚哈气田就是有天然气藏水合物矿藏。
目前冻土带水合物地质模式主要以下几种模式块状水合物、层状或片状水合物、分散状水合物、透镜体水合物。
块状水合物形成之前有一个较完好的天然气气藏,底部有充足的水源,在地壳急剧上升过程中,由于构造运动使原气藏遭到破坏,当达到水合物稳定性成所需的温度压力,并且破坏的气藏又有新的的圈闭形成,此时水与天然气将形成块状的天然气水合物。
如图所示:
图2-6冻土带水合物块状水合物形成
在地壳运动的过程中伴随着沉积,由于下部天然气渗流到上部,和岩石进入上部岩石的空隙中,如果上部隆升和沉积持续进行,且此时温度压力满满足水合物形成温度压力条件时,形成片状或层状水合物。
如图所示:
图2-7层状或片状水合物形成
有完整的层状或块状水合物形成是,由于地质构造运动,褶皱活断层剧烈时,水合物矿藏将受到破坏,水合物矿藏发生断裂破坏,此时温度压力改变水合物将分解断裂,当有新的适宜压力和温度条件时,水合物不再分解,重新稳定,形成分散状的各种形态不一的水合物结构。
图2-8分散状天然气水合物形成
透镜体水合物是是低温冰冻模式形成的水合物,它是有地下已存在的气体矿床或饱和气体的水在上升过程中受到冷冻发生相变而形成中间厚扩展向两边透镜体状的水合物。
如图所示,开始在岩石裂缝或大的颗粒空隙中存在天然气和水等流体,在温度压力达到水合物形成条件后形成透镜体状水合物。
图2-9透镜体天然气水合物形成
近日在祁连山取出水合物岩心结构如图2-10所示,通过取心表明发现水合物赋存状态复杂,研究认为祁连山天然气水合物的产出方式可能有二种,一是以薄层状、片状、团块状赋存于粉砂岩、泥岩、油页岩的裂隙面中,单层仅仅数毫米厚,肉眼能观测到的均是这类水合物。
一是以浸染状赋存于细粉砂岩的孔隙中,这类水合物肉眼难以观测到,但可通过岩心中不断冒出的气泡和水珠来间接推测存在这类水合物,也可通过红外测温中的分散状低温异常证实有这类水合物。
天然气水合物分解后将导致岩心温度降低,从A-B-C-D反映水合物含量逐渐升高[18]。
图2-10祁连山冻土区发现的天然气水合物及其燃烧照片
图2-11祁连山冻土区含天然气水合物岩心的红外热像仪照片
第3章冻土带天然气水合物注热降压开采
目前天然气水合物的开采方法常见有三种,第一种注热开采,第二种降压开采,第三种是注入化学剂[19]。
其中注热和降压开采用于麦索亚哈气田开采所取得较好的效果。
3.1降压开采原理
3.1.1降压开采机理
通过降低天然气水合物储集层的压力促使水合物脱离其稳定带,从而引起水合物的相平衡曲线的移动,使水合物分解。
为降低水合物储集层的压力,可以开采水合物储集层之下的游离天然气或是采用热激发或注入化学试剂的方法,从而使水合物不稳定而分解。
此外,在钻井过程中采用低密度钻井液也可以实现水合物储层的降压,而通过调节天然气的开采速度则可以有效控制储层压力,从而控制水合物的分解。
图3-1天然气水合物降压开采原理
不需要昂贵的连续激发是降压法的最大特点,这也使得运用该方法经济成本较低,简单易行,适合大面积开采,尤其适用于存在下伏游离气层的天然气水合物藏的开采,因而其可能成为未来水合物商业开发的有效方法之一。
降压法的缺陷在于,由于没有额外的热量注入,水合物分解只能从周围环境中吸收热量,当吸收的热量达到一定程度,周围环境温度的降低会抑制水合物的进一步分解。
因此,它对天然气水合物藏的性质有一定要求,只有当天然气水合物藏的原始温度、压力条件位于温压平衡边界附近时,降压开采法才具有较好的经济可行性。
3.1.2降压开采能量分析
水合物开采的主要能量驱动有以下几种:
天然气的膨胀能,对于底部有天然气藏的水合物层,开发前期,底部的天然气膨胀能实现天然气的能量驱动。
水合物分解后的天然气驱动,当储层压力或温度不满住水合物稳定存在条件时候,水合物开始分解,此时分解的天然气继续为地层提供能量。
其他方式的能量,在降压开采过程中往往可以结合其他的方式,使得水合物不断分解。
往往很多开发方式都是伴随着降压过程进行。
3.2注热开采原理
3.2.1注热开采机理
注热开采是利用钻探技术在天然气水合物储层中安装管道,对天然气水合物的地层进行加热,提高局部储集层温度,打破天然气水合物稳定存在的条件,从而造成天然气水合物的分解。
注热开采通常采用的蒸汽注入、热水注入、热盐水注入、火驱及电磁加热等方法。
图3-2注热开采水合物示意图
不同的注热开采方式也有不同的优缺点,蒸汽注入和火驱技术方法在开采比较薄的水合物储层时的热量损失很大,只有在储层厚度大于15m时进行开采的热量利用率才比较高。
注入热水的热损失较蒸汽和火驱小,但水合物储层内的水对于注入率又限制了该方法的应用。
热水注入法与蒸汽注入法相比,前者由于比后者消耗能量少,也不需要考虑稀释问题,从而可能有效的降低开发费用。
同时多井注热比单井注热效果在促进水合物分解上效果又更好。
注热开采法的主要不足是会造成大量的热损失,效率很低,特别是在永久冻土区,即使利用绝热管道,永冻层也会降低传递给储集层的有效热量。
图3-3水合物层单井注热与多井注热开采效果比较[20]
3.2.2注热开采能量分析
水合物分解所需要的能量主要分为五部分:
水合物分解能量、岩石骨架升温所需能量、水合物升温所需要的能量、地层中束缚水升温所需要的能量和水合物分解后分解水和天然气升温所需要的能量。
不同的水合物藏的物性参数不同,上述各部分能量所占的比例有所不同。
在理想状态下,假设岩石孔隙度为30%;水合物为甲烷水合物,水合系数为6且饱和度为80%;原始地层温度为4℃,水合物开始分解温度为12℃,地层最终温度为15℃;整个过程保持地层压力不变。
如图3-3所示,岩石升温所需要的热量,占72.04%;束缚水升温吸收的热量,占8.99%;表示水合物分解前升温吸收热量,占11.73%;水合物分解能量,占0.74%;水合物分解后分解水升温吸热,6.5%。
图3-3水合物分解各部分吸热示意图计算程序
特例当原始地层温度、水合物开始分解温度和水合物最终地层温度一致的情况下,水合物分解能量在能量比例中占到了100%。
图3-4在特例情况下水合物分解能量比例极值
可以看出,注热开采水合物的大部分能量都被非水合物升温所消耗。
但对于水合物的分解来讲,只有水合
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 完整版 冻土 水合物 有效 开采 方法 研究 毕业设计