• 제목/요약/키워드: model Optimization

검색결과 5,633건 처리시간 0.035초

반응 표면 분석법을 사용한 Rhodobacter sphaeroides PS-24 유래 carotenoid 생산 배지 최적화 (Optimization of Medium for the Carotenoid Production by Rhodobacter sphaeroides PS-24 Using Response Surface Methodology)

  • 봉기문;김공민;서민경;한지희;박인철;이철원;김평일
    • 한국유기농업학회지
    • /
    • 제25권1호
    • /
    • pp.135-148
    • /
    • 2017
  • 본 연구를 통해 논, 시설재배 밭 토양, 쓰레기장, 하천 및 호수의 퇴적 토양 등 22개소에서 분리한 총 6종의 광합성세균 중 호기 암 배양이 가능한 Rhodobacter sphaeroides PS-24를 분리하였다. 형태학적 특징으로는 그람음성의 막대모양으로, 운동성이 있었다. 분리균주의 16S rRNA 염기서열을 분석한 결과 Rhodobacter sphaeroides ATH2.4.1과 99%의 상동성을 나타내었으며, 본 연구에서 Rhodobacter spharoides PS-24로 명명하여 연구를 수행하였다. 선별균주를 modifed Van niel's yeast 배지에서 배양 후 생성된 carotenoid를 추출한 결과 $12.03{\pm}0.15mg/L$의 함량이 측정되었으며, 반응표면분석법 중 Plackett burman 분석방법과 Box-Behnken 분석방법을 통해 carotenoid 생산에 영향을 미치는 요인을 분석하고 농도를 최적화하였다. 분석결과 각각의 독립변수 yeast extract -0.4144 (1.23 g/L), $Na_2CO_3$ 0.8541 (3.71 g/L)와 $MgSO_4$ 1.00 (1.00 g/L)의 농도를 선정하였으며, 이를 바탕으로 배지 조성을 최적화한 결과 yeast extract 1.23 g, sodium acetate 1 g, $NH_4Cl$ 1.75 g, NaCl 2.5 g, $K_2HPO_4$ 2 g, $MgSO_4$ 1.0 g, mono-sodium glutamate 7.5 g, $Na_2CO_3$ 3.71 g, $NH_4Cl$ 3.5 g, $CaCl_2$ 0.01 g/ liter로 선정하였다. 최적배지를 대상으로 5 L, 50 L, 500 L scale-up을 진행한 결과 최종 carotenoid는 각각 17.98 mg/L, 18.03 mg/L, 18.11 mg/L로 조사되었다. 최적배지의 경우 modified Van niel's yeast 배지보다 약 1.5배 많은 carotenoid를 생산하였으며, 대량배양을 통한 scale-up 과정 시 carotenoid의 생산량은 크게 변화하지 않는 것으로 조사되었다. 따라서 본 연구를 바탕으로 산업적으로 다양하게 사용되고 있는 carotenoid를 생산하는 광합성세균 Rhodobacter spharoides PS-24를 개발하였으며, 본 연구를 바탕으로 유기농축산에 사용이 가능한 기능성 미생물제제를 개발하고자 한다.

