WO2014087494A1 - 計算システム及び計算方法 - Google Patents

計算システム及び計算方法 Download PDF

Info

Publication number
WO2014087494A1
WO2014087494A1 PCT/JP2012/081445 JP2012081445W WO2014087494A1 WO 2014087494 A1 WO2014087494 A1 WO 2014087494A1 JP 2012081445 W JP2012081445 W JP 2012081445W WO 2014087494 A1 WO2014087494 A1 WO 2014087494A1
Authority
WO
WIPO (PCT)
Prior art keywords
image
unit
measured
information
density length
Prior art date
Application number
PCT/JP2012/081445
Other languages
English (en)
French (fr)
Inventor
浩二 荒牧
雅成 高口
都 松井
Original Assignee
株式会社日立製作所
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 株式会社日立製作所 filed Critical 株式会社日立製作所
Priority to PCT/JP2012/081445 priority Critical patent/WO2014087494A1/ja
Priority to US14/649,924 priority patent/US20150293040A1/en
Priority to JP2014550839A priority patent/JPWO2014087494A1/ja
Publication of WO2014087494A1 publication Critical patent/WO2014087494A1/ja

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
    • G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
    • G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/10—Different kinds of radiation or particles
    • G01N2223/101—Different kinds of radiation or particles electromagnetic radiation
    • G01N2223/1013—Different kinds of radiation or particles electromagnetic radiation gamma
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/20—Sources of radiation
    • G01N2223/205—Sources of radiation natural source
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/40—Imaging
    • G01N2223/401—Imaging image processing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/60—Specific applications or type of materials
    • G01N2223/601—Specific applications or type of materials density profile
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2223/00—Investigating materials by wave or particle radiation
    • G01N2223/60—Specific applications or type of materials
    • G01N2223/631—Specific applications or type of materials large structures, walls

