Tính toán chỉ tiêu bảo vệ chống sét cho trạm biến áp trung áp 35kv

Bộ giáo dục đào tạo Trường Đại học Bách Khoa Hµ Néi - LUẬN VĂN THẠC SĨ KHOA HỌC TÍNH TỐN CHỈ TIÊU BẢO VỆ CHỐNG SÉT CHO TRẠM BIẾN ÁP TRUNG ÁP 35KV LẠI VĂN CNG Luận Văn Thạc Sĩ Khoa Học CHUYấN NGNH: H THNG IN Hà Nội - 2005 Bộ giáo dục đào tạo Trường Đại học Bách Khoa Hà Nội - LUẬN VĂN THẠC SĨ KHOA HỌC TÍNH TOÁN CHỈ TIÊU BẢO VỆ CHỐNG SÉT CHO TRẠM BIẾN P TRUNG P 35KV LI VN CNG Luận Văn Thạc SÜ Khoa Häc CHUYÊN NGÀNH: HỆ THỐNG ĐIỆN NGƯỜI HƯỚNG DẪN PGS – TS TRẦN VĂN TỚP Hµ Néi - 2005 Lời cảm ơn Tôi xin bày tỏ lòng biết ơn thầy giáo hớng dẫn PGS.TS Phạm Thợng Hàn, ngời đà tận tình hớng dẫn bảo vần đề khoa học để hoàn thành luận văn Tôi chân thành cảm ơn đồng nghiệp KS Nguyễn Quang Minh, Tiến sỹ Quản Hoàng Lâm ngời đà cung cấp tài liệu, thiết bị để thực hoàn thành luận văn Tôi chân thành cảm ơn đông nghiệp thuộc phòng công nghệ mô Viện UDCN, đồng nghiệp Ban Đo lờng Tổng cục TC-ĐL- CL đà tạo điều kiện cho thời gian học tập hoàn thành luận văn Hà Nội, ngày 22 tháng 11 năm 2005 Trang Đặt vấn đề Ung th bệnh nan y xà hội loài ngời Hàng năm, Việt nam cã tõ 50.000 –70.000 ng−êi chÕt v× ung th− Nhê phát triển y học đại, liệu pháp chuẩn đoán điều trị ung th có nhiều tiến đáng kể.Tuy nhiên tỷ lệ ung th tái phát, di cao.Các thống kê y học cho thấy phát ung th sớm, hiệu điều trị ung th da đạt tới 70 %-80% tỉ lệ di thấp Bởi vậy, song song với nghiên cứu tìm liệu pháp điều trị ung th mới, nhà khoa học tập trung nghiên cứu phơng pháp chẩn đoán ung th có hiệu cao Một phơng pháp chẩn đoán ung th dựa phân tích xử lý phổ huỳnh quang mô tế bào bị Laser kích thích theo phơng pháp Quang động học (gọi tắt phơng pháp chẩn đoán ung th Laser) Phơng pháp phận liệu pháp Quang động học đợc Dougherty TJ ứng dụng lần Mỹ vào năm 1978 Chẩn đoán ung th Laser thuộc trờng phái chẩn đoán quang phổ học hớng nghiên cứu Yhọc năm gần Phơng pháp chủ yếu dựa hình ảnh huỳnh quang đồ thị phổ mô bị kích thích Laser nguồn UV Uu điểm phơng pháp chẩn đoán đo không tiếp xúc không tác động lên đối tợng nghiên cứu Đây yêu cầu quan trọng thăm khám điều trị ung th Bên cạnh nhu cầu hệ chẩn đoán điều trị ung th− rÊt cÊp thiÕt Bëi c¸c xÐt nghiƯm ung th− tèn nhiỊu thêi gian míi cho kÕt qu¶ Trong với chẩn đoán ung th thiết bị Laser cần vài phút có kết nhanh chóng thuận tiện Vì lý tiến hành nghiên cứu đề tài Nghiên cứu phơng pháp chẩn đoán ung th Laser nhằm mục đích tạo đợc Algorithm tự động đánh giá phổ huỳnh quang tạo tiền đề xây Trang dng phần mềm chẩn đoán ung th trợ giúp bác sỹ chẩn đoán điều trị ung th Cấu trúc luận văn gồm chơng chia thành phần : + Phần I : Nghiên cứu lý thuyết gồm chơng - Chơng 1: Nghiên cứu tổng quan phơng pháp chẩn đoán ung th trình bày tóm tắt phơng pháp chẩn đoán cận lâm sàng nh : chụp X- quang, miễn dịch học - Chơng 2: trình bày đại cơng khối u ung th, phân loại, loại ung th phổi - Chơng 3: trình bày sở lý thuyết phơng pháp chẩn đoán ung th Laser + Phần II: Thực hành bao gồm - Chơng 4: Xây dựng Algorithm chẩn đoán, giải thích thuật toán làm trơn đờng cong, tính diện tích - Chơng 5: Thực nghiệm đánh giá phổ số liệu thu đợc từ thiết bị quang phổ Laser mô hinh thực nghiệm ung th chuột bệnh nhân u biểu mô tế bào đáy vành tai Trang Phần I: Các nghiên cứu lý thuyết Trang Chơng nghiên cứu Tỉng quan vỊ Khèi u ung th− BƯnh ung th− gì? Ung th l tờn chung dựng gi nhóm bệnh 200 loại khác nguồn gốc tế bào, nguyên, tiên lượng cách thức điều trị Nhưng bệnh có đặc điểm chung có phân chia khơng kiểm sốt tế bào, khả tồn phát triển quan tổ chức lạ Các ung thư thường phát triển từ tế bào ban đầu phải nhiều năm có kích thước đủ lớn để nhận thấy Qu¸ tr×nh phát triển từ tế bào thành khối ung thư trải qua nhiều giai đoạn Thông thường, tế bào lành có tuổi thọ định tuân thủ theo quy luật chung phát triển - già - chết Các tế bào chết thay tế bào Cơ thể có chế kiểm sốt quy luật cách chặt chẽ tŕi số lượng tế bào quan, tổ chức mức ổn định Bệnh ung thư bắt đầu có tế bào vượt qua chế kiểm soát thể, bắt đầu phát triển sinh sôi không ngừng nghỉ, h×nh thành đám tế bào có chung đặc điểm phát triển vô tổ chức, xâm lấn chèn ép vào quan tổ chức xung quanh Các tế bào ung thư có liên kết lỏng lẻo, dễ dàng bứt khỏi khối u mẹ, theo mạch máu mạch bạch huyết di cư đến tổ chức quan mới, bám lại tiếp tục sinh sơi nẩy nở (q tr×nh gọi “di căn”) Các ung thư chèn ép di vào quan giữ chức sống thể năo, phổi, gan, thận bệnh nhân tử vong.[1] Trang 1.1 Khái niệm khối u Không có định nghĩa u hoàn toàn xác chất nguyên ung th đến cha đợc biết rõ Đó mô tả chung để khối mô bất thờng phát triển không hài hoà với mô bình thờng kế cận, tồn mÃi mÃi sau đà ngừng kích thích sinh u Thuật ngữ ung th dùng để khối u ác tính Tất u lành ác tính có hai thành phần bản: 1, tế bào u tăng sinh tạo nên nhu mô u mô đệm nâng đỡ đợc hình thành từ mô liên kết huyết quản Mặc dù tế bào nhu mô thành phần tăng sinh quan trọng u định chất u, phát triển tiến triển u phụ thuộc nhiều vào mô đệm Việc cung cấp máu mô đệm liên kết cần thiết mô đệm nâng đỡ nghèo nàn, u mềm giống thịt Đôi tế bào u kích thích tạo mô đệm giàu sợi tạo keo gọi tạo xơ a.Một số khái niệm chung Biểu mô: loại mô đợc tạo thành nhóm tế bào nằm sát cung thực chức định Biểu mô có nguồn gốc nội bì ngoại bì.ví dụ: biểu mô giác mạc có nguồn gốc ngoại bì, tuyến phụ thuộc ống tiêu hoá gan, dày, ruột có nguồn gốc nội bì Tuỳ thuộc vị trí biểu mô thể biểu mô thực chức năng: - Bao phủ mặt thể (biểu bì da) - Lớp mặt tạng rỗng (dạ dày, rt, tư cung ) khoang miƯng, mịi - HÊp thơ tiết thể qua biểu mô - Bảo vệ thể chống lại va chạm học, chống bốc làm độ ẩm da - Lớp thu nhận cảm giác quan cảm giác Trang U ác tính biểu mô (carcinoma) đợc xếp loại dựa tế bào nguồn gốc chúng biểu mô Ung th biểu mô đợc phân chia nguồn gốc : loại có cấu trúc tuyến gọi biểu mô tuyến nh ung th biểu mô tế bào thận Loại có cấu trúc vảy xuất phát từ biểu mô thể đợc gọi biểu mô tế bào vảy U ác tính tuyến (adenoma ) đợc áp dụng cho biểu mô lành tính tạo thành cấu trúc tun cịng nh− u cã ngn gèc tõ c¸c tun nhng không thiết tạo thành từ hình thái tuyến Ví dụ u biểu mô lành tính phát sinh từ tế bào biểu mô ống thận dới dạng nhiều tuyến nhỏ xếp sát đợc gọi u tuyến U ác tính sarcom phát sinh từ mô trung mô nh thịt u có mô liên kết nên giống thịt Ví dụ sacôm xơ, sacôm mỡ, sacôm trơn với ung th xơ , mỡ Trong thực hành ngời ta thờng thêm tên quan nguồn gốc dụ ung th biểu mô tế bào thận, ung th biểu mô tế bào vảy phế quản b Phân loại khối ác tính Có nhiều cách phân loại khác nhng phổ biến theo nguồn gốc tế bào Nguồn gốc mô ác tính 1.Hình thành từ loại tế bào biểu mô trung mô Tế bào vảy U biểu mô tế bào vảy hay thợng bì Các tế bào đáy da U biểu mô tế bào đáy phần phụ Thần kinh ngoại bì U hắc tố ác tính Đờng hô hấp U biểu mô tuyến phế quản Trang Biểu mô thận U biểu mô tế bào thận Tế bào gan U biểu mô tế bào gan Mô liên kết Sarcom xơ, mỡ sụn, xơng 2.Tế bào máu tế bào liên quan Các tế bào tạo máu Bệnh bạch cầu Mô lympho U lymphô ác tính Cơ trơn Sarcom trơn Cơ vân Sarcom vân 3.Nội mô mô liên quan Mạch máu Sarcom huyết quản Bạch mạch Sarcom bạch mạch Màng khớp Sarcom màng khớp Màng nÃo U màng nÃo xâm nhập Các u hỗn hợp có nhiều loại tế bào u Các tuyến nớc bọt U hỗn hợp ác tính tuyến nớc bọt Vú Sarcom nang dạng ác tính Thận U Wilms 1.2 Đặc tính u ác tính a/ Biệt hoá biệt hoá: Thuật ngữ áp dụng cho tế bào nhu mô.Sự biệt hoá để mức độ phát triển tế bào nhu mô giống tế bào bình thờng hình thái chức U biệt hoá cao u gồm tế bào thục sinh u, không chuyên biệt Nói chung u lành tính có tính biệt hoá cao Tế bào u u trơn giống tế bào bình thờng làm cho phân biệt tế bào mức độ Chỉ tụ tập tế bào thành cục thấy rõ chất tổn thơng Trang10 5.3.3.Tìm quan hệ Ds Dk Với kết qủa đạt đợc phần 4.2.2b sử dụng đoạn chơng trinh ta xác định đợc đỉnh phổ bớc sóng 690nm là: I λ = 690 = 780.083 vµ: I λ = 690 = 343 Đại lợng chuẩn đoán Dk tính theo c«ng thøc Dk = I 1λ =690 780.083 = = 2.2743 343 I =690 Đối với phần diện tích nằm bên dới đờng cong (1) (2) Ta tính diện tích theo phơng pháp Monter Calor với thuật toán nh sau: Từ Algorithm tổng quát ta xây dựng thuật toán tính diện tích đờng cong nh sau: Đặt toàn đờng cong phổ hình vu«ng cã diƯn tÝch S = b x b víi số điểm nằm trục bớc sóng = 2048 giá trị Chia nhỏ hình vuông nói thành N ô vuông với kích thớc k = đơn vÞ Tỉng sè sÏ cã M = 1024 x 1024 ô vuông tơng ứng với số phép thử M = Đếm số lợng ô vuông N1 nằm dới đờng f1 () đờng mẫu Giá trị f() có đợc từ bảng số liệu đo b Tính S1 theo c«ng thøc S1 = ∫ f ( x)dx = N1 * S M 5, Lặp lại bớc Tính S2 với số lợng ô vuông N2 b S = ∫ f ( x)dx = N1 * S M Sau tÝnh to¸n ta thu đợc tỷ lệ Ds = S2/S1=2.1406 Trang58 Hình 4.7: Đồ thị phổ huỳnh quang cuả mô nghi ác tính (2) mô lành tính (1) Nh kết tính đợc tỷ lệ Dk =Peak1(680)/Peak2(680)=2.2833 Sau chơng trình tính diện tích S1 S2 viết ngôn ngữ Matblab: hàm dientich1(mau) hàm tính diƯn tÝch S1; hµm dientich2(mau) lµ hµm tÝnh diƯn tÝch S2 Toàn kết đợc thực đoạn chơng trinh viết ngôn ngữ Matlab nh sau: function S1= dientich1(mau) A=2048^2; j=1024; M=j^2; N1=0; k=2; for x=1:j for i=1:j y=k*i; if y lần Nh so sánh cờng độ quang phổ mô nghiên cứu với mô bình thờng đối tợng ta thu đợc kết cờng độ phổ lớn lần Theo tiêu chuẩn phân loại phổ đà trình bày chơng 1, tỷ số cho thấy đối tợng nghiên cứu đà mắc u ác tính Kết hoàn toàn phù hợp với mô hình thực nghiệm SARCOM 180 kết xét nghiệm tế bào đợc thực Bộ môn phôi Học viện quân Y 103 mà chúng đà thực thời gian qua Trang61 5.4Kết luận Trong Chơng đà trình bày số kết thực Algorithm chuẩn đoán ung th sở số liệu quang phổ huỳnh quang tổ chức tế bào nằm đối tợng sống Phơng pháp nghiên cứu cấy tế bào ung th lên đối tợng cần nghiên cứu Sau dùng thiết bị quang phổ kế để ghi lại huỳnh quang khối u mô lành tính lân cận khối u Kết đo sau đợc xử lý lọc bỏ nhiễu tính toán đề đa kết luận chuẩn đoán Bằng phơng pháp so sánh với mô hinh thực nghiệm nhận thấy kết tính toán giựa Algorithm mà đa đáng tin cậy Trang62 Kết luận Luận văn đà nghiên cứu hoàn thành số nội dung Vể lý thuyết : nghiên cứu tổng quan phơng pháp chẩn đoán ung th lựa chọn phơng pháp chẩn đoán ung th thiết bi đo quang phổ Laser để thực hành chẩn đoán khối u ung th - Nghiên cứu ung th thể khối u ung th phổi - Nghiên cứu phơng pháp làm trơn đồ thị đờng cong phổ SavitzkyGola, phơng pháp dùng phổ biến phân tích xử lý đờng đồ thị phổ nghiên cứu hoá sinh học - Nghiên cứu phơng pháp Monter Calor tính diện tích - Nghiên cứu đợc nguyên lý thu nhận phổ huỳnh quang thiết bị Laser Về mặt thực hành: đà có số kết định - Xâydựng đợc Algorithm chẩn đoán ung th biểu mô thể tế bào ®¸y cã da - Thùc nghiƯm kiĨm chøng Algorithm với số liệu thu đợc từ thực nghiệm So sánh kết rhực nghiệm với lý thuyết đạt kết tốt - Khắc phục đợc tợng nhiễu thăng giáng đờng comg phổ thiết bị Laser 2.Những hạn chế đề tài Do khó khăn điều kiện thực hành thí nghiệm nh: thuốc nhạy quang Photogem, điều kiện thực hành ngời bênhh Nên luận văn phải sử dụng lại kết nghiên cứu cũ nhóm nghiên cứu thuộc Trung tâm Công nghệ Laser năm 1999 - Số lợng thí nghiệm có thí nghiệm nên việc đánh giá Algorithm hạn chế Trang63 - Việc xây dựng phần mềm đánh giá phổ dừng lại thuật toán chạy thử Matlab 3.Hớng nghiên cứu thời gian tới Nhằm khắc phục hạn chế đề tài : thời gian tới tiếp tục nghiên cứu áp dụng Algorithm xây dựng phần mềm chẩn đoán tù ®éng kÕt nèi trùc tiÕp víi Card ®o cđa thiết bị đo quang phổ Laser để thiết bị trở thành hệ chẩn đoán ung th - Kết hợp với Bộ môn Phôi Học viên Quân Y 103, phòng khám Laser Bệnh viện Việt Đức thử nghiệm nhiều mẫu bệnh phẩm để xây dựng th viện mẫu phổ loại ung th khác - Nghiên cứu sử dung mạng Nơron để phân loại ung th Việc giải toán đòi hỏi nhiều công sức, tài ủng hộ nhiều đồng nghiệp Tuy nhiên nghiên cứu tìm phơng pháp chẩn ®o¸n ung th− míi ®Ị ®−a c¸c chÈn ®o¸n sớm bệnh ung th đấu tranh lâu dài ngời Do công việc lâu dài hứa hẹn nhiều khó khăn Trang64 Tài liệu tham khảo Tiếng Việt Lê Đình Roanh, 2001, Bệnh học khói u, Nhà xuất Y học Quản Hoàng Lâm, 2001,Nghiên cứu tác dụng liệu pháp quang động học biến đổi hình thái chỗ xung quanh mô ung th thực nghiệm SARCOM 180 , luận văn tiến sỹ Y học, Học viện Quân y Dơng Chạm Uyên, Kiểu Đình Hùng,2000, Nhận xét bơc đầu áp dụng quang động học điều trị u nÃo Bệnh viện Việt Đức, Tạp chí thông tin y dợc chuyên đề ung th Nguyễn Quang Minh, Tạ Văn Tuân công ,1999,Nghiên cứuứng dụng phát triển thiết bị Laser chẩn đoán điều trị ung th phơng pháp quang động học, Báo cáo đề tài KHCN cấp nhà nớc Đỗ Kin Sơn, Nguyễn Bửu Triều, Dơng Chạm Uyên,2000.Thông báo số kết bớc đầu ứng dụng Laser Y học Bênh viện Việt Đức Tập thể tác giả Bộ môn mô học phôi thai, 2004,Mô học, Nhà xuất Y học Bộ môn ung th, Bài giảng ung th học, 2001, Nhà xuất Yhäc TiÕng n−íc ngoµi V.B.Loschenov,V.I Konov, A.M.Prokhorov, 2000,”Photodynamic Therapy and Fluorescence Diagnostics”,Laser Physic, Vol 10, No 6, pp 1188-1207 William H Press, Saul A Teukolsky, William T Vetterling,1991,” Numerical Recipe in C: The Art of Scientific Computing “, Cambridge University Press 10 Г И Назаренко, А Б Канючевский, С Н Петров, Б М Касаткин, А.Н.Димитров,http://www.medlinks.ru/ 11 V.B Loschenov and A.A Stratonnikov, 1993,” Spectra line separation Trang65 method for sophisticated data analysis of biological tissue optical spectra”,SPIE Vol 2081 Optical Biopsy 12.Savitzky A., and Golay, M.J.E 1964, Analytical Chemistry, vol 36, pp 1627–1639 13 Jianwen Luo , Kui Ying , Ping He , Jing Bai, 2005, Digital Signal Processing 15,pp 122136 Trang66 Phụ lục chơng trình % % Chơng trinh danh gia huynh quanng cua te bµo ung th− function CT=CT clf; clc; mau=importdata ('c:\matlab6\work\norm.dat'); test3=importdata ('c:\matlab6\work\test3.dat') ; F= load('c:\matlab6\work\data\40.dat','-ASCII'); nn=length(F); h=(700-660)/(length(F)-1);% Doi ve truc toa buoc song xf=660:h:700; %subplot('position',[.1.1 80]); plot (xf,F,'b'); axis([660 700 1000]); title('Pho huynh quang cua to chuc binh thuong ','FontSize',13); xlabel ('\lambda = nm','FontSize',12); ylabel ('Relative Intensive (a.u)','FontSize',12); grid on hold on %plot(xf,mau,'k','LineWidth',3); hold off waitforbuttonpress; thamso = inputdlg({'Chieu dai bo loc(Num) :','Bac da thuc(M):'},'Thiet dat tham so loc',2,{'40','2'}); Trang67 num=str2num(char(thamso(1))); M=str2num(char(thamso(2))); [nl,nr,Y]= thamsoloc(F,nn,num,M); Wfilter=nl+nr+1; %sprintf('Kich thuoc cua so Wf=%d',Wfilter); %loc tin hieu voi cac he so nr, nl, M tinh duoc g=SG(F,nn,nl,nr,M); %g=SavGol_1(F,nl,nr,M); %xg=xf(1+nl:nn-nr); xg=xf; plot (xf,F,'b'); axis([660 700 1000]); hold on; g(1:nl)=10; g(nn-nr:nn)=0; plot(xg,g,'r','LineWidth',3); grid on title ('Pho tin hieu sau loc Savitzky-Golay','FontSize',13); xlabel ('\lambda = nm','FontSize',12); ylabel ('Relative Intensive (a.u)','FontSize',12); legend('Tin hieu co nhieu','Tin hieu sau loc',['Chieu dai cua so loc W = ',int2str(num)]); hold off waitforbuttonpress; % ve tin hieu sau loc voi tin hieu thu plot(xg,g,'r','LineWidth',3); axis([660 700 1000]); Trang68 title ('Pho tin hieu sau loc Savitzky-Golay','FontSize',13); legend('Tin hieu sau loc',['Chieu dai cua so loc W = ',int2str(num)]); grid on hold on waitforbuttonpress; axis([660 700 1000]); plot(xg,g,'r-','LineWidth',3); grid on hold on plot(xf,mau,'k','LineWidth',3); [Dk,Peak]=SS(mau,g,h); [Dk,Peak,nPeak,nPeak1,dientich,dientich1]=SS(mau,g,h); grid on hold on title ('So voi mau','FontSize',13); xlabel('\lambda = nm','FontSize',12); ylabel ('Relative Intensive (a.u)','FontSize',12); legend('Pho tin hieu duoc(S2)','Pho tin hieu mau(S1)'); S2=dt(mau); S1=dt1(g); k=(S1/S2)/4 waitforbuttonpress; legend('Pho tin hieu duoc(S2)','Pho tin hieu mau(S1)',['Ty le Ds=S2/S1=',num2str(k)],['Ty le Dk=Ifc(680)/Ifn(680)=',num2str(Peak)]); waitforbuttonpress; if(Peak>2)warndlg('Ket qua chuan doan duong tinh','Chuan doan ung thu'); hold off; end Trang69 % Savitzky-Golay Smoothing Filter function g=SG(F,nn,nl,nr,M) A=zeros(nl+nr+1,M+1); for i=-nl:nl; for j=0:M; A(i+nl+1,j+1)=i^j; end end h=zeros(M+1,1); h(1)=1; b=inv(A'*A)*h; c=zeros(nl+nr+1,1); for k=-nl:nr nm=k.^[0:M]; c(k+nl+1)=nm*b; end for i=1+nl:nn-nr, g(i) = c' * F (i-nl:i+nr); end end % ham tinh tham so boc loc function [nl,nr,Y] = thamsoloc(F,nn,num,M) for ii=32:num nl=ii; nr=ii; g=SG(F,nn,nl,nr,M); end chi_g=chig(32:num); Trang70 [Y I]=min(chi_g); nl=I+31; nr=I+31; end %Tinh dien tich function [Dk,Peak,nPeak,nPeak1,dientich,dientich1]=SS(mau,F,h) g=F; nPeak=g(680); nPeak1=mau(680); Peak=g(690)/mau(680); end function I=dt(mau) A=2048^2;%% Area of enclosing rectangle j=1024; N=j^2; %% Set the number of trials s=0; %% initialize success counter k=2; for x=1:j for i=1:j %% Begin Trials y=k*i; if y