A hybrid algorithm for the synthesis of computer-generated holograms

  • Nguyen The Anh;An Jun Won;Choe Jae Gwang;Kim Nam
    • 한국광학회:학술대회논문집
    • /
    • 한국광학회 2003년도 하계학술발표회
    • /
    • pp.60-61
    • /
    • 2003
  • A new approach to reduce the computation time of genetic algorithm (GA) for making binary phase holograms is described. Synthesized holograms having diffraction efficiency of 75.8% and uniformity of 5.8% are proven in computer simulation and experimentally demonstrated. Recently, computer-generated holograms (CGHs) having high diffraction efficiency and flexibility of design have been widely developed in many applications such as optical information processing, optical computing, optical interconnection, etc. Among proposed optimization methods, GA has become popular due to its capability of reaching nearly global. However, there exits a drawback to consider when we use the genetic algorithm. It is the large amount of computation time to construct desired holograms. One of the major reasons that the GA' s operation may be time intensive results from the expense of computing the cost function that must Fourier transform the parameters encoded on the hologram into the fitness value. In trying to remedy this drawback, Artificial Neural Network (ANN) has been put forward, allowing CGHs to be created easily and quickly (1), but the quality of reconstructed images is not high enough to use in applications of high preciseness. For that, we are in attempt to find a new approach of combiningthe good properties and performance of both the GA and ANN to make CGHs of high diffraction efficiency in a short time. The optimization of CGH using the genetic algorithm is merely a process of iteration, including selection, crossover, and mutation operators [2]. It is worth noting that the evaluation of the cost function with the aim of selecting better holograms plays an important role in the implementation of the GA. However, this evaluation process wastes much time for Fourier transforming the encoded parameters on the hologram into the value to be solved. Depending on the speed of computer, this process can even last up to ten minutes. It will be more effective if instead of merely generating random holograms in the initial process, a set of approximately desired holograms is employed. By doing so, the initial population will contain less trial holograms equivalent to the reduction of the computation time of GA's. Accordingly, a hybrid algorithm that utilizes a trained neural network to initiate the GA's procedure is proposed. Consequently, the initial population contains less random holograms and is compensated by approximately desired holograms. Figure 1 is the flowchart of the hybrid algorithm in comparison with the classical GA. The procedure of synthesizing a hologram on computer is divided into two steps. First the simulation of holograms based on ANN method [1] to acquire approximately desired holograms is carried. With a teaching data set of 9 characters obtained from the classical GA, the number of layer is 3, the number of hidden node is 100, learning rate is 0.3, and momentum is 0.5, the artificial neural network trained enables us to attain the approximately desired holograms, which are fairly good agreement with what we suggested in the theory. The second step, effect of several parameters on the operation of the hybrid algorithm is investigated. In principle, the operation of the hybrid algorithm and GA are the same except the modification of the initial step. Hence, the verified results in Ref [2] of the parameters such as the probability of crossover and mutation, the tournament size, and the crossover block size are remained unchanged, beside of the reduced population size. The reconstructed image of 76.4% diffraction efficiency and 5.4% uniformity is achieved when the population size is 30, the iteration number is 2000, the probability of crossover is 0.75, and the probability of mutation is 0.001. A comparison between the hybrid algorithm and GA in term of diffraction efficiency and computation time is also evaluated as shown in Fig. 2. With a 66.7% reduction in computation time and a 2% increase in diffraction efficiency compared to the GA method, the hybrid algorithm demonstrates its efficient performance. In the optical experiment, the phase holograms were displayed on a programmable phase modulator (model XGA). Figures 3 are pictures of diffracted patterns of the letter "0" from the holograms generated using the hybrid algorithm. Diffraction efficiency of 75.8% and uniformity of 5.8% are measured. We see that the simulation and experiment results are fairly good agreement with each other. In this paper, Genetic Algorithm and Neural Network have been successfully combined in designing CGHs. This method gives a significant reduction in computation time compared to the GA method while still allowing holograms of high diffraction efficiency and uniformity to be achieved. This work was supported by No.mOl-2001-000-00324-0 (2002)) from the Korea Science & Engineering Foundation.

  • PDF