Definitions

  • the present invention relates to a calculation system and a calculation method for calculating the structure of a structure.
  • the non-destructive test has an advantage that the possibility of a decrease in soundness due to the inspection is low as compared with the destructive test. Moreover, it has the advantage that a restoration process after the test is not required, and is applied to various structures. For example, as an element of maintenance inspection of metal pipes conforming to industrial standards, there is a method to clarify the degree of wear inside the pipe by measuring the thickness of the object part to be measured using the reflection of ultrasonic pulses. Used. However, since this method directly contacts the ultrasonic transducer with the object to be measured, it is necessary to remove the pipe covering before the measurement and restore it after the measurement. There is a disadvantage that it takes.
  • Patent Document 1 a technique is also used in which the process of removing and restoring the clothing material of the pipe is omitted by performing radiographic imaging.
  • an artificial radiation source and a detector are installed at positions that sandwich a measurement target object.
  • Patent Document 2 uses a muon of high energy secondary cosmic rays, which is a kind of natural radiation, to capture the degree of wear of the refractory in the furnace by transmitting through the furnace wall and bottom of the blast furnace for iron making. A method to clarify the above is shown.
  • Patent Document 3 discloses a technique for clarifying the degree of change in the material of the measurement object by utilizing the property of the secondary cosmic ray muon passing through the measurement object.
  • Patent Document 1 Although the above-described covering material removal / restoration step can be omitted, since the ability to transmit artificial radiation is weak, it is necessary to arrange a radiation source in the vicinity of the measurement target object. In a large-scale plant, it is impractical to install an artificial radiation source near all objects to be measured, and it is necessary to install and remove the radiation source and detector from each object to be measured at each measurement. There is a drawback that it requires a great deal of cost before and after the actual measurement.
  • the density length of an object that the flight path penetrates is provided as a contrast image using particles with extremely high transparency and straightness. There is no need to place a detector near the object. However, in many cases, the captured transmission image has a defect that the density length including not only the measurement target object but also other objects through which the particle flight path passes is reflected in the contrast image.
  • An object of the present invention is to provide an apparatus or a method for obtaining a contrast image reflecting an actual density length of only a desired measurement target object in a nondestructive test for inspection.
  • a calculation system includes an imaging unit that captures structural information of an arbitrary section by cosmic rays that pass through a structure to be measured and another structure, and a structure having structure information of the other structure An information acquisition unit, a density length calculation unit that determines whether to include the structure information of the other structure as a measurement target, calculates a density length of a structure determined as a non-target, an image captured by the imaging unit, and the And an image calculation unit that calculates structure information of the structure to be measured based on the density length of the structure determined as a non-object.
  • the calculation system includes a track acquisition unit that acquires a track of a cosmic ray that passes through a measurement target and another structure, an imaging unit that captures an image based on the track of the cosmic ray, and the space
  • a first storage unit for storing information on track of the line a second storage unit having structure information of the measurement object and the other structure, and a first storage unit based on the structure information of the other structure.
  • a density length calculating unit that calculates a density image of density length and calculates a density image of the second density length based on the structure information of the measurement object and the structure information of the other structure.
  • the calculation method includes a first step of imaging structure information of an arbitrary section with a cosmic ray that passes through a structure to be measured and another structure, and the other structure.
  • FIG. 1 is a drawing examined by the inventors of the present application prior to the invention.
  • the situation where various non-measurement target objects exist is schematically illustrated in the route of the arrow M1 which is the flight route of the muon.
  • the object to be measured 91 such as piping constituting the large-scale facility 90
  • the muon detector 99 the wall of the building 92 surrounding the object to be measured 91 or the object to be measured 91 other than the object to be measured 91
  • the flight path M1 also penetrates the pipe 93, the mountain 94 outside the building 92, and the other building 95.
  • the flight path includes non-materials and shapes of various materials, such as the walls surrounding the building to be measured and the walls and pillars that make up the external building, or earth and sand that make up the topography of mountains and hills. It is conceivable that an object to be measured exists.
  • the obtained contrast image includes information on the non-measurement target object in addition to the information caused by the object to be originally measured.
  • FIG. 2 is a diagram showing a configuration of a large-scale facility state monitoring system according to an embodiment of the present invention.
  • the large-scale facility state monitoring system 10 of this embodiment includes a track acquisition unit 20 and an information processing unit 30 as components.
  • the track acquisition unit 20 will be described below.
  • the track acquisition unit 20 includes a power supply unit 21, a detection unit 22, and a coincidence counting unit 23 as components.
  • the power supply unit 21 is a power supply unit that supplies necessary power to the detection unit 22 and the coincidence counting unit 23.
  • the detection unit 22 includes a first detector 221, an iron block 222, and a second detector 223 as components.
  • Both the first detector 221 and the second detector 223 are position sensitive detectors.
  • the position-sensitive detector here is one that can detect the coordinate value of the passing point of the detector when charged particles pass through.
  • the arrangement relationship between the first detector 221 and the second detector 223 is preferably a relationship in which the above-described planar detectors are arranged in parallel to each other. If the plane of the second detector is inclined with respect to the plane of the first detector, charged particles that have passed through the measurement object and the first detector 221 cannot be efficiently captured, and this is prevented. In order to achieve this, it is conceivable that a larger plane than the first detector 221 is required. Therefore, the detection efficiency of charged particles can be increased as the respective planes change from the inclined relationship to the parallel relationship.
  • charged particles are incident on the detection unit 22 at irregular timings from various azimuth and zenith angles.
  • These charged particles include not only high-energy muons but also soft components such as electrons.
  • the soft component is extremely inferior in the straightness in the object than the muon, and in order to obtain a contrast image based on the exact density length of the object to be measured, it is necessary to selectively detect only the muon with high straightness. .
  • the second detector 223 can transmit a single coordinate value to the coincidence counting unit 23. This is because the muon passes through the iron block 222 as it is and enters the second detector 223.
  • the second detector 223 transmits a plurality of coordinate values to the coincidence counting unit 23 almost simultaneously. . This is because the particles generate a large number of particles in the process of passing through the iron block 222, and the traveling directions of these particles are various.
  • the coincidence counting unit 23 compares the times transmitted from the first detector 221 and the second detector 223, and the same muon passes only when the difference falls within a certain range. judge.
  • the determined two sets of coordinate values are transmitted to the information processing unit 30. As will be described later, the muon trajectory can be obtained from these two sets of coordinate values determined to be simultaneous.
  • the information processing unit 30 includes a power supply unit 31, a user interface 32, a storage unit 33, and a calculation unit 34 as components.
  • the power supply unit 31 is a power supply unit that supplies necessary power to the user interface 32, the storage unit 33, and the calculation unit 34.
  • the user interface 32 includes an input unit 321 such as a keyboard and a mouse and an output unit 322 such as a display device and a printer as constituent elements.
  • the operator can perform the following operations via the user interface 32. For example, an operation for selecting one or more desired items from the options presented by the calculation unit 34, an operation for selecting a partial range of the geometric shape illustrated by the calculation unit 34, and a finite range presented by the calculation unit 34 For example, an operation for selecting a specific numerical value from the numerical values, and an operation for inputting a value to be substituted into a variable presented by the calculation unit 34.
  • the storage unit 33 includes a track information storage unit 331, a terrain information storage unit 332, and a facility information storage unit 333 as components.
  • the track information storage unit 331 combines the time information with the two sets of muon passage coordinate data transmitted from the coincidence counting unit 23 each time the muon passes and makes track (trajectory) information, and newly stores it. .
  • the terrain information storage unit 332 holds geometric data indicating the terrain around the large-scale facility 90 having the measurement target object as a component, and data indicating the material or density of the terrain object constituting the terrain.
  • the topography includes not only naturally formed mountains and hills, but also components such as artificially formed structures such as dams, tunnels, bridges, and buildings.
  • the equipment information storage unit 333 has a position and a shape within a certain range such as a door, a hoist (for example, an installed crane), in addition to a static building member such as a wall, a pillar, and a beam constituting the large-scale equipment 90. It holds geometric shape data and data indicating the material or density of an object that can store an amount of fluid objects that can change within a certain range, such as an object that can change or a tank. In addition, the facility information storage unit 332 holds various amount values that can change within the certain range, data indicating the change range of the amount, and data indicating the material or density of the fluid object. .
  • the calculation unit 34 includes an imaging unit 341, a density length calculation unit 342, and an image calculation unit 343 as components.
  • the imaging unit 341 selects track information having time information corresponding to a period specified by the operator through the user interface 32 from among a large number of track information held in the track information storage unit 331, and the muon passage coordinates thereof.
  • a contrast image of the measurement object is constructed based on the value data.
  • this contrast image is referred to as an actual measurement image.
  • the imaging unit 341 transmits real image data to the image calculation unit 343.
  • the density length calculation unit 342 constructs a virtual contrast image by a method described later based on various data stored in the topographic information storage unit 332 and the facility information storage unit 333.
  • this contrast image is referred to as a virtual image.
  • the image calculation unit 343 performs an operation that the operator designates via the user interface 32 on the actual measurement image transmitted from the imaging unit 341 and the virtual image transmitted from the density length calculation unit 342, and the result is calculated by the user. Transmit to the interface 32. This calculation includes the one that directly handles the transmitted image as the calculation result.
  • FIG. 3 is a diagram showing an example of the screen of the user interface 32.
  • the screen 500 includes a mouse pointer 501, an image display unit 502, an object list display unit 503, a display mode selection unit 504, a virtual image mode selection unit 505, an actual image time designation unit 506, and a variable state quantity designation unit 507.
  • Various operations by an operator are realized in cooperation with an input unit 321 such as a mouse.
  • the display mode selection unit 504 shows four choices “target selection”, “measured image display”, “virtual image display”, and “comparison display” as display modes, and the operator can select one of them.
  • the image display unit 502 displays a drawing of the large-scale facility 90.
  • This drawing only needs to be a drawing that can determine the state of equipment such as a plan view and a three-dimensional view.
  • a total of four types of drawings, a plan view and a three-dimensional view viewed from three directions, are displayed side by side.
  • a plurality of drawings may be combined.
  • the three-dimensional view may be displayed in the form of a cross-sectional view cut out by an arbitrary plane.
  • a measurement target flag is also held in the facility information storage unit 333 as attribute information of each object constituting the large-scale facility 90.
  • This measurement target flag is a binary state variable and indicates whether the object is a measurement target or a non-measurement target.
  • the object to be measured is highlighted in the form of drawing in a specific color or blinking on the drawing, so that the operator can easily recognize it.
  • identification names of individual objects constituting the large-scale facility 90 are displayed in the form of a list, and each object has a symbol added according to the measurement target flag, or the identification name is color-coded. Thus, the operator can easily recognize whether or not the measurement target.
  • the measurement target flag of the object can be changed.
  • the measurement target flag can also be changed by operating the object identification name with the mouse. This change result is registered in the facility information storage unit 333.
  • the image display unit 502 displays the actual image or virtual image transmitted from the image calculation unit 343, respectively.
  • the image display unit 502 displays a comparison image transmitted from the image calculation unit 343.
  • the comparison image is, for example, a color image obtained by superimposing the actual image as a blue contrast image and the virtual image as a red contrast image, or representing the difference between the actual image and the virtual image. It is an image.
  • the calculation between the actually measured image and the virtual image is performed by the image calculation unit 343.
  • the virtual image mode selection unit 505 shows three options “only target”, “only non-target”, and “all objects” as the virtual image mode, and the operator can select one of them.
  • the density length calculation unit 342 targets only the objects whose measurement target flag value is the measurement target among the objects held in the facility information storage unit 333. In addition, a line segment in which a virtual linear muon track corresponding to each pixel is cut out by each object is calculated.
  • the product of the density of the object is the density length of the object, and the sum of the density lengths of all the objects to be processed is the pixel value of the pixel. In this way, the density length is obtained for all the pixels, and a contrast image constructed by normalizing it is used as a virtual image.
  • the correspondence relationship between the virtual linear muon track and the pixel matches the correspondence relationship between the muon track and the pixel that are actually detected when the imaging unit 341 constructs the actual measurement image. , Predetermined. For an object whose position and shape can change within a certain range, or an object that can store an amount of fluid objects that can change within a certain range, such as displacement held in the equipment information storage unit 333 It is used when constructing the virtual image after referring to the value and calculating its shape as appropriate.
  • the density length calculation unit 342 has a measurement target flag value that is a non-measurement target among the objects held in the facility information storage unit 333.
  • a virtual image is constructed for all objects and all objects held in the terrain information storage unit 332. The specific construction method is the same as when the virtual image mode is “target only”.
  • the density length calculation unit 342 targets all objects held in the facility information storage unit 333 and all objects held in the terrain information storage unit 332. Then, build a virtual image.
  • the specific construction method is the same as when the virtual image mode is “target only”.
  • the virtual image mode By setting the virtual image mode to “non-target only” and the display mode to “comparison display”, the influence of the non-measurement target object is eliminated, and a contrast image of only the measurement target object can be provided.
  • contrast information on the difference between the terrain information and facility information stored in the storage unit and the actual object state is contrasted. It can be provided in the form of an image, which can be used to detect construction defects and secular changes.
  • the measured image time designation unit 506 has an input field for designating the period of track information used when the imaging unit 341 constructs the measured image. For example, the operator operates the keyboard using the start date and the end date. Etc. can be input.
  • the variable state quantity designation unit 507 has a column for designating values such as displacement of each object, which is used when the density length calculation unit 342 constructs a virtual image.
  • a sliding door named “sliding door 89D3”
  • the amount of “displacement” can be specified to a value of “4.54 m” by keyboard operation or the like.
  • a non-destructive test for the purpose of soundness inspection of piping or the like constituting a large structure such as a power plant is performed without removing / recovering the covering material of the object to be measured. It is possible to obtain a contrast image reflecting the actual density length of only a desired measurement target object without installing or removing the source.
  • the present invention it is possible to obtain an image showing a difference between an actual state of a desired measurement target object and an ideal design state. As a result, the soundness of the plant can be constantly monitored during operation, and a plant with a high operation rate and high safety can be provided.

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

