正定天气预报-正定天气预报一周天气
石家庄并非热门旅游目的地,十一人少的旅游地:
1. 赵州桥;
赵州桥坐落在河北省赵县洨河上。建于隋代(公元581-618年)大业年间(公元605-618年),由著名匠师李春设计和建造,距今已有约1400年的历史,是当今世界上现存最早、保存最完善的古代敞肩石拱桥;
2. 苍岩山;
苍岩山所在地方圆几十里,皆为荒山秃岭,唯独此山林木苍翠,古树葱笼,峭壁岩,处处绿色,故名“苍岩”;
苍岩山位于河北省石家庄市西南50公里,地处井陉县境内,总面积63平方公里。为中国历史文化名山、国家重点风景名胜区、国家4A级旅游区、苍岩山福庆寺做为苍岩山的核心景区被列为中国重点文物保护单位;
3. 天桂山;
天桂山,因山中多奇峰异洞,自然风光颇具桂林山水之神,人称"天桂山";
位于河北省会石家庄市西北90公里处的平山县境内,总面积60平方公里,是河北省著名的山岳古刹型风景名胜区;
4. 抱犊寨;
抱犊寨,旧名抱犊山,古名萆山,位于河北省石家庄鹿泉市西郊,距省会石家庄16公里,距首都北京288公里,是一外集历史人文和自然风光为一体的名山古寨。
5. 西柏坡纪念馆;
西柏坡纪念馆,位于河北省平山县西柏坡村;
1958年修建岗南水库时,西柏坡旧址被水库淹没。1971年在离水库不远的同一座山上兴建了新的西柏坡村,大体保持原貌。1977年在旧址西南又新建了西柏坡纪念馆。
6. 嶂石岩;
嶂石岩为位于石家庄西南的赞皇县境内(到赞皇县城52公里),距河北省省会石家庄市区约110公里,是太行山森林公园精华所在;
此景点内峰峦叠嶂,奇峰突兀,怪石林立,幽洞深深,泉瀑淙淙,险峻中叠现奇伟,幽清中不失秀美,是一处自然风光绝佳的旅游圣地,称得起一方绝胜。
描写过年时情景的句子
1、柘城县天气预报2、商丘市有几个县3、中国地名有哪些生僻字4、商丘电视台今天晚上天气预报5、墨迹天气的支持城市柘城县天气预报
柘城天气实况
今天
2016年6月22日
周三晴转多云
38~27℃
无持续风向微风15m/s
降雨量:0mm日出日落
今天:05:07|19:33
明天:05:08|19:34
降水概率:54%体感温度:41℃
湿度:86%
气压:1003hPa
能见度:7km
商丘市有几个县
截止2017年,商丘市下辖2个市辖区,1个县级市,6个县。
2个市辖区:睢阳区、梁园区
1个县级市:永城市
6个县:民权县、睢县、宁陵县、柘城县、虞城县、夏邑县
扩展资料:
商丘,简称“商”或“宋”,河南省地级市。位于豫鲁苏皖四省辐凑之地,素有"豫东门户"之称。商丘自古一直是兵家必争之地和商贾云集中心,隋唐大运河横贯全境,陇海铁路与京九铁路的普铁高铁、连霍高速与济广高速均在此交会,是规划建设的全国性综合交通枢纽。
商丘是国家历史文化名城,是华夏文明和中华民族的重要发祥地、中国重要古都。燧人氏、炎帝朱襄氏、颛顼、帝喾等三皇五帝以及夏朝、商朝、周朝宋国、汉晋梁国等均在此建都。商丘享有“三商之源·华商之都”、“汉兴之地”、“两宋龙潜之地”之称。
商丘自古为长寿之乡,是国家森林城市,有5处国家湿地公园。商丘位居全国18个百亿斤产粮大市,被誉为“豫东粮仓”。地下煤炭资源储量居全省之最。
参考资料:
中国地名有哪些生僻字
古语云:北不识盱眙,南不识__。
(或:西不识盱眙,东不识__。
说的是自古以来人们
就困扰于地名中的生僻字,
但由于地广人多地名就多,
一旦把这些生僻地名取消,
又增加了
重名的麻烦,而且也会抹杀了很多历史遗产。
作为地理爱好者,
我从小就一直遭受着地名的折磨,
小学时不认识邯郸,
觉得这地方太格路
了,难怪别人学不会他们走步。还有山东莒
jv3
县,读《上下五千年》齐桓公那一段,常读
作:齐桓
heng2
公,莒
lv3
国。中学受到《小石潭记》影响,一直把天津宝坻
di3
读作
chi2
另外对于浙江台
tai1
州也是偶然听到天气预报才知道正确发音。
这里简单罗列一下中国县级以上的难念地名,
很多是我查了又忘,
忘了又查的,
也会有一些
错误,希望大家给予指正:
安徽:砀
dang4
山,歙
she4
县,黟
yi1
县,谯
qiao2
城(南谯同
,枞
zong1
阳(多音字:
枞
cong1
树
,琅琊
ya2
(来源于滁州附近的琅琊山,不过古代山东东南的琅琊山更有名,
还有一个琅琊郡,容易和这个混淆
福建:柘
zhe4
荣(河南柘城同
甘肃:宕
dang4
昌,舟曲
qv1
(多音字,还有很多藏语地名带这个字
广东:番
pan1
禺(其实广东还有一些生僻字地名如:深圳,东莞,都已经为人熟知
海南:儋
dan1
州
河北:藁
gao3
城,井陉
xing2
,蠡
li3
县,蔚
yu4
县,任
ren2
丘(多音字,任县同,不同于
山东济宁任
ren4
城区
黑龙江:阿
a1
城(不同于山东东阿
e1
,讷
ne4
河
河南:_
chan2
河,郏
jia2
县,武陟
zhi4
,浚
xun4
县(多音字,
jun4
,召
shao4
陵(不同
于南召
zhao4
,泌
bi4
阳(多音字
,济
ji3
源(古时与山东济南一水相连
湖北:
_
qiao2
口,
黄陂
pi2
郧
yun2
县
(郧西同
蕲
qi2
春,
_
xiao1
亭,
曾
zeng1
都
(多
音字
湖南:醴
li3
陵,郴
chen1
州(初中学
“
郴江幸自绕郴山
”
时才会读
,酃
ling2
县(
94
年改为
炎陵县
江苏:
盱
xv1
眙
yi2
六
lu4
合
(安徽六安同
邗
han2
江
(出自邗沟,
隋代大运河的一段
江西:鄱
po
阳(鄱阳波阳改来改去,不过这个对于多数人应该不成问题
,铅
yan2
山
吉林:绿
lu4
园,珲
hun2
春(因为是长春人,这两个都没问题,但见过其他人读错
辽宁:岫
xiu4
岩,鲅
ba4
鱼圈
陕西:__(
56
年改周至
,_县(
64
年改户县
,柞
zha4
水,华
hua4
县(其他地名一般
为
hua2
,如:华宁,华坪
山东:芝罘
fu2
,临朐
qv2
,汶
wen4
上
(四川汶川同
,郯
tan2
城,茌
chi2
平,鄄
juan4
城,
郓
yun4
城(读过水浒的人都知道
,单
shan4
县(多音字
,牟
mu4
平
商丘电视台今天晚上天气预报
商丘天气预报
(未来15天
切换城市
相关地区:睢县天气民权天气虞城天气柘城天气宁陵天气永城天气
更多
10月03日(周一休
晴转多云17~29℃
东南风1级
空气良
10月04日(周二休
晴转多云18~28℃
无持续风向微风
空气良
10月05日(周三休
多云转阴16~26℃
东北风3-4级
空气良
10月06日(周四休
多云转阴15~21℃
无持续风向微风
空气良
10月07日(周五休
多云转阴15~24℃
无持续风向微风
空气良
10月08日(周六班
多云14~23℃
无持续风向微风
空气良
10月09日(周日班
多云11~23℃
东北风3-4级
空气良
墨迹天气的支持城市
直辖市
北京北京朝阳通州昌平顺义怀柔大兴平谷密云延庆丰台*石景山*房山*门头沟*八达岭*佛爷顶*汤河口*密云上甸子*斋堂*霞云岭*海淀
天津天津宝坻武清宁河静海蓟县*东丽*西青*北辰汉沽*津南塘沽*大港
上海上海嘉定浦东金山松江青浦南汇奉贤崇明*宝山*徐家汇*闵行
重庆重庆黔江北碚巴南长寿合川永川江津南川綦江潼南铜梁大足荣昌璧山垫江武隆丰都城口梁平开县巫溪巫山奉节云阳忠县石柱秀山酉阳彭水*万盛*渝北*涪陵*沙坪坝*万州天城*万州龙宝
其他省份及城市
辽宁:
沈阳沈阳苏家屯新民辽中法库康平*新城子
大连大连旅顺金州瓦房店普兰店庄河长海*皮口
鞍山鞍山海城台安岫岩
抚顺抚顺清原*章党
本溪本溪本溪县桓仁*草河口
丹东丹东东港凤城宽甸*东沟
锦州锦州凌海北镇黑山义县*北宁
营口营口盖州大石桥
阜新阜新彰武
辽阳辽阳灯塔辽阳县
盘锦盘锦大洼盘山
铁岭铁岭开原西丰昌图
朝阳朝阳北票凌源建平喀左*羊山*建平县
葫芦岛葫芦岛兴城绥中建昌
吉林:
长春长春双阳九台榆树德惠农安
吉林吉林舒兰桦甸蛟河磐石永吉*烟筒山
四平四平公主岭双辽梨树伊通*孤家子
辽源辽源东丰
通化通化梅河口集安辉南通化县柳河
白山白山靖宇长白*东岗临江
白城白城洮南大安镇赉通榆
松原松原扶余长岭乾安前郭
延边延边延吉图们敦化龙井珲春和龙安图汪清*二道*松江*罗子沟
黑龙江:
哈尔滨哈尔滨呼兰双城尚志五常阿城依兰方正宾县巴彦木兰通河延寿
齐齐哈尔齐齐哈尔讷河龙江依安泰来甘南富裕克山克东拜泉
鸡西鸡西虎林密山鸡东
鹤岗鹤岗萝北绥滨
双鸭山双鸭山集贤宝清饶河
大庆大庆肇源肇州林甸*杜蒙
伊春伊春铁力嘉荫*乌伊岭*五营
佳木斯佳木斯同江富锦桦南桦川汤原抚远
七台河七台河勃利
牡丹江牡丹江穆棱绥芬河海林宁安东宁林口
黑河黑河北安五大连池嫩江逊克孙吴
绥化绥化安达肇东海伦望奎兰西青冈庆安明水绥棱
大兴安岭大兴安岭加格达奇呼玛塔河漠河*呼中*新林
河南:
郑州郑州新郑荥阳新密巩义登封中牟
开封开封杞县通许尉氏兰考
洛阳洛阳偃师孟津新安栾川嵩县汝阳宜阳洛宁伊川
平顶山平顶山舞钢汝州宝丰叶县鲁山郏县
焦作焦作沁阳孟州修武博爱武陟温县
鹤壁鹤壁浚县淇县
新乡新乡卫辉辉县获嘉原阳延津长垣封丘
安阳安阳林州汤阴内黄滑县
濮阳濮阳清丰南乐范县台前
许昌许昌禹州长葛鄢陵襄城
漯河漯河舞阳临颍
三门峡三门峡灵宝渑池卢氏
南阳南阳邓州南召方城镇平社旗唐河新野内乡西峡淅川桐柏
商丘商丘永城民权睢县宁陵柘城虞城夏邑
信阳信阳罗山光山新县商城固始潢川淮滨息县
周口周口项城扶沟西华商水沈丘郸城淮阳太康鹿邑
驻马店驻马店西平上蔡平舆正阳确山泌阳汝南遂平新蔡
济源济源
湖北:
武汉武汉黄陂新洲*江夏*蔡甸
黄石黄石大冶阳新
襄樊襄樊襄阳枣阳老河口宜城谷城保康南漳
十堰十堰郧县郧西竹山竹溪房县*丹江口
荆州石首洪湖松滋公安监利江陵荆州
宜昌宜昌宜都当阳枝江远安兴山秭归长阳五峰三峡*宜昌县*夷陵
荆门荆门京山钟祥
孝感孝感汉川安陆应城大悟云梦
黄冈黄冈麻城武穴红安罗田英山浠水蕲春黄梅
咸宁咸宁赤壁嘉鱼通城崇阳通山
随州随州广水
恩施恩施利川建始巴东宣恩咸丰来凤鹤峰
仙桃仙桃
潜江潜江
天门天门
神农架神农架
鄂州鄂州
湖南:
长沙长沙浏阳宁乡*马坡岭
株洲株洲醴陵攸县茶陵炎陵
湘潭湘潭湘乡*韶山
衡阳衡阳南岳耒阳常宁衡阳县衡南衡山衡东祁东
邵阳邵阳武冈邵东新邵邵阳县隆回洞口绥宁新宁城步
岳阳岳阳汨罗临湘华容湘阴平江
常德常德安乡汉寿临澧澧县桃源石门
张家界张家界慈利桑植
益阳益阳沅江南县桃江安化*赫山区
郴州郴州资兴桂阳宜章永兴嘉禾临武汝城桂东安仁
永州永州祁阳东安双牌道县江永宁远蓝山新田江华
怀化怀化沅陵辰溪溆浦会同麻阳新晃芷江靖州通道
娄底娄底冷水江涟源双峰新化*冷水滩
吉首吉首泸溪凤凰花垣保靖古丈永顺龙山
广西:
南宁南宁邕宁武鸣隆安马山上林宾阳横县*南宁城区
柳州柳州柳江柳城鹿寨融安融水三江
桂林桂林阳朔临桂灵川平乐兴安灌阳荔浦资源永福龙胜恭城*全州
梧州梧州岑溪苍梧藤县蒙山
北海北海合浦*涠洲岛
防城港防城上思*防城港*东兴
贺州贺县钟山昭平富川
来宾来宾象州武宣忻城金秀
贵港贵港桂平平南
玉林玉林北流容县陆川博白兴业
百色百色凌云平果西林乐业德保田林田阳靖西田东那坡隆林
河池河池宜州天峨凤山南丹东兰都安罗城巴马环江
钦州钦州灵山浦北
崇左崇左宁明扶绥龙州大新天等*凭祥
广东:
广州广州番禺花都增城从化
深圳深圳
珠海珠海斗门*黄茅洲
汕头汕头潮阳澄海南澳
韶关韶关曲江乐昌南雄始兴仁化翁源新丰乳源
佛山佛山顺德三水*南海
江门江门新会台山开平鹤山恩平上川岛
茂名茂名高州化州信宜电白
湛江湛江廉江雷州吴川遂溪徐闻
肇庆肇庆高要四会广宁怀集封开德庆
惠州惠州惠阳博罗惠东龙门
梅州梅州兴宁梅江区梅县区大埔丰顺五华平远蕉岭
汕尾汕尾陆丰海丰
河源河源紫金龙川连平和平
阳江阳江阳春
清远清远英德连州佛冈阳山连山连南
东莞东莞
中山中山
潮州潮州饶平
揭阳揭阳普宁揭西惠来
云浮云浮罗定新兴郁南
海南:
海口海口
琼山琼山
三亚三亚
文昌文昌
琼海琼海
万宁万宁
东方东方
五指山五指山
儋州儋州
临高临高
澄迈澄迈
定安定安
屯昌屯昌
昌江昌江
白沙白沙
琼中琼中
陵水陵水
保亭保亭
乐东乐东
西沙*西沙
南沙岛南沙岛
江苏:
南京南京浦口江宁六合溧水高淳*江浦
无锡无锡江阴宜兴
徐州徐州新沂邳州丰县沛县睢宁
常州常州金坛溧阳
苏州苏州常熟张家港昆山吴江太仓*吴县东山*吴县
南通南通启东如皋通州海门海安如东吕泗
连云港连云港赣榆东海灌云灌南*燕尾港*西连岛
淮安淮安楚州淮阴涟水洪泽盱眙金湖*淮阴县
盐城盐城盐都东台大丰响水滨海阜宁射阳建湖
扬州扬州仪征高邮江都宝应
镇江镇江丹徒丹阳扬中句容
泰州泰州兴化靖江泰兴姜堰
宿迁宿迁沭阳泗阳泗洪
浙江:
杭州杭州萧山余杭建德富阳临安桐庐淳安
台州台州椒江黄岩路桥临海温岭玉环天台仙居*三门*大陈*洪家
宁波宁波鄞州镇海余姚慈溪奉化宁海象山*北仑石浦*鄞县
温州温州瑞安乐清永嘉文成平阳泰顺洞头苍南
嘉兴嘉兴海宁平湖桐乡嘉善海盐
湖州湖州德清安吉长兴
绍兴绍兴诸暨上虞嵊州新昌
金华金华兰溪义乌东阳永康武义浦江磐安
衢州衢州江山常山开化龙游
舟山舟山定海普陀岱山嵊泗
丽水丽水龙泉青田缙云遂昌松阳云和庆元*景宁
江西:
南昌南昌南昌县新建安义进贤*莲塘
景德镇乐平景德镇
萍乡萍乡莲花
九江九江武宁修水永修德安星子都昌湖口彭泽庐山*瑞昌
新余新余分宜
鹰潭鹰潭贵溪余江
赣州赣州瑞金南康信丰大余上犹崇义安远龙南定南全南宁都于都兴国会昌寻乌石城
吉安吉安吉水峡江新干永丰遂川泰和万安安福永新*宁冈*井冈山*吉安县
宜春宜春丰城樟树高安奉新万载上高宜丰靖安铜鼓
抚州南城黎川南丰崇仁乐安宜黄金溪资溪东乡广昌*抚州
上饶上饶德兴上饶县广丰玉山铅山横峰弋阳余干波阳万年婺源*鄱阳
山东:
济南济南长清章丘平阴济阳商河
青岛青岛胶州即墨平度胶南莱西*崂山
淄博淄博周村博山淄川桓台高青沂源*临淄
枣庄枣庄薛城峄城台儿庄滕州
东营东营利津垦利广饶*河口
烟台烟台福山牟平龙口莱阳莱州蓬莱招远栖霞海阳长岛
潍坊潍坊青州诸城寿光安丘高密昌邑昌乐临朐
威海威海文登乳山荣成*石岛成山头
济宁济宁曲阜兖州邹城鱼台金乡嘉祥微山汶上泗水梁山
泰安泰安新泰肥城宁阳东平泰山
日照日照五莲莒县
莱芜莱芜
临沂临沂沂南郯城沂水苍山费县平邑莒南蒙阴临沭
德州德州乐陵禹城陵县宁津庆云临邑齐河平原夏津武城
聊城聊城临清阳谷莘县茌平东阿冠县高唐*朝城
滨州滨州惠民阳信无棣沾化博兴邹平
菏泽菏泽曹县单县成武巨野郓城鄄城定陶东明
安徽:
合肥合肥长丰肥西肥东
芜湖芜湖芜湖县南陵繁昌
蚌埠蚌埠怀远固镇五河
淮南淮南凤台
马鞍山马鞍山当涂
淮北淮北濉溪
铜陵铜陵
安庆安庆桐城宿松枞阳太湖怀宁岳西望江潜山
黄山黄山市*黄山区*黄山(景区屯溪休宁歙县祁门黟县
滁州滁州天长明光全椒来安定远凤阳
阜阳阜阳界首临泉颍上阜南太和
宿州宿州萧县泗县砀山灵璧
巢湖巢湖*和县*含山*庐江*无为
六安六安寿县霍山霍邱舒城金寨
亳州亳州利辛涡阳蒙城
池州池州东至石台青阳*九华山
宣城宣城宁国广德郎溪泾县旌德绩溪
福建:
福州福州福清长乐闽侯连江罗源闽清永泰平潭*福州郊区
厦门厦门同安
莆田莆田仙游*秀屿港
泉州泉州晋江南安安溪永春德化*崇武九仙山
漳州漳州龙海云霄漳浦诏安长泰东山南靖平和华安
南平南平邵武武夷山建瓯建阳顺昌浦城光泽松溪政和
龙岩龙岩漳平武平长汀永定连城上杭
三明三明永安明溪清流宁化大田尤溪沙县将乐泰宁建宁
宁德宁德福安福鼎霞浦古田屏南寿宁周宁柘荣
河北:
石家庄石家庄藁城辛集晋州新乐平山井陉栾城正定行唐灵寿高邑赵县赞皇深泽无极元氏
唐山唐山丰润丰南遵化迁安迁西滦南玉田唐海乐亭滦县
秦皇岛秦皇岛昌黎卢龙抚宁青龙*北戴河
邯郸邯郸永年曲周馆陶魏县成安大名涉县鸡泽邱县广平肥乡临漳磁县*武安*峰峰
邢台邢台南宫沙河柏乡任县清河宁晋威县隆尧临城广宗临西内邱平乡巨鹿新河南和
保定保定涿州定州安国高碑店满城阜平徐水唐县高阳容城涞源望都安新易县曲阳蠡县雄县*顺平
张家口张家口宣化康保张北阳原赤城沽源怀安怀来崇礼尚义蔚县涿鹿万全
承德承德承德县兴隆隆化平泉滦平丰宁围场宽城
沧州沧州泊头任丘黄骅河间青县献县东光海兴盐山肃宁南皮吴桥孟村*曹妃甸
廊坊廊坊霸州三河固安永清香河大城文安大厂
衡水衡水冀州深州饶阳枣强故城阜城安平武邑景县武强
山西:
太原太原太原古交区阳曲清徐娄烦*太原南郊*太原北郊
大同大同大同县天镇灵邱阳高左云广灵浑源
阳泉阳泉平定盂县
长治长治襄垣屯留平顺黎城壶关长子武乡沁县沁源*潞城
晋城晋城高平沁水陵川阳城
朔州朔州山阴应县右玉怀仁*平鲁
晋中晋中榆次介休榆社左权和顺昔阳寿阳太谷祁县平遥灵石
运城运城永济河津芮城临猗万荣新绛稷山闻喜夏县绛县平陆垣曲
忻州忻州原平定襄五台山代县繁峙宁武静乐神池五寨岢岚河曲保德偏关*五台县豆村
临汾临汾侯马霍州曲沃翼城襄汾洪洞古县安泽浮山吉县乡宁蒲县大宁永和隰县汾西
吕梁吕梁离石孝义汾阳文水中阳兴县临县方山柳林岚县交口交城石楼
内蒙古:
呼和浩特呼和浩特托克托武川和林格尔清水河土默特左旗*呼和浩特市郊区
包头包头固阳达尔罕茂明安联合旗土默特右旗*白云鄂博满都拉
乌海乌海
赤峰赤峰宁城林西喀喇沁旗巴林右旗敖汉旗阿鲁科尔沁旗翁牛特旗克什克腾旗巴林左旗*浩尔吐*宝过图*岗子*八里罕
通辽通辽开鲁科尔沁左翼中旗科尔沁左翼后旗库伦旗奈曼旗扎鲁特旗*青龙山*巴雅尔吐胡硕*高力板*舍伯吐
鄂尔多斯鄂尔多斯东胜准格尔旗乌审旗伊金霍洛旗鄂托克旗鄂托克前旗杭锦旗达拉特旗河南*乌审召*伊克乌素
呼伦贝尔呼伦贝尔海拉尔满洲里牙克石扎兰屯根河额尔古纳陈巴尔虎旗阿荣旗新巴尔虎左旗新巴尔虎右旗鄂伦春旗莫力达瓦旗鄂温克旗小二沟图里河
巴彦淖尔巴彦淖尔五原磴口杭锦后旗乌拉特中旗乌拉特前旗乌拉特后旗*海力素那仁宝力格*大佘太
乌兰察布乌兰察布丰镇兴和卓资商都凉城化德四子王旗察哈尔右翼前旗察哈尔右翼中旗察哈尔右翼后旗
锡林郭勒盟锡林浩特二连浩特多伦阿巴嘎旗西乌珠穆沁旗东乌珠穆沁旗苏尼特左旗苏尼特右旗太仆寺旗正镶白旗正兰旗镶黄旗朱日和博克图*乌拉盖
兴安盟乌兰浩特阿尔山突泉扎赉特旗科尔沁右翼中旗*霍林郭勒*索伦*胡尔勒
阿拉善盟阿拉善左旗阿拉善右旗额济纳旗拐子湖吉兰太*孪井滩*头道湖*中泉子*巴彦诺尔贡*雅布赖*乌斯太*锡林高勒
陕西:
西安西安临潼长安高陵蓝田户县周至*杨凌
铜川铜川耀县宜君
宝鸡宝鸡陈仓凤翔岐山扶风眉县陇县千阳麟游凤县太白*宝鸡县
咸阳咸阳兴平三原泾阳乾县礼泉永寿彬县长武旬邑淳化武功*杨凌
渭南渭南华阴韩城华县潼关大荔蒲城澄城白水合阳富平华山
延安延安延长延川子长安塞志丹吴起甘泉富县洛川宜川黄龙黄陵
汉中汉中南郑城固洋县西乡勉县宁强略阳镇巴留坝佛坪
榆林榆林神木府谷横山靖边定边绥德米脂佳县吴堡清涧子洲
安康安康汉阴石泉宁陕紫阳岚皋平利镇坪旬阳白河
商洛商洛洛南丹凤商南山阳镇安柞水
宁夏:
银川银川灵武永宁贺兰
石嘴山石嘴山大武口惠农平罗*陶乐石炭井
吴忠吴忠青铜峡同心盐池*韦州*麻黄山
固原固原西吉彭阳泾源隆德*六盘山
中卫中卫海原中宁*兴仁堡
甘肃:
兰州兰州永登榆中皋兰
金昌金昌永昌
白银白银靖远景泰会宁华家岭
天水天水武山甘谷清水秦安张家川*麦积*北道区
嘉峪关嘉峪关
武威武威民勤古浪天祝乌鞘岭
张掖张掖民乐山丹临泽高台肃南
平凉平凉灵台静宁崇信华亭泾川庄浪*崆峒
酒泉酒泉玉门敦煌金塔肃北*瓜州*鼎新马鬃山
庆阳庆阳西峰庆城镇原合水华池环县宁县正宁
定西定西安定通渭临洮漳县岷县渭源陇西
武都武都成县宕昌康县文县西和礼县两当徽县
临夏临夏康乐永靖广河和政东乡
合作合作临潭卓尼舟曲迭部玛曲碌曲夏河
青海:
西宁西宁湟源湟中大通
海东*海东平安乐都民和互助化隆循化冷湖
海北*海北海晏祁连刚察门源*托勒
黄南*黄南尖扎泽库河南
海南*海南共和同德贵德兴海贵南*海南
果洛*果洛玛沁班玛甘德达日久治玛多*清水河
玉树玉树杂多治多囊谦曲麻莱托托河
海西*海西格尔木德令哈乌兰都兰天峻*野牛沟*五道梁*小灶火*诺木洪茫崖大柴旦*茶卡
新疆:
乌鲁木齐乌鲁木齐*小渠子巴仑台*白杨沟*十三间房气象站*天山大西沟*乌鲁木齐牧试站*天池*蔡家湖*达坂城
克拉玛依克拉玛依
石河子*石河子*莫索湾*乌兰乌苏*炮台
阿拉尔阿拉尔
博乐*博乐精河*阿拉山口*温泉
吐鲁番吐鲁番鄯善托克逊*红柳河*吐鲁番东坎
哈密哈密伊吾巴里坤*淖毛湖
阿克苏阿克苏温宿沙雅拜城阿瓦提库车柯坪新和乌什
和田和田墨玉皮山洛浦策勒于田民丰
喀什喀什巴楚泽普伽师叶城岳普湖麦盖提英吉沙莎车塔什库尔干
阿图什阿图什阿克陶阿合奇乌恰*吐尔尕特
库尔勒库尔勒轮台尉犁若羌且末焉耆和静和硕博湖*塔中*库米什巴音布鲁克铁干里克
昌吉昌吉阜康米泉呼图壁玛纳斯奇台吉木萨尔木垒北塔山
伊宁伊宁霍城巩留新源昭苏特克斯尼勒克察布查尔*伊宁县*霍尔果斯
塔城塔城乌苏额敏沙湾托里裕民和布克赛尔*和丰
阿勒泰阿勒泰布尔津富蕴福海哈巴河青河吉木乃
四川:
振幅处理及提高信噪比、分辨率的处理方法
吃的用的东西应有尽有,招待客人的东西数不胜数,糖果瓜子不在话下,什么开心果、腰果,榛子……
对于过年,总有一种说不出的欢悦、欣喜和久久的期待。
欢声笑语透着和谐和吉祥。在灯笼里、丰富多彩的庙会里、在声声的祝福里,迎来了新的一年。
家人下楼点燃鞭炮,顿时满城近带鞭炮声,震耳欲聋,花炮升腾五彩斑斓,整个城市沉浸在烟花爆竹声中。
空气里弥漫着过年的香甜气。愈来愈浓。
相信在寒冷的冬天,没有我的祝福您也一样的快乐,有了我的问候您会更加幸福!虽然相隔千里,我依然能够看到您的微笑,亲爱的朋友,祝您新年快乐。。
新年天气预报:傍晚下人民币,西北方有支票,局部有金块,半夜有零星美元转欧元,请带好钢盔,备好麻袋,准备发财!祝您新年提钱快乐!
新年之际祝您:东行吉祥,南走顺利,西出平安,北走无虑,中有健康;左逢源,右发达,前有福星,后有菩萨;内积千金,外行好运!祝您新年快乐!
眼看着春节就要来临了,这几天走在街上到处都是拎着大包小包购买年货的行人,每个人的脸上,身上都是一片喜气洋洋,仿佛过年的快乐从人们的心里溢了出来流淌到了全身。
一家人围在堂屋里看春晚。当电视里新年钟声敲响时,爸爸已经在院坝里点燃了接年的爆竹,噼噼啪啪的鞭炮声,在寂静的山野里是那样的清脆,响亮。孩子们在院坝里放起了烟花,一朵朵美丽的礼花在天空中闪烁、爆炸,小孩在院坝里欢呼着,跳跃着,新年就在这欢乐,祥和的氛围中来到了。
真的过年了,在忙忙碌碌行人的脚步中,在手拎肩扛物件的人群里,在一声声问候的话语间,在震动山河的鞭炮声中,在济济一堂座无虚席的年夜饭的壮观场面里,在熙熙攘攘的庙会人群里。
盼望着,盼望着,春节终于在我们的期待中,风风火火地赶来了。
天气虽然寒冷,却挡不住新年的喜气:树梢上的红气球,屋檐下的红灯笼,地面火红的鞭炮屑,还有家家户户门口贴的红春联,以及小孩子手里拿着的红包,都见证了春节给我们带来的喜悦。
除夕之夜,千家万户都溢出了欢声笑语,大街小巷响起“噼噼啪啪”的爆竹声,天空中不时升起五彩缤纷的烟花。夜空仿佛一个偌大的电视荧屏,正在播放万家庆新春的精彩节目,真是“灯火辉煌庆佳节,欢声笑语迎新年”啊!
五彩缤纷的礼花,有的像钢花飞溅,有的似麦浪滚滚,有的好比金雨阵阵,有的仿佛银花朵朵,有的恰似红燕飞舞,有的宛如孔雀开屏。
走进庙会,人们比肩接踵,一个个面带微笑,一派喜庆祥和的景象。看,一串串大糖葫芦笑弯了腰;听,嘎嘎作响的小风车笑得合不拢嘴;这边捏面人儿的捏了个孙悟空;那边吹糖人的吹了一个小羊羔儿。各种玩具各样商品,真是琳琅满目。
“春节”是我们中华民族的传统佳节。每到这时候,家家户户都要聚在一起吃团圆饭。外来打工者都急急忙忙地收拾行李,准备回家,拿着一年的工资,心中暖洋洋的;我们小孩子,取得了好成绩,心中美滋滋的,等待着家长们的夸奖……大街上、小区里、超市中,到处都洋溢着喜气,一派热闹的场面。
今天是大年三十。早晨,金色的阳光撒满大地。我和妈妈兴高采烈地去街上买东西,街道上一夜之间仿佛变了个样,大街上张灯结彩,热闹非凡,家家户户贴着“倒福”、对联,还挂着红灯笼。人们穿着节日的盛装,纷纷脸带微笑地走向各个市场,有的购物,有的观赏,一片欢乐详和的景象。
到了晚上8点钟,《春节联欢晚会》开始了,我们一家坐在电视机前看《春节联欢晚会》,晚会的节目很精彩:有优美动听的歌曲,有幽默风趣的小品,还有滑稽逗人的相声…...节目丰富多彩,看得我眼花缭乱。
“当当当……”随着新年钟声的敲响,烟花和爆竹便齐鸣起来。我和爸爸拿出了早已准备好的烟花,点燃了烟花筒,只听见“轰”的一声巨响,一颗烟花弹升到了空中,在一瞬间爆炸了,那爆炸了的烟花仿佛是一朵美丽的莲花在空中展开了花瓣,这时,一颗颗烟花又从烟花筒中喷发,像无数明亮而璀璨的流星,在天空中一闪而过……天空中还闪烁着其他的烟火,它们有的如一串串珍珠,有的如一颗颗流星,有的如一朵朵菊花,有的如一条条瀑布……让我目不暇接。一声声响亮的爆竹声,一朵朵光芒四射、灿烂无比的烟花,在天空中飞舞,各式各样,五光十色,把天空点缀成一朵鲜花, 夜空顿时变得光彩夺目。此时的场面,被笑声、鞭炮声、喊声、乐曲声合成的旋律包围着,真是热闹极了。
到了初一,就是逛公园的时候了。早晨,我们洗漱过后,早点吃的是汤圆,汤圆是用白酒、鸡蛋和红糖煮的,因为家乡的大年初一早上有“不杀生”、“不吃荤”的风俗。吃过早点,我们兴高采烈地出发了。刚到公园门口,我就被丰富的小吃吸引住了,当时我还真有些眼花缭乱呢。进到公园里面,我带着小表弟划船、座小飞机、骑小马、座猴子抬轿,爸爸还给我们照了好多照片。
今年这个春节我过的真是太有意义了,我真是快乐,让我不禁想起了王安石的《元日》:爆竹声中一岁除,春风送暖入屠苏。千门万户瞳瞳日,总把新桃换旧符。
新的一年又到了,在新的一年里,一切都将发生变化,自己又长大了一岁,我希望这一年好运一直陪伴着我,我也会更加努力的,争取在这一年中,能取得优异的成绩,长大后,为家乡的建设作出巨大的贡献
在地震资料处理中,高度保持地震波的真振幅特征,尽量提高地震记录的信噪比和分辨率,即称为“三高”处理,这一直是地震资料处理人员追求的目标。因为“三高”处理的质量直接影响到岩性参数提取以及地震勘探的精度和效果。
10.3.1 真振幅恢复
保持地震波的真振幅特征(简称保幅处理),从广义讲应包含两大方面内容:即真振幅恢复(或称振幅补偿)和其他各项处理中的振幅保持问题。本节主要讨论真振幅恢复的方法,而对其他各项处理中凡要影响到振幅特征的处理方法,则要采取相应的措施,尽可能的使振幅的相对关系保持不变。
地震记录经增益恢复处理后,其振幅特征已与地表检波器所接收到的地震波振幅特征一致。这种振幅仍不称为真振幅,我们所谓的真振幅是指由地层波阻差而产生的反射波振幅,即能反映地层岩性变化的振幅。在地表所接收到的振幅除有地层波阻抗的变化因素外,还有球面扩散因素以及非弹性衰减的因素,因此需要消除球面扩散和非弹性衰减的影响,恢复地震波的真振幅特征。
球面扩散是当波离开震源传播时由于波前扩展造成的振幅衰减。这样的振幅衰减(A)与传播距离r成反比
勘查技术工程学
其中v是界面上覆介质的平均速度;t是反射的记录时间。对球面扩散作校正需要用时变函数vt乘以数据。
非弹性衰减是弹性波能量在岩石中传播时,由于内摩擦而耗散为热被地层吸收的结果。原理部分已说明这种衰减是频率和传播距离的指数形式的函数
勘查技术工程学
其中α为非弹性衰减系数(吸收系数)
勘查技术工程学
所以,用函数eαvt乘以数据就可校正非弹性衰减。至此,真振幅恢复处理完成。
系数α可从增益恢复及球面扩散校正后的振幅-时间函数来测定。为了得到α的较好统计估计,要用一组地震道测定能量来求得衰减曲线。
还有另一种真振幅恢复的方法,这时不需要速度信息。在增益恢复之后,假设振幅衰减是指数函数。因此,按照最小平方法,用指数函数拟合增益校正后的记录,就得到真振幅校正函数(即包括球面扩散和非弹性衰减校正两者)。
前述已知,波前发散因子K为
勘查技术工程学
式中r和t分别为波的传播距离和传播时间,C和a为与地层速度有关的常数。
吸收衰减因子是
勘查技术工程学
式中α是吸收系数;b是待定的常数。波前发散和吸收衰减总的影响是
勘查技术工程学
求取a和b的方法如下。
从地震记录上读取反射波的振幅极值(波峰或波谷),以(10.3-5)为回归方程,得
勘查技术工程学
式中:ut=lnAi-lnti;Ai,ti为振幅极值及其对应的时间;N为振幅极值的点数。校正函数是a-1tebt。
为了获得有代表性的真振幅恢复参数,所选的地震道应是没有多次波及有较高的信噪比。对地质条件稳定地区,一组参数就可代表全区。在工区内地质条件有较大变化时,这些参数要重新计算。
10.3.2 提高信噪比的数字滤波处理
在地震勘探中,用于解决地质任务的地震波称为有效波,而其他波统称为干扰波。压制干扰,提高信噪比是一项贯穿地震勘探全过程的任务。除在野外数据采集中采取相应措施压制干扰外,在地震资料数字处理中数字滤波也是一项非常重要的提高信噪比的措施。
数字滤波方法是利用有效波和干扰之间频率和视速度方面的差异来压制干扰的,分别称为频率滤波和视速度滤波。又因频率滤波只需对单道数据进行运算,故称为一维频率滤波。实现视速度滤波需同时处理多道数据,故称为二维视速度滤波。本节主要介绍这两种滤波方法。
10.3.2.1 一维频率滤波
所谓一维数字滤波是指用计算机实现对单变量信号的滤波,该单变量可以是时间或频率,也可以是空间或波数。以时间或频率为例讨论一维数字滤波,其他原理相同。
(1)一维数字滤波原理
设地震记录x(t)是由有效波s(t)和干扰波n(t)组成,即
勘查技术工程学
其频谱为
勘查技术工程学
式中:X(f)为 x(t)的频谱;S(f)、N(f)分别为 s(t)、n(t)的频谱。如果 X(f)的振幅谱|X(f)|可用图10-6表示。说明有效波的振幅谱|S(f)|处在低频段,而干扰波的振幅谱处于高频段。
图10-6 有效波和干扰波频谱分布示意图
若设计一频率域函数 H(f)的振幅谱为|H(f)|,
勘查技术工程学
其图形为图10-7(a)所示。
令
勘查技术工程学
及
勘查技术工程学
在时间域有(利用傅里叶变换的褶积定理)
勘查技术工程学
称 H(f)为一维滤波器频率响应,(10.3-9)式为频率域滤波方程,h(t)为 H(f)的时间域函数,称为一维滤波器滤波因子(图10-7(b))。(10.3-11)为时间域滤波方程,y(t)和 Y(f)分别为滤波后仅存在有效波的地震记录及频谱,φx(f)、φy(f)、φh(f)分别为滤波前、滤波后地震记录及滤波器的相位谱,以上滤波主要是利用了有效波和干扰波的频率差异消除干扰波,故也称为频率滤波。
图10-7 滤波频率响应及滤波因子
以上所述的滤波器称为理想低通滤波,根据有效波和干扰波的频段分布不同,还可将滤波器分为理想带通滤波器、理想高通滤波器等。所谓理想是指滤波器的频率响应是一个矩形门,门内的有效波无畸变地通过,称为通频带,而门外的干扰波全部消除。在数字滤波中这一点实际是做不到的。因为数字滤波时所能处理的滤波因子只能是有限长,而由间断函数组成的理想滤波器的滤波因子是无限长的。实际应用中只能截断为有限长,截断后就会出现截断效应,即截断后的滤波因子所对应的频率响应不再是一个理想的矩形门,而是一条接近矩形门,但有振幅波动的曲线,这种现象称为吉普斯现象。
由于频率响应曲线在通频带内是波动的曲线,滤波后有效波必定会发生畸变。另外,在通频带外也是波动的曲线,必定不能有效地压制干扰。为了避免吉普斯现象,可采用若干方法,其中之一是镶边法。
10.3.2.2 二维视速度滤波
(1)二维视速度滤波的提出
在地震勘探中,有时有效波和干扰波的频谱成分十分接近甚至重合,这时无法利用频率滤波压制干扰,需要利用有效波和干扰波在其他方面的差异来进行滤波。如果有效波和干扰波在视速度分布方面有差异,则可进行视速度滤波。这种滤波要同时对若干道进行计算才能得到输出,因此是一种二维滤波。
地表接收的地震波动实际上是时间和空间的二维函数g(t,x),即是振动图和波剖面的组合,二者之间通过
勘查技术工程学
发生内在联系。式中k为空间波数,表示单位长度上波长的个数,f为频率,描述单位时间内振动次数,v为波速。
实际地震勘探总是沿地面测线进行观测,上述波数和速度应以波数分量kx和视速度v*代入。则有
勘查技术工程学
既然地震波动是空间变量x和时间变量t的二维函数,且空间和时间存在着密切关系,无论单独进行哪一维滤波都会引起另一维特性的变化(例如单独进行频率滤波会改变波剖面形状,单独进行波数滤波会影响振动图形,产生频率畸变),产生不良效果。那么只有根据二者的内在联系组成时间空间域(或频率波数域)滤波,才能达到压制干扰,突出有效波的目的。因此,应该进行二维滤波。
(2)二维视速度滤波的原理
二维滤波原理是建立在二维傅里叶变换基础上的。沿地面直测线观测到的地震波动g(t,x)是一个随时间和空间变化的波,通过二维正、反傅里叶变换得到其频率波数谱G(ω,kx)和时空函数。
勘查技术工程学
上式说明,g(t,x)是由无数个角频率为ω=2πf、波数为kx的平面简谐波所组成,它们沿测线以视速度v*传播。
如果有效波和干扰波的平面简谐波成分有差异,有效波的平面谐波成分以与干扰波的平面谐波成分不同的视速度传播如图10-8,则可用二维视速度滤波将它们分开,达到压制干扰,提高信噪比的目的。
(3)二维滤波的计算
图10-8 有效波和干扰波以不同成分平面简谐波的传播
二维线性滤波器的性质由其空间-时间特性h(t,x)或频率-波数特性H(ω,kx)所确定。同一维滤波一样,在时-空域中,二维滤波由输入信号g(t,x)与滤波
算子h(t,x)的二维褶积运算实现,在频率-波数域中,由输入信号的谱G(ω,kx)与滤波器的频率波数特性H(ω,kx)相乘来完成。
勘查技术工程学
由于地震观测的离散性和排列长度的有限性,必须用有限个(N个)记录道的求和来代替对空间坐标的积分。
勘查技术工程学
式中,n为原始道号;m为结果道号。
由(10.3-15)式可见,二维褶积可归结为对一维褶积的结果再求和。故测线上任一点处二维滤波的结果可由N个地震道的一维滤波结果相加得到。这时每一道用各自的滤波器处理,其时间特性hm-n(t)取决于该道与输出道之间的距离。沿测线依次计算,可以得到全测线上的二维滤波结果(图10-9)。
与理想一维滤波一样,理想二维滤波也要求在通放带内频率-波数响应的振幅谱为1,在通放带外为0,相位谱亦为0,即零相位滤波。因此,二维理想滤波器的频率-波数响应是正实对称函数(二维对称,即对两个参量均对称),空间时间因子必为实对称函数。二维滤波同样存在伪门现象和吉普斯现象,也可采用镶边法和乘因子法解决。因是二维函数,情况复杂得多,通常只采用减小采样间隔(包括时间采样间隔Δt 和频率采样间隔Δf)和增大计算点数(包括时、空二方向上的点数 M 和N)的方法。
图10-9 二维滤波计算示意图(N=5)
(4)扇形滤波
最常用的二维滤波是扇形滤波。它能滤去低视速度和高频的干扰。其频率波数响应为
勘查技术工程学
图10-10 扇形滤波器的频率波数响应
通放带在f-kx平面上构成由坐标原点出发,以f轴和kx轴为对称的扇形区域(图10-10)。因此这种滤波器称为扇形滤波器。
利用傅里叶反变换可求出其因子为
勘查技术工程学
当在计算机上实现运算时,需要离散化。对时间采样:t=nΔ,n=0,±1,±2,……,Δ为时间采样间隔,Δ=1/2fc。空间采样间隔即输入道的道间距Δx。
由标准扇形滤波器可以组构出既压制高视速度干扰,又压制低视速度干扰的切饼式滤波器,进而还可组构出同时压制高、低频干扰的带通扇形滤波器和带通切饼式滤波器。
在叠加前应用扇形滤波,压制的目标可以是面波、散射波、折射波或电缆振动产生的波。至于在叠加后的应用,则可压制从倾斜界面上产生的多次反射或侧面波。
10.3.3 提高纵向分辨率的反滤波处理
由地震波的传播理论可知,在介质中地震波是以地震子波的形式在地下传播。地面接收到的反射波地震记录是地层反射系数与地震子波的褶积。因此,地层相当一个滤波器,使反射系数序列变成了由子波组成的地震记录,降低了地震勘探的纵向分辨率。反滤波的目的就是要设计一个反滤波器,再对地震记录滤波,消除地层滤波的作用,提高地震记录的纵向分辨率。
由前所述,地震记录是地层反射系数序列r(t)与地震子波b(t)的褶积,即
勘查技术工程学
由于子波的问题,使高分辨率的反射系数脉冲序列变成了低分辨的地震记录,b(t)就相当地层滤波因子。为提高分辨率,可设计一个反滤波器,设反滤波因子为a(t),并要求a(t)与b(t)满足以下关系
勘查技术工程学
用a(t)对地震记录x(t)反滤波
勘查技术工程学
其结果为反射系数序列。以上即为反滤波的基本原理。
反滤波在具体实现时,核心是确定反滤波因子a(t)。由于地震子波的不确定性以及地震记录中噪音干扰的存在,实际中要确定精确的a(t)是非常困难的,甚至是不可能的。为此在不同的近似假设条件下,相继研究了很多种确定反滤波因子a(t)的方法,这些方法基本可以分为两大类:一类是先求取地震子波b(t),再根据b(t)求a(t);另一大类是直接从地震记录中求a(t)。每一类中又有很多不同的方法(就仅反滤波方法之多,说明了反滤波处理的难度)。下面就反滤波方法中具有代表性的几种反滤波进行讨论。
10.3.3.1 地层反滤波
地层反滤波属于先求子波b(t),再求a(t)的方法。该方法要求有测井资料以及较好的井旁地震记录道。首先由声波测井资料转换与井旁地震记录道x(t)相匹配的地层反射系数序列r(t),对r(t)及x(t)求其频谱可得频率域方程为
勘查技术工程学
即有
勘查技术工程学
式中B(ω)为子波b(t)的谱,再由子波与反滤波因子的关系有
勘查技术工程学
经反傅里叶变换得 a(t)。式中 A(ω)为反滤波因子的频谱。写成 z 变换,为 A(z)=,可见A(z)是一个有理分式,要使A(z)具有稳定性,分母多项式B(z)的根必须在单位圆外,即要求子波b(t)为最小相位。
利用测井和井旁地震道求取子波及反滤波因子,即可用该反滤波因子对测线的其他道进行反滤波。
10.3.3.2 最小平方反滤波
最小平方反滤波是最小平方滤波(或称维纳滤波、最佳滤波)在反滤波领域中的应用。
最小平方反滤波的基本思想在于设计一个滤波算子,用它把已知的输入信号转换为与给定的期望输出信号在最小平方误差的意义下是最佳接近的输出。
设输入信号为x(t),它与待求的滤波因子h(t)相褶积得到实际输出y(t),即y(t)=x(t)*h(t)。由于种种原因,实际输出y(t)不可能与预先给定的期望出(t)完全一样,只能要求二者最佳地接近。判断是否最佳接近的标准很多,最小平方误差准则是其中之一,即当二者的误差平方和为最小时,则意味着二者为最佳地接近。在这个意义下求出滤波因子h(t)所进行的滤波即为最小平方滤波。
若待求的滤波因子是反滤波因子a(t),对输入子波b(t)反滤波后的期望输出为d(t),实际输出为y(t),按最小平方原理,使二者的误差平方和为最小时求得的反滤波因子称为最小平方反滤波因子。用它对地震记录x(t)进行的反滤波为最小平方反滤波。
设输入离散信号为地震子波b(n)={b(0),b(1),…,b(m)},待求的反滤波因子a(n)={a(m0),a(m0+1),a(m0+2),……,a(m0+m)},m0为a(t)的起始时间,(m+1)为a(t)的延续长度,b(n)与a(n)的褶积为实际输出y(n),即
勘查技术工程学
为地震子波与期望输出的互相关函数。
根据最小平方原理,经推导即可得到最小平方反滤波的基本方程:
勘查技术工程学
式中,
勘查技术工程学
为地震子波的自相关函数,
勘查技术工程学
为地震子波与期望输出的互相关函数。
(10.3-24)式是一个线性方程组,写成矩阵形式为
勘查技术工程学
式中利用了自相关函数的对称性。该方程中,系数矩阵为一种特殊的正定矩阵,称为一般的托布里兹矩阵,该矩阵方程可用莱文森递推算法快速求解。
式(10.3-27)适应子波b(n)为最小相位、最大相位和混合相位。式中反滤波因子a(n)的起始时间m0与子波的相位有关,其取值规则由子波及反滤波因子的z变换确定。
10.3.3.3 预测反滤波
预测问题是对某一物理量的未来值进行估计,利用已知的该物理量的过去值和现在值得到它在未来某一时刻的估计值(预测值)的问题。它是科学技术中十分重要的问题。天气预报、地震预报、反导弹的自动跟踪等都属于这类问题。预测实质上也是一种滤波,称为预测滤波。
(1)预测反滤波原理
根据预测理论,若将地震记录x(t)看成一个平稳的时间序列,地震子波b(t)为物理可实现的最小相位信号,反射系数r(t)为互不相关的白噪声,由地震记录的褶积模型,在(t+α)时的地震记录x(t+α)为
勘查技术工程学
分析(10.3-28)式的第一项
勘查技术工程学
可见这一项是由反射系数r(t)的将来值决定的。若令第二项为
勘查技术工程学
^x(t+α)是 t 和t 以前时刻的r(t)值决定的,也就是说(t+α)可由现在和过去的资料预测,称(t+α)为预测值。求 x(t+α)与(t+α)的差值为
勘查技术工程学
ε(t+α)称为预测误差,或称为新记录。比较(10.3-28)及(10.3-29)两式,当预测值已知时,从原记录x(t+α)中减去预测值(t+α)后形成的新记录ε(t+α)中比原记录中涉及的反射系数少,与子波褶积后波形的干涉程度轻,波形易分辨,即分辨率提高了。
在上式中α称为预测距或预测步长。当α=1时,
勘查技术工程学
即有
勘查技术工程学
这时(t+1)时刻的预测误差与反射系数之间仅差一个常数b(0)。
因此,选预测距α=1,预测误差为反射系数,达到了反滤波的目的,此时称为预测反滤波。
当α>1时,预测误差为预测滤波结果,预测滤波主要用于消除多次波,尤其是消除海上鸣震。
(2)计算预测值(t+α)的方法
在预测滤波及预测反滤波中,关键是计算预测值(t+α),其方法如下。
由反滤波方程
勘查技术工程学
代入预测值(t+α)的表达式
勘查技术工程学
式中令τ=s-j,c(s)=b(j+α)a(s-j)称为预测因子。a(t)为反滤波因子。预测值(t+α)为预测因子 c(s)与地震记录的褶积。
现在需设计一个最佳预测因子c(s),使求取的预测值(t+α)与x(t+α)最接近,即使预测误差的平方和(误差能量)
勘查技术工程学
为最小。根据最小平方原理,可得线性方程组
勘查技术工程学
式中Rxx(τ)为地震记录的自相关函数
勘查技术工程学
T为相关时窗长度,m+1是预测因子长度。将(10.3-34)写成矩阵形式为
勘查技术工程学
解此方程组即可求得预测滤波因子c(t),用它对地震记录x(t)褶积可以求出未来时刻(t+α)时的最佳预测值(t+α)。
声明:本站所有文章资源内容,如无特殊说明或标注,均为采集网络资源。如若本站内容侵犯了原著者的合法权益,可联系本站删除。