반응표면분석법을 이용한 Leuconostoc mesenteroides SRCM201425의 만니톨 생산배지 최적화 (Optimization of Medium Components using Response Surface Methodology for Cost-effective Mannitol Production by Leuconostoc mesenteroides SRCM201425)

  • 하광수;신수진;정성엽;양호연;임수아;허주희;양희종;정도연
    • 생명과학회지
    • /
    • 제29권8호
    • /
    • pp.861-870
    • /
    • 2019
  • 본 연구에서는 경제성 있는 생물학적 만니톨 고생산을 위해 선별균주의 만니톨 생산 최적화 배지조성을 RSM 방법을 이용하여 확립하고자 하였다. 먼저 김치로부터 분리된 10균주의 만니톨 생산량과 과당으로부터 만니톨 전환율 분석을 통하여 SRCM201425 균주를 선발하였으며, 선발균주는 16S rRNA 유전자 염기서열과 당 발효 분석을 통하여 Leuconostoc mesenteroides로 동정하였다. Plackett-Burman design (PBD)을 이용하여 만니톨 생산에 영향을 주는 배지 인자를 선별하기 위해 총 11개의 탄소원, 질소원, 무기원소의 영향을 조사하였으며, 통계학적 분석을 통하여 최종적으로 fructose와 sucrose, peptone을 선정하였다. 만니톨 생산을 위한 선별된 각 변수의 최적 농도를 결정하기 위한 방법으로 central composite design (CCD)과 반응표면 분석법을 이용하였으며, 최종적으로 CCD를 통해 만니톨 생산을 위한 배지 조성의 최적 농도는 fructose 38.68 g/l, sucrose 30 g/l, peptone 39.67 g/l으로 예측되었으며, 통계학적 분석을 통해 실험모델의 적합성을 확인하였다. 최종적으로 MRS 배지에서의 생산량의 약 20배의 만니톨을 생산할 수 있었으며, 100 g/l의 fructose가 포함된 MRS 배지 대비 97.46%의 만니톨을 생산하면서, 산업화시 기존 배지 대비 생산비용을 크게 절감할 수 있을 것으로 예상되었다. 본 연구를 통하여 만니톨 생산을 위한 배지 조성의 최적화를 확립하였으며, 만니톨을 생산하기 위한 방법으로 주로 사용되고 있는 고비용의 촉매환원 방법을 대체할 방법을 제시할 수 있는 기초자료로 활용될 수 있을 것으로 기대된다.

전자빔 공정에서 실험계획법을 이용한 살균제 Benomyl의 제거특성 및 독성평가 (Decomposition Characteristics of Fungicides(Benomyl) using a Design of Experiment(DOE) in an E-beam Process and Acute Toxicity Assessment)

  • 유승호;조일형;장순웅;이시진;천석영;김한래
    • 대한환경공학회지
    • /
    • 제30권9호
    • /
    • pp.955-960
    • /
    • 2008
  • 본 연구는 전자빔 공정에서 실험계획법(design of experiment: DOE) 중 일반요인배치법(general factorial design)을 이용하여 2개 인자(X$_1$: benomyl concentration(mg/L), X$_2$: E-beam irradiation(Gy))를 토대로 요인(X$_1$: benomyl concentration) 1에서 3개 수준(3 level: 0.5, 1 및 1.5 mg/L)와 요인(X$_2$: E-beam irradiation) 2에서 6개 수준(6 level: 100, 800, 600, 400, 200 및 100 Gy)으로 구성된 3블록(block) 실험조합에 따라 Benomyl의 분해(Y$_1$: the % of decomposition), 무기화(Y$_2$: the % of materialization) 및 독성평가(acute toxicity assessment)를 수행하였다. 우선 HPCL 분석에 의한 Benomyl에 분해특성은 처리조합(treatment combination) 3 블록(block)의 17 및 18번을 제외한 모든 실험조건에서 100% 분해되었고 등분산(equal variance) 조건에서 일원분산분석(one-way ANOVA)결과 수준 간 유의한 차이가 없었다(p > 0.05). 전자빔 조사에 의한 Benomyl에 무기화(materialization) 특성은 각 3개의 처리조합에서 평균 46%, 36.7% 및 22%의 제거효율을 나타났고 각 조합에서 처리수준 간 예측식은 block 1(Y$_1$ = 0.024X$_1$ + 34.1(R$^2$ = 0.929)), block 2(Y$_2$ = 0.026X$_2$ + 23.1(R$^2$ = 0.976)) 및 block 3(Y$_3$ = 0.034X$_3$ + 6.2(R$^2$ = 0.98)) 등의 1차 선형 회귀식을 만족하였다. 또한 Benomyl에 무기화(materialization)에 대한 Anderson-Darling 검정을 이용한 정규성(normality)을 만족하였다(p > 0.05). 또한 무기화에 대한 반응에 대한 선형 및 비선형을 포함한 다중회귀분석(multi regression analysis)을 도출한 결과 다음과 같은 예측식 Y = 39.96 - 9.36X$_1$ + 0.03X$_2$ - 10.67X$_1{^2}$ - 0.001X$_2{^2}$ + 0.011X$_1$X$_2$(R$^2$ = 96.3%, Adjusted R$^2$ = 94.8%)을 도출하였다. 2가지 반응변수(X$_1$: benomyl concentration(mg/L), X$_2$: E-beam irradiation(Gy))에 의한 2차 반응표면 모형식 추정으로부터 정준분석을 통해 최적조건을 도출한 결과 Benomyl 초기농도(X$_1$) 0.55 mg/L, 전자빔 조사량 950 Gy에서 TOC 제거율 57.3%으로 나타났다. 마지막으로 V. fischeri를 이용한 MicrotoxTM modified 81.9% test을 이용하여 전자빔에 의한 Benomyl에 대한 급성 독성을 평가한 결과 전자빔 조사전 block 1의 0.5 mg/L에서 10.25%, block 2의 1 mg/L에서 20.14% 및 block 3의 1.5 mg/L에서 26.2%의 생물학적 방해(inhibition)작용이 발생하였으나 전자빔 조사 후 모든 조건에서 생물학적 방해영향을 나타나지 않았다.