　宇宙線による透過撮影を利用した大規模設備監視において、測定対象のみの濃淡画像を算出できる計算システムを提供する。宇宙線の飛跡に基づく濃淡画像のほかに、測定対象でない構造物の構造情報に基づく密度長の濃淡画像を構築し、宇宙線の飛跡に基づく濃淡画像の補正に利用する。

Description

計算システム及び計算方法
　本発明は、構造物の構造を計算する計算システム及び計算方法に関するものである。
　従来、発電プラント等の大型構造物の健全性を検査する目的で、その構造物を構成する配管や支持部材の破壊試験および非破壊試験が実施されている。このうち、非破壊試験は、破壊試験と比較して、検査をすることによる健全性の低下の可能性が低いとう長所を有する。また、試験実施後の復元工程が不要であるという長所を備えており、様々な構造物に適用されている。例えば、工業規格に適合した金属配管の保守検査の一要素として、超音波パルスの反射を利用し、測定対象物体部位の厚さを測定することで配管内部の損耗の度合いを明らかにする手法が用いられる。だたしこの手法は、測定対象物体に超音波振動子を直接接触するため、配管の被覆材を測定前に除去し、測定後に復元する必要があり、実際の測定の前後に多大のコストを要するという欠点がある。
　そこで、特許文献１に示されるように、放射線透過撮影を行うことで、配管の被服材を除去ならびに復元する工程を省略する手法も用いられている。この手法では、人工放射線源と検出器とを、測定対象物体を挟む位置にそれぞれ設置する。
　また、特許文献２には、自然放射線の一種である高エネルギー二次宇宙線のミュオンを利用し、製鉄用高炉の炉壁や炉底を透過撮影することで、炉の耐火物の損耗の度合いを明らかにする手法が示されている。
　さらに、特許文献３には、二次宇宙線ミュオンが測定対象を透過する際の性質を利用して、測定対象の材質の変化の度合いを明らかにする手法が示されている。
特許第４８１４９１８号 特開２００７－１２１２０２号公報 特開２０１１－１２３０４８号公報
　特許文献１に示される手法では、前記の被覆材除去復元工程を省略できるものの、人工放射線の透過能力が弱いため、測定対象物体の近傍に放射線源を配置する必要がある。大規模なプラントにおいては、全ての測定対象物体近傍に人工線源を常設することは非現実的であり、測定の都度、個々の測定対象物体に線源と検出器の設置と撤去を行う必要があり、実際の測定の前後に多大のコストを要するという欠点がある。
　特許文献２および特許文献３に示される手法では、極めて透過性・直進性の高い粒子を用いてその飛行経路が貫く物体の密度長をコントラスト像として提供おり、特許文献１のように必ずしも測定対象物体近傍に検出器を配置する必要は無い。しかし多くの場合、撮影された透過画像は、測定対象物体のみならず、粒子の飛行経路が貫く他の物体をも合わせた密度長がコントラスト像に反映されてしまう欠点がある。
　以上のような従来技術の欠点（測定対象物体の被覆材を除去・復元する、人工放射線源を設置・撤去する）を考慮し、本発明は、大規模構造物またはその一部の健全性の検査を目的とする非破壊試験において、所望の測定対象物体のみの実際の密度長を反映するコントラスト像を得る装置または方法を提供することを目的とする。
　上記目的を達成するため、本願で開示する発明の概要を説明すれば以下の通りである。
　本願発明に記載の計算システムは、測定対象の構造物と他の構造物とを透過する宇宙線により任意の区間の構造情報を撮像する撮像部と、前記他の構造物の構造情報を有する構造情報取得部と、前記他の構造物の構造情報を測定対象として含めるか判断し、非対象として判断した構造の密度長を算出する密度長算出部と、前記撮像部により撮像された画像と前記非対象として判断した構造の密度長とに基づき測定対象の構造物の構造情報を計算する画像演算部と、を有する。
　また、本願発明に記載の計算システムは、測定対象と他の構造物とを透過する宇宙線の飛跡を取得する飛跡取得部と、前記宇宙線の飛跡に基づいて撮像する撮像部と、前記宇宙線の飛跡の情報を保持する第１の記憶部と、前記測定対象ならびに前記他の構造物の構造情報を有する第２の記憶部と、前記他の構造物の構造情報に基づいて第１の密度長の濃淡像を算出し、かつ前記測定対象の構造情報ならびに前記他の構造物の構造情報に基づいて第２の密度長の濃淡像を算出する密度長算出部と、を有する。
　また、本願発明に記載の計算方法は、測定対象の構造物と他の構造物とを透過する宇宙線により任意の区間の構造情報を撮像する撮像する第１のステップと、前記他の構造物の構造情報を取得する第２のステップと、前記他の構造物の構造情報を測定対象として含めるか判断する第３のステップと、前記第３のステップにて非対象として判断した構造の密度長を算出する第４のステップと、前記第１のステップにより撮像された画像と前記第４のステップにて算出した密度長とに基づき測定対象の構造物の構造情報を計算する第５のステップと、を有する。
　本発明によれば、所望の測定対象物体に対する好適な非破壊検査が可能となる。
本発明においてミュオンの飛行経路上に存在する様々な物体を例示する概念図である。 本発明の一実施形態の大規模設備状態監視システムの構成を示す図である。 本発明の一実施形態のユーザーインターフェース画面の例を示す図である。
　図１は、発明に先立ち、本願の発明者が検討した図面である。ミュオンの飛行経路である矢印Ｍ１の経路には、様々な非測定対象物体が存在する状況を模式的に例示している。例えば、大規模設備９０を構成する配管等の測定対象物体９１をミュオン検出器９９により透過撮影することを目的とする場合に、測定対象物体９１を囲む建屋９２の壁や、測定対象物体９１以外の配管９３、建屋９２の外部の山９４、他の建築物９５も、飛行経路Ｍ１が貫いていることがわかる。
　言い換えると、二次宇宙線ミュオンは地球大気上層という遠方で生成され、検出器に入射するまでにミュオンが飛行する距離は１０ｋｍを上回ることも珍しくない。よって、その飛行経路には、測定対象物体を囲む建屋ならびにその外部の建築物を構成する壁や柱など、もしくは、山や丘などの地形を構成する土砂などの、様々な材質と形状の非測定対象物体が存在することが考えられる。
　しかし、現実にはこれらの非測定対象物体が測定時においてどのような変化をするか、一つ一つを把握することは難しい。従って、求めたコントラスト像は、本来測定すべき物体に起因する情報以外にも、非測定対象物体の情報を含んでしまうことになる。
　そこで、以下、本発明の実施の形態について図面を用いて説明する。
　図２は、本発明の一実施形態の大規模設備状態監視システムの構成を示す図である。
　本実施形態の大規模設備状態監視システム１０は、飛跡取得部２０と、情報処理部３０を構成要素として含む。
　まず、飛跡取得部２０を以下に説明する。
　飛跡取得部２０は、電源部２１と、検出部２２と、同時計数部２３を構成要素として含む。
　電源部２１は、検出部２２と同時計数部２３に、必要な電力を供給する電源部である。
　検出部２２は、第１検出器２２１、鉄塊２２２と、第２検出器２２３を構成要素として含む。
　第１検出器２２１ならびに第２検出器２２３はいずれも位置敏感型検出器である。ここでいう位置敏感型検出器とは、荷電粒子が通過した際に、その検出器のおける通過点の座標値が検出可能なものである。例えば、複数の検出素子(シンチレータ)をマトリクスになるように平面状に配置し、特定の検出素子において検出が行われたことを時間と共に検出する検出器が挙げられる。
　また、第１検出器２２１ならびに第２検出器２２３は各々を荷電粒子が通過した際に、一定の極めて短い遅延時間の後にその通過点の座標値を同時計数部２３に送信する。第１検出器２２１ならびに第２検出器２２３の配置関係は、上述した平面状の検出器が互いに平行に配置される関係であるのが好適である。第1の検出器の平面に対して第２の検出器の平面が傾斜関係にあると、測定対象物及び第１検出器２２１を透過した荷電粒子を効率よく捉えることができなくなり、それを防止するために第１検出器２２１よりも大きな平面を必要とすることも考えられる。従って、それぞれの平面が傾斜関係から平行関係になるほど荷電粒子の検出効率を上げることが出来る。
　また、検出部２２には、様々な方位角・天頂角の方向から、不規則なタイミングで荷電粒子が入射する。これら荷電粒子は高エネルギーのミュオンのみならず電子等の軟成分も含まれる。軟成分はミュオンよりも物体中における直進性が極めて劣っており、測定対象物体の正確な密度長に基づくコントラスト像を得るためには、直進性の高いミュオンのみを選択的に検出する必要がある。
　よって、ミュオンのみを選択的に検出する方法を以下に述べる。まず、第１検出器２２１を通過し鉄塊２２２に入射した荷電粒子がミュオンである場合は、第２検出器２２３は単一の座標値を同時計数部２３に送信することができる。これは、そのミュオンがそのまま鉄塊２２２を通過し、第２検出器２２３に入射するためである。
　それに対し、第１検出器２２１を通過し鉄塊２２２に入射した荷電粒子が軟成分である場合は、第２検出器２２３がほぼ同時に複数の座標値を同時計数部２３に送信することになる。これは、その粒子は鉄塊２２２を通過する過程で多数の粒子を発生させ、それらの粒子の進行方向は様々であるためである。
　以上のように、第２検出器２２３が単一の座標値を検出する場合のみを選択的に採用することにより、ミュオンのみを選択的に検出し軟成分の影響を排除できる。
　同時計数部２３は、第１検出器２２１ならびに第２検出器２２３から送信されたそれぞれの時刻を比較し、その差がある一定の範囲内に収まる場合にのみ、同一のミュオンが通過したものと判定する。そして、判定されたこれら２組の座標値を、情報処理部３０に送信する。後述するように、同時と判定されたこれら２組の座標値から、ミュオンの軌跡を求めることができる。
　次に、情報処理部３０を以下に説明する。
　情報処理部３０は、電源部３１と、ユーザーインターフェース３２と、記憶部３３と、計算部３４を構成要素として含む。
　電源部３１は、ユーザーインターフェース３２と、記憶部３３と、計算部３４に、必要な電力を供給する電源部である。
　ユーザーインターフェース３２は、キーボードやマウス等の入力部３２１と、ディスプレイ装置やプリンター等の出力部３２２を構成要素として含む。ユーザーインターフェース３２を介し、以下の操作を操作者は行うことができる。例えば、計算部３４が提示する選択肢の中から所望の１個以上の項目を選択する操作、計算部３４が図示する幾何形状の一部範囲を選択する操作、計算部３４が提示する有限の範囲の数値の中から特定の数値を選択する操作、計算部３４が提示する変数に代入すべき値を入力する操作、などが挙げられる。
　記憶部３３は、飛跡情報記憶部３３１と、地形情報記憶部３３２と、設備情報記憶部３３３を構成要素として含む。
　飛跡情報記憶部３３１は、ミュオン通過の度に同時計数部２３から送信される２組のミュオン通過座標データに、時刻情報を組み合わせて飛跡（軌跡）情報とした上で、それを新たに記憶する。
　地形情報記憶部３３２は、測定対象物体を構成要素とする大規模設備９０周辺の地形を示す幾何学的データと、その地形を構成する地形物の材質または密度を示すデータを保持する。ここで地形物とは、天然に形成された山や丘陵等のみならず、人工的に形成されたダム、隧道、橋梁、ビル等の構造物といった構成要素を含む。
　設備情報記憶部３３３は、大規模設備９０を構成する壁、柱、梁等の静的な建築部材のほか、扉や起重機（例えば設置型クレーン）等の、一定の範囲内で位置や形状が変化しうる物体や、タンク等の、一定の範囲内で変化しうる分量の流動性物体を蓄えうる物体の、幾何学的形状データと、材質または密度を示すデータを保持する。また、設備情報記憶部３３２は、前記の一定の範囲内で変化しうる各種の量の値ならびにその量の変化範囲を示すデータと、前記の流動性物体の材質または密度を示すデータを保持する。
　計算部３４は、撮像部３４１と、密度長算出部３４２と、画像演算部３４３を構成要素として含む。
　撮像部３４１は、飛跡情報記憶部３３１に保持された多数の飛跡情報のうち、操作者がユーザーインターフェース３２を介して指定する期間に該当する時刻情報を持つ飛跡情報を選び出し、それらのミュオン通過座標値データを基に測定対象物のコントラスト像を構築する。以後、このコントラスト像を実測像と表記する。撮像部３４１は実像データを画像演算部３４３に送信する。
　密度長算出部３４２は、地形情報記憶部３３２ならびに設備情報記憶部３３３に保持された各種データを基に、後述の方法により仮想的なコントラスト像を構築する。以後、このコントラスト像を仮想像と表記する。
　画像演算部３４３は、撮像部３４１から送信される実測像と、密度長算出部３４２から送信される仮想像に対し、操作者がユーザーインターフェース３２を介して指定する演算を行い、その結果をユーザーインターフェース３２に送信する。この演算には、送信された像をそのまま演算結果として扱うものも含まれる。
　図３は、ユーザーインターフェース３２の画面の例を示す図である。画面５００は、マウスポインタ５０１、画像表示部５０２、物体一覧表示部５０３、表示モード選択部５０４、仮想像モード選択部５０５、実測像時間指定部５０６、可変状態量指定部５０７を備え、キーボードやマウス等の入力部３２１と連携して操作者による各種操作を実現する。
　表示モード選択部５０４には表示モードとして４つの選択肢「対象選択」、「実測像表示」、「仮想像表示」、「比較表示」が示され、そのうち１つを操作者は選択できる。
　選択されている表示モードが「対象選択」の場合、画像表示部５０２には、大規模設備９０の図面が表示される。この図面は、平面図、立体図等、設備の状況が判別できる図面であればよく、また、例えば３方向から見た平面図と立体図の計４種類の図面が並べて表示されるように、複数にの図面を組み合わせてもよい。さらに、立体図については、任意の平面により切り取られた断面図の形で表示されてもよい。
　これらの図面データは、設備情報記憶部３３３に保持されている。大規模設備９０を構成する個々の物体の属性情報として、測定対象フラグも、設備情報記憶部３３３に保持されている。この測定対象フラグは、二値の状態変数であり、その物体が測定対象であるか、非測定対象であるかを示す。
　測定対象である物体は、図面上で特定の色に描画されたり、点滅表示されるなどの形で、ハイライト表示され、操作者は容易に認識できる。物体一覧表示部５０３には、大規模設備９０を構成する個々の物体の識別名称が一覧表の形で表示され、各物体が測定対象フラグに応じて、シンボルが付記されたり、識別名称が色分けされ、操作者は容易に測定対象か否かを認識できる。
　そして、操作者がマウス操作により、図面に記載された配管などの物体を選択することで、その物体の測定対象フラグを変更できる。また、物体の識別名称をマウス操作することによっても、測定対象フラグを変更できる。この変更結果は、設備情報記憶部３３３に登録される。
　選択されている表示モードが「実測像表示」または「仮想像表示」の場合、画像表示部５０２には、それぞれ、画像演算部３４３から送信される実測像または仮想像が表示される。
　選択されている表示モードが「比較表示」の場合、画像表示部５０２には、画像演算部３４３から送信される比較画像が表示される。ここで比較画像とは、例えば、実測像を青色のコントラスト像とし、仮想像を赤色のコントラスト像とした上で、両者を重ね合わせたカラー画像であったり、実測像と仮想像の差分を表す画像である。これら実測像と仮想像の間の演算は、画像演算部３４３が行う。
　仮想像モード選択部５０５には、仮想像モードとして３つの選択肢「対象のみ」、「非対象のみ」、「全物体」が示され、そのうち１つを操作者は選択できる。
　選択されている仮想像モードが「対象のみ」の場合、密度長算出部３４２は、設備情報記憶部３３３に保持されている物体のうち、測定対象フラグの値が測定対象である物体のみを対象に、各画素に対応する仮想的な直線状のミュオン飛跡が個々の物体により切り取られる線分を算出する。そしてその物体の密度との積を、その物体の密度長とし、対象とされる全ての物体の密度長の和を、その画素の画素値とする。こうして全画素について密度長を求め、それを例えば正規化することで構築されたコントラスト像を仮想像とする。
　ここで、仮想的な直線状のミュオン飛跡は、それと画素との対応関係が、撮像部３４１が実測像を構築する際の実際に検出されたミュオン飛跡と画素との対応関係に一致するように、予め定める。なお、一定の範囲内で位置や形状が変化しうる物体や、一定の範囲内で変化しうる分量の流動性物体を蓄えうる物体については、設備情報記憶部３３３に保持されている変位などの値を参照し、適宜その形状を算出した上で、上記の仮想像構築の際に用いられる。
　また、選択されている仮想像モードが「非対象のみ」の場合、密度長算出部３４２は、設備情報記憶部３３３に保持されている物体のうち、測定対象フラグの値が非測定対象である物体すべてと、地形情報記憶部３３２に保持されている物体すべてを対象に、仮想像を構築する。具体的な構築方法は、仮想像モードが「対象のみ」の場合と同じである。
　選択されている仮想像モードが「全物体」の場合、密度長算出部３４２は、設備情報記憶部３３３に保持されている物体すべてと、地形情報記憶部３３２に保持されている物体すべてを対象に、仮想像を構築する。具体的な構築方法は、仮想像モードが「対象のみ」の場合と同じである。
　仮想像モードを「非対象のみ」とし、表示モードを「比較表示」とすることで、非測定対象物体の影響が排除され、測定対象の物体のみのコントラスト像を提供できる。また、仮想像モードを「全物体」とし、表示モードを「比較表示」とすることで、記憶部に保持されている地形情報ならびに設備情報と、実際の物体の状態との差異に関する情報をコントラスト像の形で提供でき、施工不良や経年変化などの発見に役立てることができる。
　実測像時間指定部５０６には、撮像部３４１が実測像を構築する際に用いる飛跡情報の期間を指定するための入力欄があり、例えば開始日と終了日の形で、操作者がキーボード操作等により入力できる。
　可変状態量指定部５０７には、密度長算出部３４２が仮想像を構築する際に用いる、各物体の変位等の値を指定するための欄があり、例えば、「引戸89D3」という名称の引戸の「変位」という量を、「4.54m」という値に、キーボード操作等により指定できる。
　最後に、本発明によれば、発電プラント等の大型構造物を構成する配管等の健全性検査を目的とする非破壊試験を、測定対象物体の被覆材を除去・復元することなく、人工放射線源を設置・撤去することなく、かつ、所望の測定対象物体のみの実際の密度長を反映するコントラスト像を得ることができる。
　また、本発明によれば、所望の測定対象物体の実際の状態と、設計上の理想的な状態との差異を示す像を得ることができる。その結果、プラントの健全性を操業時に常時監視することが可能となり、稼働率と安全性の高いプラントを提供できるようになる。
１０　大規模設備状態監視システム
２０　飛跡取得部
３０　情報処理部
９０　大規模設備
９１　測定対象物体