과점산업(寡占産業)에서의 진입제한가격(進入制限價格) (Limit Pricing by Noncooperative Oligopolists)

  • 남일총
    • KDI Journal of Economic Policy
    • /
    • 제12권1호
    • /
    • pp.127-148
    • /
    • 1990
  • 이 논문의 기본목표는 Bain 이래 논란의 대상이 되어 온 기존기업들에 의한 진입제한가격(進入制限價格)의 형성이 일반적으로 가능한가, 가능하다면 어떠한 메커니즘을 통해 형성되는 가를 밝히고, 진입제한가격이론(進入制限價格理論)이 한국경제에 갖는 의의를 찾아보는데 있다, 이 논문에서 밝혀질 주요결과는 다음과 같다. 첫째, 다수의 기존기업(旣存企業)이 각자의 이윤극대화(利潤極大化)를 추구하며 카르텔을 형성하지 않는 때에도 기존기업(旣存企業)과 잠재적(潛在的) 신규기업간(新規企業間)에 정보의 불균형이 존재하는 경우 진입제한가격(進入制限價格)이 채택될 가능성이 있다. 둘째, 이러한 과점기업(寡占企業)들에 의한 진입제한가격형성(進入制限價格形成)은 암묵적 담합의 새로운 형태로 해석할 수 있다. 셋째, 진입제한가격형성(進入制限價格形成)은 각종 회계자료(會計資料)가 공표되지 않을 경우에 가능하다. 넷째, 기존기업(旣存企業)의 수(數)가 증가하여 산업(産業)이 완전경쟁산업(完全競爭産業)에 접근해 감에 따라 진입제한가격(進入制限價格)이 형성될 가능성은 사라지게 된다.

  • PDF