Claims (9)

  1. 　測定対象の構造物と他の構造物とを透過する宇宙線により任意の区間の構造情報を撮像する撮像部と、
    　前記他の構造物の構造情報を有する構造情報取得部と、
    　前記他の構造物の構造情報を測定対象として含めるか判断し、非対象として判断した構造の密度長を算出する密度長算出部と、
    　前記撮像部により撮像された画像と前記非対象として判断した構造の密度長とに基づき測定対象の構造物の構造情報を計算する画像演算部と、を有する計算システム。
  2. 　請求項１において、
    　前記測定対象の構造物の構造情報は、前記測定対象の構造物の密度長であることを特徴とする計算システム。
  3. 　請求項１において、
    　前記画像演算部は前記測定対象の構造物の構造情報により、前記測定対象の構造物の損耗を評価することを特徴とする計算システム。
  4. 　測定対象と他の構造物とを透過する宇宙線の飛跡を取得する飛跡取得部と、
    　前記宇宙線の飛跡に基づいて撮像する撮像部と、
    　前記宇宙線の飛跡の情報を保持する第１の記憶部と、前記測定対象ならびに前記他の構造物の構造情報を有する第２の記憶部と、
    　前記他の構造物の構造情報に基づいて第１の密度長の濃淡像を算出し、かつ前記測定対象の構造情報ならびに前記他の構造物の構造情報に基づいて第２の密度長の濃淡像を算出する密度長算出部と、を有する計算システム。
  5. 　請求項４において、
    　前記計算システムは、前記撮像部により撮像した画像と前記第１の密度長の濃淡像とに基づいて前記撮像した画像を補正し、前記撮像した画像と前記第２の密度長の濃淡像との差異を示す濃淡像を演算する画像演算部と、
    　を有する計算システム。
  6. 　請求項４において、
    　前記飛跡取得部は測定対象を挟むように配置された第１乃至第２検出器と、前記第１検出器ならびに前記第２検出器に接続され検出された宇宙線の強度を係数する同時計数部と、を有する計算システムであって、
    　前記同時計数部は、前記宇宙線がミュオンであるか判断し、ミュオンと判断した宇宙線の飛跡のみを取得することを特徴とする計算システム。
  7. 　測定対象の構造物と他の構造物とを透過する宇宙線により任意の区間の構造情報を撮像する撮像する第１のステップと、
    　前記他の構造物の構造情報を取得する第２のステップと、
    　前記他の構造物の構造情報を測定対象として含めるか判断する第３のステップと、
    　前記第３のステップにて非対象として判断した構造の密度長を算出する第４のステップと、
    　前記第１のステップにより撮像された画像と前記第４のステップにて算出した密度長とに基づき測定対象の構造物の構造情報を計算する第５のステップと、を有する計算方法。
  8. 　請求項７において、
    　前記測定対象の構造物の構造情報は、前記測定対象の構造物の密度長であることを特徴とする計算方法。
  9. 　請求項７において、
    　前記測定対象の構造物の構造情報により、前記測定対象の構造物の損耗を評価する第６のステップと、を有する計算方法。