국내 일사량 추정을 위한 Angstrom-Prescott계수의 평가 (Assessment of the Angstrom-Prescott Coefficients for Estimation of Solar Radiation in Korea)

  • 현신우;김광수
    • 한국농림기상학회지
    • /
    • 제18권4호
    • /
    • pp.221-232
    • /
    • 2016
  • 농업 생태계의 생산성을 예측하기 위한 모델의 필수 입력변수인 일사량은 비교적 적은 수의 기상관측소에서만 관측되고 있어, 이들 관측값을 대신하기 위해 일사량을 추정하는 모델들이 사용되고 있다. 특히, 간단한 계수를 사용하여 일조시간을 이용하는 Angstrom-Prescott(AP) 모델이 일사량 추정을 위해 가장 널리 쓰이고 있다. 국내에서 보편적으로 적용가능한 AP모델의 계수값을 탐색하기 위해 국내 20개 기상관측소의 30년간의 일단위 관측자료를 입력자료로 사용하여 경험적으로 얻어진 계수와 Frere and Popov(1979)가 제시한 계수($AP_{Frere}$)를 이용한 일사량을 추정하고, 이들의 신뢰도를 분석하였다. 또한, 전역최적화 과정을 통해 시공간적으로 신뢰도가 높은 일사량을 얻을 수 있는 AP계수의 범위를 탐색하였다. 분석을 위해 월별, 년도별, 지역별로 추정값과 측정값 사이의 일치도를 계산하였다. $AP_{Frere}$를 사용한 결과 작물 생산성 예측을 위한 조건에서 일치도가 높게 나타났다. $AP_{Frere}$를 사용하였을 때 전역최적화를 통해 추정한 AP계수($AP_{max}$)를 사용하였을 경우 보다 일치도가 낮았으나 경험적으로 얻어진 계수($AP_{Choi}$)보다는 일치도가 높은 일사량 추정이 가능하였다. 전역최적화를 통해 일사량 추정치의 신뢰도를 분석한 결과, 신뢰도가 높은 일사량을 얻을 수 있는 AP계수의 범위는 년도별로는 좁게 분포하였으나 월별, 지역별로는 넓게 분포하였다. 그 중에서도 변이가 작은 범위는 지역별 일치도가 월별 일치도보다 넓게 분포하였다. $AP_{Frere}$는 각각의 경우에 대해 일치도가 높고 변이가 작은 범위에 속해 국내 조건에서 $AP_{Frere}$를 적용할 경우, 신뢰도 높은 일사량 추정치를 얻을 수 있을 것으로 예상되었다.

ANN 및 SVR 알고리즘을 활용한 최적 터널지보패턴 선정에 관한 연구 (A study on the optimization of tunnel support patterns using ANN and SVR algorithms)

  • 이제겸;김양균;이승원
    • 한국터널지하공간학회 논문집
    • /
    • 제24권6호
    • /
    • pp.617-628
    • /
    • 2022
  • 터널 건설 시 암반 등급에 따라 다양한 지보재를 적절히 병용하여 지보패턴을 결정하고 시공이 이루어진다. 이 과정에서 시공 경험이 풍부한 전문가의 기술적 판단이 필요한데, 터널 설계의 초기 단계인 타당성 조사 및 기본설계 단계에는 상대적으로 짧은 수행기간과 부족한 자료 및 예산으로 인해 설계에 많은 어려움이 존재한다. 터널 건설의 급증과 함께 축적된 설계 데이터와 머신러닝을 활용한다면, 지보패턴 설계를 보다 신속하고 신뢰도 있게 수행할 수 있다. 따라서 본 연구에서는 암반등급 판정 및 해당 암반등급에 적합한 지보패턴 설계를 자동화하고자 국내 48개 도로터널의 설계자료 및 지반조사 자료를 수집하였으며, 8개의 입력항목(암종, 전기비저항, 심도, 터널연장, 터널연장에 따른 방재등급, 위험도지수에 따른 방재등급, 터널 종류, 터널 단면적)과 11개의 출력항목(암반등급, 숏크리트 제원 2개 항목, 록볼트 제원 3개 항목, 강지보 제원 3개 항목, 콘크리트 라이닝 2개 항목)에 대한 데이터를 정리하였다. 이와 같이 정리된 데이터를 활용하여 2가지 머신러닝 알고리즘(SVM, ANN)을 활용하여 3가지 머신러닝 모델(S1, A1, A2)을 개발하였으며, 세 가지 모델의 성능을 비교해본 결과 출력값의 데이터 형식에 따라 서로 다른 손실함수를 적용한 ANN 기반의 A2 모델이 가장 뛰어난 성능을 보였다. 본 연구를 통해 머신러닝을 활용한 지보패턴 설계의 가능성을 확인할 수 있었으며, 향후 지속적으로 실제 설계에 사용함으로써 단점을 보완하고 적용성을 개선해 나간다면 설계에 보다 큰 도움을 줄 수 있는 지보패턴 설계 모델을 개발할 수 있을 것으로 기대된다.

Word2Vec을 활용한 제품군별 시장규모 추정 방법에 관한 연구 (A Study on Market Size Estimation Method by Product Group Using Word2Vec Algorithm)

  • 정예림;김지희;유형선
    • 지능정보연구
    • /
    • 제26권1호
    • /
    • pp.1-21
    • /
    • 2020
  • 인공지능 기술의 급속한 발전과 함께 빅데이터의 상당 부분을 차지하는 비정형 텍스트 데이터로부터 의미있는 정보를 추출하기 위한 다양한 연구들이 활발히 진행되고 있다. 비즈니스 인텔리전스 분야에서도 새로운 시장기회를 발굴하거나 기술사업화 주체의 합리적 의사결정을 돕기 위한 많은 연구들이 이뤄져 왔다. 본 연구에서는 기업의 성공적인 사업 추진을 위해 핵심적인 정보 중의 하나인 시장규모 정보를 도출함에 있어 기존에 제공되던 범위보다 세부적인 수준의 제품군별 시장규모 추정이 가능하고 자동화된 방법론을 제안하고자 한다. 이를 위해 신경망 기반의 시멘틱 단어 임베딩 모델인 Word2Vec 알고리즘을 적용하여 개별 기업의 생산제품에 대한 텍스트 데이터를 벡터 공간으로 임베딩하고, 제품명 간 코사인 거리(유사도)를 계산함으로써 특정한 제품명과 유사한 제품들을 추출한 뒤, 이들의 매출액 정보를 연산하여 자동으로 해당 제품군의 시장규모를 산출하는 알고리즘을 구현하였다. 실험 데이터로서 통계청의 경제총조사 마이크로데이터(약 34만 5천 건)를 이용하여 제품명 텍스트 데이터를 벡터화 하고, 한국표준산업분류 해설서의 산업분류 색인어를 기준으로 활용하여 코사인 거리 기반으로 유사한 제품명을 추출하였다. 이후 개별 기업의 제품 데이터에 연결된 매출액 정보를 기초로 추출된 제품들의 매출액을 합산함으로써 11,654개의 상세한 제품군별 시장규모를 추정하였다. 성능 검증을 위해 실제 집계된 통계청의 품목별 시장규모 수치와 비교한 결과 피어슨 상관계수가 0.513 수준으로 나타났다. 본 연구에서 제시한 모형은 의미 기반 임베딩 모델의 정확성 향상 및 제품군 추출 방식의 개선이 필요하나, 표본조사 또는 다수의 가정을 기반으로 하는 전통적인 시장규모 추정 방법의 한계를 뛰어넘어 텍스트 마이닝 및 기계학습 기법을 최초로 적용하여 시장규모 추정 방식을 지능화하였다는 점, 시장규모 산출범위를 사용 목적에 따라 쉽고 빠르게 조절할 수 있다는 점, 이를 통해 다양한 분야에서 수요가 높은 세부적인 제품군별 시장정보 도출이 가능하여 실무적인 활용성이 높다는 점에서 의의가 있다.