PCT/JP2012/081445 2012-12-05 2012-12-05 計算システム及び計算方法 WO2014087494A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
PCT/JP2012/081445 WO2014087494A1 (ja) 2012-12-05 2012-12-05 計算システム及び計算方法
US14/649,924 US20150293040A1 (en) 2012-12-05 2012-12-05 Calculation system and calculation method
JP2014550839A JPWO2014087494A1 (ja) 2012-12-05 2012-12-05 計算システム

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2012/081445 WO2014087494A1 (ja) 2012-12-05 2012-12-05 計算システム及び計算方法

Publications (1)

Publication Number Publication Date
WO2014087494A1 true WO2014087494A1 (ja) 2014-06-12

Family

ID=50882945

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2012/081445 WO2014087494A1 (ja) 2012-12-05 2012-12-05 計算システム及び計算方法

Country Status (3)

Country Link
US (1) US20150293040A1 (ja)
JP (1) JPWO2014087494A1 (ja)
WO (1) WO2014087494A1 (ja)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017015413A (ja) * 2015-06-28 2017-01-19 株式会社サイエンスインパクト 放射線計算装置、放射線計算方法、放射線計算プログラム
JP2017044642A (ja) * 2015-08-28 2017-03-02 株式会社東芝 構造物の検査装置及びその検査方法
JP2022047821A (ja) * 2020-09-14 2022-03-25 株式会社東芝 非破壊物質組成識別装置および非破壊物質組成識別方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI639102B (zh) * 2016-08-10 2018-10-21 張雅如 一種指標顯示裝置、指標控制裝置、指標控制系統及其相關方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08166359A (ja) * 1994-12-15 1996-06-25 Hitachi Medical Corp ごみ処理施設における処理不適物検査システム
JP2011257431A (ja) * 1997-12-18 2011-12-22 Siemens Ag モジュール構造のイメージング装置

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7633062B2 (en) * 2006-10-27 2009-12-15 Los Alamos National Security, Llc Radiation portal monitor system and method
US8288721B2 (en) * 2007-04-23 2012-10-16 Decision Sciences International Corporation Imaging and sensing based on muon tomography
WO2009107575A1 (ja) * 2008-02-25 2009-09-03 大学共同利用機関法人高エネルギー加速器研究機構 複合構造物の非破壊検査装置及び非破壊検査方法
RU2503953C2 (ru) * 2008-08-27 2014-01-10 ЛОС АЛАМОС НЭШНЛ СЕКЬЮРИТИ, ЭлЭлСи Построение изображения по заряженным частицам, создаваемым космическими лучами
WO2012167381A1 (en) * 2011-06-07 2012-12-13 Atomic Energy Of Canada Limited Detecting high atomic number materials with cosmic ray muon tomography

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08166359A (ja) * 1994-12-15 1996-06-25 Hitachi Medical Corp ごみ処理施設における処理不適物検査システム
JP2011257431A (ja) * 1997-12-18 2011-12-22 Siemens Ag モジュール構造のイメージング装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
AKIHIKO SHINOTAKE: "Observation for the Inner Structure of Blast Furnace by Cosmic-Ray Muon Radiography", BULLETIN OF THE IRON AND STEEL INSTITUTE OF JAPAN (FERRUM), vol. 14, no. 11, 1 November 2009 (2009-11-01), pages 700 - 705 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2017015413A (ja) * 2015-06-28 2017-01-19 株式会社サイエンスインパクト 放射線計算装置、放射線計算方法、放射線計算プログラム
JP2017044642A (ja) * 2015-08-28 2017-03-02 株式会社東芝 構造物の検査装置及びその検査方法
JP2022047821A (ja) * 2020-09-14 2022-03-25 株式会社東芝 非破壊物質組成識別装置および非破壊物質組成識別方法
JP7476058B2 (ja) 2020-09-14 2024-04-30 株式会社東芝 非破壊物質組成識別装置および非破壊物質組成識別方法

Also Published As

Publication number Publication date
US20150293040A1 (en) 2015-10-15
JPWO2014087494A1 (ja) 2017-01-05

Similar Documents

Publication Publication Date Title
US10551280B2 (en) Arrangement and method for inspecting an object, in particular a building
Farahani et al. A coupled 3D laser scanning and digital image correlation system for geometry acquisition and deformation monitoring of a railway tunnel
US20170052150A1 (en) System and Method for Progressive Damage Monitoring and Failure Event Prediction in a Composite Structure
US20170122909A1 (en) Non-destructive system and method for detecting structural defects
AU2015258535B2 (en) Electromagnetic radiation detection apparatus for hidden fire danger in mine, and method therefor
US10190992B2 (en) Structure status determination device, status determination system, and status determination method
US10571404B2 (en) Defect inspection apparatus, method, and program
Sabato et al. Noncontact sensing techniques for AI-aided structural health monitoring: a systematic review
JP2008232998A (ja) 構造物の応力変動分布の測定方法およびその測定装置、ならびに構造物の欠陥検出方法および構造物の危険性把握方法
JP2020510198A (ja) パイプ検査のための拡張現実可視化
CN103075641A (zh) 非接触式管道磁检测方法
WO2014087494A1 (ja) 計算システム及び計算方法
JP2019194605A (ja) 対象物、特に建造物を検査するための構成および方法
Chen et al. Computer vision application programming for settlement monitoring in a drainage tunnel
CN1898698B (zh) 控制核反应堆元件的非破坏性方法
JP2017219386A (ja) 放射線配管診断システム
Benton et al. Close-range photogrammetry in underground mining ground control
KR100967084B1 (ko) 균열 모니터링 시스템, 균열 모니터링 방법 및 균열모니터링 프로그램을 기록한 컴퓨터로 읽을 수 있는 매체
Minichan et al. H-canyon air exhaust tunnel inspection vehicle development
Yasuda et al. Fatigue crack detection system based on IoT and statistical analysis
JP2013213748A (ja) 内部構造検査装置および内部構造監視システム
JP7312864B2 (ja) 対象物、特に建造物を検査するための構成および方法
Shi et al. Leveraging point cloud data for detecting building Façade deteriorations caused by neighboring construction
Fahle et al. Evaluation of SLAM-based mobile lidar workflows and data quality for mine-wide underground geotechnical monitoring
JP2019152497A (ja) ミュオグラフィを用いる構造物の非破壊検査方法及び装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 12889681

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2014550839

Country of ref document: JP

Kind code of ref document: A

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 14649924

Country of ref document: US

122 Ep: pct application non-entry in european phase

Ref document number: 12889681

Country of ref document: EP

Kind code of ref document: A1