응력 구속조건을 고려한 동하중의 등가정하중으로의 변환 (Transformation of Dynamic Loads into Equivalent Static Load based on the Stress Constraint Conditions)

  • 김현기;김의영;조맹효
    • 한국전산구조공학회논문집
    • /
    • 제26권2호
    • /
    • pp.165-171
    • /
    • 2013
  • 동하중을 고려하는 구조해석은 전산자원과 시간측면에서 상당한 어려움이 따르기 때문에 외력을 이상적인 정하중으로 가정하는 것이 일반적이다. 그러나 정하중 조건으로 해석된 결과는 구조물의 안전설계 측면에서 충분한 신뢰를 주기 어렵다. 최근에는, 동하중의 영향을 받는 구조물의 효과적인 구조해석을 위해 동하중을 등가정하중으로 변환하는 기법이 제안되어 왔다. 이 기법은 최적화를 통해 구속조건을 만족하는 최소의 등가정하중을 구하는데, 구속조건은 임계시간의 변위를 사용하고, 등가정하중 분포 자유도는 경험적으로 선정하여 왔다. 그러나 안전설계 관점에서는 응력 구속조건을 적용하는 것이 타당하며, 경험적 자유도 선정은 몇 개의 자유도에 과도한 하중이 부과되거나 구조물의 거동에 영향력이 없는 자유도들이 선정될 가능성이 있다. 본 연구에서는 등가응력 구속조건을 고려하는 등가정하중 최적화 방법을 제안하고, 축소시스템 개념을 도입한 주자유도, 구속조건 요소 자유도, 외부하중 자유도로 구성되는 등가정하중 분포 자유도의 구성방법을 제안한다. 수치예제에서는 제안된 방법으로 구해진 등가정하중을 사용하여 등가응력을 구하고 동하중 해석 결과와 비교함으로써 제안된 방법을 통한 구조해석 방법이 구조안전성 측면에서 타당함을 보인다.

Box-Behnken법을 이용한 Geobacter lovleyi의 TCE 탈염소화 공정 최적화 연구 (The Statistical Optimization of TCE Dechlorination by Geobacter lovleyi Using Box-Behnken Design)

  • 차재훈;안상우;천석영;박재우;장순웅
    • 한국지반환경공학회 논문집
    • /
    • 제13권11호
    • /
    • pp.37-42
    • /
    • 2012
  • 본 연구는 TBOS를 기질로 적용한 Geobacter lovleyi의 TCE 탈염소화 공정에 관하여 조사하였다. Geobacter lovleiy를 이용한 TCE 탈염소화 반응은 TCE 초기농도와 Geobacter lovleyi 주입량, 초기 TBOS 농도의 독립변수를 수학적으로 표시하였고, 반응표면법(RSM)을 활용하여 모델화하였다. 이 실험들은 Box-Behnken Design(BBD)을 통해 계획된 15개의 실험을 통해 이루어졌다. RSM을 통하여 TCE 제거효율과 독립변수들의 모델식이 도출되었다 : $Y_1$= -11.50(initial concentration of TCE) + 4.25(단백질 양, 주입된 Geobacter lovleyi의 양) - 4.75(initial concentration of TBOS) - ${6.58X_1}^2$ - ${8.583X_3}^2$ + 93.67, $Y_2=-10.92X_1+5.06X_2-4.89X_3-{4.93X_3}^2-2.19X_1X_2+2.54X_1X_3-2.19X_2X_3+16.71$. 도출된 반응모델은 수정결정계수는 각 0.975, 0.934로 1에 가깝게 나타났으며, 모델의 기여율이 높은 것으로 나타났다. 또한, 통계적 분석결과 TCE 탈염소화 효율에 미치는 영향은 TCE 초기농도 > TBOS 초기농도 > 단백질 양의 순으로 나타났으며, 상호항의 영향은 나타나지 않았다.